我想计算正态分布的均值μ的最大似然估计,方差= 36。
但以下程序显示错误:
set.seed(36)
x <- rnorm(50,mean=20,sd=6)
loglikelikelihood<- function(mu,x){
L = (-50)*log(sqrt(2*pi)*6)-((1/72)*(sum(x^2)-2*mu*sum(x)+50*mu^2))
}
optimize(f=loglik,X=x,interval=c(0,100),maximum=T)
答案 0 :(得分:4)
loglikelihood
而不是loglik
和x
而不是X
,它会起作用。
> optimize(f=loglikelihood,x=x,interval=c(0,100),maximum=T)
$maximum
[1] 20.86679
$objective
[1] -157.6814