Mam nadzieję, że ktoś może pomóc w tym, co uważam za stosunkowo proste pytanie, i myślę, że znam odpowiedź, ale bez potwierdzenia stała się ona czymś, czego po prostu nie mogę być pewien.
Mam dane zliczania jako zmienną odpowiedzi i chcę zmierzyć, jak ta zmienna zmienia się wraz z proporcjonalną obecnością czegoś.
Bardziej szczegółowo, zmienną odpowiedzi są zliczenia obecności gatunku owada w wielu miejscach, więc na przykład z tego miejsca pobiera się próbki 10 razy, a gatunek ten może wystąpić 4 razy.
Chcę sprawdzić, czy koreluje to z proporcjonalną obecnością grupy gatunków roślin w ogólnej społeczności roślin w tych miejscach.
Oznacza to, że moje dane wyglądają następująco (to tylko przykład)
Site, insectCount, NumberOfInsectSamples, ProportionalPlantGroupPresence
1, 5, 10, 0.5
2, 3, 10, 0.3
3, 7, 9, 0.6
4, 0, 9, 0.1
Dane obejmują również losowy wpływ na lokalizację.
Myślałem o dwóch metodach, jedną byłby model liniowy ( lmer
) z owadami zamienionymi na proporcje np
lmer.model<-lmer(insectCount/NumberOfInsectSamples~
ProportionalPlantGroupPresence+(1|Location),data=Data)
Drugi to dwumianowy GLMM ( glmer
) np
glmer.model <- glmer(cbind(insectCount,NumberOfInsectSamples-insectCount)~
ProportionalPlantGroupPresence+(1|Location),
data=Data,family="binomial")
Uważam, że dwumianowy błysk jest poprawną metodą, jednak dają one całkiem inne wyniki. Nie mogę znaleźć ostatecznej odpowiedzi w sieci, nie czując się trochę niepewnie, i chcę się upewnić, że się nie mylę.
Każda pomoc lub wgląd w alternatywne metody byłyby bardzo mile widziane.