#=========================================
#  histogram of 10000 male heights (18 yr old
#  males, Hong Kong, early 1990s).
#=========================================
HW=read.table("heightweight.txt",header=TRUE)
attach(HW)
hist(height,freq=FALSE)

mu=mean(height)
sig=sd(height)

x=seq(60,75,.01)
y=dnorm(x,mu,sig)
points(x,y,type="l")

windows()
qqnorm(height)

detach(HW)