1Demographic inference using the site frequency spectrum (SFS) is a common way to understand 2 historical events affecting genetic variation. However, most methods for estimating demography 3 from the SFS assume random mating within populations, precluding these types of analyses in in-4 bred populations. To address this issue, we developed a model for the expected SFS that includes 5 inbreeding by parameterizing individual genotypes using beta-binomial distributions. We then 6 take the convolution of these genotype probabilities to calculate the expected frequency of bial-7 lelic variants in the population. Using simulations, we evaluated the model's ability to co-estimate 8 demography and inbreeding using one-and two-population models across a range of inbreeding 9 levels. We also applied our method to two empirical examples, American pumas (Puma concolor) 10 and domesticated cabbage (Brassica oleracea var. capitata), inferring models both with and with-11 out inbreeding to compare parameter estimates and model fit. Our simulations showed that we 12 are able to accurately co-estimate demographic parameters and inbreeding even for highly inbred 13 populations (F = 0.9). In contrast, failing to include inbreeding generally resulted in inaccurate 14 parameter estimates in simulated data and led to poor model fit in our empirical analyses. These 15 results show that inbreeding can have a strong effect on demographic inference, a pattern that was 16 especially noticeable for parameters involving changes in population size. Given the importance 17 of these estimates for informing practices in conservation, agriculture, and elsewhere, our method 18 provides an important advancement for accurately estimating the demographic histories of these 19 species. 20 Estimating the demographic history of closely related populations or species is an important first 23 step in understanding the interplay of the evolutionary forces shaping genetic variation. Diver-24 gence, migration, changes in population size, and other historical events all contribute to pop-25 ulation allele frequency dynamics over time, a process that can be modeled using a variety of 26 approaches. Connecting the expectations from these models with observed genomic data is of-27 ten achieved using the site frequency spectrum (SFS), a genome-wide summary of genetic poly-28 morphism within and between populations (Sawyer and Hartl 1992; Adams and Hudson 2004; 29 Caicedo et al. 2007; Gutenkunst et al. 2009; Nielsen et al. 2009). The ease and affordability of col-30 lecting genomic SNP data make inferences of demography using the SFS especially appealing, 31 highlighting their importance in gaining insights into the historical factors affecting neutral varia-32 tion in populations. Several recent analyses have also applied SFS-based methods to infer the fit-33 ness effects of mutations (Kim et al. 2017; Tataru et al. 2017; Fortier et al. 2019), allowing researchers 34 to model patterns of selection while simultaneously controlling for demography (Willia...