2

I'm looking at Springtail community analysis and want to see what, if any, affect Litter Type (LT), Site Type (ST), and Sampling Day (SD) have on springtail abundance. LT and ST both have 3 factors each while SD has 2.

I'm still learning the ropes with R, but initially thought the best way to do this was along the lines of:

model1 <- glm(Abundance ~ LT * ST * SD, data=mydata, family=poisson())

summary(model1)
anova(model1, test="Chisq")

Does this seem logical? Am I on the right track? Further more I'd like to perform some sort of posthoc TUKEYHSD test to get more information on any variables which DO significantly affect Springtail Abundance. Is this possible?

Ferdi
  • 5,179
O.Cowan
  • 21
  • 1
  • Could you explain what "abundance" means, scientifically, and how that is measured by your response variable Abundance? – whuber Nov 05 '18 at 14:01
  • 1
    "Abundance" here is the total number of Springtails recorded in each litter trap. Each litter trap (n=102) has an abundance value of the total number of individuals counted in the trap (values range from 1 to 350) – O.Cowan Nov 05 '18 at 15:53
  • Have a look at my answers here and here. Maybe that will help you to get further. – Stefan Nov 06 '18 at 02:08

0 Answers0