##Hierarchical sampling and variance components analysis 
library(lmerTest)
ex.data<-read.table("W:/teaching/stat579/data/hre.txt",header=TRUE)
str(ex.data)
head(ex.data)
names(ex.data)
attach(ex.data)
myfit<-lmer(subject~1+(1|town/district/street/family/gender))
summary(myfit)
#standard error
sds<-c(0.9245321,1.150604,1.131932,1.489864,1.923191,3.917264)
#variance
vars<-sds^2
100*vars/sum(vars)