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]]

