Simulate Z-scores

gistfile1.r

#####SIMULATE & PLOT Z-SCORES#########
 
nSims <- 100000 #number of simulated experiments
p <-numeric(nSims) #set up empty container for all simulated p-values
z <-numeric(nSims) #set up empty container for z-scores
 
for(i in 1:nSims){ #for each simulated experiment
x<-rnorm(n = 23, mean = 100, sd = 20) #produce simulated participants
y<-rnorm(n = 23, mean = 110, sd = 20) #produce simulated participants
t<-t.test(x,y) #perform the t-test
p[i]<-t$p.value #get the p-value and store it
z[i]<-qnorm(1-(t$p.value/2)) #convert to the z-score and store it
}
 
#now plot the histogram NOTE if you get error some 'x' not counted; maybe 'breaks' do not span range of 'x' increase max scale end from 8 to higher value
hist(p, main="Histogram of p-values", xlab=("Observed p-value"), breaks =20)
hist(z, main="Histogram of z-scores", xlab=("Observed z-score"), seq(0,8,by=0.245))