Hi everyone, could anybody help me with the code to make a MANOVA with a split-plot design (any help will be very welcomed). The design was as follows:
-there were three blocks, every block was divided in to three treatments, and in every treatment there were different species. I'm thinking species and treatments as factors (fixed both of them), treatments as the plots, species as the subplot...for every specie were some (4) individuals from which some measures were taken. I don't understad how to do the analysis in the nlme package (lme function), I've been traying to reconstruct in that function the analysis of the paper manufacturer from Montgomery's "Desing and analysis of experiments" book, but I understad how to do it. Could someone tell me how to do it in lme function... or any other suggestion... ################################################# #### The code for the analysis in aov function is: #the data: a<-c(30,35,37,36,34,41,38,42,29,26,33,36,28,32,40,41,31,36,42,40,31,30,32,40,31,37,41,40,35,40,39,44,32,34,39,45) temp<-c(rep(c(200,225,250,275),9));temp<-as.factor(temp) pulpa<-c(rep(c(rep(1,4),rep(2,4),rep(3,4)),3));pulpa<-as.factor(pulpa) B<-c(rep(1,12),rep(2,12),rep(3,12));B<-as.factor(B) #aov function that gives in fact same output as shown in the book summary(aov(a~temp*pulpa+Error(B/pulpa+B/temp))) ################################################## here yo can read the paper manufacturer example: https://onlinecourses.science.psu.edu/stat503/node/71 thanks for your attention, and have a nice day Daniel [[alternative HTML version deleted]] _______________________________________________ R-sig-ecology mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
