SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
EDGE-pro into DESeq, DESeq run error? epistatic Bioinformatics 5 09-02-2015 04:32 AM
javascript still not working merkutinaa RNA Sequencing 2 03-13-2014 11:00 AM
Working in MiReap njsphd Bioinformatics 1 01-03-2014 09:13 AM
JIGSAW not working volavii Bioinformatics 2 10-22-2013 07:42 AM
working in polyscan awalisa Bioinformatics 0 02-23-2012 08:01 PM

Reply
 
Thread Tools
Old 02-15-2016, 11:11 PM   #1
jay2008
Member
 
Location: Australia

Join Date: Sep 2010
Posts: 44
Default DESeq baseVarFunc is not working

can you help me to make function "baseVarFunc" below (from supplementary 2 of DESeq paper) work?
I got error:
Error in baseVarFunc(cdsFly, "A")(10^xg) : could not find function "rvf"

here is the code (file fly_RNA_counts.tsv was downloaded from supplementary 3).




baseVarFunc <- function( cds, cond ) {
rvf <- rawVarFunc( cds, cond )
sf <- sizeFactors(cds)[ conditions(cds) == cond ]
xim <- sum(1/sf) / length(sf)
function( q ) rvf( q ) + xim * q
}

rawVarFunc <- function( cds, condOrName ) {
stopifnot( is( cds, "CountDataSet" ) )
res <- cds@fitInfo[[ as.character(condOrName) ]]
if( is.null(res) ) {
res <- cds@fitInfo[[ cds@dispTable[ as.character(condOrName) ] ]]
if( is.null(res) )
stop( sprintf( "No base variance function found for condition or with name '%s'.", condOrName ) )
}
res
}


countsTableFly <- read.delim( "fly_RNA_counts.tsv" )
condsFly <- c( "A", "A", "B", "B" )
rownames( countsTableFly ) <- paste( "Gene", 1:nrow(countsTableFly), sep="_" )
cdsFly <- newCountDataSet( countsTableFly, condsFly )
cdsFly <- estimateSizeFactors( cdsFly )
cdsFly <- estimateDispersions( cdsFly )
gg<-log10( baseVarFunc(cdsFly,"A")( 10^xg ) ) ## this sentence got error


> sessionInfo()
R version 3.2.3 (2015-12-10)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

locale:
[1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
[5] LC_TIME=English_United States.1252

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

other attached packages:
[1] DESeq_1.22.1 lattice_0.20-33 locfit_1.5-9.1 Biobase_2.30.0 BiocGenerics_0.16.1

loaded via a namespace (and not attached):
[1] IRanges_2.4.7 XML_3.98-1.3 grid_3.2.3 xtable_1.8-2 DBI_0.3.1
[6] stats4_3.2.3 RSQLite_1.0.0 genefilter_1.52.1 annotate_1.48.0 S4Vectors_0.8.11
[11] splines_3.2.3 RColorBrewer_1.1-2 geneplotter_1.48.0 survival_2.38-3 AnnotationDbi_1.32.3
jay2008 is offline   Reply With Quote
Old 02-16-2016, 02:26 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Why are you trying to use DESeq? It's out dated and shouldn't be used for new projects (that's what DESeq2 is for).
dpryan is offline   Reply With Quote
Old 02-16-2016, 02:58 AM   #3
jay2008
Member
 
Location: Australia

Join Date: Sep 2010
Posts: 44
Default

I am going to investigate the difference between DESeq and DESeq2. the reason is that I can better understand the advantages of DESeq2 by reading, testing original DESeq.
you are right, as a user, who should simply pick up better tool, but I want to know more about DESeq than a normal user. if you can solve the problem, it will be great appreciated. thanks
Jay
jay2008 is offline   Reply With Quote
Old 02-16-2016, 03:55 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

It looks like that part was in error and I'm not familiar enough with how the 2010 version of DESeq worked internally to be able to reverse engineer things quickly. I would suggest just skipping that and following the vignette and user guide of the most recent version.
dpryan is offline   Reply With Quote
Old 02-16-2016, 12:27 PM   #5
jay2008
Member
 
Location: Australia

Join Date: Sep 2010
Posts: 44
Default

thanks, the code is from supplementary 2. actually, i try to do it work in current DESeq version.
can anyone else help this issue?
jay2008 is offline   Reply With Quote
Old 02-18-2016, 11:32 AM   #6
Wolfgang Huber
Senior Member
 
Location: Heidelberg, Germany

Join Date: Aug 2009
Posts: 109
Default

The original post was cross-posted, and is answered here: https://support.bioconductor.org/p/78427/#78550
__________________
Wolfgang Huber
EMBL
Wolfgang Huber is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:31 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO