# HG changeset patch # User fubar # Date 1420167187 18000 # Node ID a9ffdf5f7d2a1e7a87d71caeece869b4f2999926 # Parent a95f5496d3306032ec995b5b21ca62eb40dd4cba with tests diff -r a95f5496d330 -r a9ffdf5f7d2a test-data/edgeRtest1out.html --- a/test-data/edgeRtest1out.html Fri Nov 28 06:45:33 2014 -0500 +++ b/test-data/edgeRtest1out.html Thu Jan 01 21:53:07 2015 -0500 @@ -8,26 +8,26 @@
-KernSmooth 2.23 loaded +[1] @@@ Quantiles for non-zero row counts: + + 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% + + 1.0 1.0 2.0 3.0 4.0 8.0 13.0 24.0 86.6 753.0 2325567.0 + +[1] Read 3242 contigs. Removed 1494 with no reads. After filtering at count quantile =0.3, there are 1141 contigs + +[1] "@@ using genecards substitution for urls" + +[1] @@ Total low count contigs per sample = 1,1,1,1,1,1,1,1 + +[1] Using samples: liver_X11706Liv_CAAAAG_L003_R1_001_trimmed.fastq_bwa.sam.bam,liver_X11700Liv_ATTCCT_L003_R1_001_trimmed.fastq_bwa.sam.bam,liver_X11698Liv_ACTGAT_L003_R1_001_trimmed.fastq_bwa.sam.bam,liver_X11699Liv_ATGAGC_L003_R1_001_trimmed.fastq_bwa.sam.bam,heart_X11706He_AGTTCC_L001_R1_001_trimmed.fastq_bwa.sam.bam,heart_X11699He_GGCTAC_L001_R1_001_trimmed.fastq_bwa.sam.bam,heart_X11698He_TAGCTT_L001_R1_001_trimmed.fastq_bwa.sam.bam,heart_X11700He_CTTGTA_L001_R1_001_trimmed.fastq_bwa.sam.bam + +[1] Using design matrix: + + (Intercept) groupliver + +1 1 1 + +2 1 1 + +3 1 1 + +4 1 1 -Copyright M. P. Wand 1997-2009 +5 1 0 + +6 1 0 + +7 1 0 + +8 1 0 + +attr(,"assign") + +[1] 0 1 + +attr(,"contrasts") + +attr(,"contrasts")$group + +[1] contr.treatment + +[1] "Raw sample read totals 2443751,1644652,1682104,1806045,1440960,1341813,2888924,1428365" + +[1] heart liver + + ++ +
Attaching package: ‘gplots’ @@ -157,9 +221,9 @@ lowess -Loading required package: limma +Loading required package: methods -Loading required package: methods +Loading required package: limma Loading required package: splines @@ -187,7 +251,7 @@ The following objects are masked from ‘package:base’: - anyDuplicated, append, as.data.frame, as.vector, cbind, colnames, do.call, duplicated, eval, evalq, Filter, Find, get, intersect, is.unsorted, lapply, Map, mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rep.int, rownames, sapply, setdiff, sort, table, tapply, union, unique, unlist + anyDuplicated, append, as.data.frame, as.vector, cbind, colnames, duplicated, eval, evalq, Filter, Find, get, intersect, is.unsorted, lapply, Map, mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rep.int, rownames, sapply, setdiff, sort, table, tapply, union, unique, unlist Loading required package: IRanges @@ -197,7 +261,7 @@ space -Loading required package: GenomeInfoDb +Loading required package: XVector Loading required package: Rcpp @@ -209,21 +273,21 @@ final dispersion estimates -Warning messages: +you had estimated gene-wise dispersions, removing these -1: In bplt(at[i], wid = width[i], stats = z$stats[, i], out = z$out[z$group == : +you had estimated fitted dispersions, removing these - Outlier (-Inf) in boxplot 1 is not drawn +you had estimated gene-wise dispersions, removing these -2: In bplt(at[i], wid = width[i], stats = z$stats[, i], out = z$out[z$group == : +you had estimated fitted dispersions, removing these - Outlier (-Inf) in boxplot 3 is not drawn +Warning message: -3: In bxp(list(stats = c(-0.430723026286372, -0.127900608036896, 0.474159383291067, : +closing unused connection 4 (edgeR.log) - some notches went outside hinges ('box'): maybe set notch=FALSE +Warning message: -4: In par(defpar) : calling par(new=TRUE) with no plot +In sink() : no sink to remove@@ -232,11 +296,11 @@ (Click on a thumbnail image to download the corresponding original PDF image)
-[1] Common Dispersion = 0.228651459675839 CV = 0.478175134940995 getPriorN = 3.33333333333333 +[1] Common Dispersion = 0.228651460998105 CV = 0.478175136323613 getPriorN = 3.33333333333333 [1] "No GLM fit outlier genes found\n" -[1] Heart Liver - -[1] "Raw sample read totals 2443751,1644652,1682104,1806045,1440960,1341813,2888924,1428365" - -[1] raw contig counts by sample: - - Liver_X11706Liv_CAAA Liver_X11700Liv_ATTC Liver_X11698Liv_ACTG Liver_X11699Liv_ATGA Heart_X11706He_AGTTC Heart_X11699He_GGCTA Heart_X11698He_TAGCT Heart_X11700He_CTTGT - - Min. :0.0043 Min. :0.0043 Min. :0.0043 Min. :0.0043 Min. :0.0043 Min. :0.0043 Min. :0.0043 Min. :0.0043 - - 1st Qu.:0.3032 1st Qu.:0.0043 1st Qu.:0.3032 1st Qu.:0.0043 1st Qu.:0.0043 1st Qu.:0.0043 1st Qu.:0.0043 1st Qu.:0.0043 - - Median :0.7789 Median :0.6031 Median :0.6998 Median :0.6031 Median :0.4786 Median :0.4786 Median :0.4786 Median :0.4786 - - Mean :1.0519 Mean :1.0343 Mean :0.9855 Mean :0.9966 Mean :0.9210 Mean :0.9428 Mean :1.0205 Mean :0.9753 - - 3rd Qu.:1.5410 3rd Qu.:1.6335 3rd Qu.:1.4473 3rd Qu.:1.5534 3rd Qu.:1.5770 3rd Qu.:1.5855 3rd Qu.:1.7161 3rd Qu.:1.6022 - - Max. :5.8209 Max. :5.6905 Max. :5.7999 Max. :5.7215 Max. :5.3609 Max. :5.3589 Max. :5.6967 Max. :5.3702 - - NA's :650 NA's :969 NA's :664 NA's :886 NA's :902 NA's :957 NA's :821 NA's :950 - -[1] normalised contig counts by sample: - - Liver_X11706Liv_CAAA Liver_X11700Liv_ATTC Liver_X11698Liv_ACTG Liver_X11699Liv_ATGA Heart_X11706He_AGTTC Heart_X11699He_GGCTA Heart_X11698He_TAGCT Heart_X11700He_CTTGT - - Min. :-Inf Min. :-Inf Min. :-Inf Min. :-Inf Min. :-Inf Min. :-Inf Min. :-Inf Min. :-Inf - - 1st Qu.: 0 1st Qu.:-Inf 1st Qu.: 0 1st Qu.:-Inf 1st Qu.:-Inf 1st Qu.:-Inf 1st Qu.:-Inf 1st Qu.:-Inf - - Median : 0 Median : 0 Median : 0 Median : 0 Median : 0 Median : 0 Median : 0 Median : 0 - - Mean :-Inf Mean :-Inf Mean :-Inf Mean :-Inf Mean :-Inf Mean :-Inf Mean :-Inf Mean :-Inf - - 3rd Qu.: 1 3rd Qu.: 1 3rd Qu.: 1 3rd Qu.: 1 3rd Qu.: 1 3rd Qu.: 1 3rd Qu.: 1 3rd Qu.: 1 - - Max. : 6 Max. : 6 Max. : 6 Max. : 5 Max. : 5 Max. : 5 Max. : 6 Max. : 5 - -[1] Using ncol rawrs= 8 - -[1] # edgeR Top tags\n +[1] @@ edgeR Top tags\n Name logFC logCPM LR PValue adj.p.value Dispersion totreads URL @@ -438,19 +458,19 @@ Mir149 Mir149 -7.008984 8.861767 484.30321 2.473909e-107 1.411365e-104 0.04959937 6164 Mir149 -Mir208b Mir208b -13.291635 9.905945 417.69758 7.737466e-93 2.942816e-90 0.10508096 14756 Mir208b +Mir208b Mir208b -13.291635 9.905945 417.69758 7.737463e-93 2.942815e-90 0.10508096 14756 Mir208b -Mir122a Mir122a 10.514683 12.478088 415.17429 2.740528e-92 7.817355e-90 0.10803882 90428 Mir122a +Mir122a Mir122a 10.514683 12.478088 415.17429 2.740525e-92 7.817349e-90 0.10803882 90428 Mir122a -Mir204 Mir204 -7.498162 7.634507 341.30678 3.313429e-76 7.561245e-74 0.06907958 2601 Mir204 +Mir204 Mir204 -7.498162 7.634507 341.30678 3.313430e-76 7.561247e-74 0.06907958 2601 Mir204 -Mir499 Mir499 -13.577454 8.700078 325.79199 7.930762e-73 1.508166e-70 0.12042284 6527 Mir499 +Mir499 Mir499 -13.577454 8.700078 325.79199 7.930755e-73 1.508165e-70 0.12042284 6527 Mir499 -Mir490 Mir490 -8.534394 6.991023 303.17184 6.710365e-68 1.093790e-65 0.07949711 1741 Mir490 +Mir490 Mir490 -8.534394 6.991023 303.17184 6.710366e-68 1.093790e-65 0.07949711 1741 Mir490 -Mir192 Mir192 6.953853 17.169364 217.22866 3.638312e-49 5.189142e-47 0.12700995 2325567 Mir192 +Mir192 Mir192 6.953853 17.169364 217.22867 3.638307e-49 5.189135e-47 0.12700995 2325567 Mir192 -Mir802 Mir802 11.440805 6.593380 212.88058 3.231647e-48 4.097010e-46 0.12273671 1514 Mir802 +Mir802 Mir802 11.440805 6.593380 212.88059 3.231644e-48 4.097007e-46 0.12273671 1514 Mir802 Mir1948 Mir1948 7.418142 7.252734 195.66958 1.840248e-44 2.099723e-42 0.12060221 2404 Mir1948 @@ -458,21 +478,21 @@ Mir23a Mir23a -3.153807 9.529402 177.53185 1.676248e-40 1.593833e-38 0.04442763 10118 Mir23a -Mir181c Mir181c -3.767686 10.639598 169.87390 7.883296e-39 6.919108e-37 0.06368883 23605 Mir181c +Mir181c Mir181c -3.767686 10.639598 169.87390 7.883295e-39 6.919107e-37 0.06368883 23605 Mir181c -Mir3073 Mir3073 10.686337 5.859950 164.86740 9.778601e-38 7.969560e-36 0.14069249 904 Mir3073 +Mir3073 Mir3073 10.686337 5.859950 164.86740 9.778593e-38 7.969554e-36 0.14069249 904 Mir3073 Mir181d Mir181d -3.643963 7.300371 162.18591 3.767663e-37 2.865936e-35 0.05729574 2139 Mir181d Mir195 Mir195 -3.203683 8.215089 150.20548 1.563314e-34 1.114838e-32 0.05235020 3962 Mir195 -Mir10b Mir10b -5.182616 13.946466 147.24792 6.926822e-34 4.649120e-32 0.12268790 197340 Mir10b +Mir10b Mir10b -5.182616 13.946466 147.24793 6.926819e-34 4.649118e-32 0.12268790 197340 Mir10b -Mir101b Mir101b 3.759962 11.863187 136.31359 1.703813e-31 1.080028e-29 0.07961343 59019 Mir101b +Mir101b Mir101b 3.759962 11.863187 136.31359 1.703812e-31 1.080028e-29 0.07961343 59019 Mir101b Mir378 Mir378 -3.115599 8.119617 126.76408 2.092233e-29 1.256441e-27 0.05942391 4075 Mir378 -Mir27a Mir27a -3.064687 10.642480 124.98911 5.117477e-29 2.919521e-27 0.06113852 21886 Mir27a +Mir27a Mir27a -3.064687 10.642480 124.98911 5.117477e-29 2.919520e-27 0.06113852 21886 Mir27a Mir182 Mir182 5.057509 8.846381 123.17765 1.275060e-28 6.927826e-27 0.13653707 7189 Mir182 @@ -480,7 +500,7 @@ Mir199b Mir199b -5.520119 4.792610 102.10724 5.259607e-24 2.609223e-22 0.13417024 370 Mir199b -Mir181a-2 Mir181a-2 -3.000177 7.637692 101.38361 7.578821e-24 3.603098e-22 0.06896653 2817 Mir181a-2 +Mir181a-2 Mir181a-2 -3.000177 7.637692 101.38361 7.578821e-24 3.603098e-22 0.06896654 2817 Mir181a-2 Mir125b-2 Mir125b-2 -2.987759 8.144514 91.72544 9.957640e-22 4.488356e-20 0.07737381 3837 Mir125b-2 @@ -490,43 +510,43 @@ Mir215 Mir215 -3.058208 6.447966 84.35278 4.139167e-20 1.686711e-18 0.08138517 1182 Mir215 -Mir133b Mir133b -8.383611 3.584760 83.96681 5.031521e-20 1.960318e-18 0.17482280 159 Mir133b +Mir133b Mir133b -8.383611 3.584760 83.96681 5.031517e-20 1.960318e-18 0.17482280 159 Mir133b -Mir150 Mir150 -2.883446 8.307765 83.91918 5.154209e-20 1.960318e-18 0.08008123 4229 Mir150 +Mir150 Mir150 -2.883446 8.307765 83.91918 5.154210e-20 1.960318e-18 0.08008123 4229 Mir150 Mir3074-2 Mir3074-2 -2.778308 7.935651 83.74839 5.619282e-20 2.040616e-18 0.07424646 3470 Mir3074-2 -Mir24-2 Mir24-2 -2.778307 7.935651 83.71222 5.723023e-20 2.040616e-18 0.07427992 3470 Mir24-2 +Mir24-2 Mir24-2 -2.778307 7.935651 83.71222 5.723024e-20 2.040616e-18 0.07427992 3470 Mir24-2 Mir193 Mir193 5.176579 4.801090 83.19222 7.445011e-20 2.574169e-18 0.14794861 421 Mir193 Scarna17 Scarna17 2.182159 9.244479 81.91330 1.421894e-19 4.771710e-18 0.04982909 9224 Scarna17 -Mir214 Mir214 -3.271172 6.271755 80.43948 2.997457e-19 9.771711e-18 0.09566584 1048 Mir214 +Mir214 Mir214 -3.271172 6.271755 80.43948 2.997458e-19 9.771712e-18 0.09566584 1048 Mir214 -Snord104 Snord104 2.330488 11.053611 79.50529 4.809370e-19 1.524303e-17 0.05915990 33458 Snord104 +Snord104 Snord104 2.330488 11.053611 79.50529 4.809369e-19 1.524303e-17 0.05915990 33458 Snord104 Mir200a Mir200a 7.201555 4.139422 77.35503 1.428304e-18 4.365755e-17 0.19287764 264 Mir200a Mir200b Mir200b 6.525423 5.752604 77.31985 1.453976e-18 4.365755e-17 0.26237966 888 Mir200b -Mir21 Mir21 2.923147 13.825255 75.51798 3.620939e-18 1.059357e-16 0.09395834 229120 Mir21 +Mir21 Mir21 2.923147 13.825255 75.51798 3.620938e-18 1.059357e-16 0.09395834 229120 Mir21 Mir203 Mir203 1.956427 8.767610 75.17870 4.299815e-18 1.226522e-16 0.04381710 6739 Mir203 Mir155 Mir155 -3.886731 5.068563 73.81316 8.587210e-18 2.389758e-16 0.12522673 463 Mir155 -Cyp3a25 Cyp3a25 8.681501 3.972085 72.29680 1.851472e-17 5.029834e-16 0.23125383 226 Cyp3a25 +Cyp3a25 Cyp3a25 8.681501 3.972085 72.29680 1.851471e-17 5.029829e-16 0.23125383 226 Cyp3a25 -Rabggtb Rabggtb 1.934093 10.298211 72.02043 2.129809e-17 5.651423e-16 0.04596646 19535 Rabggtb +Rabggtb Rabggtb 1.934093 10.298211 72.02043 2.129809e-17 5.651422e-16 0.04596646 19535 Rabggtb -Mir23b Mir23b -2.100584 10.184110 71.44225 2.854936e-17 7.403367e-16 0.05416378 16387 Mir23b +Mir23b Mir23b -2.100584 10.184110 71.44225 2.854935e-17 7.403367e-16 0.05416378 16387 Mir23b -Snord52 Snord52 2.207491 10.217554 71.27974 3.100028e-17 7.860292e-16 0.05941483 18059 Snord52 +Snord52 Snord52 2.207491 10.217554 71.27974 3.100027e-17 7.860292e-16 0.05941483 18059 Snord52 -Gm5441 Gm5441 -6.881248 3.538457 70.05615 5.764005e-17 1.429724e-15 0.20097284 142 Gm5441 +Gm5441 Gm5441 -6.881248 3.538457 70.05615 5.764004e-17 1.429724e-15 0.20097284 142 Gm5441 -Mir153 Mir153 -6.857671 3.517446 69.37600 8.137283e-17 1.975455e-15 0.20158808 140 Mir153 +Mir153 Mir153 -6.857671 3.517446 69.37600 8.137282e-17 1.975455e-15 0.20158808 140 Mir153 Mir132 Mir132 -2.858294 5.938312 64.52507 9.531204e-16 2.265647e-14 0.09274248 857 Mir132 @@ -534,28 +554,11 @@ Snord91a Snord91a 2.654072 6.557504 62.40549 2.795431e-15 6.379174e-14 0.08637410 1437 Snord91a -[1] # 416 tags significant at adj p= 0.05 +[1] @@ 416 tags significant at adj p= 0.05-
@@ -564,59 +567,9 @@ Got TCols=1 5 6 7; CCols=2 3 4 8 -[1] # Quantiles for non-zero row counts: - - 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% - - 1.0 1.0 2.0 3.0 4.0 8.0 13.0 24.0 86.6 753.0 2325567.0 - -[1] Read 3242 contigs. Removed 1494 with no reads. After filtering at count quantile =0.3, there are 1141 contigs - -[1] "@@ using genecards substitution for urls" - -[1] # urls sample - -[1] 0610005C13Rik 0610007N19Rik 0610008F07Rik 0610009L18Rik 0610012G03Rik - -[6] 0610031O16Rik - -[1] # Total low count contigs per sample = 1 # Total low count contigs per sample = 1 # Total low count contigs per sample = 1 # Total low count contigs per sample = 1 # Total low count contigs per sample = 1 # Total low count contigs per sample = 1 # Total low count contigs per sample = 1 # Total low count contigs per sample = 1 - -[1] Using samples: Liver_X11706Liv_CAAAAG_L003_R1_001_trimmed.fastq_bwa.sam.bam,Liver_X11700Liv_ATTCCT_L003_R1_001_trimmed.fastq_bwa.sam.bam,Liver_X11698Liv_ACTGAT_L003_R1_001_trimmed.fastq_bwa.sam.bam,Liver_X11699Liv_ATGAGC_L003_R1_001_trimmed.fastq_bwa.sam.bam,Heart_X11706He_AGTTCC_L001_R1_001_trimmed.fastq_bwa.sam.bam,Heart_X11699He_GGCTAC_L001_R1_001_trimmed.fastq_bwa.sam.bam,Heart_X11698He_TAGCTT_L001_R1_001_trimmed.fastq_bwa.sam.bam,Heart_X11700He_CTTGTA_L001_R1_001_trimmed.fastq_bwa.sam.bam - -[1] Using design matrix: - - (Intercept) groupLiver - -1 1 1 +R version 3.0.2 (2013-09-25) -2 1 1 - -3 1 1 - -4 1 1 - -5 1 0 - -6 1 0 - -7 1 0 - -8 1 0 - -attr(,"assign") - -[1] 0 1 - -attr(,"contrasts") - -attr(,"contrasts")$group - -[1] contr.treatment - -R version 3.1.0 (2014-04-10) - -Platform: x86_64-unknown-linux-gnu (64-bit) +Platform: x86_64-pc-linux-gnu (64-bit) locale: @@ -628,11 +581,11 @@ other attached packages: - [1] RColorBrewer_1.0-5 DESeq2_1.4.5 RcppArmadillo_0.4.400.0 Rcpp_0.11.2 GenomicRanges_1.16.4 GenomeInfoDb_1.0.2 IRanges_1.22.10 BiocGenerics_0.10.0 edgeR_3.6.8 limma_3.20.9 gplots_2.14.1 stringr_0.6.2 + [1] RColorBrewer_1.0-5 DESeq2_1.2.10 RcppArmadillo_0.4.400.0 Rcpp_0.11.2 GenomicRanges_1.14.4 XVector_0.2.0 IRanges_1.20.7 BiocGenerics_0.8.0 edgeR_3.4.2 limma_3.18.13 gplots_2.15.0 stringr_0.6.2 loaded via a namespace (and not attached): - [1] annotate_1.42.1 AnnotationDbi_1.26.0 Biobase_2.24.0 bitops_1.0-6 caTools_1.17 DBI_0.2-7 gdata_2.13.3 genefilter_1.46.1 geneplotter_1.42.0 grid_3.1.0 gtools_3.4.1 KernSmooth_2.23-12 lattice_0.20-29 locfit_1.5-9.1 RSQLite_0.11.4 stats4_3.1.0 survival_2.37-7 XML_3.98-1.1 xtable_1.7-3 XVector_0.4.0 + [1] annotate_1.40.1 AnnotationDbi_1.24.0 Biobase_2.22.0 bitops_1.0-6 caTools_1.17.1 DBI_0.3.0 gdata_2.13.3 genefilter_1.44.0 grid_3.0.2 gtools_3.4.1 KernSmooth_2.23-13 lattice_0.20-29 locfit_1.5-9.1 RSQLite_0.11.4 stats4_3.0.2 survival_2.37-7 XML_3.98-1.1 xtable_1.7-4@@ -642,28 +595,27 @@