Как установить уровень лечения для Даннетта в R

Question

Таким образом, код работает просто отлично, но он использует первый уровень лечения - ABA;Тем не менее, мне нужно установить его в качестве моей переменной Control.Есть ли способ сделать это, или я должен просто переименовать мои переменные, чтобы сначала получить мой элемент управления?

summary(plant) 
library(multcomp) 
grow<-glht(mod4,linfct = mcp(Type = "Dunnett"), 
alternative = "greater") 
summary(grow)

Answer 1

Решение с функцией contrMat:

library(multcomp)
fit <- lm(Sepal.Width ~ Species, data = iris)
groups <- table(iris$Species)
# setosa versicolor  virginica 
#     50         50         50 
contrastMatrix <- contrMat(groups, base = 2) # 2 <=> versicolor
glht(fit, linfct = mcp(Species=contrastMatrix), alternative='less')

---

> glht(fit, linfct = mcp(Species=contrastMatrix), alternative='less')

     General Linear Hypotheses

Multiple Comparisons of Means: User-defined Contrasts

Linear Hypotheses:
                            Estimate
setosa - versicolor >= 0       0.658
virginica - versicolor >= 0    0.204