public void EMpar(Database db, double prior) {
int j, np;
double x, y, tmp, tmpx, tmpy, tmpsx, tmpsy, tmpsxy;
np = db.nPoints();
tmpsx = tmpsy = tmpsxy = kmx = kmy = 0;
for (j = 0; j < np; j++) {
x = db.xVal(j);
y = db.yVal(j);
kmx += probs[j] * x;
kmy += probs[j] * y;
tmpsx += probs[j] * x * x;
tmpsy += probs[j] * y * y;
tmpsxy += probs[j] * x * y;
}
tmp = np * weight;
kmx /= tmp;
kmy /= tmp;
ksx = Math.sqrt(tmpsx / tmp - kmx * kmx);
ksy = Math.sqrt(tmpsy / tmp - kmy * kmy);
ksxy = tmpsxy / tmp - kmx * kmy;
if (ksx < mins) ksx = mins;
if (ksy < mins) ksy = mins;
weight = 0.9 * weight + 0.1 * prior;
}
