transpose matrix for clr

biobakery · Nov 27, 2024 · 940b0bb · 940b0bb
1 parent b139d5f
commit 940b0bb
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/metaphlan/utils/calculate_diversity.R b/metaphlan/utils/calculate_diversity.R
@@ -224,8 +224,8 @@ if (opt$diversity == "beta"){
   # CLR or Aitchison
   if (opt$metric == "clr" || opt$metric == "aitchison"){
     # Centered Log-Ratio
-    ait_norm_mpa_table <- compositions::clr(mpa_table)
-    mat <- as.matrix(compositions::as.data.frame.rmult(ait_norm_mpa_table))
+    ait_norm_mpa_table <- compositions::clr(t(mpa_table))
+    mat <- t(as.matrix(compositions::as.data.frame.rmult(ait_norm_mpa_table)))
 
     if (opt$metric == "aitchison"){
       # Aitchison