Fitting the copCAR regression model for discrete areal data