text: user4595: I have RNA-seq from two sequencing batches; Lab technician says that he has run the RNA expression quantification two times in bathes 1 and 2 for example tumor 1 in batch 1 and tumor 1 in batch 2 , normal 2 in batch1 and normal 2 in batch 2. PCA shows however shows two experimental runs are too close to each other and can be merged
  [![enter image description here][1]][1]
  Now my question is can I compare only a tumour sample with its own matched normal sample in one patient to get differentially expressed genes for this patient by using experimental runs (technical replicates )? Design would be Tumour1_batch1 Tumour1_batch2 versus Normal1_batch1 Normal1_bacth2 
  > head(mycols)
                Condition
  Normal1_batch1  N
  Normal1_batch2  N     
  Tumour1_batch1  T     
  Tumour1_batch2  T     
  
  
  user235: No. If the only difference between Tumour_batch1 and Tumour_batch2 was the library prep/sequencing run, then all this will tell you is if the particular sample of RNA taken from the tumour is different from the particular sample of RNA taken from normal. It almost certainly is.
  Indeed if you did technical replicates of two RNA samples, both taken from normal, you would almost certainly find differences. 
  What you want to find is the differences in RNA level between tumour in general and normal in general. To do that you need to make an estimate of the mean and variance of RNA levels in the tumour and in the normal. To do that you need multiple samples from each. 
  If you need differences within a single patient, it would be more honest just to compute the fold changes without any implication of statistical certainty. You would then need to think about alternate schemes for measuring the confidence in inferences drawn from the data when you consider the actual biological questions.
  
original_text: 4595: <p>I have RNA-seq from two sequencing batches; Lab technician says that he has run the RNA expression quantification two times in bathes 1 and 2 for example tumor 1 in batch 1 and tumor 1 in batch 2 , normal 2 in batch1 and normal 2 in batch 2. PCA shows however shows two experimental runs are too close to each other and can be merged</p>
  
  <p>[![enter image description here][1]][1]</p>
  
  <p>Now my question is can I compare only <code>a tumour sample with its own matched normal sample in one patient</code> to get differentially expressed genes for this <code>patient</code> by using experimental runs (technical replicates )? Design would be <code>Tumour1_batch1</code> <code>Tumour1_batch2</code> <strong><em>versus</em></strong> <code>Normal1_batch1</code> <code>Normal1_bacth2</code> </p>
  
  <pre><code>&gt; head(mycols)
                Condition
  Normal1_batch1  N
  Normal1_batch2  N     
  Tumour1_batch1  T     
  Tumour1_batch2  T     
  </code></pre>
   235: <p>No. If the only difference between Tumour_batch1 and Tumour_batch2 was the library prep/sequencing run, then all this will tell you is if the particular sample of RNA taken from the tumour is different from the particular sample of RNA taken from normal. It almost certainly is.</p>
  
  <p>Indeed if you did technical replicates of two RNA samples, both taken from normal, you would almost certainly find differences. </p>
  
  <p>What you want to find is the differences in RNA level between tumour in general and normal in general. To do that you need to make an estimate of the mean and variance of RNA levels in the tumour and in the normal. To do that you need multiple samples from each. </p>
  
  <p>If you need differences within a single patient, it would be more honest just to compute the fold changes without any implication of statistical certainty. You would then need to think about alternate schemes for measuring the confidence in inferences drawn from the data when you consider the actual biological questions.</p>
  