dat<-read.csv("http://www.matsuka.info/data_folder/datWA01.txt") # Male mean.M <-mean(dat$h[dat$gender=="M"]) stddev = 10 n.M = length(dat$h[dat$gender=="M"]) z.value=(mean.M-171)/(sqrt(stddev^2/n.M)) (1-pnorm(abs(z.value)))*2 # Female mean.F <-mean(dat$h[dat$gender=="F"]) stddev = 5 n.F = length(dat$h[dat$gender=="F"]) z.value=(mean.F-158)/(sqrt(stddev^2/n.F)) (1-pnorm(abs(z.value)))*2