R 中的无限 AIC - 编程之家

如何解决R 中的无限 AIC

我正在尝试使用 R 的 SAEmix 包执行非线性混合效应模型，但我不断收到 -2LL、AIC 和 BIC 的 INF。我试图模拟“k”对“p”的影响，数据集按“id”分组。模型函数是a*exp(-k/b)，是连续的。但是，随机参数 a 和 b 的估计是合理的，所以我不确定出了什么问题。代码如下所示：

library(SAEmix)
SAEmix.data <- SAEmixData(name.data = dataset,name.group = "id",name.predictors = "k",name.response = "p")
model1 <- function(psi,id,x){
  k <- x[,1]
  a <- psi[id,1]
  b <- psi[id,2]
  
  ypred <- a*exp(-k/b)
  return (ypred)
}
SAEmix.model <- SAEmixModel(model = model1,psi0  = c(a=100,b=50000))
SAEmix.options <- list(map=TRUE,fim=TRUE,ll.is=FALSE,displayProgress=FALSE,seed=632545)
SAEmix.fit1    <- SAEmix(SAEmix.model,SAEmix.data,SAEmix.options)
SAEmix.fit1@results