早上好,我应该指定一个重复测量模型,其中我对每个患者进行两次访问(患者数 = 100)。 我在问:
1) 下面的这个模型是否正确指定为重复测量模型?你能写下描述它的数学公式吗?
mmrm <- gls(Outcome~ Treatment+Visit+ TreatmentVisit, data=sim, corr=corSymm(form= ~ 1 |patient_id), weights = varIdent(form = ~ 1 | Visit) )*
另外,如果我应用公式 resid(mmrm) 我得到一个分割输出,前 200 个残差(100 名患者,每次访问两次)都是不同的,而后 200 个都是每次访问相同。为了测试正态性,我会说我只需要考虑输出的前 200 个残差,所以 my_resid<-head(resid(mmrm),100)。
2) 你能不能告诉我这些 my_resid 到底是什么,或者写下数学公式?
非常感谢