Wyodrębnianie elementów lsmobj z pakietu lsmeans w R

Question

Zastanawiam się, jak uzyskać element lsmobj z pakietu lsmeans w 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 

chcę wyodrębnić lsmeans i SE kolumny fiber.lsm. Jakakolwiek pomoc będzie doceniona. Z góry dziękuje za twoją pomoc.

Answer 1

fiber.lsm jest obiektem S4, którego nie można poddać podsekcji. Ale możesz użyć obiektu summary. Zwraca ramkę danych wyników modelu.

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