library(ggplot2)

dirs<- list.dirs(".", recursive=F)

dirs<- dirs[grep("_K9me3", dirs)]


data_len<- data.frame(matrix(0, nrow=0, ncol=2));
colnames(data_len)<- c("tag", "len");

ref_data<- read.table("/home/Resource/Genome/mm10/mm10.chrom.sizes", header=F)
ref_size<- sum(as.numeric(as.character(ref_data$V2)))

data_percent<- data.frame(matrix(0, nrow=length(dirs), ncol=2));
colnames(data_percent)<- c("tag", "percent");


for(i in 1:length(dirs)){

    tag<- gsub("_K9me3", "", strsplit(dirs[i], "/")[[1]][2])
    data<- read.table(paste0(dirs[i], "/K9me3_LearnModel/mm10_dense.bedGraph"), header=F)

    len<- data$V3 -data$V2
    part_len<- cbind.data.frame(tag, len)
    data_len<- rbind.data.frame(data_len, part_len);
    percent<- 100*sum(len)/ref_size
    data_percent[i,]<- c(tag, percent);
   }

data_len$len<- as.numeric(as.character(data_len$len))

p<- ggplot(data_len, aes(x=log2(len), color=tag)) + geom_density();

p<- p+ ggtitle("K9me3 log2(peak length)");

pdf("peak_len_dist_K9me3.pdf");
print(p);
dev.off();
############3

data_percent$percent<- as.numeric(as.character(data_percent$percent));

q<- ggplot(data_percent, aes(x=tag, y=percent, fill=tag)) + geom_bar(stat="identity")
q<- q+ ggtitle("K9me3 peak percent of genome");

pdf("peak_len_percent_K9me3.pdf");
print(q);
dev.off();