[BioC] DESeq without replicates, Problem trying to obtain significant hits
Melissa.Martin at lshtm.ac.uk
Melissa.Martin at lshtm.ac.uk
Wed Mar 14 15:50:56 CET 2012
I am trying to use the DEseq package without replicates. When I look for significant hits all of them are False. Could this be a real result? Below are the commands that I enter.
> countsTable <-read.delim("spore_readcounts.txt",header=TRUE,stringsAsFactors=TRUE)
> rownames( countsTable ) <- countsTable$gene_id
> countsTable <- countsTable[ , -1 ]
> head(countsTable)
> cds <- newCountDataSet( countsTable, conds )
> cds <- estimateSizeFactors( cds )
> sizeFactors( cds )
> head( counts( cds, normalized=TRUE ) )
> cds <- cds[ ,c( "1", "2" ) ]
> cds <- estimateDispersions( cds, method="blind", sharingMode="fit-only" )
> res <- nbinomTest( cds, "1", "2" )
> plotDE( res )
Error: could not find function "plotDE"
> addmargins( table( res_sig = res$padj < .1, res_sig = res$padj < .1 ) )
res_sig
res_sig FALSE Sum
FALSE 4960 4960
Sum 4960 4960
Many Thanks
Melissa
More information about the Bioconductor
mailing list