### ecdf for meltpt of beeswax - Example A p.379

library(Hmisc)

bees = matrix(scan("C:\\Users\\Eric A. Suess\\Desktop\\beeswax.dat"),ncol=2,byrow=T)
bees

meltpt = bees[,1]
phydcarb = bees[,2]

x.bar = mean(meltpt)
x.sd = sd(meltpt)		

z = (meltpt - x.bar)/x.sd		# stansardized values

# The ecdf command plots the ECDF.

X11()
plot(ecdf(z))

# The goodness of fit tests H_0: ecdf = a specified distribution.

shapiro.test(z)

ks.test(z, "pnorm", 0, 1)