Changes to performance function: Added NA counts for Phi and Zeta, changed variance bias estimate to mean variance, enabled na.rm option on many calculations because NAs are being tracked elsewhere