На этот вопрос уже есть ответ здесь:

Я пытаюсь бежатьanova() в R и столкнуться с некоторыми трудностями. Это то, что я сделал до сих пор, чтобы пролить свет на мой вопрос.

Здесьstr() моих данных к этому моменту.

 str(mhw)
'data.frame':   500 obs. of  5 variables:
$ r    : int  1 2 3 4 5 6 7 8 9 10 ...
$ c    : int  1 1 1 1 1 1 1 1 1 1 ...
$ grain: num  3.63 4.07 4.51 3.9 3.63 3.16 3.18 3.42 3.97 3.4 ...
$ straw: num  6.37 6.24 7.05 6.91 5.93 5.59 5.32 5.52 6.03 5.66 ...
$ Quad : Factor w/ 4 levels "NE","NW","SE",..: 2 2 2 2 2 2 2 2 2 2 ...

Столбец r - это числовое значение, указывающее, в какой строке в поле находится отдельный график. Столбец c - это числовое значение, указывающее, в каком столбце расположен отдельный график.
Столбец Quad соответствует географическому положению в поле, к которому относится каждый участок

Quad <- ifelse(mhw$c > 13 & mhw$r < 11, "NE",ifelse(mhw$c < 13 & mhw$r < 11,"NW", ifelse(mhw$c < 13 & mhw$r >= 11, "SW","SE")))
mhw <- cbind(mhw, Quad)

Я подошелlm() следующее

nov.model <-lm(mhw$grain ~ mhw$straw)
anova(nov.model)

Этоanova() для всего поля, которое проверяет урожай зерна против урожайности соломы для каждого участка в наборе данных.

Моя проблема в том, что я хочу управлять человекомanova() для столбца Quad моих данных, чтобы проверить урожай зерна и урожай соломы в каждом квадранте.

возможноwith() может это исправить Я никогда не использовал его раньше, и я нахожусь в процессе изучения R в настоящее время. Любая помощь будет принята с благодарностью.