Limma Ebayes on SILAC data
1
0
Entering edit mode
@reventropy2003-14321
Last seen 2.9 years ago

Is there any need to set up a design or contrast matrix using limma ebayes on ratiometric replicates?  I have three replicate ratios that are log2 transformed.  I'm running ebayes:

 

data = read.csv(paste(inpath,fname,sep=""),header=T,stringsAsFactors=F,na.strings="NaN",sep="\t")

fit = lmFit(data[,c("rep1","rep2","rep3")]))
fit2 = eBayes(fit)
tt = topTable(fit2, adjust="BH",number=Inf,confint=T,sort.by='none')  

 

Where "rep 1-3" are the columns containing the ratio measurements.  Is this correct?  

 

Thanks,

Jeremy

ebayes limma bioconductor • 752 views
ADD COMMENT
0
Entering edit mode

Hi Jeremy,

I meet the same problem, I also want to use limma package to deal with my SILAC data, but I really confused what codes should I use. Should I build a design matrix?

If you could reply me, so appreciated it!

Thanks so much! Best, Peiru

ADD REPLY
1
Entering edit mode
@gordon-smyth
Last seen 22 minutes ago
WEHI, Melbourne, Australia

I am not familiar with SILAC ratiometric data but, if it produces log-expression ratios, then you could analyse the results like those from two colour microarrays. Your data appears to be similar to the design considered in Section 11.2.1 of the limma User's Guide. In that case, you don't need to create a design matrix -- the default linear model with just an intercept will do the job -- so your code seems correct.

ADD COMMENT

Login before adding your answer.

Traffic: 232 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