Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- m1<-gam(Abundance~s(Temperature)+s(Temperature, by=Site,m=1)+s(Site,k=11,bs="re")+s(Year), method="REML", family = binomial(link="logit"),weights = Trap,data=c2)
- m2<-gamm(Abundance~s(Temperature)+s(Temperature, by=Site,m=1)+s(Site,k=11,bs="re")+s(Year),correlation = corAR1(form = ~ Year|Site), method="REML", family = binomial(link="logit"),weights = Trap,data=c2)
- m3<-gamm(Abundance~s(Temperature)+s(Temperature, by=Site,m=1)+s(Site,k=11,bs="re"),correlation = corAR1(form = ~ Year|Site), method="REML", family = binomial(link="logit"),weights = Trap,data=c2)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement