code: def disasters(value=disasters_array,
                early_mean=early_mean,
                late_mean=late_mean,
                switchpoint=switchpoint):
      
      return pm.poisson_like(
          value[:switchpoint], early_mean) + pm.poisson_like(value[switchpoint:], late_mean)
docs: Annual occurences of coal mining disasters.