如何使用mgcv拟合重复测量GAM模型？

问题描述

我们希望得到一些关于使用 mgcv R 包的 gam 模型的指导。我们想知道随时间测量的变量是否会影响我们的结果变量。换句话说，如果变量“X”在控制时间后发生变化。

我们有以下变量：

1. ID - 主题识别
2. 年 - 连续协变量、重复测量
3. Var X - 连续变量
4. Var Y - 连续变量
5. 结果 - 连续变量

我们不完全确定以下模型是否能回答我们的问题：

m1 <- gam(Outcome ~ s(VarX) + s(vary) + s(year) +
      s(year,ID,bs = "fs"),family = gaussian,data = dat,method = 'REML')

或者，如果下面的模型是考虑到每个变量随时间为每个 ID 进行测量的模型：

m2 <- gam(Outcome ~ s(VarX) + s(vary) + s(year) +
     s(year,bs = "fs",by = VarX) +
     s(year,by = vary),method = 'REML')

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）