我最近开始学习使用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代码在帖子中正确显示吗?
由于