The statistical thermodynamics of adsorption in caged zeolites is developed by treating the zeolite as an ensemble of M identical cages or subsystems. Within each cage adsorption is assumed to occur onto a lattice of n identical sites. Expressions for the average occupancy per cage are obtained by minimizing the Helmholtz free energy in the canonical ensemble subject to the constraints of constant M and constant number of adsorbates N. Adsorbate-adsorbate interactions in the Brag-Williams or mean field approximation are treated in two ways. The local mean field approximation ͑LMFA͒ is based on the local cage occupancy and the global mean field approximation ͑GMFA͒ is based on the average coverage of the ensemble. The GMFA is shown to be equivalent in formulation to treating the zeolite as a collection of interacting single site subsystems. In contrast, the treatment in the LMFA retains the description of the zeolite as an ensemble of identical cages, whose thermodynamic properties are conveniently derived in the grand canonical ensemble. For a z coordinated lattice within the zeolite cage, with ⑀ aa as the adsorbate-adsorbate interaction parameter, the comparisons for different values of ⑀ aa * ϭ⑀ aa z/2kT, and number of sites per cage, n, illustrate that for Ϫ1Ͻ⑀ aa * Ͻ0 and nу10, the adsorption isotherms and heats of adsorption predicted with the two approaches are similar. In general, the deviation between the LMFA and GMFA is greater for smaller n and less sensitive to n for ⑀ aa * Ͼ0. We compare the isotherms predicted with the LMFA with previous GMFA predictions ͓K. G. Ayappa, C. R. Kamala, and T. A. Abinandanan, J. Chem. Phys. 110, 8714 ͑1999͔͒ ͑which incorporates both the site volume reduction and a coverage-dependent ⑀ aa ͒ for xenon and methane in zeolite NaA. In all cases the predicted isotherms are very similar, with the exception of a small steplike feature present in the LMFA for xenon at higher coverages.