R/splot.mol.mcia.R

In omicade4: Multiple co-inertia analysis of omics datasets

splot.mol.mcia <-
function(x, axis1, axis2, col, pch, gene.nlab) {
  ndata <- length(x$coa)
  co <- x$mcoa$Tco[, c(axis1, axis2)]
  c <- as.numeric(x$mcoa$TC$"T")
  if (length(col) == 1)
    col <- rep(col, ndata)
  if (length(pch) == 1)
    pch <- rep(pch, ndata)
  
  plotgenes(co, colpoints=0, nlab=gene.nlab, pch=".", axis1=1, axis2=2, sub="variable space")
  
  for (i in 1:ndata) {
    points(co[c %in% i, ], col=col[i], pch=pch[i])
  }
}