Don't like ads? PRO users don't see any ads ;-)
Public Pastes
Guest

distribution of means

By: organofcorti on Feb 29th, 2012  |  syntax: R  |  size: 0.74 KB  |  hits: 83  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. library(actuar)
  2. library(bigmemory)
  3. mat.length<-10^6
  4. mat1<-matrix(0, mat.length,4)
  5. write.csv(mat1,"maged.answer1.csv")
  6. r.mat<-read.big.matrix("maged.answer1.csv",type="double",
  7. header=F, backingfile="maged.answer1.bin",descriptorfile="maged.answer1.desc")
  8. r.mat<-r.mat[,-1];r.mat<-r.mat[-1,]
  9. colnames(r.mat)<-c('loop#','rlnorm','rllogis','rgamma')
  10. for(i in 1: mat.length){
  11.         r.mat[i,1]<-i
  12.         r.mat[i,2]<-mean(rlnorm(330,-0.0627,0.4085))
  13.         r.mat[i,3]<-mean(rllogis(330,4.3110,1.0653))
  14.         r.mat[i,4]<-mean(rgamma(330,6.080,5.949))
  15.         print(i)
  16. }
  17.  
  18. length(subset(r.mat$rlnorm,r.mat$rlnorm<=1))/length(r.mat$rlnorm)
  19. length(subset(r.mat$rllogis,r.mat$rllogis <=1))/length(r.mat$rllogis)
  20. length(subset(r.mat$rgamma,r.mat$rgamma <=1))/length(r.mat$rgamma)
create a new version of this paste RAW Paste Data