如何在R中的lsmeans函数的输出中保存/引用估计值和SE

时间:2014-11-07 16:05:30

标签: r mixed-models random-effects

我正在使用lsmeans函数调查数据中的时间依赖性:

lme=lme(attraction~factor(time),random=~1|id, data=na.exclude(subject))  
lme.lms=lsmeans(lme, "time")
summary(lme.lms)

time    lsmean         SE df asymp.LCL asymp.UCL  
   1 0.5823399 0.01805961 NA 0.5469394 0.6177403  
   2 0.5662435 0.01805961 NA 0.5308430 0.6016439  
   3 0.5225464 0.01805961 NA 0.4871459 0.5579468  
   4 0.4938745 0.01805961 NA 0.4584740 0.5292750  
   5 0.4884408 0.01805961 NA 0.4530403 0.5238412  
   6 0.5079754 0.01805961 NA 0.4725749 0.5433758  
   7 0.4521263 0.01805961 NA 0.4167258 0.4875268  
   8 0.4604106 0.01808727 NA 0.4249559 0.4958653  

 ![plot(lme.lms)][1]

我需要移调此图的X轴和Y轴,以便lsmean位于Y轴上,时间位于X轴上。我不知道如何做到这一点,而没有找到一种方法来保存lsmean估计和SE输出作为他们自己的对象。但是,我无法弄清楚如何做到这一点。

对于其他功能,我可能会执行lme.lms$lsmeanlme.lms$SE之类的操作,但是当我尝试此操作时,会出现以下错误:

Error in subject.lme2.lms$lsmean : $ operator not defined for this S4 class

我已经在线查看并发现lsmeans输出属于“lsmobj”类,但我不知道如何操纵该类,似乎无法通过我读过的内容弄明白。

1 个答案:

答案 0 :(得分:1)

如果你这样做

summ = summary(lme.lsm)

然后summ继承自data.frame,您可以根据需要绘制它。