Mercurial > repos > iuc > aldex2
diff aldex2.R @ 1:75214276e2b7 draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/aldex2 commit f9acc7c5710d6eb99a883cb83b7967e6af972f1b
author | iuc |
---|---|
date | Fri, 14 Jun 2024 21:21:29 +0000 |
parents | f4d0bd4b4d6d |
children |
line wrap: on
line diff
--- a/aldex2.R Wed Jun 29 07:36:45 2022 +0000 +++ b/aldex2.R Fri Jun 14 21:21:29 2024 +0000 @@ -69,14 +69,16 @@ iterate <- get_boolean_value(opt$iterate) if (opt$analysis_type == "aldex") { - aldex_obj <- aldex(reads = reads_df, - conditions_vector, - mc.samples = opt$num_mc_samples, - test = opt$aldex_test, - effect = effect, - include.sample.summary = include_sample_summary, - denom = opt$denom, - iterate = iterate) + aldex_obj <- aldex( + reads = reads_df, + conditions_vector, + mc.samples = opt$num_mc_samples, + test = opt$aldex_test, + effect = effect, + include.sample.summary = include_sample_summary, + denom = opt$denom, + iterate = iterate + ) } else { # Generate Monte Carlo samples of the Dirichlet distribution for each sample. Convert each # instance using a log-ratio transform. This is the input for all further analyses. @@ -117,22 +119,26 @@ } else if (opt$analysis_type == "aldex_kw") { aldex_obj <- aldex.kw(aldex_clr_obj) } else if (opt$analysis_type == "aldex_plot") { - aldex_obj <- aldex(reads = reads_df, - conditions_vector, - mc.samples = opt$num_mc_samples, - test = opt$aldex_test, - effect = effect, - include.sample.summary = include_sample_summary, - denom = opt$denom, - iterate = iterate) + aldex_obj <- aldex( + reads = reads_df, + conditions_vector, + mc.samples = opt$num_mc_samples, + test = opt$aldex_test, + effect = effect, + include.sample.summary = include_sample_summary, + denom = opt$denom, + iterate = iterate + ) png(filename = opt$output) - aldex.plot(x = aldex_obj, - type = opt$plot_type, - test = opt$plot_test, - cutoff.pval = opt$cutoff_pval, - cutoff.effect = opt$cutoff_effect, - xlab = opt$xlab, - ylab = opt$ylab) + aldex.plot( + x = aldex_obj, + type = opt$plot_type, + test = opt$plot_test, + cutoff.pval = opt$cutoff_pval, + cutoff.effect = opt$cutoff_effect, + xlab = opt$xlab, + ylab = opt$ylab + ) dev.off() } else if (opt$analysis_type == "aldex_plot_feature") { png(filename = opt$output)