I have a data set representing the abundance of a protein in a population of cells. Based on our understanding of the biology behind this, I expect there to be two subpopulations - one in which this protein is abundant and one in which it is not. The observed data seem to fit this hypothesis well, and certainly look bimodal (I'll try to attach an example below). If we assume that any subpopulations are themselves normally-distributed, is there a way to quantify the probability that this data is multimodal as opposed to being unimodal?
Looking through similar questions elsewhere on this exchange, a mixture model seems to be designed to do this, but I have no idea how to implement that. Is this something that could realistically be accomplished by a statistics novice?
Many thanks, O