Analyzing binary outcome data with small clusters: A simulation study