Posted on by

データ解析基礎論a W05


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