ACalculate the matrix of deviance residuals
getDevMat(X, thetaMat, mu)
X
the data matrix
thetaMat
the matrix of dispersions
mu
the matrix of means
The matrix of deviance residuals