getwd()
setwd("/users/rg/gbogu/projects/repeats/correcting_Gv19_gene_structures_using_stringtie_models/stringTieModels")

library(ggplot2)

ThreePrime_ExonLength_GV19genes <- read.table("55755_3pMostexons_gv19.plot", head=T) 
ggplot(ThreePrime_ExonLength_GV19genes, aes(x=length)) + geom_density() + theme(text = element_text(size=20))
summary(ThreePrime_ExonLength_GV19genes$length)

ThreePrime_ExonLength_StringTieTranscripts <- read.table("first_53_tissues.combined.3pMostExons_noartifacts.plots", head=T)
summary(ThreePrime_ExonLength_StringTieTranscripts$length)
ggplot(ThreePrime_ExonLength_StringTieTranscripts, aes(x=length)) + geom_density() + theme(text = element_text(size=20))

ThreePrime_ExonLength_Overlapped_Transcripts <- read.table("10909_improved_3pMostExons.plots", head=T)
summary(ThreePrime_ExonLength_Overlapped_Transcripts$length)
ggplot(ThreePrime_ExonLength_Overlapped_Transcripts, aes(x=length)) + geom_density() + theme(text = element_text(size=20))

library(reshape2)
comparisions <- read.table("comprision_in_lengths.plots", head=T)
summary(comparisions$stringTieExonLength)
m<-melt(comparisions)
ggplot(m, aes(x=value, colour=variable)) + geom_density()+ theme(text = element_text(size=20))+xlab("length")+xlim(0,12000)
ggplot(m, aes(x=variable, y=value, fill=variable)) + geom_boxplot() +guides(fill=FALSE)+ theme(text = element_text(size=20))