README.mkdn 1.7 KB

123456789101112131415161718192021222324252627282930313233343536
  1. # Human globin protocol stats
  2. This is a pair of spreadsheets summarizing the globin-reducing
  3. properties of an experimental RNA-seq protocol.
  4. - [method-select.xlsx](method-select.xlsx) shows the results for
  5. several different methods of globin reduction. The most important
  6. column is "Non-HB", which represents the fraction of total reads
  7. that map to non-globin genes.
  8. - [concentration-select.xlsx](concentration-select.xlsx) shows the
  9. results of selecting the best globin blocking method and optimizing
  10. the concentration of blockers as well as the number of
  11. hybridizations.
  12. # Cyno globin plots
  13. This is a series of example plots for evaluation of a similar globin
  14. reduction protocol that was designed for cynomolgus monkeys.
  15. - [BCVplots.pdf](BCVplots.pdf) and [corrplot.pdf](corrplot.pdf) show,
  16. respectively, that the biological coefficient of variation is not
  17. increased, and the sample-to-sample correlation is not decreased, by
  18. the globin reduction protocol.
  19. - [cyno-vs-hg19.pdf](cyno-vs-hg19.pdf) shows excellent correlation
  20. between total read counts on the cyno and human genomes, indicating
  21. that the cyno annotation is reasonably complete.
  22. - [maplot.pdf](maplot.pdf) shows an MA plot of the effects of globin
  23. reduction. Notably, there are real and significant effects on
  24. non-globin genes, but they are consistent and predictable.
  25. - [pval-comparisons.pdf](pval-comparisons.pdf) shows the comparison
  26. between p-values from edgeR, limma-voom, and DESeq2 on the exact
  27. same differential expression test. Surprisingly, despite the
  28. significant algorithmic differences, edgeR and limma are in quite
  29. close agreement. DESeq2 id overly liberal because it does not
  30. account for the negative estimation bias in negative-binomial
  31. dispersions.