[问题类型]: 程式谘询(我想用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 --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 39.13.194.176 (台湾) ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1589567739.A.25C.html ※ 编辑: getsimple (39.13.194.176 台湾), 05/16/2020 02:37:59 ※ 编辑: getsimple (39.13.194.176 台湾), 05/16/2020 02:42:42 ※ 编辑: getsimple (39.13.194.176 台湾), 05/16/2020 02:43:57 ※ 编辑: getsimple (39.13.194.176 台湾), 05/16/2020 02:45:04