library(nortest)

a=c(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9)
x=matrix(0,ncol=length(a),nrow=100)


for (i in 1:length(a))  x[,i]=(1-a[i])*rnorm(100)+a[i]*rexp(100)

library(car)
library(nortest)
par(mfrow=c(3,3))

for (i in 1:length(a)) qqPlot(x[,i])

for (i in 1:length(a)) qqnorm(x[,i])

par(mfrow=c(1,1))
# moze faktycznie lepiej na jednej klatce i przewijac strzalka (lewy gorny rog nad rysunkiem)
for (i in 1:length(a)) qqnorm(x[,i]) 


ad=c()
for (i in 1: length(a))   ad[i]=ad.test(x[,i])[2]
ad<-as.numeric(ad)
plot(ad) # pierwsze piec rozkladow zostaje odrzucona watpliwosci co do szostego


li=c()
for (i in 1: length(a))  li[i]=lillie.test(x[,i])[2]
li<-as.numeric(li)
plot(li) # pierwsze piec rozkladow zostaje odrzucona


s=c()
for (i in 1: length(a))  s[i]=shapiro.test(x[,i])[2]
s<-as.numeric(s)
plot(s) # pierwsze piec rozkladow zostaje odrzucona

#itd.