![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
DESeq2 error : the model matrix is not full rank | NicoBxl | Bioinformatics | 9 | 03-14-2017 06:10 AM |
Wired DESeq2 error | sindrle | Bioinformatics | 5 | 10-25-2013 05:29 AM |
Error Message in nbinomLRT in DESeq2 | ToddB | Bioinformatics | 13 | 09-05-2013 06:22 AM |
Tophat Error: Error: segment-based junction search failed with err =-6 | sjnewhouse | RNA Sequencing | 8 | 03-19-2013 04:14 AM |
Tophat 2 error: Parse error at line 297359: missing colon in auxiliary data | magbju | RNA Sequencing | 2 | 01-15-2013 03:57 AM |
![]() |
|
Thread Tools |
![]() |
#1 |
Member
Location: Norway Join Date: Jan 2010
Posts: 83
|
![]()
Hi,
I am simply following the DESeq2 tutorial, but I get an error during the varianceStabilizingTransformation: Code:
> vsd <- varianceStabilizingTransformation(dds, blind=TRUE) Error in varianceStabilizingTransformation(dds, blind = TRUE) : unused argument (blind = TRUE) > |
![]() |
![]() |
![]() |
#2 |
Devon Ryan
Location: Freiburg, Germany Join Date: Jul 2011
Posts: 3,478
|
![]()
What's the output of:
Code:
sessionInfo() |
![]() |
![]() |
![]() |
#3 |
Member
Location: Norway Join Date: Jan 2010
Posts: 83
|
![]() Code:
> sessionInfo() R version 3.0.2 (2013-09-25) Platform: x86_64-apple-darwin10.8.0 (64-bit) locale: [1] C attached base packages: [1] grid parallel stats graphics grDevices utils datasets methods base other attached packages: [1] pasilla_0.2.18 DESeq_1.14.0 lattice_0.20-24 locfit_1.5-9.1 DEXSeq_1.8.0 [6] Biobase_2.22.0 parathyroidSE_1.0.4 DESeq2_1.2.5 RcppArmadillo_0.3.920.1 Rcpp_0.10.6 [11] BiocInstaller_1.12.0 cummeRbund_2.4.0 Gviz_1.6.0 rtracklayer_1.22.0 GenomicRanges_1.14.3 [16] XVector_0.2.0 IRanges_1.20.5 fastcluster_1.1.11 reshape2_1.2.2 ggplot2_0.9.3.1 [21] RSQLite_0.11.4 DBI_0.2-7 BiocGenerics_0.8.0 loaded via a namespace (and not attached): [1] AnnotationDbi_1.24.0 BSgenome_1.30.0 Biostrings_2.30.1 GenomicFeatures_1.14.2 Hmisc_3.12-2 [6] MASS_7.3-29 RColorBrewer_1.0-5 RCurl_1.95-4.1 Rsamtools_1.14.1 XML_3.95-0.2 [11] annotate_1.40.0 biomaRt_2.18.0 biovizBase_1.10.3 bitops_1.0-6 cluster_1.14.4 [16] colorspace_1.2-4 dichromat_2.0-0 digest_0.6.3 genefilter_1.44.0 geneplotter_1.40.0 [21] gtable_0.1.2 hwriter_1.3 labeling_0.2 latticeExtra_0.6-26 munsell_0.4.2 [26] plyr_1.8 proto_0.3-10 rpart_4.1-3 scales_0.2.3 splines_3.0.2 [31] statmod_1.4.18 stats4_3.0.2 stringr_0.6.2 survival_2.37-4 tools_3.0.2 [36] xtable_1.7-1 zlibbioc_1.8.0 |
![]() |
![]() |
![]() |
#4 |
Devon Ryan
Location: Freiburg, Germany Join Date: Jul 2011
Posts: 3,478
|
![]()
It's because you have both DESeq and DESeq2 loaded. The version of that function in DESeq has a different syntax. That's why when you loaded DESeq, you got a message like:
Code:
The following objects are masked from ‘package:DESeq2’: estimateSizeFactorsForMatrix, getVarianceStabilizedData, plotPCA, varianceStabilizingTransformation |
![]() |
![]() |
![]() |
#5 |
Member
Location: Norway Join Date: Jan 2010
Posts: 83
|
![]()
Ok thanks!
I see that loading the Pasilla package from the DESeq2 tutorial forces DESeq to be loaded also. Do you know if there is any way to use the Pasilla dataset without loading DESeq? |
![]() |
![]() |
![]() |
#6 |
Devon Ryan
Location: Freiburg, Germany Join Date: Jul 2011
Posts: 3,478
|
![]()
Hmm, I never noticed that! You can just load pasilla first and then load DESeq2. That'll mask the DESeq functions. The pasilla library also requires and loads DEXSeq, which seems a bit silly as well.
|
![]() |
![]() |
![]() |
#7 |
Devon Ryan
Location: Freiburg, Germany Join Date: Jul 2011
Posts: 3,478
|
![]()
I went ahead and shot Alejandro Reyes (the maintainer for pasilla) an email. Hopefully he can just change the pasilla package so that no one else gets bit by this.
|
![]() |
![]() |
![]() |
#8 |
Senior Member
Location: Heidelberg Join Date: Aug 2010
Posts: 165
|
![]()
Thanks dpryan for your e-mail. After a short discussion with the authors of DESeq2, we agreed that the most convenient solution is that pasilla does not load DESeq or DEXSeq anymore. These changes are in the latest version of pasilla (0.2.19), which should be available in the next couple of days via biocLite.
In the meantime, you could load the objects from the pasilla vignette without loading the package like this: data("pasillaGenes", package="pasilla") data("pasillaExons", package="pasilla") Thanks for noticing this issue! Alejandro |
![]() |
![]() |
![]() |
Tags |
deseq2 |
Thread Tools | |
|
|