#out<-read.table("/Users/reese/teaching/651/data/thetai.out",header=T)
#ttt<-read.table("golfdataR.dat",header=F)

nbin<-29
aa<-(1:nbin)/(nbin+1)
nplayer<-504
bincnts<-rep(0,nbin)
BX2<-rep(0,niter)
for(i in 1:niter){
  bincnts<-0
  for(j in 1:nplayer){
    quants<-qnorm(aa,res[i,j],sqrt(res[i,nplayer+3]))
    quants<-c(0,quants,100)
    bincnts<-bincnts+hist(c(golfdat[pid==j,2],golfdat[pid==j,3],golfdat[pid==j,4],golfdat[pid==j,4],golfdat[pid==j,5]),plot=F,breaks=quants,na.rm=T)$count   
  }
  n<-sum(bincnts)
  BX2[i]<-sum((bincnts-n/nbin)^2/(n/nbin))
}