How to design limma for two populations from each tumor, multiple tumors from several of the patients
0
0
Entering edit mode
JMINCHAN • 0
@390ece06
Last seen 4 months ago
United States

I have an experiment where I perform pseudobulk differential expression using limma-voom. In particular, I obtain expression profiles of cell types A and B from two subtypes of tumor samples X and Y. However, some of the tumor samples come from the same patient. I want to compare the difference between X and Y while adjusting for cell type and controlling for the fact that there are pairs of cell types coming from each sample and some samples come from the same patient. I am wondering what the best approach would be? So far, I use the model

~ cell_type + subtype

and use duplicateCorrelation with block = patient (treating patient as a random effect). However, I feel like this ignores the hierarchical nature of the paired cell types per sample, and some samples matched to the same patient. Is there any way to use duplicateCorrelation to capture this hierarchical design?

limma • 187 views
ADD COMMENT
0
Entering edit mode

What differential expression comparison are you planning to make? A vs B for each subtype? Or X vs Y for each cell type?

ADD REPLY

Login before adding your answer.

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