# # hbstats.R # h<-unlist(read.table("nhb.dat")) summary(h) sd(h)