Weighted least squares estimation for exchangeable binary data

Parametric models of discrete data with exchangeable dependence structure present substantial computational challenges for maximum likelihood estimation. Coordinate descent algorithms such as the Newton’s method are usually unstable, becoming a hit or miss adventure on initialization with a good starting value. We propose a method for computing maximum likelihood estimates of parametric models for finitely exchangeable binary data, formalized as an iterative weighted least squares algorithm.