我正在使用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$lsmean
和lme.lms$SE
之类的操作,但是当我尝试此操作时,会出现以下错误:
Error in subject.lme2.lms$lsmean : $ operator not defined for this S4 class
我已经在线查看并发现lsmeans输出属于“lsmobj”类,但我不知道如何操纵该类,似乎无法通过我读过的内容弄明白。
答案 0 :(得分:1)
如果你这样做
summ = summary(lme.lsm)
然后summ
继承自data.frame
,您可以根据需要绘制它。