1# Smooth plot of wcFst, pFst or abba-baba
2#
3#usage:  nohup R --vanilla < plotPfst --args smoothedpFst.txt wcFst|pFst|abba-baba
4
5cmd_args <- commandArgs(trailingOnly = TRUE)
6
7plotPfst<-function(x){
8	require("ggplot2")
9	dat<-read.table( x[1], header=FALSE )
10	dat$V2<-1:length(dat$V2)
11	pngName<-paste(c(x, format(Sys.time(), "%a%b%d_%H_%M_%S.png")), collapse="_")
12
13	theplot<-NULL
14
15	if(x[2] == "pFst"){
16		theplot<-ggplot(dat, aes(x=V2, y=-log10(V5), colour=V4))+geom_point()+theme_grey(15)+labs(x="SNP index", y="-log10(smoothed pFst)")+scale_colour_continuous(low="grey", high="red", name="variants in window")
17	}
18	if(x[2] == "wcFst"){
19		theplot<-ggplot(dat, aes(x=V2, y=V5, colour=V4))+geom_point()+theme_grey(15)+labs(x="SNP index", y="smoothed wcFst")+scale_colour_continuous(low="grey", high="red", name="variants in window")
20	}
21	if(x[2] == "xpEHH"){
22		theplot<-ggplot(dat, aes(x=V2, y=V5, colour=V4))+geom_point()+theme_grey(15)+labs(x="SNP index", y="smoothed xpEHH")+scale_colour_continuous(low="grey", high="red", name="variants in window")
23	}
24	if(x[2] == "abba-baba"){
25		theplot<-ggplot(dat, aes(x=V3, y=V5, colour=V4))+geom_point()+theme_grey(15)+labs(x="SNP index", y="D-statistic")+scale_colour_continuous(low="grey", high="red", name="variants in window")
26	}
27	ggsave(filename=pngName, width=20, height=4, units="in", theplot)
28}
29
30plotPfst(cmd_args)
31