# Achenbach data in Fig 3.3 data <- read.table("C:/Documents and Settings/David Howell/My Documents/Webs/methods/DataFiles/Fig3-3.dat" , header = T) attach(data) fulldata <- rep(Score, Freq) qqnorm(fulldata, main = "Normal Q-Q of Achenbach's Index", xlab = "Theoretical Quantiles", ylab = "Sample Quantiles", plot.it = TRUE, datax = T) qqline(fulldata, datax = T, col = 2)