Notably, popular contemporary polymer chemistry textbooks either do not discuss, or only briefly mention linearized methods derived from simplified copolymer equations, in the context of characterizing equilibrium copolymerization. Moreover, some authors acknowledge the presence of reversible propagation in their experimental system yet choose to ignore the inclusion of these effects when characterizing copolymerization reactivity, likely due to misconceptions about the accuracy and useability of preexisting methods. Finally, despite the wide array of methods and simplifications employed in all the above cases, qualitative agreements between the models and experimental data were generally observed, further complicating choice of method, as even simplified or linearized methods can provide suitable fits to data under conditions where their key assumptions are violated (vide infra). ultimately present a simple method to fit equilibrium copolymerization data to extract meaningful copolymerization reactivity information, including reactivity ratios. Our motivation to resolve this knowledge gap is threefold. First, polymers synthesized through mechanisms with accessible ceiling temperatures can produce circular, chemically recyclable materials, and we envision copolymerization of two or more suitable comonomers to be an enabling tactic within such systems to access circular materials with performance advantaged properties. For example, butyrolactone can be copolymerized with caprolactone or cyclic acetals to modify crystallinity and mechanical properties, and myriad bio-derived methacrylate-type monomers can be copolymerized to tune glass transition temperature and tensile properties. Second, the recent development of cleavable comonomers that polymerize under equilibrium conditions have shown great promise as additives to imbue degradability to both thermoplastics and thermosets. In the context of ring-opening metathesis polymerization (ROMP), low-strain cyclic silyl ethers and enol ethers were recently reported by our groups to copolymerize with traditional norbornene-based monomers to yield deconstructable linear, graft, and thermoset copolymers. Similarly, cyclic disulfides, isocyanides, and thioamides have recently emerged as competent cleavable comonomers with broad classes of vinyl monomers via radical polymerization. In this context, the distribution of these comonomers within a polymer backbone (i.e., the comonomer sequence) dictates degradation properties; thus, accurate understanding of comonomer reactivity is vital for predicting and engineering the end-of-life properties of these materials. Finally, the emerging use of data driven research practices in combination with machine learning methods relies on large bodies of accurate and high-quality data to develop broadly useful tools for property and reactivity predictions. In this work, we investigate three distinct methods to characterize equilibrium copolymerization under batch conditions and extract reactivity information: (1) stochastic simulations, (2) kinetic and population balance modeling, and (3) direct numerical integration of an equilibrium copolymer equation. Stochastic simulations can incorporate arbitrary reactivity complexity with minimal practical issues, allowing us to adapt previously reported scripts and methods applied to irreversible copolymerization systems; however, the use of such methods for the characterization of copolymerization data is scarce and, to our knowledge, no rigorous comparisons to canonical methods for irreversible systems (e.g., the Meyer-Lowry equation) have been conducted. On the other hand, population balance modeling can be viewed as a kinetically accurate way to describe monomer consumption over the course of a polymerization reaction, as has been shown previously for irreversible copolymerization. Unfortunately, previous attempts to model equilibrium copolymerization with population balance equations have noted an apparent need to account for every possible monomer sequence when both monomers undergo reversible polymerization, which necessitates the use of a computationally infeasible number of equations; thus, simplifications have been widely adopted. Here, to make kinetic and population balance modeling of equilibrium copolymerization tractable, we introduce a novel technique that preserves the relative rates of monomer consumption reactions while significantly reducing the number of equations required to accurately describe the system. Finally, we reevaluate three general equilibrium copolymer equations from the literature, rectify critical errors present in their initial reports and show that they are broadly equivalent despite their disparate forms.