我有一个基于数据集EM6的lmer模型modelstatus。我根据模型中的预测值在EM6中创建了一个新变量：
modelstatus＆lt; -lmer（MS_TOT~acstatus2 + ExamStage + acstatus2 * ExamStage + AIS.1 + Level.1 + MS_TOT.1 + I +（ExamStage | id），EM6
所以我的结果变量是运动评分（MS_TOT，通常会随着时间的推移而改善）。 ExamStage是时间（以月为单位）。我感兴趣的治疗方法是acstatus2（= 0,1或2）。我还控制了他们的基线运动评分（MS_TOT.1）。然后我找到了这个分数的预测值：
2016-11-07 16:08:05 DEBUG Service - Request: Started
2016-11-07 16:08:05 DEBUG Service - ##########
2016-11-07 16:08:05 DEBUG Service - Response: Completed at Mon 2016.11.07 at 04:08:05 PM EST
2016-11-07 16:08:05 DEBUG Service - IAM-548792. Internal Error
然后，我希望根据治疗组（acstatus2）和时间变量（ExamStage，几个月）看到这些预测值，所以我做了：
EM6$predictedvalues<-predict(modelstatus, EM6)
View(EM6$predictedvalues)
然而，正如您所看到的，三个治疗组之间的基线预测值（在ExamStage = 1）差别很大（48.78,49.81,62.07）。我想要做的是这种精确比较，btu具有相等的基线值（即如果保持不变，治疗0,1或2的影响是什么？）。这可能吗？
谢谢！