Hi,
Thanks for mixOmics
. I have been finding it very useful at my end for horizontal integration. I have a query with regards to the outputs from circosPlot
.
I see in my plot there are multiple correlations between variables from protein with variables from other blocks like mRNA and metabolites. Is it possible to extract the information or network from the object that shows how one variable or analyte from protein has +ve or -ve correlation with multiple analytes of other blocks (e.g. mRNA and metabolites)? In other words, I am interested to pull those analytes, group them, and characterize what biological process or molecular functions they enrich for? I did extract the corr.mat
using below snippet
and plotted the correlation matrix but that does not let me make such inference that I can from circosPlot
. Any assistance will be appreciated.
circosPlot(MyResult.diablo.tune.P, cutoff = 0.8, line = TRUE,
color.blocks= c('darkorchid', 'brown1', 'lightgreen'),
color.cor = c("chocolate3","grey20"), size.labels = 1, size.variables = 0.6, size.legend = 0.8)
corMat.P <- circosPlot(MyResult.diablo.tune.P, cutoff = 0.8, ncol.legend = 2, size.legend = 0.8)
corMat
corrplot::corrplot(corMat.P, order="hclust",
col=brewer.pal(n=8, name="RdYlBu"),addrect = 4,tl.cex = 0.55)