remove(list=ls())
setwd("/Users/michael/Dropbox/Angewandte Statistik SoSe11/Uebungen/Blatt04/R")
#setwd("/Volumes/MichaelsHD/Dropbox/Angewandte Statistik SoSe11/Uebungen/Blatt04")
set.seed(23514)

#250 Realisierungen von X
X <- rchisq(250,df=5)
X

#250 Realisierungen von Y
Y <- rchisq(250,df=13)
Y

#berechne Q
Q <- X/Y*13/5
Q

#250 Realisierungen von Z
Z <- rf(250,df1=5,df2=13)
Z


#plotte direkt in .pdf
pdf("./QQPlots.pdf")
op <- par(mfrow=c(2,2))

#QQ-Plot von Q gegen N(0,1)
qqnorm(Q)
#Gerade hinzufügen
qqline(Q)

#QQ-Plot von Q gegen Z
qqplot(Z,Q,xlab="Quantile von Z",ylab="Quantile von Q",main="QQ-Plot von Z und Q")

par(op)
#schliesse .pdf
dev.off()