#Generate some fake data for the analysis#
#alpha was generated from N(0, 5), beta=10, random error
#was generated from N(0, 2)

sigma<-sqrt(2)
tao<-sqrt(5)
beta<-10
l<-7
m<-5
alpha<-rep(NA, l)
alphasq<-rep(NA, l)
y<-rep(NA, l*m)
er<-rep(NA, l*m)
dim(y)<-c(l,m)
dim(er)<-c(l,m)
ybar<-rep(NA, l)
betah<-0
sumer<-0
alphasq<-0
for (i in (1:l)) {
  ybar<-0
  alpha[i]<-signif(rnorm(1, 0, tao),4)
  alphasq<-(alphasq+alpha[i]*alpha[i])
    for (j in (1:m)) {
       er[i,j]<-signif(rnorm(1, 0, sigma), 4)
       y[i,j]<-beta+alpha[i]+er[i,j]
       sumer<-sumer+er[i,j]*er[i,j]
       }
  ybar<-(sum(y[i,])/m)
  betah<-(betah+sum(y[i,]-alpha[i]))
  write((y[i,][!is.na(y[i,])]), file="obs4.dat",ncolumns=length(y[i,][!is.na(y[i,])]),append=TRUE)
  }