SHARE
TWEET

Untitled

a guest Aug 22nd, 2019 73 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. library("mixOmics")
  2. gtex.plsda <- plsda(X, Y, ncomp = 2)
  3. background = background.predict(gtex.plsda, comp.predicted = 2, dist = "max.dist")
  4. plotIndiv(gtex.plsda , comp = 1:2, group = Y, ind.names = FALSE, ellipse = TRUE,
  5.           legend = TRUE, title = 'PLSDA on GTEX Skeletal Muscles', background = background)
  6. plotLoadings(gtex.plsda, comp = 1, title = 'Loadings on comp 1', contrib = 'max',
  7.              method = 'median', ndisplay = 10, size.name = 0.6)
  8. plotLoadings(gtex.plsda, comp = 2, title = 'Loadings on comp 2', contrib = 'max',
  9.              method = 'median', ndisplay = 10, size.name = 0.6)
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top