GAMM4 时间变量平滑样条
GAMM4 smoothing spline for time variable
我正在构建一个 GAMM 模型（这是第一次）来比较双相情感障碍 (BD) 样本与对照 (HC) 样本的认知表现的纵向斜率。该研究设计被称为“加速纵向研究”，其中参与者年龄跨度很大，年龄在 25-60 岁之间，随访时间为 2 年（HC 组）和 4 年（BD 组）。
假设 (1) BD 组处理速度的年变化率总体上高于健康对照组，表明 BD 组的认知能力下降速度比 HC 组快。
这是我的 R 代码公式，我认为它有点偏离：
'''RUN2 <- gamm4(BACS_SC_R ~ group + s(VISITMONTH, bs = "cc") +
s(VISITMONTH, bs = "cc", by=group), random=~(1|SUBNUM), data=Df, REML = TRUE)'''
visitmonth 变量编码为“从第一次访问开始的月数”，访问 1 将等于 0，随后的访问（每年 3 次）编码为从访问 1 开始经过的月数。在这种情况下，循环平滑是否正确？
我计划在模型中添加额外的变量（即外周炎症）以预测 BD 中认知轨迹的个体斜率。
如有其他建议，将不胜感激。谢谢！
如果 VISITMONTH 超过几年（即对于 BD 观察，我们将在 {0, 1, 2, ..., 48} 中得到 VISITMONTH（四年）），那么不，你不想要循环平滑，除非有一些 4 年的周期性，这意味着 0 和 11 应该被限制为相同。
默认的薄板样条 bs = 'tp' 应该足够了。
我还假设 VISITMONTH 有很多可能的值，因为不是每个人都在相同的每月间隔进行随访？否则你不会有很多自由度可用于时间平滑。
此处组编码为有序因素吗？如果是这样，那就太好了； by smooth 将对参考水平（一定要将 HC 设置为参考水平）和其他水平之间的差异进行编码，以便您可以直接在摘要中看到 BD 差异的测试组.
不清楚您是如何处理 HC 的跟进时间比 BD 组少的事实。看起来该模型有 VISITMONTH 代表研究的全部时间，而不仅仅是 winthin-year 任期。那么，在没有观察到 HC 组的 2 年中，您打算如何比较 BD 组和 HC 组？
我正在构建一个 GAMM 模型（这是第一次）来比较双相情感障碍 (BD) 样本与对照 (HC) 样本的认知表现的纵向斜率。该研究设计被称为“加速纵向研究”，其中参与者年龄跨度很大，年龄在 25-60 岁之间，随访时间为 2 年（HC 组）和 4 年（BD 组）。
假设 (1) BD 组处理速度的年变化率总体上高于健康对照组，表明 BD 组的认知能力下降速度比 HC 组快。
这是我的 R 代码公式，我认为它有点偏离： '''RUN2 <- gamm4(BACS_SC_R ~ group + s(VISITMONTH, bs = "cc") + s(VISITMONTH, bs = "cc", by=group), random=~(1|SUBNUM), data=Df, REML = TRUE)'''
visitmonth 变量编码为“从第一次访问开始的月数”，访问 1 将等于 0，随后的访问（每年 3 次）编码为从访问 1 开始经过的月数。在这种情况下，循环平滑是否正确？
我计划在模型中添加额外的变量（即外周炎症）以预测 BD 中认知轨迹的个体斜率。
如有其他建议，将不胜感激。谢谢！
如果 VISITMONTH 超过几年（即对于 BD 观察，我们将在 {0, 1, 2, ..., 48} 中得到 VISITMONTH（四年）），那么不，你不想要循环平滑，除非有一些 4 年的周期性，这意味着 0 和 11 应该被限制为相同。
默认的薄板样条 bs = 'tp' 应该足够了。
我还假设 VISITMONTH 有很多可能的值，因为不是每个人都在相同的每月间隔进行随访？否则你不会有很多自由度可用于时间平滑。
此处组编码为有序因素吗？如果是这样，那就太好了； by smooth 将对参考水平（一定要将 HC 设置为参考水平）和其他水平之间的差异进行编码，以便您可以直接在摘要中看到 BD 差异的测试组.
不清楚您是如何处理 HC 的跟进时间比 BD 组少的事实。看起来该模型有 VISITMONTH 代表研究的全部时间，而不仅仅是 winthin-year 任期。那么，在没有观察到 HC 组的 2 年中，您打算如何比较 BD 组和 HC 组？