# 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 @@