[問題] 對數常數mle及信賴區間

看板R_Language作者getsimple (getsimple)時間5年前 (2020/05/16 02:35)推噓3(3推 0噓 4→)

留言7則, 3人參與討論串1/1

[問題類型]: 程式諮詢(我想用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://www.ptt.cc/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

推

tnlomtt

05/16 12:54, 5年前 , 1^F

05/16 12:54, 1^F

推

empireisme

05/17 23:59, 5年前 , 2^F

05/17 23:59, 2^F

→

getsimple

05/18 13:51, 5年前 , 3^F

05/18 13:51, 3^F