Doubly robust pseudo-likelihood for incomplete hierarchical binary data. (February 2020)