#!/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()