1
#usage: nohup R --vanilla < plotPfst --args pFst.txt
3
cmd_args <- commandArgs(trailingOnly = TRUE)
7
dat<-read.table( x, header=FALSE )
10
theplot<-ggplot(dat, aes(x=V2, y=V9))+geom_point()+geom_segment(aes(x=V2, xend=V2, y=V10, yend=V11))+labs(x="KB position", y="Fst")+theme_grey(15)
11
pngName<-paste(c(x, format(Sys.time(), "%a%b%d_%H_%M_%S.png")), collapse="_")
12
ggsave(filename=pngName, width=20, height=4, units="in", theplot)
b'\\ No newline at end of file'