Hi,
Thanks for this simple and great package.
I’m now using PLS-DA with transcriptome and Proteome expression data.
In my data, each feature is divided into a total of five groups. For this, we performed plsda to find out which feature can explain the well-divided of these five groups and obtained a vip score using the vip function.
In my case, plsda was performed with comp=4, and the result of vip seems to give a score corresponding to each component. After that, I want to find out which group the feature with high vip score contributes to.
As a result of the function plotloadings(), the contents of the column ‘GroupContrib’ in the result datatable determine which group each feature contributes to. I wonder if it is possible to interpret the result by integrating the result of vip().
Or, I would like to get advice on how to select features with a VIP score and find out if it corresponds to five groups.