IN="TEST2.bdg.valpeak.tsv"
FIN=read.table(IN,header=FALSE)
print(head(FIN))
DV<-data.frame(FIN[which(FIN$V3 == -1),])
DP<-data.frame(FIN[which(FIN$V3 == 1),])
png(paste0(IN,".valley.png"),600,600)
boxplot(DV$V2)
dev.off()
png(paste0(IN,".peak.png"),600,600)
boxplot(DP$V2)
dev.off()
~
相关文章：