def _get_mswd(self):
self.valid_mswd=False
if self._degree==1:
# a = self.intercept
# b = self.slope
coeffs=self._calculate_coefficients()
if not len(coeffs):
self.calculate()
coeffs = self._calculate_coefficients()
if len(coeffs):
x = self.xs
y = self.ys
sx = self.xserr
sy = self.yserr
if not len(sx):
sx=zeros(self.n)
if not len(sy):
sy=zeros(self.n)
x=self._clean_array(x)
y=self._clean_array(y)
sx=self._clean_array(sx)
sy=self._clean_array(sy)
m=calculate_mswd2(x, y, sx, sy, coeffs[1], coeffs[0])
self.valid_mswd=validate_mswd(m, len(ys), k=2)
return m
else:
return 0
else:
return super(OLSRegressor, self)._get_mswd()
def _get_mswd(self):
self.valid_mswd = False
if self._degree == 1:
# a = self.intercept
# b = self.slope
coeffs = self._calculate_coefficients()
if not len(coeffs):
self.calculate()
coeffs = self._calculate_coefficients()
if len(coeffs):
x = self.xs
y = self.ys
sx = self.xserr
sy = self.yserr
if not len(sx):
sx = zeros(self.n)
if not len(sy):
sy = zeros(self.n)
x = self._clean_array(x)
y = self._clean_array(y)
sx = self._clean_array(sx)
sy = self._clean_array(sy)
m = calculate_mswd2(x, y, sx, sy, coeffs[1], coeffs[0])
self.valid_mswd = validate_mswd(m, len(ys), k=2)
return m
else:
return 0
else:
return super(OLSRegressor, self)._get_mswd()
def test_chi_squared(self):
x, y, wx, wy = pearson()
b = -0.4807
a = 5.4806
ex = 1 / wx ** 0.5
ey = 1 / wy ** 0.5
v = 1.4833
m = calculate_mswd2(x, y, ex, ey, a, b, correlated_errors=True)
self.assertAlmostEqual(m, v, 4)
v = 2.2045
m = calculate_mswd2(x, y, ex, ey, a, b)
self.assertAlmostEqual(m, v, 4)
def test_chi_squared(self):
x, y, wx, wy = pearson()
b = -0.4807
a = 5.4806
ex = 1 / wx**0.5
ey = 1 / wy**0.5
v = 1.4833
m = calculate_mswd2(x, y, ex, ey, a, b, correlated_errors=True)
self.assertAlmostEqual(m, v, 4)
v = 2.2045
m = calculate_mswd2(x, y, ex, ey, a, b)
self.assertAlmostEqual(m, v, 4)
def _get_mswd(self):
if not self._slope:
self.calculate()
a = self.intercept
b = self.slope
x, y, sx, sy = self.clean_xs, self.clean_ys, self.clean_xserr, self.clean_yserr
v = calculate_mswd2(x, y, sx, sy, a, b, corrcoeffs=self._calculate_correlation_coefficients())
self.valid_mswd = validate_mswd(v, len(x), k=2)
return v
def _get_mswd(self):
if not self._slope:
self.calculate()
a = self.intercept
b = self.slope
x, y, sx, sy = self.clean_xs, self.clean_ys, self.clean_xserr, self.clean_yserr
v = calculate_mswd2(x, y, sx, sy, a, b,
corrcoeffs=self.calculate_correlation_coefficients())
self.valid_mswd = validate_mswd(v, len(x), k=2)
return v
def _get_mswd(self):
if not self._slope:
self.calculate()
a = self.intercept
b = self.slope
x = self.xs
y = self.ys
sx = self.xserr
sy = self.yserr
#print a
#print b
#print x
#print y
#print sx
#print sy
#print '----'
return calculate_mswd2(x, y, sx, sy, a, b)
def _get_mswd(self):
if not self._slope:
self.calculate()
a = self.intercept
b = self.slope
x = self.xs
y = self.ys
sx = self.xserr
sy = self.yserr
# print a
# print b
# print x
# print y
# print sx
# print sy
# print '----'
return calculate_mswd2(x, y, sx, sy, a, b)
