具有匹配对和重复测量的模型

Question

我会删除这是否太松散编程但我的搜索已经显示NULL所以我希望有人可以提供帮助。

我的设计具有案例/控制匹配对设计，并具有重复测量。在R

中寻找模型/功能/包

我2采取time=1衡量标准，2采取time=2衡量标准。我的案例/控件状态为Group（2个级别），匹配的对ID为match_id，并希望估算Group，time和speed上的互动效果1}}，一个连续变量。

我想做这样的事情：

（reg_id是实际参与者ID）

speed_model <- geese(speed ~ time*Group, id = c(reg_id,match_id),
        data=dataforGEE, corstr="exchangeable", family=gaussian)

我希望通过reg_id对人员中的自相关进行建模，但也可以通过match_id

对匹配的对中的自相关进行建模

但我明白了：

Error in model.frame.default(formula = speed ~ time * Group, data = dataFullGEE,  : 
variable lengths differ (found for '(id)')

geese或GEE一般不会处理2组ID的聚类吗？有没有办法做到这一点？我确定有。

感谢您提供任何帮助。

Answer 1

对于Cross Validated来说，这绝对是一个更好的问题，但由于每个主题只有2个观察点，我会考虑ANCOVA模型：

鹅（speed_at_time_2~speed_at_time_1 *组，id = c（match_id），         data = dataforGEE，corstr =＆＃34; exchangeable＆＃34;，family = gaussian）

关于ANCOVA的使用，您可能会发现this reference有用。