Question: JunctionSeq Can't produce Summary Plots
gravatar for rachelrodgers12
2.6 years ago by
rachelrodgers120 wrote:


I have run JunctionSeq on a large RNA-seq data set, and have managed to produce the comprehensive set of plots using buildAllPlots() following the JunctionSeq user manual.  However, when I try to run plotDispEsts(jscs), a plot is produced but seems corrupted and won't open.  Here is the output to the console which looks like things are running ok:

> plotDispEsts(jscs)
     abundance ranges from 0.010150733280822 to 288462.038735504
     dispersion ranges from 1.00626741453259e-06 to 104
     Plotting dispersions from 0.00112604526834281 to 104

Next when I try to plotMA(jscs, FDR.threshold=0.05) I receive this error:

Error in plotMA(jscs, FDR.threshold = 0.05) :
  Error: you must specify the name of the column containing the desired fold change variable.

I have tried consulting the manuals to see what argument I am missing, but I can't find it.

Any help is appreciated,



Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu precise (12.04.5 LTS)

 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C
 [9] LC_ADDRESS=C               LC_TELEPHONE=C

attached base packages:
[1] parallel  stats4    stats     graphics  grDevices utils     datasets
[8] methods   base

other attached packages:
 [1] Cairo_1.5-9                BiocInstaller_1.24.0
 [3] JunctionSeq_1.4.0          SummarizedExperiment_1.4.0
 [5] Biobase_2.34.0             GenomicRanges_1.26.3
 [7] GenomeInfoDb_1.10.3        IRanges_2.8.1
 [9] S4Vectors_0.12.1           BiocGenerics_0.20.0

loaded via a namespace (and not attached):
 [1] genefilter_1.56.0    statmod_1.4.29       locfit_1.5-9.1
 [4] splines_3.3.3        lattice_0.20-35      colorspace_1.3-2
 [7] htmltools_0.3.5      base64enc_0.1-3      survival_2.41-3
[10] XML_3.98-1.5         DBI_0.6              foreign_0.8-67
[13] BiocParallel_1.8.1   RColorBrewer_1.1-2   plyr_1.8.4
[16] stringr_1.2.0        zlibbioc_1.20.0      munsell_0.4.3
[19] gtable_0.2.0         DESeq2_1.14.1        htmlwidgets_0.8
[22] memoise_1.0.0        latticeExtra_0.6-28  knitr_1.15.1
[25] geneplotter_1.52.0   AnnotationDbi_1.36.2 htmlTable_1.9
[28] Rcpp_0.12.9          acepack_1.4.1        xtable_1.8-2
[31] scales_0.4.1         backports_1.0.5      checkmate_1.8.2
[34] plotrix_3.6-4        Hmisc_4.0-2          annotate_1.52.1
[37] XVector_0.14.0       gridExtra_2.2.1      ggplot2_2.2.1
[40] digest_0.6.12        stringi_1.1.2        grid_3.3.3
[43] tools_3.3.3          bitops_1.0-6         magrittr_1.5
[46] RSQLite_1.1-2        lazyeval_0.2.0       RCurl_1.95-4.8
[49] tibble_1.2           Formula_1.2-1        cluster_2.0.6
[52] Matrix_1.2-8         data.table_1.10.4    assertthat_0.1
[55] rpart_4.1-10         nnet_7.3-12



junctionseq • 329 views
ADD COMMENTlink written 2.6 years ago by rachelrodgers120

Update:  found dispersion and MA plots within the directory with all other plots produced by buildAllPlots().  I assume this is what I am looking for.


ADD REPLYlink written 2.6 years ago by rachelrodgers120
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 172 users visited in the last hour