search for: lmm1

...uot;Region", y="Recruit Density"), units = list(y="dm^2")) 1: my first question is this code correctly creating a grouped data object for a nested design with one nested random factor? I than used lme to build my full model with the code; (i had to transform my response) lmm1 <- lme(asinh(density_recruit) ~ region, data=datag, random=list(site=~1)) 2: Is this code correct, i have tried other variations such as lmm1 <- lme(asinh(density_recruit) ~ region, data=datag, random=~1| site) but i essentially get the same results. I can assess significance of my fixed...

...ughing as I type) with a simple answer. Is the multcomp package appropriate for using with a lme built under the nlme package. I know i can get it to work (i.e report p-values for a tukey test) but i am not unsure if this is appropriate for a linear mixed effects model. Below is the code i used lmm1 <- lme(asinh(density_recruit) ~ region, data=datag, random=list(site=~1)) mc<-glht(lmm1, linfct = mcp(region = "Tukey")) Simultaneous Tests for General Linear Hypotheses Multiple Comparisons of Means: Tukey Contrasts Fit: lme.formula(fixed = asinh(density_recruit) ~ regi...