40 #include "EST_multistats.h"
47 for (i = 0; i < v.
n(); ++i)
58 for (i = 0; i < v.
n(); ++i)
92 v[j] += pow(m(i, j) - u(j), float(2.0));
102 v = sample_variance(m);
105 for (j = 0; j < v.
n(); ++j)
124 c(j, k) += (m(i, j) - u(j)) * (m(i, k) - u(k));
143 r(j, k) = c(j, k)/(s(j) * s(k));
164 z(i, j) = (m(i, j) - sub(j)) / div(j);
179 cout <<
"pop mean " << gu;
181 for (i = 0; i < n; ++i)
182 for (j = 0; j < n; ++j)
185 for (k = 0; k < p; ++k)
186 P(i, j) += pow(gu(i, k) - gu(j, k),
float(2.0)) / gv(k);
213 for (i = 0; i < n; ++i)
214 for (j = 0; j < n; ++j)
216 transpose(row(gu, i),ui);
217 transpose(row(gu, j),uj);
218 P(i, j) = single_mahal(ui, uj, invv);
230 for (k = 0; k < n; ++k)
231 P += pow(ui(k) - uj(k),
float(2.0)) / v(k);
243 for (i = 0; i < num_pop; ++i)
257 for (i = 0; i < num_pop; ++i)
258 l += in[i].num_rows();
263 for (k = l = 0; k < num_pop; ++k)
264 for (j =0; j < n; ++j, ++l)
265 for (i = 0; i < in[i].
num_rows(); ++i)
266 msum(l, j) = in[k](i, j);