2016-07-07 6 views
0

私はlsmeansで提供されたorangesのデータを使用しています。LSMEANSの信頼限界を抽出するには?

library(lsmeans) 
oranges.rg1<-lm(sales1 ~ price1 + price2 + day + store, data = oranges) 
days.lsm <- lsmeans(oranges.rg1, "day") 
days_contr.lsm <- contrast(days.lsm, "trt.vs.ctrl", ref = c(5,6)) 

信頼区間はplot(contrast(days.lsm, "trt.vs.ctrl", ref = c(5,6)))をplotingによって可視化することができるが、それらは私がdata.frameに信頼区間を抽出することができますどのようにdays_contr.lsm

> days_contr.lsm 
contrast  estimate  SE df t.ratio p.value 
1 - avg(5,6) -7.8538769 2.194243 23 -3.579 0.0058 
2 - avg(5,6) -6.9234858 2.127341 23 -3.255 0.0125 
3 - avg(5,6) 0.2462789 2.155529 23 0.114 0.9979 
4 - avg(5,6) -4.6760034 2.110761 23 -2.215 0.1184 

で示したされていませんか?

> days_contr.lsm 
contrast  estimate  SE df t.ratio p.value lower.CL upper.CL 
1 - avg(5,6) -7.8538769 2.194243 23 -3.579 0.0058 ?  ? 
2 - avg(5,6) -6.9234858 2.127341 23 -3.255 0.0125 ?  ? 
3 - avg(5,6) 0.2462789 2.155529 23 0.114 0.9979 ?  ? 
4 - avg(5,6) -4.6760034 2.110761 23 -2.215 0.1184 ?  ? 

答えて

2

confint(contrast(days.lsm, "trt.vs.ctrl", ref = c(5,6)))細かい

を働きました
関連する問題