Add nlogtheta and nlogbeta to stats, add newton estimate.