Issue #28 resolved
Kim-Anh Le Cao repo owner created an issue

plotindiv with sgccda, pch option (from workshop material)

data(nutrimouse) X1 = nutrimouse$gene[1:50] X2 = nutrimouse$gene[51:120] X3 = nutrimouse$lipid

Y = nutrimouse$diet data = list(Gene1 = X1, Gene2 = X2, Lipid = X3)

design = matrix(c(0,1,1,1,0,1,1,1,0), ncol = 3, nrow = 3, byrow = TRUE, dimnames=c(list(names(data), names(data))))

nutrimouse.sgccda <- wrapper.sgccda(blocks = data, Y = Y, design = design, keep = list(X1 = c(20,20), X2 = c(20, 20), X3 = c(10,10)), ncomp = c(2, 2, 2), scheme = "centroid" )

this is not outputting the pch

plotIndiv(nutrimouse.sgccda, block = c(1,2,3), pch = 16, group = nutrimouse$diet, plot.ellipse = TRUE, add.legend =TRUE)

Comments (18)

  1. Log in to comment