c=dim(controldata1)[2]
r=dim(controldata1)[1]
for(i in 1:r){
controldata1[i,c+1]<-sum(controldata1[i,1:c])
}
var=c(0)
for(i in 1:c+1){
(var[i]<-var(controldata1[,i]))
}
a<-(c/(c-1))*(1-sum(var[1:c])/var[c+1])