Fitting lineares Modell / ANOVA nach Gruppe [duplizieren]

Diese Frage hat hier bereits eine Antwort:

Lineare Regression und Gruppieren nach in R 10 Antworten

Ich versuche zu rennenanova() in R und auf Schwierigkeiten stoßen. Dies ist, was ich bis jetzt getan habe, um etwas Licht in meine Frage zu bringen.

Hier ist dasstr() meiner Daten bis zu diesem Punkt.

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

Spalte r ist ein numerischer Wert, der angibt, in welcher Zeile des Feldes sich ein einzelnes Diagramm befindet. Spalte c ist ein numerischer Wert, der angibt, in welcher Spalte sich ein einzelnes Diagramm befindet.
Column Quad entspricht der geografischen Position in dem Feld, in dem sich die einzelnen Diagramme befinden.

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)

Ich habe ein @ gepaslm() wie folg

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

Das ist einanova() für das gesamte Feld, das den Getreideertrag mit dem Strohertrag für jede Parzelle im Datensatz vergleicht.

Mein Problem ist, dass ich eine Person ausführen möchteanova() für die Quad-Spalte meiner Daten zum Testen des Getreideertrags und des Strohertrags in jedem Quadranten.

Vielleicht einwith() könnte das beheben. Ich habe es noch nie benutzt und bin gerade dabei, R zu lernen. Jede Hilfe wäre sehr dankbar.

Antworten auf die Frage(2)

Ihre Antwort auf die Frage