lsmeans和连续时间

时间:2018-03-28 13:07:39

标签: mixed-models continuous lsmeans

我想帮助解决这个问题。 我希望通过调整临床因素来表达生活质量评分随时间的演变。 我在一个超过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

$对比 对比估计SE df z.ratio p.value (没有)nonEst NA NA NA NA

结果的平均值为:prono1,sex1

非常感谢

1 个答案:

答案 0 :(得分:0)

尝试

lsmeans(mod_mix2, specs ="delai",
    at = list(delai = c(3,6,8,12)))