Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Surrogate Variable Analysis for Complex Experimental Design
1
gravatar for tucanj
3.7 years ago by
tucanj • 70
Canada
tucanj • 70 wrote:

I am conducting differential expression analysis in limma for a time course microarray experiment. I have 3 time points (1,2,3) and two conditions (control,disease). My contrast is formulated as:

pheno<-factor(c("disease.1","disease.1","disease.2","disease.2","disease.3","disease.3","control.1","control.1","control.2","control.2","control.3","control.3")
​model<-model.matrix(~0+pheno)
contrast<-makeContrasts(disease.1-control.1, disease.2-control.2, disease.3-control.3,levels=model)
fit<-lmFit(expression,model)
fit2<-contrasts.fit(fit,contrast)
fit2<-eBayes(fit2)

topTableF(fit2)

How can I apply SVA with this experimental setup? Or is there an alternative? The challenge I am facing is SVA seems designed for two class experiments, because the null model requires leaving out the variable of interest, however in this case there isn't one variable but many.

Thanks!
limma sva R • 1.8k views
ADD COMMENT • link •
modified 3.7 years ago by andrew.j.skelton73 ♦ 5.3k • written 3.7 years ago by tucanj • 70
1
gravatar for andrew.j.skelton73
3.7 years ago by
andrew.j.skelton73 ♦ 5.3k
London
andrew.j.skelton73 ♦ 5.3k wrote:

You apply the sva function to a normalised matrix of expression values. I'm assuming that you're just looking for potential surrogate variables and not trying to correct for a known batch effect?

Edit: Also, I'd recommend you check out the time course section of the limma users guide http://www.bioconductor.org/packages/release/bioc/vignettes/limma/inst/doc/usersguide.pdf
ADD COMMENT • link modified 3.7 years ago • written 3.7 years ago by andrew.j.skelton73 ♦ 5.3k

Thanks for the response! Yes not adjusting for any known batch effects. I didn't realize that SVA did not require a null model, but I just tested it an it works (another thing I had to change was to include an intercept in my model because sva did not work without). As well, I purposefully chose this model but thanks for referring me to limma user guide.

ADD REPLY • link written 3.7 years ago by tucanj • 70
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 775 users visited in the last hour