Abstract The Kepler mission observed thousands of transiting exoplanet candidates around hundreds of thousands of FGK dwarf stars. He et al. applied forward modeling to infer the distribution of intrinsic architectures of planetary systems, developed a clustered Poisson point process model for exoplanetary systems ( SysSim ) to reproduce the marginal distributions of the observed Kepler population, and they showed that orbital periods and planet radii are clustered within a given planetary system. Here, we extend the clustered model to explore correlations between planetary systems and their host-star properties. We split the sample of Kepler FGK dwarfs into two halves and model the fraction of stars with planets (0.5–10 R ⊕ and 3–300 days), f swpa , as a linear function of the Gaia DR2 color. We confirm previous findings that the occurrence of these planetary systems rises significantly toward later-type (redder) stars. The fraction of stars with planets increases from for F2V dwarfs to for mid-K dwarfs. About half ( ) of all solar-type (G2V) dwarfs harbor a planetary system between 3 and 300 days. This simple model can closely match the observed multiplicity distributions of both the bluer and redder halves in our sample, suggesting that the architectures of planetary systems around stars of different spectral types may be similar aside from a shift in the overall fraction of planet-hosting stars.