double IndependentMvnSuf::sample_var(int i)const{ double ni = n_[i]; if(ni - 1 < std::numeric_limits<double>::epsilon()){ return 0; } double ybari = ybar(i); double ss = sumsq_[i] - ni * ybari * ybari; return ss / (ni - 1); }
double NeRegSuf::SST()const{ return sumsqy - n()*pow(ybar(),2); }
SpdMatrix WMS::center_sumsq()const{ return center_sumsq(ybar());}
double QrRegSuf::SST()const{ // if(!current) refresh_qr(); return sumsqy - n()*pow(ybar(),2); }