75 | | CHECK_CLOSE_EX(mean, num_mean2(hepdf.get(), xb, yb, nbins(0), nbins(1)), 0.1); |
76 | | CHECK_CLOSE_EX(1.0, normcoef(hepdf.get(), xb, yb, nbins(0), nbins(1)), 0.1); |
| 81 | CHECK_CLOSE_EX(mean, num_mean2(hepdf.get(), xb, yb, nbins(0), nbins(1)), tolerance); |
| 82 | CHECK_CLOSE_EX(1.0, normcoef(hepdf.get(), xb, yb, nbins(0), nbins(1)), tolerance); |
| 93 | |
| 94 | if (mrv.get()) { |
| 95 | RV crv = hepdf->_rv().subt(*mrv); |
| 96 | shared_ptr<epdf> m = hepdf->marginal(*mrv); |
| 97 | shared_ptr<mpdf> c = hepdf->condition(crv); |
| 98 | mepdf mm(m); |
| 99 | |
| 100 | Array<mpdf *> aa(2); |
| 101 | aa(0) = c.get(); |
| 102 | aa(1) = &mm; |
| 103 | mprod mEp(aa); |
| 104 | |
| 105 | int n = nsamples; |
| 106 | mat smp = mEp.samplecond(vec(0), n); |
| 107 | vec emu = sum(smp, 2) / n; |
| 108 | CHECK_CLOSE_EX(mean, emu, tolerance); |
| 109 | |
| 110 | if (R.rows() > 0) { |
| 111 | mat er = (smp * smp.T()) / n - outer_product(emu, emu); |
| 112 | CHECK_CLOSE_EX(R, er, tolerance); |
| 113 | } |
| 114 | |
| 115 | // test of pdflog at zero |
| 116 | vec zero(0); |
| 117 | vec zeron(hepdf->dimension()); |
| 118 | for (int i = 0; i < zeron.size(); ++i) { |
| 119 | zeron(i) = 0; |
| 120 | } |
| 121 | |
| 122 | CHECK_CLOSE_EX(hepdf->evallog(zeron), mEp.evallogcond(zeron, zero), tolerance); |
| 123 | } |