tanu/LSHTM_analysis
1
0
Fork 0
Code Issues Pull requests Releases Wiki Activity

added ns prefix to SNPs for unambiguity

Browse source
This commit is contained in:
Tanushree Tunstall 2020-10-13 13:37:22 +01:00
parent 59911687c8
commit 1f9ea3f789
3 changed files with 7 additions and 31 deletions
Download patch file Download diff file Expand all files Collapse all files

4
scripts/plotting/basic_barplots_LIG.R
View file

@ -94,7 +94,7 @@ OutPlot_lig_count = g + geom_bar(aes(fill = ligand_outcome)
, legend.title = element_text(size=my_als) , legend.title = element_text(size=my_als)
, plot.title = element_blank()) + , plot.title = element_blank()) +
labs(title = "" labs(title = ""
, y = "Number of SNPs" , y = "Number of nsSNPs"
#, fill="ligand_outcome" #, fill="ligand_outcome"
) + ) +
scale_fill_discrete(name = "Ligand Outcome" scale_fill_discrete(name = "Ligand Outcome"
@ -182,7 +182,7 @@ OutPlot_lig_pos_count = g + geom_bar(aes (alpha = 0.5)
, axis.title.y = element_text(size = my_als) , axis.title.y = element_text(size = my_als)
, plot.title = element_blank()) + , plot.title = element_blank()) +
labs(x = "Number of SNPs" labs(x = "Number of nsSNPs"
, y = "Number of Sites") , y = "Number of Sites")
print(OutPlot_lig_pos_count) print(OutPlot_lig_pos_count)

4
scripts/plotting/basic_barplots_PS.R
View file

@ -92,7 +92,7 @@ OutPlot_count = g + geom_bar(aes(fill = duet_outcome)
, legend.title = element_text(size=my_als) , legend.title = element_text(size=my_als)
, plot.title = element_blank()) + , plot.title = element_blank()) +
labs(title = "" labs(title = ""
, y = "Number of SNPs" , y = "Number of nsSNPs"
#, fill="DUET Outcome" #, fill="DUET Outcome"
) + ) +
scale_fill_discrete(name = "DUET Outcome" scale_fill_discrete(name = "DUET Outcome"
@ -180,7 +180,7 @@ OutPlot_pos_count = g + geom_bar(aes (alpha = 0.5)
, axis.title.y = element_text(size = my_als) , axis.title.y = element_text(size = my_als)
, plot.title = element_blank()) + , plot.title = element_blank()) +
labs(x = "Number of SNPs" labs(x = "Number of nsSNPs"
, y = "Number of Sites") , y = "Number of Sites")
print(OutPlot_pos_count) print(OutPlot_pos_count)

30
scripts/plotting/basic_barplots_foldx.R
View file

@ -71,7 +71,7 @@ theme_set(theme_grey())
# start plot 1 # start plot 1
#-------------- #--------------
g = ggplot(df, aes(x = foldx_outcome)) g = ggplot(df, aes(x = foldx_outcome))
OutPlot_count = g + geom_bar(aes(fill = foldx_outcome) foldx_outcome_count = g + geom_bar(aes(fill = foldx_outcome)
, show.legend = TRUE) + , show.legend = TRUE) +
geom_label(stat = "count" geom_label(stat = "count"
, aes(label = ..count..) , aes(label = ..count..)
@ -87,39 +87,15 @@ OutPlot_count = g + geom_bar(aes(fill = foldx_outcome)
, legend.title = element_text(size=my_als) , legend.title = element_text(size=my_als)
, plot.title = element_blank()) + , plot.title = element_blank()) +
labs(title = "" labs(title = ""
, y = "Number of SNPs" , y = "Number of nsSNPs"
#, fill="FoldX Outcome" #, fill="FoldX Outcome"
) + ) +
scale_fill_discrete(name = "FoldX Outcome" scale_fill_discrete(name = "FoldX Outcome"
, labels = c("Destabilising", "Stabilising")) , labels = c("Destabilising", "Stabilising"))
print(OutPlot_count) print(foldx_outcome_count)
dev.off() dev.off()
table(df$foldx_outcome) table(df$foldx_outcome)
#======================================================================= #=======================================================================
g = ggplot(df, aes(x = foldx_outcome))
OutPlot_count = g + geom_bar(aes(show.legend = TRUE)) +
geom_label(stat = "count"
, aes(label = ..count..)
, color = "black"
, show.legend = FALSE
, size = 10) +
theme(axis.text.x = element_blank()
, axis.title.x = element_blank()
, axis.title.y = element_text(size=my_als)
, axis.text.y = element_text(size = my_ats)
, legend.position = c(0.73,0.8)
, legend.text = element_text(size=my_als-2)
, legend.title = element_text(size=my_als)
, plot.title = element_blank()) +
labs(title = ""
, y = "Number of SNPs"
#, fill="FoldX Outcome"
) + scale_fill_manual("legend", values = c("Destabilising" = "black", "Stabilising" = "orange")) +
scale_fill_discrete(name = "FoldX Outcome"
, labels = c("Destabilising", "Stabilising"))
print(OutPlot_count)