Hi Davide,

I looked at the code for the functions zinbAIC and zinbBIC. It seems like the rows and columns are switched in the if function:

if ((nSamples(model) != nrow(x))|(nFeatures(model) != ncol(x))) {
stop("x and model should have the same dimensions!")
}

Shouldn't the rows of the x data matrix correspond to the features (genes) and columns to the samples (cells)?

Thanks,
-Reuben