Given a long molecular dynamics trajectory which consists of hundreds of association and dissociation events, the theoretical formulas to calculate the affinity and dissociation rate constant are presented. The derivation is based on the survival function of the associated complex, and it emphasizes the nearest-neighbor distribution function. The applicability of this brute-force approach is demonstrated by the simulation of methane association in water. (C) 2003 American Institute of Physics.