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

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)