Which normal Q-Q plot is of the real data?

qqfind.R

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)
}