Факторы повторного заказа llm

Question

Я пытаюсь переупорядочить факторы в простом lm / aov, поэтому я могу изменить точку пересечения.

Я пробовал изменять порядок различных функций, но lm всегда по умолчанию использует тот же самый перехват . Чтобы дать контекст, выращивайте растения с 3 обработками: с высоким, средним и низким содержанием питательных веществ.

DF$treat <- factor( DF$treat , ordered = FALSE )

переупорядочивание

arrange(DF,desc(treat))

и, наконец, создание совершенно нового .csv

Но обработка High всегда является перехватом независимо от порядка. Я хочу посмотреть на эффекты лечения относительно низкого содержания питательных веществ. Я знаю, что могу сделать post ho c с Tukeys, чтобы получить взаимодействия, но мне интересно просто использовать lm-сравнения. Ясно, что это связано с lm / aov (я сделал и то, и другое). Хотите знать, может ли кто-нибудь помочь указать, что мне здесь не хватает? Я знаю, явно что-то очень фундаментальное .....

Answer 1

relevel - ваш друг.

lm(Sepal.Length ~  Petal.Width + relevel(Species, ref="virginica"), iris)
# Call:
#   lm(formula = Sepal.Length ~ Petal.Width + relevel(Species, ref = "virginica"), 
#      data = iris)
# 
# Coefficients:
#                               (Intercept)                                    Petal.Width  
#                                   4.73036                                        0.91690  
# relevel(Species, ref = "virginica")setosa  relevel(Species, ref = "virginica")versicolor  
#                                   0.05009                                       -0.01017  

aov(Sepal.Length ~  Petal.Width + relevel(Species, ref="virginica"), iris)
# Call:
#   aov(formula = Sepal.Length ~ Petal.Width + relevel(Species, ref = "virginica"), 
#       data = iris)
# 
# Terms:
#                 Petal.Width relevel(Species, ref = "virginica") Residuals
# Sum of Squares     68.35344                             0.03460  33.78029
# Deg. of Freedom           1                                   2       146
# 
# Residual standard error: 0.4810113
# Estimated effects may be unbalanced

Конечно, вы также можете использовать его заранее, чтобы получить более чистый результат.

iris <- transform(iris, Species.r=relevel(Species, ref="virginica"))