We present a general theory for the calculation of various characteristic properties of the beam of atoms emerging from a resonator in one-atom-maser experiments. The beam is described in terms of the statistics of the detector clicks. The evolution of the state of the maser photons between clicks is governed by a nonlinear master equation. The nonlinearity originates in the necessity to account for the atoms that escape detection. Despite the permanent reductions of the photon state, resulting from the detections, the steady state of the conventional linear master equation determines the statistics of the detector clicks. The whole process is ergodic, in the sense that a single run of the experiment contains all reproducible statistical data, provided the duration of the run is much longer than all relevant correlation times. The nonlinear master equation is used to calculate the distribution of waiting times between detector clicks. Other statistical properties of the clicks that are derived include correlation functions and variances of the counting statistics. The formalism is applied to standard one-atom-maser experiments and to parity measurements on both unpumped and pumped cavities. We find that, for the standard one-atom-maser operation, none of the said statistical properties is a simple immediate indicator for a sub-Poissonian variance of the photon number. For example, the detector clicks may be antibunched although the photon distribution has a super-Poissonian variance.