我的分配遇到问题,该分配会生成100个n = 100的样本进行指数分布。然后,我必须使用最大似然估计(MLE)来估计参数。我已经将问题编码如下,但是遇到了一个错误,指出“子分配类型修复中的类型不兼容(从闭包到双精度)”。我能知道这是什么意思吗?有什么解决的办法吗?
您的回复将不胜感激。
nsam=10^2
minuslogL=rep(1,nsam)
MaxLikeEst=rep(1,nsam)
for(i in 1:nsam){
exp=rexp(100,5)
minuslogL[i]=function(lambda){-sum(dexp(exp,lambda[i],log = TRUE))}
library(stats4)
MaxLikeEst[i]=mle(minuslogL[i],start = list(lambda=5))
summary(MaxLikeEst)
i=i+1
}
Error in minuslogL[i] <- function(lambda) { :
incompatible types (from closure to double) in subassignment type fix