Quantitative microbial risk assessment (QMRA) is a computational science leveraged to optimize infectious disease controls at both population and individual levels. Often, diverse populations will have different health risks based on a population's susceptibility or outcome severity due to heterogeneity within the host. Unfortunately, due to a host homogeneity assumption in the microbial dose-response models' derivation, the currentQMRAmethod of modeling exposure volume heterogeneity is not an accurate method for pathogens such as Legionella pneumophila. Therefore, a newmethod to model within-group heterogeneity is needed. The method developed in this research uses USA national incidence rates from the Centers for Disease Control and Prevention (CDC) to calculate proxies for the morbidity ratio that are descriptive of the within-group variability. From these proxies, an example QMRA model is developed to demonstrate their use. This method makes the QMRA results more representative of clinical outcomes and increases population-specific precision. Further, the risks estimated demonstrate a significant difference between demographic groups known to have heterogeneous health outcomes after infection. The method both improves fidelity to the real health impacts resulting from L. pneumophila infection and allows for the estimation of severe disability-adjusted life years (DALYs) for Legionnaires' disease, moderate DALYs for Pontiac fever, and post-acute DALYs for sequela after recovering from Legionnaires' disease.