Thank you kevin for your kind reply. I did the mean of each sample so now i have 18 samples instead of 54. i am trying to plot heat map correlation between species and chemical data.( i have 2 matrices, Species data and chemical data for each sample) What i want to do it is to show with a star the sig difference between correlation. The problem that padj value is 1 in almot all the correlations. Package ( microbiome)
y <- log10(abund_table)
x <- as.matrix(meta_table)
correlations <- associate(x, y, method = "pearson", mode = "matrix", p.adj.threshold = 1, n.signif = 1)
p <- heat(correlation.table, "X1", "X2", fill = "Correlation", star = "p.adj", p.adj.threshold = 1)
print(p)