[BioC] SeqGSEA estiGeneNBstat()

[SeqGSEA-user [guest]]

Hi I'm having trouble running the estiGeneNBstat() function.

Here is the error message I receive:

Error in estiGeneNBstat(counts) : Please run estiExonNBstat first.

But I have already run the estiExonNBstat() first. I basically just plugged in my data as show in the example in the vignette. Below is my code. I just feed in the counts from estiExonNBstat() into estiGeneNBstat(), so I'm not sure why I'm getting that error. Any insight will be appreciated. Thank you!

CODE:

#######################################################Step 1: DS analysis
# load exon read count data 
#RCS <- loadExonCountData(case.files, control.files) 
counts <- loadExonCountData(ps.files, cb.files) 

# remove genes with low exprssion 
#RCS <- exonTestability(RCS, cutoff=5) 
counts <- exonTestability(counts, cutoff=5) 

geneTestable <- geneTestability(counts) 

counts <- subsetByGenes(counts, unique(geneID(counts))[ geneTestable ]) 

# get gene IDs, which will be used in initialization of gene set 
#geneIDs <- unique(geneID(RCS)) 
geneIDs <- unique(geneID(counts))

# calculate DS NB statistics 
counts <- estiExonNBstat(counts) 	
counts <- estiGeneNBstat(counts) 


 -- output of sessionInfo(): 

> sessionInfo()
R version 3.0.0 (2013-04-03)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

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

other attached packages:
[1] SeqGSEA_1.0.2      foreach_1.4.1      biomaRt_2.16.0     DESeq_1.12.0       lattice_0.20-15    locfit_1.5-9.1     Biobase_2.20.0    
[8] BiocGenerics_0.6.0

loaded via a namespace (and not attached):
 [1] annotate_1.38.0      AnnotationDbi_1.22.6 codetools_0.2-8      compiler_3.0.0       DBI_0.2-7            doParallel_1.0.3    
 [7] genefilter_1.42.0    geneplotter_1.38.0   grid_3.0.0           IRanges_1.18.1       iterators_1.0.6      RColorBrewer_1.0-5  
[13] RCurl_1.95-4.1       RSQLite_0.11.4       splines_3.0.0        stats4_3.0.0         survival_2.37-4      tools_3.0.0         
[19] XML_3.95-0.2         xtable_1.7-1        
> 

--
Sent via the guest posting facility at bioconductor.org.