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