Question: Running model separately for each cell line or as the same model?
gravatar for Ahdee
4 months ago by
United States
Ahdee40 wrote:

Hi, suppose I have an experiment where I compare Drug A to control however I do this experiment with 3 cell lines. Lets say I run an RNA.seq on this and PCA analysis  shows a clear separation of cell lines but within each cell line, Drug and control separates out nicely.  The question is should I run the the analysis separately for each cell line and take the intersection OR should I run it all in the same model but correcting for cell line?  For example using limma it would look something like this. 

design <-  model.matrix(~0+ key$group + key$cell  )

the header would look something like this: 

 control  druga  cell1 cell2 

and the makeContrasts(pn=drug1-control  , levels=design)

if I do take this method how would I define my contrast to see the DGE for all three cell lines seperately? 


Or should I just do the analysis separately for each cell line and then take the intersection of the DGE from each cellline? 

so in summary its a two part question. 

1. should I run the model together or separately?

2. and if I run it together how do I define my contrast so that I find the DGE for each cell separately.  Say I want to find DGE only in cell3? 

thanks in advance! 



ADD COMMENTlink modified 4 months ago by k.vitting.seerup20 • written 4 months ago by Ahdee40
gravatar for k.vitting.seerup
4 months ago by
European Union
k.vitting.seerup20 wrote:

You kinda answered that yourself - if you are interested in the general effects of the drug (independent of cell type) you should use the combined model. 

If you are interested in the cell-type specific effects you would need to either include an interaction term in the combined model (model.matrix(~0+ key$group * key$cell  ) or do 3 separate models.

I would probably go with one combined model as you will utilize mode of the available data to estimate the variance.

ADD COMMENTlink written 4 months ago by k.vitting.seerup20
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: 110 users visited in the last hour