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