Differential expression at pathway level
1
0
Entering edit mode
8 months ago
Chris ▴ 340

Hi Biostar,

I would like to do differential expression at pathway level by following this tutorial. I have microarray data. Would you please tell me how to get an object like leukimia_eset? Thank you so much!

data(leukemia)
leukemia_eset

https://bioconductor.org/packages/devel/bioc/vignettes/GSVA/inst/doc/GSVA.html

GSVA • 635 views
ADD COMMENT
2
Entering edit mode
8 months ago
AquaDeath ▴ 80

Hi!

The dataset on the tutorial is an ExpressionSet object. There is a detailed manual for building this type of Bioconductor object here, but esentially you need the assay info (intensities measured in microarray experiment), the probes info and the samples info as basic elements of ExpressionSet.

ADD COMMENT
1
Entering edit mode

Thank you so much for the info!

ADD REPLY
0
Entering edit mode

Hi,

microarray <- data.matrix(microarray, rownames.force = NA)
ExpressionSet_object <- ExpressionSet(assayData=microarray)

I can create an expressionset object but I don't have subtype to run the next step:

mod <- model.matrix(~ factor(leukemia_es$subtype))

My microarray data can classify samples into 2 or 3 sub groups. Do you have any suggestion to add subtype into my expressionset object:

ExpressionSet (storageMode: lockedEnvironment)
assayData: 17873 features, 5626 samples 
  element names: exprs 
protocolData: none
phenoData: none
featureData: none
experimentData: use 'experimentData(object)'
Annotation:  
ADD REPLY
2
Entering edit mode

The ExpressionSet object you are building still lacks two key "slots": featureData and phenoData. In this context, features refer to information regarding microarray probes. You should have some type of file with probe IDs (depending of the array used) that can be mapped later to genes with packages such as AnnotationDbi.

Then, phenotypic data depends on the experimental design of your run. For what you're telling me you can classify your samples in several groups depending on a factor or biological condition you want to study. You should build a data.frame whose rows correspond to the columns of the microarray matrix with the samples info (I'm assuming your microarray matrix is in the standard -omics form: rows are different probes and columns are intensities measured for each sample), and add that as pData. With this, you'd be able to use that factor the same way as the tutorial uses subtype to build a model matrix for applying limma. This last one is crucial for differential expression analysis in microarrays, so definitely worth to study it carefully before performing any enrichment analysis.

There are good limma examples in the vignettes of the package, but I'll also leave you some practical example pages:

Good luck with your analysis!

ADD REPLY

Login before adding your answer.

Traffic: 2087 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6