Es geht um vegetationskundliche Aufnahmen. Dabei wurden auf verschiedenen Plots (20 Plots) pro Jahr (2011 - 2014) die Deckungen gemessen. Die Plots unterscheiden sich nach Behandlung. Ungünstigerweise (Planung durch andere -.-) haben wir pro Behandlungsweise unterschiedlich viele Plots, also Behandlung 1: 13 Plots, Behandlung 2: 5 Plots, Behandlung 3: 2 Plots. Für meine Analysen benutze ich R!
Variablen:
- Code: Alles auswählen
$ Row.names :Class 'AsIs' chr [1:83] "1a1" "1a2" "1a3" "1a4" ...
$ shannondiv: num 0.506 1.201 1.379 1.297 0.443 ...
$ Jahr : Factor w/ 4 levels "2011","2012",..: 1 2 3 4 1 2 3 4 1 2 ...
$ Behandlung: Factor w/ 5 levels "Entfernen","Entfernen und Belassen",..: 2 1 1 1 5 5 5 5 3 3 ...
$ Sedspur : num 55 25 30 2 52 1 3 3 37 0.3 ...
usw. (wie sedspur noch viele andere Pflanzenvariablen, das sind Deckungsgrade, interessieren mich aber hier nicht)
$ PlotID : Factor w/ 24 levels "1a","1b","1c",..: 1 1 1 1 2 2 2 2 3 3 ...
1a1 bezeichnet die Fläche in 2011, 1a2 in 2012, 1a3 in 2013, usw.
Ich möchte testen, ob shannondiv durch die Behandlung beeinflusst wird. Meine Random effects stellen denke ich PlotID/Jahr dar.
Ich rechne eine Anova mit den angegebenen Random effects Error(PlotID/Jahr)
Das Modell sieht wie folgt aus:
- Code: Alles auswählen
modelbehandlung <- aov(shannondiv ~ Behandlung+Error(PlotID/Jahr), data=totalanalyse)
Der Fehler sieht folgendermaßen aus:
- Code: Alles auswählen
Warning message:
In aov(shannondiv ~ Behandlung + Error(PlotID/Jahr), data = totalanalyse) :
Error() Modell ist singulär
Habe ich da etwas falsch gerechnet? die Flächen sind unterteilt in Jahre und unterteilt in Plots, es gibt also zum Jahr 2011 einen Plot 5, zu Jahr 2012 einen Plot 5 und zu Jahr 2013 einen Plot 5 (usw.).
Ich gehe davon aus, dass die unterschiedliche Anzahl Plots je Behandlungen erstmal garkeine statistischen Tests rechtfertigen, richtig? ich kann also mit meinen Daten nichts weiter anfangen als zum Beispiel eine NMDS oder DCA durchzuführen, richtig?
Bin mir gerade zu unsicher in diesem Fall. Danke für eure Antworten und wenn ihr zusätzliche Informationen benötigt fragt mich einfach.
beste Grüße,
Lukki
