Analysis of genome-scale count data in Bioconductor

More documents

Recommendations

Info

Consequences• Looks like dispersion is much larger for lowerabundance tags• Including this in the model would decreaseability to call low abundance tags DE (butfurther increase power for high abundancetags; is perhaps more correct)• DESeq has been designed to deal with this• edgeR will soon also include an option forallowing dispersion to vary with abundance
Concluding remarks• Must understand and account forbiological variability (overdispersion) inRNA-seq data• Negative binomial model, sharinginformation between genes• Exact and multiple testing for accurate p-values
Page 1 and 2:
Analysis of genome-scalecount data
Page 3 and 4:
Sequencing experiments used for:Seq
Page 5 and 6:
Example:RNA-seq(or similar)for gene
Page 7 and 8:
Summarization
Page 9 and 10: What does genome-scale count• e.g
Page 11 and 12: Count data• Count data (e.g. RNA-
Page 13 and 14: Poisson arises naturally frommultin
Page 15 and 16: Technical replicationDNApopulationL
Page 17 and 18: Statistical models• For count dat
Page 19 and 20: Poisson describes technicalvariance
Page 21 and 22: Overdispersion: extra-Poissonvariat
Page 23 and 24: Sources of variation: technical and
Page 25 and 26: Coefficient of variation• Dispers
Page 27 and 28: Common dispersion model• One appr
Page 29 and 30: High read number is relevant for RN
Page 31 and 32: “Composition” ofsampled DNA can
Page 33 and 34: Outline1. Applications2. Summarizat
Page 35 and 36: Extending the common dispersionmode
Page 37 and 38: Weighted Likelihood• WL is the in
Page 39 and 40: Statistical testing for countdata
Page 41 and 42: Significance testing• Simple hypo
Page 43 and 44: Further considerations• RNA-seq e
Page 45 and 46: Binomial exact testingobserved• P
Page 47 and 48: Effect of dispersion> d.tuch$counts
Page 49 and 50: Limitations of exact testing• E.g
Page 51 and 52: More complicated experiments• We
Page 53 and 54: Example: Cancer dataset• RNA-seq
Page 55 and 56: GLM results> glm.res.com[o1[1:10],]
Page 57 and 58: Mean-dispersion relationship• The
Page 59: Also seems true for more datasets
Page 63 and 64: R Practical
Page 65 and 66: Reading in data• Read the data in
Page 67 and 68: Multidimensionalscaling plot• Use
Page 69 and 70: Exact test in edgeR: common disp> d
Page 71: Exact test in edgeR: tagwise disp>
show all

Analysis of genome-scale count data in Bioconductor

Create successful ePaper yourself

Delete template?

Save as template?