Which normal Q-Q plot is of the real data?
qqfind <- function(x) { | |
# Generates the normal Q-Q plot for the given data | |
# as well as for three samples of the same size | |
# from a normal distribution with matching | |
# mean and variance. | |
# | |
# If it's obvious which one is the real data, | |
# it might not be normally distributed. | |
# | |
# Return value is the index of the plot | |
# of the real data. | |
mean <- mean(x) | |
sd <- sd(x) | |
par(mfrow=c(2, 2)) | |
real <- sample(1:4, 1) | |
for (i in 1:4) { | |
if (i == real) { | |
qqnorm(x) | |
qqline(x) | |
} else { | |
y <- rnorm(length(x), mean, sd) | |
qqnorm(y) | |
qqline(y) | |
} | |
} | |
par(mfrow=c(1, 1)) | |
return(real) | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment