Traditional mass action laws that govern defect concentrations are based on the assumptions of negligible interactions and random distributions, and are not applicable in the high-defect-concentration regime. General equations applicable in this regime can be obtained using a statistical thermodynamic formulation that explicitly incorporates long-range Coulombic interactions and the generalized exclusion effects. This is demonstrated using calcium-doped cerium oxide as a model system.