我想在这个问题上得到你的帮助。我希望通过调整临床因素来表达生活质量评分随时间的演变。我在一个由 1000 多名患者组成的队列中工作,每 3 个月提供一次生活质量表。时间被认为是连续的,因为患者没有完全同时填写表格。
这是我的问题,我使用了具有固定因素(性别和预后)的线性混合模型。我想从这个模型中获得 3、6、9 和 12 个月的平均分数。我使用了 lsmeans 函数,但获得的分数对应于平均时间。如何在这 4 个不同时间从我的模型中获得平均分数?您将在我的代码下方找到我使用 lsmeans 得到的结果
mod_mix2 <- lme(score_utilite~delai+prono1+sex1,random = ~ delai|numero_patient, data = qdv,na.action=na.omit,method="ML") lsmeans(mod_mix2, specs="delai")
$lsmeans delai lsmean SE df lower.CL upper.CL 10.21976 0.8145542 0.005835597 1016 0.803103 0.8260054
结果在以下水平上取平均值:prono1、sex1 使用的置信水平:0.95
$contrasts 对比估计 SE df z.ratio p.value (nothing) nonEst NA NA NA NA
结果在以下水平上平均:prono1,sex1
非常感谢