add code with number counts from bethermin et. al. 2011