Displaying 1 result from an estimated 1 matches for "treatcod".
Did you mean:
treatcode
2013 May 02
1
multivariate, hierarchical model
...ated so i've decided to include that in the final model. Positively skewed
I'm thinking of doing a mixed effects model with random intercepts as the treatment code and zones within the treatment cose and random slopes as the seasons.
mod2 <- lme(tthm ~ cl2free, random= ~ seasons| treatcode/loc_code)
but that doesn't work.
these seems to work good:
mod2 <- lme(tthm ~ cl2free, random= ~ 1| loc_code, data=new.data, method="ML")
mod3 <- lme(tthm ~ cl2free, random= ~ 1| treatcode/loc_code, data=new.data, method ="ML")
mod2 has a lower...