Esta pregunta ya tiene una respuesta aquí:

Estoy tratando de correranova() en R y se encuentra con alguna dificultad. Esto es lo que he hecho hasta ahora para ayudar a arrojar algo de luz sobre mi pregunta.

Aquí está elstr() de mis datos a este punto.

 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 ...

La columna r es un valor numérico que indica qué fila en el campo reside un gráfico individual. La columna c es un valor numérico que indica qué columna reside un gráfico individual.
Columna Quad corresponde a la ubicación geográfica en el campo en el que reside cada parcela

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)

Me he quedado unlm() como sigue

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

Esto es unanova() para todo el campo, que prueba el rendimiento de grano contra el rendimiento de paja para cada parcela en el conjunto de datos.

Mi problema es que quiero dirigir un individuoanova() para la columna Quad de mis datos para probar el rendimiento de grano y el rendimiento de paja en cada cuadrante.

tal vez unwith() podría arreglar eso. Nunca lo he usado antes y estoy en el proceso de aprender R actualmente. Cualquier ayuda sería muy apreciada.