add parameterization = "theta" (phase 1); add WRMR to list of fit statistics