1
样条混合效果模型
我正在用样条曲线项拟合混合效果模型,该应用程序的随时间变化的趋势是曲线线性的。但是,我要评估的是曲线趋势是否是由于个体偏离线性而出现的,还是在组水平上使组水平拟合显得曲线线性的一种影响。我给出了一个可重现的示例,该示例使JM包中的数据集变得无聊。 library(nlme) library(JM) data(pbc2) fitLME1 <- lme(log(serBilir) ~ ns(year, 2), random = ~ year | id, data = pbc2) fitLME2 <- lme(log(serBilir) ~ year, random = ~ ns(year, 2) | id, data = pbc2) 本质上,我想知道其中哪一个更适合我的数据。但是通过比较anova给了我一个不祥的警告: Model df AIC BIC logLik Test L.Ratio p-value fitLME1 1 7 3063.364 3102.364 -1524.682 fitLME2 …