最近使用开发的过程中出现了一个小问题,顺便记录一下原因和方法--估计实例
同济 习题 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
估计和实例
---------------------------------