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?