analyse.R 522 B

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187
  1. #!/usr/bin/Rscript
  2. args <- commandArgs(TRUE)
  3. dist <- read.table(args[1])
  4. #png("fig.png", width = 800, height = 400)
  5. pdf("fig.pdf")
  6. #x11()
  7. #options(devices="X11")
  8. vdist = dist$V1[1:100]
  9. v = ks.test(vdist,pnorm,mean(vdist),sd(vdist))
  10. print(v)
  11. hist(vdist,col = "blue",breaks = 50,freq = FALSE)
  12. den <- density(vdist)
  13. lines(den, col = "red")
  14. dev.off()