##### Nested Designs - Nested ANCOVA design


### import data - dataset 2 (Nested ANCOVA)


datum=read.csv(file.choose())
head(datum)


###Analyze data without random effect of individual


results=lm(Size~Age+Sex,data=datum)
summary(results)
anova(results)
### notice inappropriate den degree of freedom
summary(results)


###Analyze data: block by individual, which is nested within sex


results=lme(Size~Age+Sex,data=datum,random=~1|Individual)
summary(results)
anova(results)
### notice appropriate den degree of freedom
### accounts for pseudoreplicated samples (pseudoreplicated by sex)

### Notice that for 'Sex' Individual is the sample
### For 'Age', each measurement in a sample/replicate