Observations of exoplanetary systems provide clues about the intrinsic distribution of planetary systems, their architectures, and how they formed. We develop a forward modeling framework for generating populations of planetary systems and "observed" catalogs by simulating the Kepler detection pipeline (SysSim). We compare our simulated catalogs to the Kepler DR25 catalog of planet candidates, updated to include revised stellar radii from Gaia DR2. We constrain our model based on the observed 1-D marginal distributions of orbital periods, period ratios, transit depths, transit depth ratios, transit durations, transit duration ratios, and transit multiplicities. Models assuming planets with independent periods and sizes do not adequately account for the properties of the multi-planet systems. Instead, a clustered point process model for exoplanet periods and sizes provides a significantly better description of the Kepler population, particularly the observed multiplicity and period ratio distributions. We find that 0.56 +0.18 −0.15 of FGK stars have at least one planet larger than 0.5R ⊕ between 3 and 300 days. Most of these planetary systems (∼ 98%) consist of one or two clusters with a median of three planets per cluster. We find that the Kepler dichotomy is evidence for a population of highlyinclined planetary systems and is unlikely to be solely due to a population of intrinsically single planet systems. We provide a large ensemble of simulated physical and observed catalogs of planetary systems from our models, as well as publicly available code for generating similar catalogs given user-defined parameters.