Marginal modeling of spatial clustered binary data