[问题] 使用R plot likelihood function

楼主: AmigoSafin   2018-11-27 11:38:34
[问题类型]:
程式咨询(我想用R 做某件事情,但是我不知道要怎么用R 写出来)
[软件熟悉度]:
入门(写过其他程式,只是对语法不熟悉)
[问题叙述]:
请简略描述你所要做的事情,或是这个程式的目的
各位大大好
想请教一下如何写R code画出laplace distribution的likelihood function?
p.d.f: f_x(x)=(1/2\theta)exp{(-1/\theta)*|x|}
现在generate a sample with n=100, theta=2
我把function写出来为:
library(VGAM)
theta
<-2;
n <-
100
x <-rlaplace(n,scale=theta)
fn <- function(theta){(1/2*theta)^n*exp((-1/theta)*sum(abs(X)))}
这部分没问题
但接下来该如何用这个function画图就困扰了我
plot(fn)是一条直线
有如我的心电图(眼神死)
还请各位高手帮忙给小鲁一点方向
谢谢~~
作者: HumuHumu (呼姆呼姆)   2018-11-27 15:59:00
作者: ksherry (K雪梨)   2018-11-29 11:39:00
fn 里面 input 没有 n
楼主: AmigoSafin   2018-11-29 21:08:00
我后来做的跟一楼差不多https://reurl.cc/AKn93 像这样 我没有用sort我晚点看看一楼大大的跑出来图是不是一样 还有我不懂二楼大
作者: andrew43 (讨厌有好心推文后删文者)   2018-11-30 00:00:00
没定义x啊原PO。
楼主: AmigoSafin   2018-11-30 00:11:00
x <-rlaplace(n,scale=theta) 这个是吗?
作者: andrew43 (讨厌有好心推文后删文者)   2018-11-30 01:46:00
好吧。尽量写在一起写完整吧。我和ksherry看到的毛病应该是一样的,养成习惯把functio自定义的function需用到的参数都写出来,包括theta,x,n而不是让R自动到外部去帮你找内部未定义的变量在哪里这不是错误,但养成好习惯以后比较不会出错,也方便延用
楼主: AmigoSafin   2018-11-30 10:30:00
我知道了!!谢谢A大~~我完食这题了 谢谢大家!已修正 抱歉太粗心了

Links booklink

Contact Us: admin [ a t ] ucptt.com