R 使用最大似然函数时如何消除optim中的错误
我试图用GEV分布拟合数据。 我的任务是模拟一个排队系统,并实现再生方法,以找到GEV分布的参数,即每个再生周期中的最大队列应符合的参数。 最大值列表命名为 最大值 我尝试在stats4包中使用mle函数 我实现了以下代码,以使数据符合GEV分布:R 使用最大似然函数时如何消除optim中的错误,r,queue,simulation,data-fitting,mle,R,Queue,Simulation,Data Fitting,Mle,我试图用GEV分布拟合数据。 我的任务是模拟一个排队系统,并实现再生方法,以找到GEV分布的参数,即每个再生周期中的最大队列应符合的参数。 最大值列表命名为 最大值 我尝试在stats4包中使用mle函数 我实现了以下代码,以使数据符合GEV分布: LL = function(mu, sigma, gamma) { if (gamma!=0 && (1+gamma*(maximum-mu)/sigma)>0 && sigma!=0) {
LL = function(mu, sigma, gamma) {
if (gamma!=0 && (1+gamma*(maximum-mu)/sigma)>0 && sigma!=0) {
return(sum(log(sigma)+(1+1/gamma)*log(1+gamma*(maximum-mu)/sigma)+(1+gamma*(maximum-mu)/sigma)^(-1/gamma)))
} else
if (gamma==0 && sigma!=0) {
return(sum(exp(-(maximum-mu)/sigma)+(maximum-mu)/sigma+log(sigma)))
} else {
return(0)
}
}
res = mle(minuslogl=LL)
我使用对数似然函数。
然后我犯了一个错误:
optim中的错误(开始,f,method=method,hessian=TRUE,…):
无限值来自“optim”
这一错误的原因是什么?或者是否有其他方法使数据符合GEV分布
提前谢谢你的回答