模型匹配对和重复措施

我会删除，如果这太松散编程但我的搜索已变成NULL所以我希望有人可以提供帮助。模型匹配对和重复措施

我有一个设计，有一个案例/控制配对设计与重复测量。在R

寻找一种模式/功能/包我在time=12措施，并在time=22措施。我有Case/Control状态为Group（2级），并且配对id为match_id，并且想估计Group,time和speed（一个连续变量）上的交互作用。

我想要做这样的事情：

（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的集群吗？有没有办法做到这一点？我确定有。

感谢您提供任何帮助。

来源

2015-04-17 Andrew Taylor

现在想想大声，但我可以只使用'id = match_id'，因为那也会在人物自相关内部捕捉到匹配之内？或者我说废话。 –

或者我只需要执行'id = reg_id'，然后将用于匹配的变量（例如'年龄习惯性体重高度'）添加到模型中。这是一些同事所推荐的。 –

这是definatly为Cross Validated一个更好的问题，而是因为你有每门课正好2周的观察，我会考虑的协方差分析模型：

鹅（speed_at_time_2〜speed_at_time_1 *集团，ID = C（match_id），数据= dataforGEE，corstr =“可交换”，family =高斯）

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

来源

2015-04-17 16:23:45 Jthorpe

那么我没有2个观察每个主题，但两个观察和两个时间点（4总观察）。这是否改变了模型/代码？谢谢！ –

而你的编码方式，我需要根据时间将数据从“长”扩展到“宽”，正确吗？ –

关于您的第二条评论，是的，您需要传播数据，以便您有一列使用'speed_at_time1'，另一列使用'speed_at_time2'。关于您的第一条评论，您是否在第一时间针对每个主题采取了两项措施，因为您的测量中存在一定数量的错误，并且您想通过多次测量或其他原因提高准确度？如果答案是肯定的，那么平均两次测量值可能是有意义的，这样您就可以对每个主题每个时间点的速度进行一次（更准确的）估计。 – Jthorpe

模型匹配对和重复措施

回答

相关问题