| 34 | Var( const VEC_TYPE& object_ ) : object(object_){} |
| 35 | |
| 36 | double get() const{ |
| 37 | const double average = mean(object).get(); |
| 38 | const R_xlen_t sample_size = object.size(); |
| 39 | double sum_squared_deviations = 0.0; |
| 40 | for (R_xlen_t i = 0; i != sample_size; ++i) |
| 41 | sum_squared_deviations += std::pow(object[i] - average, 2.0); |
| 42 | return sum_squared_deviations / (sample_size - 1); |
| 43 | } |
| 44 | |
| 45 | private: |
| 46 | const VEC_TYPE& object ; |