估计实例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
估计和实例
---------------------------------