[问题类型]:
程式咨询(我想用R 做某件事情,但是我不知道要怎么用R 写出来)
[软件熟悉度]:
新手(没写过程式,R 是我的第一次)
[问题叙述]:
各位好
我目前正在练习使用R计算log-likelihood的mle
*资料如下: (M1、M2都当作failed,None是censored资料)
https://i.imgur.com/yPquzXB.jpg
*目前想要重建这张表格的lognormal那部分
https://i.imgur.com/kJHeGc8.jpg
*遇到的问题是mle不晓得怎么用R跑
我试过optimze跟nlminb,但好像是不能同时求2个参数的mle (?)其实我也不太懂...
但他出现没有引数值
不过我就是要求参数的mle
所以不太知道要怎么解决...
像这样:
https://i.imgur.com/zDfc6Al.jpg
*也有试过想要微分求极值
但好像因为我用plnorm这个函数,所以R不能微分
像这样:
https://i.imgur.com/rKpPvqm.jpg
*主要就是想要求出mle,也就是以下这个等高线最中间那个最高点的位置(x轴是mu、y轴
是sigma)
https://i.imgur.com/0kSfOxP.jpg
*这是我的code :
https://i.imgur.com/co49GK7.jpg
*另外我想求信赖区间
但跑出来也是差了一点
*公式如下:
https://i.imgur.com/wKtn27R.jpg
我其实还没算出mle
所以是先按照表上的值代入
但发现总是差了一点
答案是 (9.91, 10.53)
我算出来的是 (9.92, 10.376)
这是我的code :
https://i.imgur.com/opV6RXd.jpg
扣除0.147是因为95%信赖区间是exp(卡方)=0.147
再用解根的方式求解
红线是0.147的位置
也就是这样:
https://i.imgur.com/SkbYLAE.jpg
上网查了很久...
但试过很多方法我都跑不太出来
拜托各位大大帮帮我
非常谢谢各位~~
手机排版,如不方便阅读请见谅
[程式范例]:
Code网址:
https://ideone.com/V6jzsf
[环境叙述]:
3.6.3
[关键字]:
Likelihood , mle , lognormal