def obs_likelihood(obs, freq):
	"""Bernoulli probability of observing binary presence / absence given circulating frequency"""
	"""obs is 0 or 1 binary indicator for absence / presence"""
	return (1-freq)**(1-obs) * freq**obs