theta_estimator function

Compute Watterson genetic diversity (Theta estimator)