估計執行個體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