powered by
make.varcov(regressor.matrix, y, subset=is.finite(y), casewt=NULL)
make.regressor.matrix
regressor.matrix
sum(subset==TRUE) - 1