估計實例R語言: 極大似然估計實例

最近使用開發的過程中出現了一個小問題，順便記錄一下原因和方法--估計實例

    同濟<概率論與數理統計> 習題 7.2

    某廠晶體管壽命屈服 E(lamda) 指數分布, Lamda未知, 且Lamda>0, 隨機抽取樣本壽命如下(小時)
518,612,713,388,434
用極大似然估計其平均壽命

    

    指數分布:

    

    

    構造似然函數

    L(lamda) 

    

f <- function(lamda){

    

   logL = n*log(lamda) - lamda*sum(x)

    

   return (logL)

    每日一道理
風，漸漸吹起，吹亂了我的發絲，也讓我的長裙有些飄動。綠葉仿佛在風中起舞，離開了樹，投向了大地，卻不知這樣會枯萎，我彎下腰，輕輕拾起一片樹葉，那非常有序的莖脈，是一種美的點綴。我有些哀嘆：綠葉啊，綠葉，你這般美麗地從樹上輕輕飄下，隨風起舞，卻不知已被人稱之為落葉！

    

}

    

    

執行  optimize

    

x= c(518,612,713,388,434)

    

n = length(x)

    

optimize(f,c(0,1),maximum = TRUE)

> a = optimize(f,c(0,1),maximum = TRUE)

$maximum

[1] 0.001878689

$objective

[1] -36.39261

平均壽命 =1/lamda

> 1/a$maximum

[1] 532.2862  

文章結束給大家分享下程序員的一些笑話語錄： 騰訊總舵主馬化騰，有人曾經戲稱如果在Z國選舉總統，馬化騰一定當選，因為只要QQ來一個彈窗”投馬總，送Q幣”即可。

--------------------------------- 原創文章 By
估計和實例
---------------------------------