如何解决如何使用mgcv拟合重复测量GAM模型?
我们希望得到一些关于使用 mgcv R 包的 gam 模型的指导。我们想知道随时间测量的变量是否会影响我们的结果变量。换句话说,如果变量“X”在控制时间后发生变化。
我们有以下变量:
- ID - 主题识别
- 年 - 连续协变量、重复测量
- Var X - 连续变量
- Var Y - 连续变量
- 结果 - 连续变量
我们不完全确定以下模型是否能回答我们的问题:
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 举报,一经查实,本站将立刻删除。