In benchdamic 1.4.0 (Bioconductor 3.16) in the function fitNB (and also in the other fit functions), should

else if (is.matrix(object)) {
        NULL
    }

be

else if (is.matrix(object)) {
        counts <- object
    }