anchor.diff:    call counts_reset()
   
   end subroutine plugin
  +
  +subroutine plugin_from_multiplicities(n, nk1, zk1, estimate)
  +  ! plugin estimator - no prior, no regularization
  +  use constants
  +  use counter
  +  implicit none
  +
  +  integer, intent(in) :: n
  +  integer, intent(in) :: nk1(n)
  +  integer, intent(in) :: zk1(n)
  +  real(float64),  intent(out) :: estimate
  +
  +  integer :: i
  +  real(float64)   :: ni
  +
  +  call counts_fit(nk1, zk1)
  +
  +  estimate = 0.0_float64
  +  do i = 1,size(nk)
  +     ni = nk(i)
  +     if (ni > 0) estimate = estimate - zk(i)*ni*log(ni*1.0_float64)
  +  end do
  +  estimate = estimate / n_data + log(n_data)
  +
  +  call counts_reset()
  +
  +end subroutine plugin_from_multiplicities
   
   subroutine pseudo(n,counts,nc,alpha,estimate)
     ! pseudocount estimator(s)
  
positive.code: 
  end subroutine pseudo
  
  subroutine ww(n, counts, nc, alpha, estimate, err_estimate)
    ! posterior mean entropy (averaged over Dirichlet distribution) given alpha
    use constants
  
negative.code: 
  end module nsb_mod
  
  subroutine plugin(n, counts, estimate)
    ! plugin estimator - no prior, no regularization
    use constants
    use counter
  