diff --git a/scripts/plotting/barplots_subcolours.R b/scripts/plotting/barplots_subcolours.R index 4e4806a..1f98bec 100755 --- a/scripts/plotting/barplots_subcolours.R +++ b/scripts/plotting/barplots_subcolours.R @@ -124,4 +124,4 @@ print(outPlot_bp_lig) dev.off() ######################################################################= # End of script -######################################################################= \ No newline at end of file +######################################################################= diff --git a/scripts/plotting/basic_barplots_combined.R b/scripts/plotting/basic_barplots_combined.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/corr_adjusted_PS_LIG.R b/scripts/plotting/corr_adjusted_PS_LIG.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/dirs.R b/scripts/plotting/dirs.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/dist_plots_check.R b/scripts/plotting/dist_plots_check.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/extreme_muts.R b/scripts/plotting/extreme_muts.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/get_plotting_dfs.R b/scripts/plotting/get_plotting_dfs.R old mode 100644 new mode 100755 index 89b477c..5876d8d --- a/scripts/plotting/get_plotting_dfs.R +++ b/scripts/plotting/get_plotting_dfs.R @@ -39,8 +39,8 @@ import_dirs(drug, gene) #--------------------------- # call: plotting_data() #--------------------------- -if (!exists("infile_params") && exists("gene")){ -#if (!is.character(infile_params) && exists("gene")){ # when running as cmd +#if (!exists("infile_params") && exists("gene")){ +if (!is.character(infile_params) && exists("gene")){ # when running as cmd #in_filename_params = paste0(tolower(gene), "_all_params.csv") #for pncA in_filename_params = paste0(tolower(gene), "_comb_afor.csv") # part combined for gid infile_params = paste0(outdir, "/", in_filename_params) @@ -67,8 +67,8 @@ cat("\nLigand distance cut off, colname:", LigDist_colname #-------------------------------- # call: combining_dfs_plotting() #-------------------------------- -if (!exists("infile_metadata") && exists("gene")){ -#if (!is.character(infile_metadata) && exists("gene")){ # when running as cmd +#if (!exists("infile_metadata") && exists("gene")){ +if (!is.character(infile_metadata) && exists("gene")){ # when running as cmd in_filename_metadata = paste0(tolower(gene), "_metadata.csv") # part combined for gid infile_metadata = paste0(outdir, "/", in_filename_metadata) cat("\nInput file for gene metadata not specified, assuming filename: ", infile_metadata, "\n") diff --git a/scripts/plotting/get_plotting_dfs_with_lig.R b/scripts/plotting/get_plotting_dfs_with_lig.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/ggcorr_all_PS_LIG.R b/scripts/plotting/ggcorr_all_PS_LIG.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/hist_af_or_base.R b/scripts/plotting/hist_af_or_base.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/hist_af_or_combined.R b/scripts/plotting/hist_af_or_combined.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/legend_adjustment.R b/scripts/plotting/legend_adjustment.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/lineage_basic_barplots_combined.R b/scripts/plotting/lineage_basic_barplots_combined.R old mode 100644 new mode 100755 index 4b63587..94bcd4a --- a/scripts/plotting/lineage_basic_barplots_combined.R +++ b/scripts/plotting/lineage_basic_barplots_combined.R @@ -3,9 +3,10 @@ getwd() setwd("~/git/LSHTM_analysis/scripts/plotting/") getwd() ######################################################### -# TASK: Basic lineage barplot showing numbers +# TASK: Basic lineage barplots showing numbers # Output: Basic barplot with lineage samples and mut count +# + SNP diversity ########################################################## # Installing and loading required packages @@ -16,6 +17,27 @@ source("../functions/bp_lineage.R") #=========== # input #=========== +#drug = 'streptomycin' +#gene = 'gid' + +spec = matrix(c( + "drug" , "d", 1, "character", + "gene" , "g", 1, "character", + "data_file1" , "fa", 2, "character", + "data_file2" , "fb", 2, "character" +), byrow = TRUE, ncol = 4) + +opt = getopt(spec) + +drug = opt$drug +gene = opt$gene +infile_params = opt$data_file1 +infile_metadata = opt$data_file2 + +if(is.null(drug)|is.null(gene)) { + stop("Missing arguments: --drug and --gene must both be specified (case-sensitive)") +} + source ('get_plotting_dfs.R') cat("Directories imported:" diff --git a/scripts/plotting/lineage_bp_data.R b/scripts/plotting/lineage_bp_data.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/lineage_dist_combined_PS.R b/scripts/plotting/lineage_dist_combined_PS.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/lineage_dist_dm_om_combined_PS.R b/scripts/plotting/lineage_dist_dm_om_combined_PS.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/opp_mcsm_muts.R b/scripts/plotting/opp_mcsm_muts.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/or_plots_combined.R b/scripts/plotting/or_plots_combined.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/other_plots_combined.R b/scripts/plotting/other_plots_combined.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/other_plots_data.R b/scripts/plotting/other_plots_data.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/output_tables.R b/scripts/plotting/output_tables.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/ps_plots_combined.R b/scripts/plotting/ps_plots_combined.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/resolving_ambiguous_muts.R b/scripts/plotting/resolving_ambiguous_muts.R old mode 100644 new mode 100755 diff --git a/scripts/plotting/running_plotting_scripts.txt b/scripts/plotting/running_plotting_scripts.txt index 332293e..74f67de 100644 --- a/scripts/plotting/running_plotting_scripts.txt +++ b/scripts/plotting/running_plotting_scripts.txt @@ -113,12 +113,28 @@ note: - fb flag has default if not supplied -#=================== -# Add LINEAGE ONE -#=================== -# Lineage_bp.R -creates Count and Diversity plot +#==================================== +# lineage_basic_barplots_combined.R +#==================================== +#----------------------------------------------------------------------- +./lineage_basic_barplots_combined.R-d streptomycin -g gid +#----------------------------------------------------------------------- +It replaces + ## lineage_basic_barplot.R +These have been moved to redundant/ + + +sources: + ## get_plotting_dfs.R + ## functions//bp_lineage.R" + + outputs: 1 svg in the plotdir + ## basic_lineage_barplots_combined.svg + +note: + - fa flag has default if not supplied + - fb flag has default if not supplied ######################################################################## # TODO