def bernoulli_lh(ho, ha, s, n):
	"""
	Returns the likelihood ratio for independently distributed bernoulli random variables.
	Parameters
	----------
	ho : float
	   null hypothesis
	ha : float
	   alternative hypothesis
	s : float or int
	   number of successes in sample
	n : float or int
	   total number of elements in sample
	Returns
	-------
	float
	   likelihood ratio of model
	"""

	null_lh = (ho ** s) * (1 - ho)**(n - s)
	alt_lh = (ha ** s) * (1 - ha)**(n - s)

	return alt_lh / null_lh