Changeset 742 for trunk/test
 Timestamp:
 Jan 13, 2007, 10:14:05 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/test/regression_test.cc
r741 r742 131 131 << std::endl; 132 132 ok=false; 133 } 134 135 // Comparing LinearWeighted and PolynomialWeighted(1) 136 *error << " comparing LinearWeighted and PolynomialWeighted(1)" 137 << std::endl; 138 linear_w.fit(x,y,w); 139 regression::PolynomialWeighted polynomial_w(1); 140 polynomial_w.fit(x,y,w); 141 if ( fabs(linear_w.beta()polynomial_w.fit_parameters()(1))>10e7 ){ 142 *error << "error: beta and fit_parameters(1) not equal" << std::endl; 143 *error << " beta = " << linear_w.beta() << std::endl; 144 *error << " fit_parameters(1) = " 145 << polynomial_w.fit_parameters()(1) << std::endl; 146 ok = false; 147 } 148 if ( fabs(polynomial_w.fit_parameters()(0)linear_w.alpha()+ 149 linear_w.beta()*1990)>0.0001){ 150 *error << "error: fit_parameters(0) = " 151 << polynomial.fit_parameters()(0)<< std::endl; 152 *error << "error: alphabeta*m_x = " 153 << linear.alpha()linear.beta()*1990 << std::endl; 154 ok = false; 155 } 156 if ( fabs(polynomial_w.s2()linear_w.s2())>0.0001){ 157 *error << "error: chisq not same in linear and polynomial(1)" 158 << std::endl; 159 ok = false; 160 } 161 if ( fabs(polynomial_w.predict(1.0)linear_w.predict(1.0))>0.0001){ 162 *error << "error: predict not same in linear and polynomial(1)" 163 << std::endl; 164 ok = false; 165 } 166 if ( fabs(polynomial_w.standard_error2(1985)linear_w.standard_error2(1985)) 167 >0.0001){ 168 *error << "error: standard_error not same in linear and polynomial(1)" 169 << "\n polynomial: " << polynomial.standard_error2(1985) 170 << "\n linear: " << linear.standard_error2(1985) 171 << std::endl; 172 ok = false; 133 173 } 134 174
Note: See TracChangeset
for help on using the changeset viewer.