+
+ if (opts->statistics & STATS_XPROD)
+ {
+ double cov = gsl_matrix_get (cv, r, col_index);
+ const double xprod_dev = cov * w;
+ cov *= w / (w - 1.0);
+
+ tab_double (t, c + heading_columns, row + 2, 0, xprod_dev, NULL);
+ tab_double (t, c + heading_columns, row + 3, 0, cov, NULL);
+ }