Estrazione di elementi di lsmobj dal pacchetto lsmeans in R

Question

Mi chiedo come ottenere l'elemento di lsmobj dal pacchetto lsmeans in R.

require(lsmeans) 
fiber.lm <- lm(strength ~ diameter + machine, data = fiber) 
fiber.lsm <- lsmeans (fiber.lm, "machine") 
fiber.lsm 
machine lsmean  SE df lower.CL upper.CL 
A  40.38241 0.7236252 11 38.78972 41.97510 
B  41.41922 0.7444169 11 39.78077 43.05767 
C  38.79836 0.7878785 11 37.06426 40.53247 

Confidence level used: 0.95 


str(fiber.lsm) 
'lsmobj' object with variables: 
    machine = A, B, C 

Voglio estrarre lsmeans e SE colonne di fiber.lsm. Qualsiasi aiuto sarà molto apprezzato. Grazie in anticipo per il vostro aiuto.

Answer 1

fiber.lsm è un oggetto S4 che non può essere sostituito. Ma puoi usare summary dell'oggetto. Restituisce un frame di dati dei risultati del modello.

s <- summary(fiber.lsm) 

class(s) 
# [1] "summary.ref.grid" "data.frame" 

s[c("lsmean", "SE")] 
# lsmean  SE 
# 40.38241 0.7236252 
# 41.41922 0.7444169 
# 38.79836 0.7878785 

c(s[c("lsmean", "SE")]) 
# $lsmean 
# [1] 40.38241 41.41922 38.79836 
# 
# $SE 
# [1] 0.7236252 0.7444169 0.7878785