## Apply the arcsine-square-root transformation to the proportions

asin_table <- asin(sqrt((t(t(counts_table) / colSums(counts_table)))))

## Get significant clusters and sort them by significance

sign_clusters <- names(which(sort(da_out2$adjp[, "adjp_BCRXLvsRef"]) < FDR_cutoff))

## Get the adjusted p-values for the significant clusters

sign_adjp <- da_out2$adjp[sign_clusters , "adjp_BCRXLvsRef", drop=FALSE]

## Normalize the transformed proportions to mean = 0 and sd = 1