二分法求最大似然估计r语言_估计实例R语言: 极大似然估计实例

最近使用开发的过程中出现了一个小问题，顺便记录一下原因和方法--估计实例

同济 习题 7.2

某厂晶体管寿命屈服 E(lamda) 指数分布, Lamda未知, 且Lamda>0, 随机抽取样本寿命如下(小时) 518,612,713,388,434 用极大似然估计其平均寿命

指数分布:

构造似然函数

L(lamda)

f

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

估计和实例

---------------------------------