protected void plotParameterPlane() {
double x1 = x.getElement(0);
double x2 = x.getElement(1);
double dx1 = Math.sqrt(cx.getElement(0, 0));
double dx2 = Math.sqrt(cx.getElement(1, 1));
double rho = (cx.getElement(1, 0)) / (dx1 * dx2);
// prepare size of plot
double xmin = x1 - 2. * dx1;
double xmax = x1 + 2. * dx1;
double ymin = x2 - 2. * dx2;
double ymax = x2 + 2. * dx2;
DatanGraphics.openWorkstation(getClass().getName(), "E3Min_2.ps");
DatanGraphics.setFormat(0., 0.);
DatanGraphics.setWindowInComputingCoordinates(xmin, xmax, ymin, ymax);
DatanGraphics.setViewportInWorldCoordinates(.2, .9, .16, .86);
DatanGraphics.setWindowInWorldCoordinates(-.414, 1., 0., 1.);
DatanGraphics.setBigClippingWindow();
DatanGraphics.chooseColor(2);
DatanGraphics.drawFrame();
DatanGraphics.drawScaleX("x_1");
DatanGraphics.drawScaleY("x_2");
DatanGraphics.drawBoundary();
DatanGraphics.chooseColor(5);
// draw data point with errors (and correlation)
DatanGraphics.drawDatapoint(1, 1., x1, x2, dx1, dx2, rho);
DatanGraphics.chooseColor(2);
DatanGraphics.drawCaption(1., caption);
// draw confidence contour
double fcont = mllg.getValue(x) + .5;
int nx = 100;
int ny = 100;
double dx = (xmax - xmin) / (int) nx;
double dy = (ymax - ymin) / (int) ny;
MinLogLikeGaussCont mllgc = new MinLogLikeGaussCont();
// System.out.println(" x = " + x.toString() + ", mllgc.getValue(x) = " +
// mllgc.getValue(x.getElement(0), x.getElement(1)));
DatanGraphics.setBigClippingWindow();
DatanGraphics.chooseColor(1);
DatanGraphics.drawContour(xmin, ymin, dx, dy, nx, ny, fcont, mllgc);
// draw asymmetric errors as horiyontal and vertical bars
DatanGraphics.chooseColor(3);
double[] xpl = new double[2];
double[] ypl = new double[2];
for (int i = 0; i < 2; i++) {
if (i == 0) xpl[0] = x1 - dxasy[0][0];
else xpl[0] = x1 + dxasy[0][1];
xpl[1] = xpl[0];
ypl[0] = ymin;
ypl[1] = ymax;
DatanGraphics.drawPolyline(xpl, ypl);
}
for (int i = 0; i < 2; i++) {
if (i == 0) ypl[0] = x2 - dxasy[1][0];
else ypl[0] = x2 + dxasy[1][1];
ypl[1] = ypl[0];
xpl[0] = xmin;
xpl[1] = xmax;
DatanGraphics.drawPolyline(xpl, ypl);
}
DatanGraphics.closeWorkstation();
}
