def __init__(self, x, y):
x = marray(x, copy=False, subok=True, dtype=float_, order='F').ravel()
y = marray(y, copy=False, subok=True, dtype=float_, order='F').ravel()
if x.size != y.size:
msg = "Incompatible size between observations (%s) and response (%s)!"
raise ValueError(msg % (x.size, y.size))
idx = x.argsort()
self._x = x[idx]
self._y = y[idx]
self._mask = mask_or(self._x._mask, self._y._mask,
copy=False, small_mask=False)
def __init__(self, n):
self._seasonal = marray(nempty((n,), float_))
self._trend = marray(nempty((n,), float_))
self._weights = marray(nempty((n,), float_))
self._residuals = marray(nempty((n,), float_))
def __init__(self, y):
self.y = marray(y, subok=True, copy=False).ravel()
self._mask = self.y._mask
if self._mask.any():
raise ValueError("Masked arrays should be filled first!")
self.y_eff = self.y.compressed()
def __init__(self, n):
self._fval = marray(nempty((n,), dtype=float_, order='F'))
self._rw = marray(nempty((n,), dtype=float_, order='F'))
self._fres = marray(nempty((n,), dtype=float_, order='F'))
#####---------------------------------------------------------------------------
#--- --- Loess ---
#####---------------------------------------------------------------------------
loess = _mloess.loess
loess_anova = _mloess.anova
################################################################################
if __name__ == '__main__':
from maskedarray.testutils import assert_almost_equal
from maskedarray import masked_values
from numpy import fromiter
import os
if 1:
NOx = marray([4.818, 2.849, 3.275, 4.691, 4.255, 5.064, 2.118, 4.602,
2.286, 0.970, 3.965, 5.344, 3.834, 1.990, 5.199, 5.283,
-9999, -9999, 3.752, 0.537, 1.640, 5.055, 4.937, 1.561])
NOx = maskedarray.masked_values(NOx, -9999)
E = marray([0.831, 1.045, 1.021, 0.970, 0.825, 0.891, 0.71, 0.801,
1.074, 1.148, 1.000, 0.928, 0.767, 0.701, 0.807, 0.902,
-9999, -9999, 0.997, 1.224, 1.089, 0.973, 0.980, 0.665])
gas_fit_E = numpy.array([0.665, 0.949, 1.224])
newdata = numpy.array([0.6650000, 0.7581667, 0.8513333, 0.9445000,
1.0376667, 1.1308333, 1.2240000])
coverage = 0.99
rfile = open(os.path.join('tests','gas_result'), 'r')
results = []
for i in range(8):
rfile.readline()
z = fromiter((float(v) for v in rfile.readline().rstrip().split()),