作者xx801114 (獅子)
看板Statistics
標題[程式] 費雪信息矩陣計算
時間Sat May 16 23:07:38 2020
[軟體程式類別]:R
[程式問題]:
計算費雪信息矩陣
[軟體熟悉度]:
新手
[問題敘述]:
可靠度的資料,採用Lognormal來計算
資料分別有Exact failure time跟censored .
把likelihood function 打出來後，想取Log後利用微分來計算費雪信息數，但是R好像沒
辦法直接對lognormal的pdf及cdf微分，想請問各位大大，有什麼方式可以計算費雪信息
矩陣？ 謝謝大家
[程式範例]:
寫出來的likelihood function:
likelihood=function(x,y){
dlnorm(6700,x,y)*(1-plnorm(6950,x,y))*(1-plnorm(7820,x,y))*(1-plnorm(8790,x,
y))*dlnorm(9120,x,y)*(1-plnorm(9660,x,y))*(1-plnorm(9820,x,y))*
(1-plnorm(11310,x,y))*(1-plnorm(11690,x,y))*(1-plnorm(11850,x,y))*(1-plnor
m(11880,x,y))*(1-plnorm(12140,x,y))*dlnorm(12200,x,y)*
(1-plnorm(12870,x,y))*dlnorm(13150,x,y)*(1-plnorm(13330,x,y))*(1-plnorm(13
470,x,y))*(1-plnorm(14040,x,y))*dlnorm(14300,x,y)*
dlnorm(17520,x,y)*(1-plnorm(17540,x,y))*(1-plnorm(17890,x,y))*(1-plnorm(18
450,x,y))*(1-plnorm(18960,x,y))*(1-plnorm(18980,x,y))*
(1-plnorm(19410,x,y))*dlnorm(20100,x,y)*(1-plnorm(20100,x,y))*(1-plnorm(20
150,x,y))*(1-plnorm(20320,x,y))*dlnorm(20900,x,y)*dlnorm(22700,x,y)*
(1-plnorm(23490,x,y))*dlnorm(26510,x,y)*(1-plnorm(27410,x,y))*dlnorm(27490
,x,y)*(1-plnorm(27890,x,y))*(1-plnorm(28100,x,y))
}
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 39.9.105.183 (臺灣)
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1589641660.A.36E.html
1F：推 ksherry: 可以考慮 numerical Hessian matrix 05/17 20:21
2F：→ xx801114: 收到，謝謝您 05/18 09:28
3F：推 cuylerLin: 對數常態的信息矩陣理論上推導得出來啊，有closed form 05/18 23:32