Marginal logistic regression for spatially clustered binary data