n <- 1000

pdf(file="qqplot_examples.pdf",width=8,height=6)
par(mfrow=c(2,3))

## 1) gaussian distribution
x <- rnorm(n)*5+10
qqnorm(x)  
## equivalent to above:
##> plot( qnorm( ((1:n)-0.5)/n ),   sort(x))

## 2) exponential distribution
x <- rexp(n)
qqnorm(x)

## 3) double exponential distribution
x <- sample(c(-1,1),n,replace=TRUE) * rexp(n)
qqnorm(x)

## 4) mixture of two gaussians
x <- c( rnorm(n/2), rnorm(n/2)*10)
qqnorm(x)

## 5) t-2 distribution
x <- rt(n,2)
qqnorm(x)

## 6) normal with outliers
x <- rnorm(n)
x[10] <- 34
qqnorm(x)


dev.off()