Mercurial > repos > iuc > aldex2

diff aldex2.R @ 1:75214276e2b7 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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)