似然方程如图所示: (参数只有一个是k, xi是样本观察值(xi有五个观察值3,1,6,8,7))
以下是我的R Codes, 但运行不出来 , 10个论坛币求答案,哭着感谢
x<-c(3,1,6,8,7)
likelifn<- function(x, k){
return(prod((k^x)*exp(-k)/prod(x)))
}
loglikelifn<-function(x,par){
return(-sum(log(likelifn(x, k = par[1]))) )
}
optimize(par, loglikelifn, x)