search for: density_recruit

...360. Thus i have as a design Response: density Factors: regions (3 levels, fixed) sites (18 levels, random) I started with creating a dataframe that included density, region (coded 1, 2, 3), and site (coded 1,2...18) I then created a grouped data object with the code; datag <- groupedData(density_recruit ~ region | site, data=data, labels = list(x="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 wi...

...h 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) ~ region, data = datag, random = l...