Hi there
I've fitted the following model:








mod15.new <- groupedData(TTHM ~ PH + CL2_FREE + treatzone| Year/seasons, data=as.data.frame(new.data2), FUN=mean, labels = list(x="Year/seasons", y= "Year, TTHM"))







fml2 <- lme(mod15.new)
And plotted it. The BLUPs returned are strange: there's clusters of predictions aggregated on either side of the slopes. Not sure why this is. Any ideas?I could attach the plot if you'd like to see. 
Many thanks 		 	   		  
	[[alternative HTML version deleted]]

