EdgeR blocking factor - missing paired sample
1
0
Entering edit mode
Lucy ▴ 30
@lucy-17014
Last seen 3 months ago
United Kingdom

Hi,

I am trying to decide on the best design for my EdgeR differential-expression analysis. I have RNA-seq data from four patients - matched blood and liver for three of the patients and only liver from one of the patients (as below).

Donor  Tissue
1      Blood
1      Liver
2      Blood
2      Liver
3      Blood
3      Liver
4      Liver

The two options I am considering for the design formula are:

  1. ~ Tissue
  2. ~ Donor + Tissue

Is it possible to use the second approach when one of the donors only has a liver sample? What is the best approach to get around this issue?

Many thanks for the advice.

Best,

Lucy

EdgeR edgeR RNA-seq differential expression analysis blocking factors • 167 views
ADD COMMENT
1
Entering edit mode
@gordon-smyth
Last seen 6 hours ago
WEHI, Melbourne, Australia

You need to use formula 2. if you wish to correct for donor baseline expression levels. It will work fine but the fourth donor will not contribute to the results.

ADD COMMENT
0
Entering edit mode

That's great thank you.

ADD REPLY

Login before adding your answer.

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