Hallo,
Da fallen mir Dinge ins Auge:
1. Die ratio ist als nominale Variable (factor) hinterlegt. Werte wie "0,5" legen nahe, dass man das auch metrisch betrachten könnte. Das hängt ein wenig von Deiner Einschätzung der Sachlage zusammen, ob Du das nominal oder metrisch haben willst, aber das musst Du aktiv entscheiden, bevor Du was rechnen willst.
2. Die Zahl der lebenden Zellen schein an Tag 0 häufig 10 zu sein. Ist 10 so etwas wie eine logische Obergrenze, können es also höchstens zehn sein, weil nur zehn anfangs eingesetzt wurden? Das würde dann nicht einem Poisson-Modell, sondern einem Binomial-Modell entsprechen. Hängt vom Versuchsaufbau ab.
3. Wenn Dich Tag 0 und Tag 3 jeweils unabhängig voneinander interessieren, dann hast Du vor dem Rechnen sicher dana gefiltert, oder? Oder sollen die Tage mit ins Modell genommen werden um eine gemeinsame Regression zu rechnen?
zu 1. das macht tatsächlich nur als nominale Variable Sinn, da ich die Verhältnisse beliebig gewählt habe und nicht variieren können.
zu 2. an Tag 1 wurden zehn Räuber in das System eingesetzt (das hatte ich nicht explizit gesagt, entschuldige bitte), daher ist der Anfangswert immer 10. Dann wurde am 5. Tag gezählt (sorry, hatte vorher Tag 3 geschrieben, das war ein anderes Experiment) und es können sowohl mehr als auch weniger sein (oder gar keine, falls alle gestorben sind).
zu 3. eigentlich interessiert mehr Tag 5 als 0, da an Tag 0 ja alles von mir festgelegt wurde. Beide zusammen hingegen will ich mir auch ansehen, um das Wachstum zu untersuchen.
"So ähnlich" ist immer schwierig. Wenn Du Ergebnisse zeigst, schreib bitte den genauen Aufruf von GLM dazu, sonst führt das früher oder später garantiert zu Missverständnissen.
ich habe durch ein Buch die Formel zusammenbekommen:
was genau genommen nicht nur ein ähnliches, sondern exakt das gleiche Ergebnis bringt (was mich etwas beruhigt, also habe ich zumindest den Grundaufbau verstanden).
Das summary Deines Modells ist sicher wesentlich umfangreicher als diese vier Zeilen.
Ich habe da tatsächlich wohl etwas abgeschnitten:
Code: Alles auswählen
Call:
glm(formula = alive ~ clone * ratio, family = poisson(), data = dat)
Deviance Residuals:
Min 1Q Median 3Q Max
-4.3205 -0.7353 -0.0524 0.9099 3.5529
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 2.25129 0.13245 16.997 < 2e-16 ***
cloneC3 0.08408 0.18350 0.458 0.64680
cloneC4 0.05129 0.18496 0.277 0.78153
cloneC7 0.21963 0.17784 1.235 0.21685
cloneS1 0.23361 0.17729 1.318 0.18761
cloneS6 0.48955 0.16822 2.910 0.00361 **
ratio0.5 -0.37949 0.20781 -1.826 0.06783 .
ratio1 -0.35417 0.20626 -1.717 0.08596 .
ratio2 -0.60906 0.22316 -2.729 0.00635 **
cloneC3:ratio0.5 0.10250 0.28388 0.361 0.71805
cloneC4:ratio0.5 0.42828 0.27518 1.556 0.11962
cloneC7:ratio0.5 0.46062 0.26506 1.738 0.08225 .
cloneS1:ratio0.5 0.11016 0.27433 0.402 0.68801
cloneS6:ratio0.5 0.28954 0.25632 1.130 0.25865
cloneC3:ratio1 -0.08408 0.28926 -0.291 0.77130
cloneC4:ratio1 0.44948 0.27269 1.648 0.09929 .
cloneC7:ratio1 0.31100 0.26708 1.164 0.24425
cloneS1:ratio1 0.29701 0.26672 1.114 0.26546
cloneS6:ratio1 0.28748 0.25454 1.129 0.25872
cloneC3:ratio2 0.09285 0.30489 0.305 0.76072
cloneC4:ratio2 0.62559 0.28786 2.173 0.02976 *
cloneC7:ratio2 0.45726 0.28334 1.614 0.10656
cloneS1:ratio2 -0.05668 0.30119 -0.188 0.85072
cloneS6:ratio2 0.10182 0.28002 0.364 0.71615
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 427.01 on 143 degrees of freedom
Residual deviance: 315.28 on 120 degrees of freedom
AIC: 918.97
Number of Fisher Scoring iterations: 5
Das ist ein biologisches Forschungsthema aber noch keine statistische Hypothese. Erst wenn die steht, lässt sich abschließend sagen, wie gut ein glm passt.
Meine Hypothesen für diesen Versuch sind die: Die Räuber der verschiedenen Algenklone reagieren unterschiedlich stark auf die Anwesenheit von Plastik. C7 sind die Algen, die am besten wachsen und damit können die Räuber damit am besten überleben. Das Plastikverhältnis 2 hat den größten negativen Einfluss auf die Räuber.
Ich hoffe, das macht es ein bisschen verständlicher.
Vielen lieben Dank!