作者AmigoSafin ()
看板R_Language
标题[问题] 使用R plot likelihood function
时间Tue Nov 27 11:38:34 2018
[问题类型]:
程式谘询(我想用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)是一条直线
有如我的心电图（眼神死）
还请各位高手帮忙给小鲁一点方向
谢谢～～
--
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 8.41.66.212
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1543289917.A.5EB.html
2F：推 ksherry: fn 里面 input 没有 n 11/29 11:39
3F：→ AmigoSafin: 我後来做的跟一楼差不多 11/29 21:08
5F：→ AmigoSafin: 看看一楼大大的跑出来图是不是一样 还有我不懂二楼大 11/29 21:10
6F：→ andrew43: 没定义x啊原PO。 11/30 00:00
7F：→ AmigoSafin: x <-rlaplace(n,scale=theta) 这个是吗？ 11/30 00:11
8F：→ andrew43: 好吧。尽量写在一起写完整吧。 11/30 01:46
9F：→ andrew43: 我和ksherry看到的毛病应该是一样的，养成习惯把functio 11/30 02:10
10F：→ andrew43: 自定义的function需用到的参数都写出来，包括theta,x,n 11/30 02:11
11F：→ andrew43: 而不是让R自动到外部去帮你找内部未定义的变数在哪里 11/30 02:13
12F：→ andrew43: 这不是错误，但养成好习惯以後比较不会出错，也方便延用 11/30 02:22
13F：→ AmigoSafin: 我知道了！！谢谢A大～～我完食这题了 谢谢大家！ 11/30 10:30
※ 编辑: AmigoSafin (129.21.70.200), 11/30/2018 10:32:02
14F：→ AmigoSafin: 已修正 抱歉太粗心了 11/30 10:32