/** * Contructs the probability matrix for mixture estimation, given a set of support points and a * set of intervals. * * @param s the set of support points * @param intervals the intervals * @return the probability matrix */ public PaceMatrix probabilityMatrix(DoubleVector s, PaceMatrix intervals) { int ns = s.size(); int nr = intervals.getRowDimension(); PaceMatrix p = new PaceMatrix(nr, ns); for (int i = 0; i < nr; i++) { for (int j = 0; j < ns; j++) { p.set( i, j, Maths.pchisq(intervals.get(i, 1), s.get(j)) - Maths.pchisq(intervals.get(i, 0), s.get(j))); } } return p; }