Question: Limma Ebayes on SILAC data
gravatar for reventropy2003
13 months ago by
reventropy200310 wrote:

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,'none')  


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




ADD COMMENTlink modified 13 months ago by Gordon Smyth35k • written 13 months ago by reventropy200310
gravatar for Gordon Smyth
13 months ago by
Gordon Smyth35k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth35k wrote:

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 COMMENTlink written 13 months ago by Gordon Smyth35k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 323 users visited in the last hour