uni/sose2021/r/w02/L02-4-exerc.R

n <- 300
u <- runif(n)
x <- qnorm(u)
y <- rnorm(n)

par(mfrow=c(1,2), mar=c(2,2,1,1))
qqplot(u, y, main="uniform vs normal")
abline(0, 1, col="gray")
qqplot(x, y, main="normal (inversion) vs normal")
abline(0, 1, col="gray")