Question

我正在尝试可视化saemix为我计算的模型。我有以下形式的数据（有更多变量，但是这些变量不重要）：

data.frame：377磅。 12个变量中：
  $ PPnumber：因子w /   126级“”，“ FB_001”，“ FB_002”，..：2 2 2 2 3 3 3 4 4 4 ...
  $ ageyrs：num 13.6 14.2 14.6 15.6 14.3 ...
  $ ABSSq：数字37.9 46.2 43.3 51 41.1 ...

我正在尝试针对每个年龄段的10个PP数字计算ABSSq的预测。这样，我应该能够使用ggplot为10位参与者绘制一条漂亮的线条。但是，我发现很难获得saemix预测。

我浏览了saemix文档，发现预测应该可以解决问题，但是我一遍又一遍地遇到相同的错误。

dat<-saemixData(name.data = mdatfull, name.group = "PPnumber",
name.predictors = "ageyrs", name.response = "ABSSq")

#this is my model
AsModel<-function(psi,PPnumber,x){
  t <- x[,1]
  a <- 70
  b <- psi[PPnumber,1]
  g1 <- psi[PPnumber,2]
  fpred <- a - b*exp(-g1*(t-12))
  return(fpred)
}

#and here is the fit
seamodel<-saemixModel(model=AsModel,psi0 = c(b=35,g1=2))
saemix.options<-list(map=TRUE, fim=TRUE, ll.is=FALSE, displayProgress=FALSE, seed=1)
saemix.fit1 <- saemix(seamodel, dat, saemix.options)

#but here is where I run into problems
pframe1SAEM<-with(mdatplotSAEM,expand.grid(PPnumber=unique(PPnumber),ageyrs=seq(12,18,0.1)))
pframe1$ABSSq<-predict(saemix.fit1,pframe1SAEM,type="ipred") ## here the error occurs

我遇到了同样的错误

Error in `[.data.frame`(newdata, "data") : undefined columns selected

我无法弄清楚我到底在做什么错。任何帮助都将不胜感激。

如何获得saemix模型的模型预测？

0 个答案: