r - R の nlme モデルに共変量を追加する

Question

私は最近、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)

ご覧のとおり、2 つの変数 asy と int には 3 つの共変量があります。共変量のすべての可能な組み合わせを生成し、最適なモデルを探したいと思います。これまでのところ、これに関連するトピックは見つかりませんでしたので、役立つアイデアをいただければ幸いです。rコードを投稿に適切に表示する方法を教えてもらえますか?

ありがとう