Неверный аргумент в модели lsmeans?

Question

Я получаю сообщение об ошибке

lsmDrugs=lsmeans(model1, -drug) Ошибка в -drug: неверный аргумент для унарного оператора "

Любая помощь, пожалуйста?

library(knitr)
library(lsmeans)

drug <- c(rep("A", 2),rep("B", 2),rep("C", 2))
drugDiff <- c(-14, -4, 5, -1, -2, 6)

data<-data.frame(Drug=as.factor(drug),Difference=drugDiff)
data

model1 <- aov(drugDiff ~ drug, data = data)

anova(model1)
kable(anova(model1), format = "markdown")

lsmDrugs = lsmeans(model1, -drug)

summary(contrast(lsmDrugs, method="pairwise", 
                 adjust = "tukey"), infer = c(T, T), level = 0.95, side = "twosided")

Answer 1

У вас есть опечатка, как вы должны использовать тильду - ~ вместо - charachter. Тильды используются для формул в R. Пожалуйста, смотрите код ниже:

library(knitr)
library(lsmeans)

drug <- c(
  rep("A", 2),
  rep("B", 2),
  rep("C", 2)
)

drugDiff <- c(-14, -4, 5, -1, -2, 6)

data <- data.frame(Drug = as.factor(drug), Difference = drugDiff)
data

model1 <- aov(drugDiff ~ drug, data = data)

anova(model1)
kable(anova(model1), format = "markdown")

lsmDrugs = lsmeans(model1, ~ drug) # tilde instead of minu, used for formula

summary(
  contrast(lsmDrugs, method = "pairwise",
           adjust = "tukey"),
  infer = c(T, T),
  level = 0.95,
  side = "twosided"
)

Выход:

 contrast      estimate       SE df lower.CL upper.CL t.ratio p.value
 A - B    -1.100000e+01 5.773503  3 -35.1261  13.1261  -1.905  0.2817
 A - C    -1.100000e+01 5.773503  3 -35.1261  13.1261  -1.905  0.2817
 B - C    -1.776357e-15 5.773503  3 -24.1261  24.1261   0.000  1.0000

Confidence level used: 0.95 
Conf-level adjustment: tukey method for comparing a family of 3 estimates 
P value adjustment: tukey method for comparing a family of 3 estimates