X <- read.table("/home/longa/www/modules/digitize/lab/snowRweb.txt", header=T) 
attach(X) 
names(X) 
y <- X$rate
z <- X$trans

hist(y,main="Untransformed Rate")
qqnorm(y,main="Untransformed Rate");
qqline(y, col = 2)

hist(z, main="Transformed Rate")
#     hist.default(x, breaks, freq = NULL, probability = !freq,
#          include.lowest = TRUE, right = TRUE,
#          col = NULL, border = par("fg"),
#          main = paste("Histogram of" , deparse(substitute(x))),
#          xlim = range(breaks),  ylim = range(counts, 0),
#          xlab = deparse(substitute(x)),  ylab,
#          axes = TRUE, plot = TRUE, labels = FALSE, ...)
 
qqnorm(z, main="Transformed Rate");
#     qqnorm(y, ylim, main = "Normal Q-Q Plot", xlab = "Theoretical Quantiles",
#            ylab = "Sample Quantiles", plot.it = TRUE, ...)

qqline(z, col = 2)