客服电话
APP下载
迪恩网络APP
随时随地掌握行业动态
官方微信
扫描二维码
关注迪恩网络微信公众号
|
最近使用开发的过程中出现了一个小问题，顺便记录一下原因和方法--估计实例 同济<概率论与数理统计> 习题 7.2 某厂晶体管寿命屈服 E(lamda) 指数分布, Lamda未知, 且Lamda>0, 随机抽取样本寿命如下(小时)
指数分布:
构造似然函数 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 |
请发表评论