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

bees <- matrix(scan("H:\\hw3\\beeswax.dat"),ncol=2,byrow=T)
bees

meltpt <- bees[,1]
phydcarb <- bees[,2]

x.bar <- mean(meltpt)
x.sd <- stdev(meltpt)		

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

# The command cdf.compare compares the ecdf to a specified distribution.

cdf.compare(z, distribution = "normal")

# The command ks.gof tests H_0: ecdf = a specified distribution.

ks.gof(z, distribution = "normal")