I=1000 women=0; men=0; propwomen=0; propmen=0 chisqsum=c(1:I) judge=array(1:3, c(7,2, I)) chisq=array(1:3, c(7, 2, I)) #Generate new observed counts for each sample judge[1,1,1:I]=rbinom(I, 354, .261) judge[1,2,]=354-judge[1,1,1:I] judge[2,1,1:I]=rbinom(I, 730, .261) judge[2,2,]=730-judge[2,1,1:I] judge[3,1,1:I]=rbinom(I, 405, .261) judge[3,2,]=405-judge[3,1,1:I] judge[4,1,1:I]=rbinom(I, 226, .261) judge[4,2,]=226-judge[4,1,1:I] judge[5,1,1:I]=rbinom(I, 111, .261) judge[5,2,]=111-judge[5,1,1:I] judge[6,1,1:I]=rbinom(I, 552, .261) judge[6,2,]=552-judge[6,1,1:I] judge[7,1,1:I]=rbinom(I, 597, .261) judge[7,2,]=597-judge[7,1,1:I] #determining gender breakdown for each repetition for (i in 1:I) { women[i]=sum(judge[,1,i]) men[i]=sum(judge[,2,i]) propwomen[i]=women[i]/(women[i]+men[i]) propmen[i]=men[i]/(women[i]+men[i]) } #determine the chisq cell contribution for each cell for (i in 1:I) { chisq[1,1,i]=(judge[1,1, i]-propwomen[i]*354)^2/(propwomen[i]*354) chisq[1,2,i]=(judge[1,2,i]-propmen[i]*354)^2/(propmen[i]*354) chisq[2,1,i]=(judge[2,1,i]-propwomen[i]*730)^2/(propwomen[i]*730) chisq[2,2,i]=(judge[2,2,i]-propmen[i]*730)^2/(propmen[i]*730) chisq[3,1,i]=(judge[3,1,i]-propwomen[i]*405)^2/(propwomen[i]*405) chisq[3,2,i]=(judge[3,2,i]-propmen[i]*405)^2/(propmen[i]*405) chisq[4,1,i]=(judge[4,1,i]-propwomen[i]*226)^2/(propwomen[i]*226) chisq[4,2,i]=(judge[4,2,i]-propmen[i]*226)^2/(propmen[i]*226) chisq[5,1,i]=(judge[5,1,i]-propwomen[i]*111)^2/(propwomen[i]*111) chisq[5,2,i]=(judge[5,2,i]-propmen[i]*111)^2/(propmen[i]*111) chisq[6,1,i]=(judge[6,1,i]-propwomen[i]*552)^2/(propwomen[i]*552) chisq[6,2,i]=(judge[6,2,i]-propmen[i]*552)^2/(propmen[i]*552) chisq[7,1,i]=(judge[7,1,i]-propwomen[i]*597)^2/(propwomen[i]*597) chisq[7,2,i]=(judge[7,2,i]-propmen[i]*597)^2/(propmen[i]*597) } #sum the chisquare contributions across all 14 cells for (i in 1:I){ chisqsum[i]=sum(chisq[,,i]) }