DEXSeq: Error in perGeneQValue with a single gene
0
0
Entering edit mode
@matmu
Last seen 12 weeks ago
Germany

I am doing a candidate gene study with DEXSeq v1.36.0 and would like to calculate test statistics for a single gene.

When running the code


dxd.cand = testForDEU(dxd[<gene_xy>, ], reducedModel = reducedModel)
dxr.cand = DEXSeqResults(dxd.cand, independentFiltering=FALSE)
qval.cand = perGeneQValue(dxr.cand)

function perGeneQValue throws an error:

Error in base::rowSums(x, na.rm = na.rm, dims = dims, ...) : 
  'x' must be an array of at least two dimensions

For two or more genes it works fine. Is this a bug or can a gene Q value not be calculated for a single gene?

DEXSeq • 268 views
ADD COMMENT

Login before adding your answer.

Traffic: 530 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6