Description Usage Arguments Details Author(s) References See Also Examples

Function to visualise correlation between components from different data sets

1 | ```
plotDiablo(x, ncomp = 1, legend = TRUE, legend.ncol, ...)
``` |

`x` |
object of class inheriting from |

`ncomp` |
Which component to plot calculated from each data set. Has to be lower than the minimum of |

`legend` |
boolean. Whether the legend should be added. Default is TRUE. |

`legend.ncol` |
Number of columns for the legend. Default to |

`...` |
not used |

The function uses a plot.data.frame to plot the component `ncomp`

calculated from each data set to visualise whether DIABLO (block.splsda) is successful at maximising the correlation between each data sets' component. The lower triangular panel indicated the Pearson's correlation coefficient, the upper triangular panel the scatter plot.

Amrit Singh

Singh A., Gautier B., Shannon C., Vacher M., Rohart F., Tebbutt S. and LĂȘ Cao K.A. (2016). DIABLO - multi omics integration for biomarker discovery. Submitted.

`block.splsda`

and http://www.mixOmics.org/mixDIABLO for more details.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | ```
data('breast.TCGA')
Y = breast.TCGA$data.train$subtype
data = list(mrna = breast.TCGA$data.train$mrna,
mirna = breast.TCGA$data.train$mirna, prot = breast.TCGA$data.train$protein)
# set number of component per data set
ncomp = 3
# set number of variables to select, per component and per data set (arbitrarily set)
list.keepX = list(mrna = rep(20, 3), mirna = rep(10,3), prot = rep(10,3))
# set up a full design where every block is connected
design = matrix(1, ncol = length(data), nrow = length(data),
dimnames = list(names(data), names(data)))
diag(design) = 0
design
BC.diablo = block.splsda(X = data, Y = Y, ncomp = ncomp, keepX = list.keepX, design = design)
plotDiablo(BC.diablo, ncomp = 1)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.