resid[i]<-mu[i]-moral[i]
llh[i]<- -2*tau *pow((moral[i]-mu[i]),2)


sumllh <-sum(llh[])
mean.resid <-sum(resid[])/N

part1[i]<- pow(mu[i]-mean(moral[]),2)
part2[i]<-pow(moral[i]-mean(moral[]),2)
se[i] <-pow( (mu[i]-moral[i]),2)

r2<-sum(part1[])/sum(part2[])
sumse<-sum(se[])
meanse<-mean(se[])