6
$\begingroup$

I was wondering if there is any way to apply classification algorithms (e.g random forest) on microarray data but not using the genes as predictors/features but the pathway they belong to.

The thing is that the expressions of the genes that belong to the same pathway should somehow, grouped together and transformed into a new value that will describe that pathway for each specific sample.

Let me give an example.

Assuming that we have g1, g2, g3, g4, g5, g6, g7, g8, g9, g10 as genes and p1, p2, p3 as pathways.

And g1, g2, g3 belong to p1 while g3, g4, g5, g6belong to p2 and g2, g6, g7, g8, g9, g10 belong to p3

One thought is to calculate the mean expression of the above gene groups per sample and then run the classification algorithm on top of these values. But I would like to hear/read your opinion on such an approach and also your suggestions on using a different method to do this.

The final goal would be the extraction of the most significant classification predictors (pathways).

$\endgroup$
6
  • $\begingroup$ Do you assume that expression of genes are similar in pathways? Because in a real experiment they are not. $\endgroup$
    – benn
    Feb 27, 2018 at 10:34
  • $\begingroup$ No, I don't want to assume something like that but I mentioned the mean in order to group the genes of each pathway. This approach might be totally wrong but that's why I posted here. To get feedback and different opinions. $\endgroup$
    – J. Doe
    Feb 27, 2018 at 10:59
  • $\begingroup$ My point is, that I haven't seen an algorithm that can transform the expression values of different genes into one value that represents a given pathway. $\endgroup$
    – benn
    Feb 27, 2018 at 11:23
  • $\begingroup$ Oh ok. I see... $\endgroup$
    – J. Doe
    Feb 27, 2018 at 11:32
  • 1
    $\begingroup$ IPA is commercial. Isn't it? Can you provide me with free R-based ones to understand what exactly you suggest? $\endgroup$
    – J. Doe
    Feb 27, 2018 at 13:35

1 Answer 1

6
$\begingroup$

Have a look at the GSVA package. It allows to convert a matrix with genes x Samples to a pathways x Samples using several methods ssgsea, gage, gsva...

Afterwards you can use that matrix as input for differential expression of pathways or classification algorithms or whatever.


However it depends on the input of the "pathways" you give it. Make sure that the pathways are correct. Also consider that the pathway topology might change with conditions.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.