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