[1]Dept.ofEcologyandEvolutionaryBiology,UniversityofColorado,Boulder,CO,USA,80309\keywordsoccupancymodel,speciesdistributionmodel,parasite,hierarchicalBayesianmodel,jointdistribution,latentstate Symbiont occurrence is often influenced by host occurrence, and vice versa. These influences lead to host-symbiont distributions that are correlated at multiple levels. For instance, interactions between co-infecting symbionts within host individuals might result in correlations in the abundance of two symbiont species across individual hosts. Similarly, interactions between symbiont transmission and host population dynamics might result in correlations between symbiont and host distributions at local levels (e.g. across habitat patches). Often, these correlations may be of biological interest, because they may reflect species interactions. Here, we present a flexible hierarchical Bayesian modeling framework that represents multi-level correlations in joint host-symbiont distributions. We demonstrate the basic method using simulated data, and provide a more realistic empirical case study using data from a multi-host multi-parasite amphibian system. Jointly modeling host and symbiont distributions facilitates a deeper understanding of patterns of co-occurrence and abundance that can be used to gain insights into host-symbiont interactions.