1

我最近开始学习使用nlme模型拟合数据。我想知道是否有一种更合乎逻辑的方法可以在下面的 nlme 模型中生成所有可能的协变量组合(对于变量asy和):int

fmzeta.nlme <- nlme(
               z ~ asy - (asy - int)*exp(-c*t), 
               fixed = list(asy ~ age + wfa + vl,
                            int ~ age + wfa + vl, c ~ 1),
               random = int + asy ~ 1,
               start = c(-1, 0, -2, 0, 0, 0, 0, 0, 0.03),
               data = DATA)

如您所见,asy 和 int 这两个变量有三个协变量。我想生成所有可能的协变量组合,然后寻找最适合的模型。到目前为止,我还没有找到任何与此相关的主题,因此我将不胜感激任何有用的想法。有人还可以告诉我如何使 r 代码在帖子中正确显示吗?

谢谢

4

0 回答 0