distribution of means

library(actuar)
library(bigmemory)
mat.length<-10^6
mat1<-matrix(0, mat.length,4)
write.csv(mat1,"maged.answer1.csv")
r.mat<-read.big.matrix("maged.answer1.csv",type="double",
header=F, backingfile="maged.answer1.bin",descriptorfile="maged.answer1.desc")
r.mat<-r.mat[,-1];r.mat<-r.mat[-1,]
colnames(r.mat)<-c('loop#','rlnorm','rllogis','rgamma')
for(i in 1: mat.length){
    r.mat[i,1]<-i
    r.mat[i,2]<-mean(rlnorm(330,-0.0627,0.4085))
    r.mat[i,3]<-mean(rllogis(330,4.3110,1.0653))
    r.mat[i,4]<-mean(rgamma(330,6.080,5.949))
    print(i)
}

length(subset(r.mat$rlnorm,r.mat$rlnorm<=1))/length(r.mat$rlnorm)
length(subset(r.mat$rllogis,r.mat$rllogis <=1))/length(r.mat$rllogis)
length(subset(r.mat$rgamma,r.mat$rgamma <=1))/length(r.mat$rgamma)