//maximum Gaussian probabilty plot correlation coefficient
double sampleR1 = max.getMaximum();
// Repeat maximization starting equidistant from the final value of lambdaOne on the opposite side from the starting estimate
start[0] = lambda1 - (start[0] - lambda1);
max.nelderMead(bcf, start, step, maxzTol);
coeff = max.getParamValues();
this.lambdaOne = coeff[0];
this.transformedSampleR = max.getMaximum();
// Choose solution with the largest Gaussian probabilty plot correlation coefficient