Search
Question: Limma Ebayes on SILAC data
0
gravatar for reventropy2003
17 days ago by
reventropy20030 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,sort.by='none')  

 

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

 

Thanks,

Jeremy

ADD COMMENTlink modified 17 days ago by Gordon Smyth32k • written 17 days ago by reventropy20030
1
gravatar for Gordon Smyth
17 days ago by
Gordon Smyth32k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth32k 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 17 days ago by Gordon Smyth32k
Please log in to add an answer.

Help
Access

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