PTT
Submit
Submit
选择语言
正體中文
简体中文
PTT
R_Language
[问题] 使用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
是这样吗
https://pastebin.com/eScbMQck
作者:
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大~~我完食这题了 谢谢大家!已修正 抱歉太粗心了
继续阅读
[问题] 有重复观测值的回归该如何写
leo790124
[问题] 能否ㄧ次汇出多个档案呢? (write.table)
ronald736
[问题] 读取多年维度为矩阵的变量值(nc档)
AndrewShi
[问题] stata 指令转成R
warpdrive
[问题] 将资料时间按月份计算(循环)
AndrewShi
[问题] 选取资料特定区域作图
AndrewShi
[问题] 原始档案为数字,R读入后变成时间格式,求转换法
paranoia5201
Re: [问题] 输出的结果,如何完整贴到文件呢?(求救)
TreeMan
[问题] 限制CPU使用上限?
Edster
[问题] 输出的结果,如何完整贴到文件呢?(求救)
ronald736
Links
booklink
Contact Us: admin [ a t ] ucptt.com