scikit-misc 0.5.1.dev3__cp313-cp313-win_amd64.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- scikit_misc-0.5.1.dev3-cp313-cp313-win_amd64.whl +0 -0
- scikit_misc-0.5.1.dev3.dist-info/LICENSE +26 -0
- scikit_misc-0.5.1.dev3.dist-info/METADATA +122 -0
- scikit_misc-0.5.1.dev3.dist-info/RECORD +21 -0
- scikit_misc-0.5.1.dev3.dist-info/WHEEL +4 -0
- skmisc/__config__.py +117 -0
- skmisc/__init__.py +43 -0
- skmisc/_distributor_init.py +10 -0
- skmisc/_version.py +4 -0
- skmisc/loess/__init__.py +64 -0
- skmisc/loess/_loess.cp313-win_amd64.lib +0 -0
- skmisc/loess/_loess.cp313-win_amd64.pyd +0 -0
- skmisc/loess/tests/__init__.py +0 -0
- skmisc/loess/tests/co2_data +87 -0
- skmisc/loess/tests/co2_results_double +7 -0
- skmisc/loess/tests/co2_results_real +7 -0
- skmisc/loess/tests/gas_data +3 -0
- skmisc/loess/tests/gas_result +10 -0
- skmisc/loess/tests/madeup_data +4 -0
- skmisc/loess/tests/madeup_result +16 -0
- skmisc/loess/tests/test_loess.py +284 -0
|
File without changes
|
|
@@ -0,0 +1,26 @@
|
|
|
1
|
+
Copyright (c) 2016, Hassan Kibirige
|
|
2
|
+
All rights reserved.
|
|
3
|
+
|
|
4
|
+
Redistribution and use in source and binary forms, with or without
|
|
5
|
+
modification, are permitted provided that the following conditions are met:
|
|
6
|
+
|
|
7
|
+
* Redistributions of source code must retain the above copyright notice,
|
|
8
|
+
this list of conditions and the following disclaimer.
|
|
9
|
+
* Redistributions in binary form must reproduce the above copyright
|
|
10
|
+
notice, this list of conditions and the following disclaimer in the
|
|
11
|
+
documentation and/or other materials provided with the distribution.
|
|
12
|
+
* Neither the name of nor the names of its contributors may be used to
|
|
13
|
+
endorse or promote products derived from this software without specific
|
|
14
|
+
prior written permission.
|
|
15
|
+
|
|
16
|
+
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
|
|
17
|
+
AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
|
18
|
+
IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
|
19
|
+
ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
|
|
20
|
+
LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
|
|
21
|
+
CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
|
|
22
|
+
SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
|
|
23
|
+
INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
|
|
24
|
+
CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
|
|
25
|
+
ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
|
|
26
|
+
POSSIBILITY OF SUCH DAMAGE.
|
|
@@ -0,0 +1,122 @@
|
|
|
1
|
+
Metadata-Version: 2.1
|
|
2
|
+
Name: scikit-misc
|
|
3
|
+
Version: 0.5.1.dev3
|
|
4
|
+
Summary: Miscellaneous tools for scientific computing.
|
|
5
|
+
Author-Email: Hassan Kibirige <has2k1@gmail.com>
|
|
6
|
+
License: Copyright (c) 2016, Hassan Kibirige
|
|
7
|
+
All rights reserved.
|
|
8
|
+
|
|
9
|
+
Redistribution and use in source and binary forms, with or without
|
|
10
|
+
modification, are permitted provided that the following conditions are met:
|
|
11
|
+
|
|
12
|
+
* Redistributions of source code must retain the above copyright notice,
|
|
13
|
+
this list of conditions and the following disclaimer.
|
|
14
|
+
* Redistributions in binary form must reproduce the above copyright
|
|
15
|
+
notice, this list of conditions and the following disclaimer in the
|
|
16
|
+
documentation and/or other materials provided with the distribution.
|
|
17
|
+
* Neither the name of nor the names of its contributors may be used to
|
|
18
|
+
endorse or promote products derived from this software without specific
|
|
19
|
+
prior written permission.
|
|
20
|
+
|
|
21
|
+
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
|
|
22
|
+
AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
|
23
|
+
IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
|
24
|
+
ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
|
|
25
|
+
LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
|
|
26
|
+
CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
|
|
27
|
+
SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
|
|
28
|
+
INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
|
|
29
|
+
CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
|
|
30
|
+
ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
|
|
31
|
+
POSSIBILITY OF SUCH DAMAGE.
|
|
32
|
+
|
|
33
|
+
Classifier: Intended Audience :: Science/Research
|
|
34
|
+
Classifier: License :: OSI Approved :: BSD License
|
|
35
|
+
Classifier: Operating System :: Unix
|
|
36
|
+
Classifier: Operating System :: MacOS
|
|
37
|
+
Classifier: Operating System :: Microsoft :: Windows
|
|
38
|
+
Classifier: Operating System :: POSIX
|
|
39
|
+
Classifier: Programming Language :: C
|
|
40
|
+
Classifier: Programming Language :: Fortran
|
|
41
|
+
Classifier: Programming Language :: Python :: 3
|
|
42
|
+
Classifier: Programming Language :: Python :: 3.10
|
|
43
|
+
Classifier: Programming Language :: Python :: 3.11
|
|
44
|
+
Classifier: Programming Language :: Python :: 3.12
|
|
45
|
+
Classifier: Programming Language :: Python :: 3.13
|
|
46
|
+
Classifier: Topic :: Scientific/Engineering
|
|
47
|
+
Project-URL: homepage, https://has2k1.github.io/scikit-misc/stable
|
|
48
|
+
Project-URL: repository, https://github.com/has2k1/scikit-misc
|
|
49
|
+
Project-URL: documentation, https://has2k1.github.io/scikit-misc/stable
|
|
50
|
+
Project-URL: changelog, https://has2k1.github.io/scikit-misc/stable/changelog.html
|
|
51
|
+
Requires-Python: >=3.10
|
|
52
|
+
Requires-Dist: numpy>=1.23.5
|
|
53
|
+
Provides-Extra: all
|
|
54
|
+
Requires-Dist: scikit-misc[doc]; extra == "all"
|
|
55
|
+
Requires-Dist: scikit-misc[build]; extra == "all"
|
|
56
|
+
Requires-Dist: scikit-misc[lint]; extra == "all"
|
|
57
|
+
Requires-Dist: scikit-misc[test]; extra == "all"
|
|
58
|
+
Requires-Dist: scikit-misc[dev]; extra == "all"
|
|
59
|
+
Provides-Extra: doc
|
|
60
|
+
Requires-Dist: numpydoc; extra == "doc"
|
|
61
|
+
Requires-Dist: sphinx; extra == "doc"
|
|
62
|
+
Provides-Extra: build
|
|
63
|
+
Requires-Dist: cython; extra == "build"
|
|
64
|
+
Requires-Dist: build; extra == "build"
|
|
65
|
+
Requires-Dist: spin; extra == "build"
|
|
66
|
+
Requires-Dist: meson; extra == "build"
|
|
67
|
+
Requires-Dist: meson-python; extra == "build"
|
|
68
|
+
Requires-Dist: ninja; extra == "build"
|
|
69
|
+
Provides-Extra: lint
|
|
70
|
+
Requires-Dist: ruff; extra == "lint"
|
|
71
|
+
Provides-Extra: test
|
|
72
|
+
Requires-Dist: pytest-cov; extra == "test"
|
|
73
|
+
Provides-Extra: dev
|
|
74
|
+
Requires-Dist: wheel; extra == "dev"
|
|
75
|
+
Requires-Dist: twine; extra == "dev"
|
|
76
|
+
Description-Content-Type: text/x-rst
|
|
77
|
+
|
|
78
|
+
###########
|
|
79
|
+
scikit-misc
|
|
80
|
+
###########
|
|
81
|
+
|
|
82
|
+
================= =================
|
|
83
|
+
Latest Release |release|_
|
|
84
|
+
License |license|_
|
|
85
|
+
Build Status |buildstatus|_
|
|
86
|
+
Coverage |coverage|_
|
|
87
|
+
Documentation Documentation_
|
|
88
|
+
================= =================
|
|
89
|
+
|
|
90
|
+
Miscellaneous tools for data analysis and scientific computing.
|
|
91
|
+
|
|
92
|
+
Installation
|
|
93
|
+
============
|
|
94
|
+
|
|
95
|
+
**Official version**
|
|
96
|
+
|
|
97
|
+
.. code-block:: console
|
|
98
|
+
|
|
99
|
+
$ pip install scikit-misc
|
|
100
|
+
|
|
101
|
+
|
|
102
|
+
**Development version**
|
|
103
|
+
|
|
104
|
+
.. code-block:: console
|
|
105
|
+
|
|
106
|
+
$ pip install git+https://github.com/has2k1/scikit-misc.git@main
|
|
107
|
+
|
|
108
|
+
|
|
109
|
+
|
|
110
|
+
.. |release| image:: https://img.shields.io/pypi/v/scikit-misc.svg
|
|
111
|
+
.. _release: https://pypi.python.org/pypi/scikit-misc
|
|
112
|
+
|
|
113
|
+
.. |license| image:: https://img.shields.io/pypi/l/scikit-misc.svg
|
|
114
|
+
.. _license: https://pypi.python.org/pypi/scikit-misc
|
|
115
|
+
|
|
116
|
+
.. |buildstatus| image:: https://github.com/has2k1/scikit-misc/workflows/build/badge.svg?branch=main
|
|
117
|
+
.. _buildstatus: https://github.com/has2k1/scikit-misc/actions?query=branch%3Amain+workflow%3A%22build%22
|
|
118
|
+
|
|
119
|
+
.. |coverage| image:: https://codecov.io /github/has2k1/scikit-misc/coverage.svg?branch=main
|
|
120
|
+
.. _coverage: https://codecov.io/github/has2k1/scikit-misc?branch=main
|
|
121
|
+
|
|
122
|
+
.. _Documentation: https://has2k1.github.io/scikit-misc
|
|
@@ -0,0 +1,21 @@
|
|
|
1
|
+
scikit_misc-0.5.1.dev3-cp313-cp313-win_amd64.whl,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
|
|
2
|
+
skmisc/__config__.py,sha256=eSJrjRHQXxhjDEKc5mJNPKr1RjYGXoTN-p3YF7U4Aaw,3247
|
|
3
|
+
skmisc/__init__.py,sha256=otcK3Qg-nH3gkQPc2hV1K486pKceOZs2DGzXTJ3OR3w,1323
|
|
4
|
+
skmisc/_distributor_init.py,sha256=BaHRy7buKOTHXYEh3CgfZ09qQNcQJDeiT4L-cceDz94,353
|
|
5
|
+
skmisc/_version.py,sha256=qcwRrLYA7T_Gr9T0z6xrPXoscV3HiQmIhjZT_KgoqaI,154
|
|
6
|
+
skmisc/loess/__init__.py,sha256=usOeOL6T9H_Rki2s00Hov9mcz10ml5thdb1mnCAU_hw,2249
|
|
7
|
+
skmisc/loess/_loess.cp313-win_amd64.lib,sha256=bSMip566LMyOjSOKsG2s_cos7ieDleg4CZqVWYR7v48,1996
|
|
8
|
+
skmisc/loess/_loess.cp313-win_amd64.pyd,sha256=mJSwx9_7GsDYkoKZNv9GB8_J71wovPkyJP2zqXPY4Ek,279040
|
|
9
|
+
skmisc/loess/tests/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
|
|
10
|
+
skmisc/loess/tests/co2_data,sha256=HhEOebV0VQoMGz7uxRlQXSZ-xUxYvmmbDXA4SqgnJFU,2571
|
|
11
|
+
skmisc/loess/tests/co2_results_double,sha256=THSQkoOE2CoN1KgRllenB67YPbHqLzRzHyuNAonZ8Xw,50138
|
|
12
|
+
skmisc/loess/tests/co2_results_real,sha256=eI-tZWWEPYFmniG6qPTsAjvIuAJHdydoo-8AzwbAqCs,31346
|
|
13
|
+
skmisc/loess/tests/gas_data,sha256=nKGvpmhCuF01MXXzgJcCHjdpu9Tos6HB8nPhewaWarE,430
|
|
14
|
+
skmisc/loess/tests/gas_result,sha256=tmM_nkv8Dmq3LEyh79A86bVLtDxB7VWWYMIN5N9nZls,816
|
|
15
|
+
skmisc/loess/tests/madeup_data,sha256=yl319uCS2fVBEV67B2HuVUMOulIcCOAOkiZqLOOji18,2920
|
|
16
|
+
skmisc/loess/tests/madeup_result,sha256=7rNBbe8nQ3dM3nIaoT66Yc_jzdyCGpKid4-f67JFSY8,4467
|
|
17
|
+
skmisc/loess/tests/test_loess.py,sha256=sZk0w0TsU84yrmmlE9sDeQMZgfdHqFuwgKELZDj4p64,10537
|
|
18
|
+
scikit_misc-0.5.1.dev3.dist-info/LICENSE,sha256=mpejZQSZn17HBnYiQv7u1TLJCGEcuCONSjoLInSnIFU,1504
|
|
19
|
+
scikit_misc-0.5.1.dev3.dist-info/METADATA,sha256=O3HoiGFwL0p4kJDRY6kAKENLLNWTZh5O-kcUSdLH21k,4997
|
|
20
|
+
scikit_misc-0.5.1.dev3.dist-info/WHEEL,sha256=suq8ARrxbiI7iLH3BgK-82uzxQ-4Hm-m8w01oCokrtA,85
|
|
21
|
+
scikit_misc-0.5.1.dev3.dist-info/RECORD,,
|
skmisc/__config__.py
ADDED
|
@@ -0,0 +1,117 @@
|
|
|
1
|
+
# This file is generated by scikit-miscs's build process
|
|
2
|
+
# It contains system_info results at the time of building this package.
|
|
3
|
+
from enum import Enum
|
|
4
|
+
|
|
5
|
+
__all__ = ["show"]
|
|
6
|
+
_built_with_meson = True
|
|
7
|
+
|
|
8
|
+
|
|
9
|
+
class DisplayModes(Enum):
|
|
10
|
+
stdout = "stdout"
|
|
11
|
+
dicts = "dicts"
|
|
12
|
+
|
|
13
|
+
|
|
14
|
+
def _cleanup(d):
|
|
15
|
+
"""
|
|
16
|
+
Removes empty values in a `dict` recursively
|
|
17
|
+
This ensures we remove values that Meson could not provide to CONFIG
|
|
18
|
+
"""
|
|
19
|
+
if isinstance(d, dict):
|
|
20
|
+
return { k: _cleanup(v) for k, v in d.items() if v != '' and _cleanup(v) != '' }
|
|
21
|
+
else:
|
|
22
|
+
return d
|
|
23
|
+
|
|
24
|
+
|
|
25
|
+
CONFIG = _cleanup(
|
|
26
|
+
{
|
|
27
|
+
"Compilers": {
|
|
28
|
+
"c": {
|
|
29
|
+
"name": "msvc",
|
|
30
|
+
"linker": "link",
|
|
31
|
+
"version": "19.29.30158",
|
|
32
|
+
"commands": "cl",
|
|
33
|
+
},
|
|
34
|
+
"cython": {
|
|
35
|
+
"name": "cython",
|
|
36
|
+
"linker": "cython",
|
|
37
|
+
"version": "3.0.12",
|
|
38
|
+
"commands": "cython",
|
|
39
|
+
},
|
|
40
|
+
"fortran": {
|
|
41
|
+
"name": "gcc",
|
|
42
|
+
"linker": "ld.bfd",
|
|
43
|
+
"version": "10.3.0",
|
|
44
|
+
"commands": "gfortran",
|
|
45
|
+
},
|
|
46
|
+
},
|
|
47
|
+
"Build Dependencies": {
|
|
48
|
+
},
|
|
49
|
+
"Machine Information": {
|
|
50
|
+
"host": {
|
|
51
|
+
"cpu": "x86_64",
|
|
52
|
+
"family": "x86_64",
|
|
53
|
+
"endian": "little",
|
|
54
|
+
"system": "windows",
|
|
55
|
+
},
|
|
56
|
+
"build": {
|
|
57
|
+
"cpu": "x86_64",
|
|
58
|
+
"family": "x86_64",
|
|
59
|
+
"endian": "little",
|
|
60
|
+
"system": "windows",
|
|
61
|
+
},
|
|
62
|
+
"cross-compiled": bool("False".lower().replace('false', '')),
|
|
63
|
+
},
|
|
64
|
+
"Python Information": {
|
|
65
|
+
"path": r"C:\Users\runneradmin\AppData\Local\Temp\build-env-xnjrc587\Scripts\python.exe",
|
|
66
|
+
"version": "3.13",
|
|
67
|
+
},
|
|
68
|
+
}
|
|
69
|
+
)
|
|
70
|
+
|
|
71
|
+
|
|
72
|
+
def _check_pyyaml():
|
|
73
|
+
import yaml
|
|
74
|
+
|
|
75
|
+
return yaml
|
|
76
|
+
|
|
77
|
+
|
|
78
|
+
def show(mode=DisplayModes.stdout.value):
|
|
79
|
+
"""
|
|
80
|
+
Show libraries and system information on which scikit-misc was built
|
|
81
|
+
and is being used
|
|
82
|
+
|
|
83
|
+
Parameters
|
|
84
|
+
----------
|
|
85
|
+
mode : {`'stdout'`, `'dicts'`}, optional.
|
|
86
|
+
Indicates how to display the config information.
|
|
87
|
+
`'stdout'` prints to console, `'dicts'` returns a dictionary
|
|
88
|
+
of the configuration.
|
|
89
|
+
|
|
90
|
+
Returns
|
|
91
|
+
-------
|
|
92
|
+
out : {`dict`, `None`}
|
|
93
|
+
If mode is `'dicts'`, a dict is returned, else None
|
|
94
|
+
|
|
95
|
+
Notes
|
|
96
|
+
-----
|
|
97
|
+
1. The `'stdout'` mode will give more readable
|
|
98
|
+
output if ``pyyaml`` is installed
|
|
99
|
+
|
|
100
|
+
"""
|
|
101
|
+
if mode == DisplayModes.stdout.value:
|
|
102
|
+
try: # Non-standard library, check import
|
|
103
|
+
yaml = _check_pyyaml()
|
|
104
|
+
|
|
105
|
+
print(yaml.dump(CONFIG))
|
|
106
|
+
except ModuleNotFoundError:
|
|
107
|
+
import warnings
|
|
108
|
+
import json
|
|
109
|
+
|
|
110
|
+
warnings.warn("Install `pyyaml` for better output", stacklevel=1)
|
|
111
|
+
print(json.dumps(CONFIG, indent=2))
|
|
112
|
+
elif mode == DisplayModes.dicts.value:
|
|
113
|
+
return CONFIG
|
|
114
|
+
else:
|
|
115
|
+
raise AttributeError(
|
|
116
|
+
f"Invalid `mode`, use one of: {', '.join([e.value for e in DisplayModes])}"
|
|
117
|
+
)
|
skmisc/__init__.py
ADDED
|
@@ -0,0 +1,43 @@
|
|
|
1
|
+
from importlib.metadata import PackageNotFoundError, version
|
|
2
|
+
|
|
3
|
+
try:
|
|
4
|
+
__version__ = version('scikit_misc')
|
|
5
|
+
except PackageNotFoundError:
|
|
6
|
+
# package is not installed
|
|
7
|
+
pass
|
|
8
|
+
finally:
|
|
9
|
+
del version
|
|
10
|
+
del PackageNotFoundError
|
|
11
|
+
|
|
12
|
+
__all__ = ['__version__']
|
|
13
|
+
|
|
14
|
+
try:
|
|
15
|
+
from skmisc.__config__ import show as show_config # noqa: F401
|
|
16
|
+
except ImportError as err:
|
|
17
|
+
msg = """Error importing skmisc: you cannot import skmisc while
|
|
18
|
+
being in skmisc source directory; please exit the skmisc source
|
|
19
|
+
tree first, and relaunch your python intepreter."""
|
|
20
|
+
raise ImportError('\n\n'.join([err.message, msg])) # type: ignore
|
|
21
|
+
else:
|
|
22
|
+
__all__.append('show_config')
|
|
23
|
+
|
|
24
|
+
|
|
25
|
+
def test(args=None, plugins=None):
|
|
26
|
+
"""
|
|
27
|
+
Run tests
|
|
28
|
+
"""
|
|
29
|
+
from pathlib import Path
|
|
30
|
+
# The doctests are not run when called from an installed
|
|
31
|
+
# package since the pytest.ini is not included in the
|
|
32
|
+
# package.
|
|
33
|
+
try:
|
|
34
|
+
import pytest
|
|
35
|
+
except ImportError:
|
|
36
|
+
msg = "To run the tests, you must install pytest"
|
|
37
|
+
raise ImportError(msg)
|
|
38
|
+
path = str(Path(__file__).parent)
|
|
39
|
+
if args is None:
|
|
40
|
+
args = [path]
|
|
41
|
+
else:
|
|
42
|
+
args.append(path)
|
|
43
|
+
return pytest.main(args=args, plugins=plugins)
|
|
@@ -0,0 +1,10 @@
|
|
|
1
|
+
""" Distributor init file
|
|
2
|
+
|
|
3
|
+
Distributors: you can add custom code here to support particular distributions
|
|
4
|
+
of scikit-misc.
|
|
5
|
+
|
|
6
|
+
For example, this is a good place to put any checks for hardware requirements.
|
|
7
|
+
|
|
8
|
+
The scikit-misc standard source distribution will not put code in this file, so
|
|
9
|
+
you can safely replace this file with your own version.
|
|
10
|
+
"""
|
skmisc/_version.py
ADDED
skmisc/loess/__init__.py
ADDED
|
@@ -0,0 +1,64 @@
|
|
|
1
|
+
"""
|
|
2
|
+
=================================================
|
|
3
|
+
Locally-weighted regression (:mod:`skmisc.loess`)
|
|
4
|
+
=================================================
|
|
5
|
+
|
|
6
|
+
Loess is a procedure for estimating a regression surface by a multivariate
|
|
7
|
+
smoothing procedure. A linear or quadratic function of the independent
|
|
8
|
+
variables is fit in a moving fashion that is analogous to how a moving
|
|
9
|
+
average is computed for a time series.
|
|
10
|
+
|
|
11
|
+
Compared to approaches that fit global parametric functions, loess
|
|
12
|
+
substantially increases the domain of surfaces that can be estimated without
|
|
13
|
+
distortion. However, analogues of the statistical procedures used
|
|
14
|
+
in parametric function fitting -- for example, ANOVA and t intervals --
|
|
15
|
+
involve statistics whose distributions are well approximated by familiar
|
|
16
|
+
distributions.
|
|
17
|
+
|
|
18
|
+
.. autosummary::
|
|
19
|
+
:toctree: generated/
|
|
20
|
+
:template: cython_class.rst
|
|
21
|
+
|
|
22
|
+
loess
|
|
23
|
+
loess_inputs
|
|
24
|
+
loess_model
|
|
25
|
+
loess_control
|
|
26
|
+
loess_outputs
|
|
27
|
+
loess_prediction
|
|
28
|
+
loess_confidence_intervals
|
|
29
|
+
loess_anova
|
|
30
|
+
|
|
31
|
+
Source
|
|
32
|
+
------
|
|
33
|
+
The original source code was written by William S. Cleveland, Eric Grosse
|
|
34
|
+
and Ming-Jen Shyu. It is available at http://www.netlib.org/a/dloess. It
|
|
35
|
+
was initially adapted to for use in Scipy by Pierre GF Gerard-Marchant.
|
|
36
|
+
|
|
37
|
+
For more see references [1]_ [2]_ and [3]_.
|
|
38
|
+
|
|
39
|
+
.. [1] W. S. Cleveland, E. Grosse, and M. J. Shyu. Local Regression Models.
|
|
40
|
+
In J. M. Chambers and T. Hastie, editors, Statistical Models in S,
|
|
41
|
+
pages 309--376. Chapman and Hall, New York, 1992.
|
|
42
|
+
|
|
43
|
+
.. [2] W. S. Cleveland, S. J. Devlin, and E. Grosse. Regression by Local
|
|
44
|
+
Fitting: Methods, Properties, and Computing. Journal of Econometrics, 37:
|
|
45
|
+
pp. 87--114. 1988.
|
|
46
|
+
|
|
47
|
+
.. [3] W. S. Cleveland. Robust Locally Weighted Regression and Smoothing
|
|
48
|
+
Scatterplots. Journal of the American Statistical Association, 74:
|
|
49
|
+
pp. 829--836. 1979.
|
|
50
|
+
"""
|
|
51
|
+
from ._loess import (
|
|
52
|
+
loess,
|
|
53
|
+
loess_anova,
|
|
54
|
+
loess_confidence_intervals,
|
|
55
|
+
loess_control,
|
|
56
|
+
loess_inputs,
|
|
57
|
+
loess_model,
|
|
58
|
+
loess_outputs,
|
|
59
|
+
loess_prediction,
|
|
60
|
+
)
|
|
61
|
+
|
|
62
|
+
__all__ = ['loess', 'loess_model', 'loess_control', 'loess_inputs',
|
|
63
|
+
'loess_model', 'loess_outputs', 'loess_prediction',
|
|
64
|
+
'loess_confidence_intervals', 'loess_anova']
|
|
Binary file
|
|
Binary file
|
|
File without changes
|
|
@@ -0,0 +1,87 @@
|
|
|
1
|
+
315.58 316.39 316.79 317.82
|
|
2
|
+
318.39 318.22 316.68 315.01
|
|
3
|
+
314.02 313.55 315.02 315.75
|
|
4
|
+
316.52 317.1 317.79 319.22
|
|
5
|
+
320.08 319.7 318.27 315.99
|
|
6
|
+
314.24 314.05 315.05 316.23
|
|
7
|
+
316.92 317.76 318.54 319.49
|
|
8
|
+
320.64 319.85 318.7 316.96
|
|
9
|
+
315.17 315.47 316.19 317.17
|
|
10
|
+
318.12 318.72 319.79 320.68
|
|
11
|
+
321.28 320.89 319.79 317.56
|
|
12
|
+
316.46 315.59 316.85 317.87
|
|
13
|
+
318.87 319.25 320.13 321.49
|
|
14
|
+
322.34 321.62 319.85 317.87
|
|
15
|
+
316.36 316.24 317.13 318.46
|
|
16
|
+
319.57 320.23 320.89 321.54
|
|
17
|
+
322.2 321.9 320.42 318.6
|
|
18
|
+
316.73 317.15 317.94 318.91
|
|
19
|
+
319.73 320.78 321.23 322.49
|
|
20
|
+
322.59 322.35 321.61 319.24
|
|
21
|
+
318.23 317.76 319.36 319.5
|
|
22
|
+
320.35 321.4 322.22 323.45
|
|
23
|
+
323.8 323.5 322.16 320.09
|
|
24
|
+
318.26 317.66 319.47 320.7
|
|
25
|
+
322.06 322.23 322.78 324.1
|
|
26
|
+
324.63 323.79 322.34 320.73
|
|
27
|
+
319 318.99 320.41 321.68
|
|
28
|
+
322.3 322.89 323.59 324.65
|
|
29
|
+
325.3 325.15 323.88 321.8
|
|
30
|
+
319.99 319.86 320.88 322.36
|
|
31
|
+
323.59 324.23 325.34 326.33
|
|
32
|
+
327.03 326.24 325.39 323.16
|
|
33
|
+
321.87 321.31 322.34 323.74
|
|
34
|
+
324.61 325.58 326.55 327.81
|
|
35
|
+
327.82 327.53 326.29 324.66
|
|
36
|
+
323.12 323.09 324.01 325.1
|
|
37
|
+
326.12 326.62 327.16 327.94
|
|
38
|
+
329.15 328.79 327.53 325.65
|
|
39
|
+
323.6 323.78 325.13 326.26
|
|
40
|
+
326.93 327.84 327.96 329.93
|
|
41
|
+
330.25 329.24 328.13 326.42
|
|
42
|
+
324.97 325.29 326.56 327.73
|
|
43
|
+
328.73 329.7 330.46 331.7
|
|
44
|
+
332.66 332.22 331.02 329.39
|
|
45
|
+
327.58 327.27 328.3 328.81
|
|
46
|
+
329.44 330.89 331.62 332.85
|
|
47
|
+
333.29 332.44 331.35 329.58
|
|
48
|
+
327.58 327.55 328.56 329.73
|
|
49
|
+
330.45 330.98 331.63 332.88
|
|
50
|
+
333.63 333.53 331.9 330.08
|
|
51
|
+
328.59 328.31 329.44 330.64
|
|
52
|
+
331.62 332.45 333.36 334.46
|
|
53
|
+
334.84 334.29 333.04 330.88
|
|
54
|
+
329.23 328.83 330.18 331.5
|
|
55
|
+
332.8 333.22 334.54 335.82
|
|
56
|
+
336.45 335.97 334.65 332.4
|
|
57
|
+
331.28 330.73 332.05 333.54
|
|
58
|
+
334.65 335.06 336.32 337.39
|
|
59
|
+
337.66 337.56 336.24 334.39
|
|
60
|
+
332.43 332.22 333.61 334.78
|
|
61
|
+
335.88 336.43 337.61 338.53
|
|
62
|
+
339.06 338.92 337.39 335.72
|
|
63
|
+
333.64 333.65 335.07 336.53
|
|
64
|
+
337.82 338.19 339.89 340.56
|
|
65
|
+
341.22 340.92 339.26 337.27
|
|
66
|
+
335.66 335.54 336.71 337.79
|
|
67
|
+
338.79 340.06 340.93 342.02
|
|
68
|
+
342.65 341.8 340.01 337.94
|
|
69
|
+
336.17 336.28 337.76 339.05
|
|
70
|
+
340.18 341.04 342.16 343.01
|
|
71
|
+
343.64 342.91 341.72 339.52
|
|
72
|
+
337.75 337.68 339.14 340.37
|
|
73
|
+
341.32 342.45 343.05 344.91
|
|
74
|
+
345.77 345.3 343.98 342.41
|
|
75
|
+
339.89 340.03 341.19 342.87
|
|
76
|
+
343.74 344.55 345.28 347
|
|
77
|
+
347.37 346.74 345.36 343.19
|
|
78
|
+
340.97 341.2 342.76 343.96
|
|
79
|
+
344.82 345.82 347.24 348.09
|
|
80
|
+
348.66 347.9 346.27 344.21
|
|
81
|
+
342.88 342.58 343.99 345.31
|
|
82
|
+
345.98 346.72 347.63 349.24
|
|
83
|
+
349.83 349.1 347.52 345.43
|
|
84
|
+
344.48 343.89 345.29 346.54
|
|
85
|
+
347.66 348.07 349.12 350.55
|
|
86
|
+
351.34 350.8 349.1 347.54
|
|
87
|
+
346.2 346.2 347.44 348.67
|
|
@@ -0,0 +1,7 @@
|
|
|
1
|
+
-0.131254643360572 0.517069428110988 1.08182067584919 2.13326637358770 2.73375546874132 2.23997776835891 0.929699515256501 -1.11444707420409 -2.67573821245580 -2.87856330488736 -1.90209785853191 -0.960288396846844 -0.127309599506098 0.519123971628999 1.09845833595895 2.15377745092389 2.74178279831138 2.24390654627401 0.927034947716549 -1.11840574881049 -2.69421608102192 -2.89832584564760 -1.91246650303411 -0.960459432315943 -0.123323614572555 0.521218070743582 1.11513416618278 2.17432369186148 2.74984228497017 2.24786402897563 0.924395632660881 -1.12234432639115 -2.71267900802084 -2.91807776589179 -1.92282884807142 -0.960626468256293 -0.119335930046407 0.523313128694124 1.13181021448592 2.19487004452363 2.75790177699876 2.25182151436214 0.921756317605213 -1.12628287349056 -2.73114187405725 -2.93782932476804 -1.93319053133537 -0.960792260330621 -0.115346419561576 0.525409317286149 1.14848669811335 2.21541273077976 2.76595350089103 2.25576409079413 0.919094952776809 -1.13025072305144 -2.74964129524388 -2.95762143862511 -1.94359676941072 -0.961001887200727 -0.111663528494415 0.527041806941144 1.16591699949743 2.23572475782599 2.77558209490346 2.26100862164176 0.916135515605680 -1.13521851311631 -2.76853269649236 -2.97819389766126 -1.95447404199440 -0.961469922313078 -0.107860654017195 0.528789664565914 1.18345805341101 2.25613699562502 2.78530035789200 2.26632885561176 0.913237815703161 -1.14013859201476 -2.78739041267661 -2.99874290502888 -1.96533809630523 -0.961929689327584 -0.104054461617129 0.530539382987096 1.20099951099456 2.27654943992924 2.79501863022095 2.27164909425196 0.910340115800643 -1.14505866933551 -2.80624812570546 -3.01929207422557 -1.97620247742963 -0.962391060045031 -0.100251149809378 0.532282176944238 1.21853000024234 2.29694352384120 2.80471115010115 2.27693520990275 0.907399922368041 -1.15002684464570 -2.82515954118310 -3.03989522748933 -1.98712112423937 -0.962900372915200 -9.657142836469074E-002 0.534498365407295 1.23663732630217 2.31763175742248 2.81524537246355 2.28255222950951 0.904114205312311 -1.15622306722147 -2.84480538227602 -3.06114244925604 -1.99798435060416 -0.963119381070332 -9.277853340694056E-002 0.536818220007927 1.25483881112409 2.33840092971464 2.82584731348562 2.28822285505067 0.900867981465680 -1.16239113573849 -2.86443440846052 -3.08237920554798 -2.00884346092780 -0.963336241276855 -8.898545859310647E-002 0.539138117293329 1.27304020145945 2.35917005259441 2.83644925016947 2.29389347842271 0.897621757619049 -1.16855899932509 -2.88406302478419 -3.10361326189795 -2.01969758122834 -0.963539280350215 -8.516973153606598E-002 0.541497365991697 1.29129764237755 2.38001852337338 2.84715383206899 2.29969629210904 0.894537269185294 -1.17453457130906 -2.90346879331549 -3.12459976725524 -2.03027944733590 -0.963462846700087 -8.318137128004088E-002 0.543878489204072 1.30556339626621 2.40086196051427 2.86349649009386 2.30662073374580 0.887792197369512 -1.18505484118622 -2.92012438706208 -3.14478006706923 -2.03869306366509 -0.962605359946822 -8.131952325850622E-002 0.546127682709739 1.31969180297595 2.42157577340269 2.87971724679271 2.31343828536589 0.880955246846406 -1.19565053554695 -2.93683895106848 -3.16500456560729 -2.04713610718260 -0.961766931458097 -7.946636457773126E-002 0.548371989282220 1.33381912516008 2.44228903304709 2.89593798152913 2.32025582600477 0.874118296323300 -1.20624598994303 -2.95355303514557 -3.18522530648506 -2.05557211530886 -0.960907763364043 -7.757876207755750E-002 0.550679154240705 1.34803772029682 2.46313591470192 2.91233468733382 2.32730163521162 0.867561911867873 -1.21650906054413 -2.96988291770040 -3.20502042283162 -2.06354107589496 -0.959566493512336 -7.730219052307574E-002 0.550357935231823 1.35866963091892 2.48275195707344 2.93421989150165 2.33628089748756 0.856814006163955 -1.23047890110198 -2.98132949467161 -3.22331217432112 -2.07013703786751 -0.958710477392881 -7.712537096548630E-002 0.549926745420134 1.36918135193231 2.50225440094597 2.95599808828023 2.34516971366972 0.845992215661742 -1.24450607549849 -2.99281685452186 -3.24163247691575 -2.07674931917128 -0.957864064655987 -7.695143884179889E-002 0.549494111512313 1.37969307218733 2.52175684968569 2.97777629555155 2.35405853509826 0.835170425159529 -1.25853325141979 -3.00430421742168 -3.25995284087020 -2.08336172014513 -0.957018509640572 -7.677910249098863E-002 0.549056980809610 1.39019739462547 2.54124599510881 2.99953529400655 2.36292081968353 0.824314769787128 -1.27259985902399 -3.01583657881712 -3.27831926911595 -2.09002125120596 -0.956215270085997 -7.659034326116144E-002 0.548303305832745 1.40189817662095 2.56059404473287 3.02216449393333 2.37249493004545 0.812882095499132 -1.28748162695347 -3.02788580296873 -3.29720423678765 -2.09686072038062 -0.955524593494270 -7.628314748519592E-002 0.547660207647820 1.41370167565417 2.58003235634854 3.04487150080557 2.38213274135591 0.801499016162785 -1.30232615353693 -3.03991013937994 -3.31607228287846 -2.10369123413391 -0.954828406701153 -7.597388672782887E-002 0.547018261938728 1.42550541465766 2.59947079083601 3.06757851345117 2.39177055555305 0.790115936826438 -1.31717061821239 -3.05193435197517 -3.33493973126345 -2.11052067629153 -0.954130057562446 -7.566137287494026E-002 0.546380413810568 1.43731409572749 2.61891435943207 3.09029085224759 2.40141467672878 0.778740145296439 -1.33200480401653 -3.06394529463410 -3.35378911923031 -2.11732726754918 -0.953403895217910 -7.602451339366684E-002 0.546239431030685 1.44876402506477 2.63888506323516 3.11434437457380 2.41097132000011 0.766666985700136 -1.34824458585656 -3.07630256711981 -3.37247512859970 -2.12367610934446 -0.952323810467218 -7.635206537097822E-002 0.546127432756093 1.46023633487122 2.65887223673730 3.13840845582895 2.42053368947150 0.754594719575008 -1.36448616515093 -3.08866432798535 -3.39116569389385 -2.13002957460944 -0.951247011186959 -7.668156481946027E-002 0.546014308367871 1.47170833992158 2.67885925398630 3.16247252933391 2.43009605506780 0.742522453449881 -1.38072769771166 -3.10102599538363 -3.40985565045688 -2.13638191587943 -0.950167084869187 -7.700593418790123E-002 0.545910597459656 1.48319404185191 2.69886726639568 3.18656489627965 2.43969768454088 0.730500421637537 -1.39690705311059 -3.11331354277110 -3.42846020971871 -2.14263758255784 -0.948981962890397
|
|
2
|
+
315.551939725323 315.644481572339 315.737023419355 315.831019342046 315.925015264737 316.017871370006 316.110727475274 316.202019778800 316.293312082325 316.398957194298 316.504602306271 316.616296186514 316.727990066757 316.807925805088 316.887861543420 316.931141044518 316.974420545617 317.009107349720 317.043794153822 317.079929437763 317.116064721703 317.145706993755 317.175349265808 317.217323426682 317.259297587556 317.333560803835 317.407824020115 317.503177240388 317.598530460661 317.693263457235 317.787996453809 317.877822068571 317.967647683333 318.050567993525 318.133488303716 318.210189243932 318.286890184147 318.358269704290 318.429649224433 318.490571336492 318.551493448551 318.604564508711 318.657635568871 318.704397175047 318.751158781223 318.799682725764 318.848206670306 318.893484228360 318.938761786415 318.971518019974 319.004274253534 319.036061188978 319.067848124423 319.111286397250 319.154724670078 319.201719446166 319.248714222255 319.284772974631 319.320831727008 319.357168859719 319.393505992430 319.440776405357 319.488046818284 319.536228132110 319.584409445937 319.624345612089 319.664281778242 319.706243792973 319.748205807704 319.801823191447 319.855440575190 319.916163786446 319.976886997702 320.047187868720 320.117488739738 320.196045022062 320.274601304387 320.342956883518 320.411312462648 320.477813606563 320.544314750478 320.619464561199 320.694614371920 320.767508281146 320.840402190372 320.894032842931 320.947663495491 321.000026728028 321.052389960566 321.121797898971 321.191205837375 321.260278580623 321.329351323871 321.385910769388 321.442470214905 321.490940957730 321.539411700554 321.589024066796 321.638636433037 321.692886430859 321.747136428681 321.807410271409 321.867684114138 321.932598715622 321.997513317105 322.066533621897 322.135553926688 322.222083283889 322.308612641090 322.397868000686 322.487123360281 322.560303090986 322.633482821692 322.714287308585 322.795091795478 322.908260930562 323.021430065647 323.150533438384 323.279636811122 323.406616081944 323.533595352765 323.665482985268 323.797370617771 323.924810654343 324.052250690915 324.157171357403 324.262092023892 324.359777947874 324.457463871855 324.556286968541 324.655110065228 324.752623165413 324.850136265599 324.960022466940 325.069908668282 325.197416497441 325.324924326600 325.446490625083 325.568056923567 325.666280376633 325.764503829700 325.848864548032 325.933225266364 326.017071987071 326.100918707778 326.175668162728 326.250417617679 326.323939678717 326.397461739755 326.484110353133 326.570758966511 326.668917429334 326.767075892157 326.861808991717 326.956542091276 327.033335451432 327.110128811587 327.191444480573 327.272760149559 327.385077493485 327.497394837411 327.629431366801 327.761467896192 327.914194748308 328.066921600425 328.248121894823 328.429322189222 328.638637270929 328.847952352636 329.060521894823 329.273091437010 329.451251567944 329.629411698878 329.765631868200 329.901852037522 329.994517707214 330.087183376907 330.156606311731 330.226029246556 330.265626534325 330.305223822094 330.323366805201 330.341509788308 330.368464497079 330.395419205849 330.433097452678 330.470775699507 330.495675736625 330.520575773744 330.540523844711 330.560471915678 330.599888688797 330.639305461916 330.698128574083 330.756951686250 330.830811600114 330.904671513979 330.997038064963 331.089404615948 331.204868521107 331.320332426266 331.430668355679 331.541004285091 331.626799906690 331.712595528288 331.781599415280 331.850603302271 331.910793078141 331.970982854011 332.037661841837 332.104340829664 332.188600437080 332.272860044496 332.384631494404 332.496402944312 332.631585554655 332.766768164998 332.918152479767 333.069536794536 333.229342023138 333.389147251741 333.545866305444 333.702585359147 333.848586388998 333.994587418849 334.126660196019 334.258732973190 334.387671725906 334.516610478623 334.648009639758 334.779408800893 334.902725346625 335.026041892356 335.139417944176 335.252793995995 335.361458254357 335.470122512718 335.573790214812 335.677457916906 335.779511353252 335.881564789598 335.989359512658 336.097154235719 336.218065610003 336.338976984288 336.476357256642 336.613737528996 336.769402814510 336.925068100023 337.091013012734 337.256957925444 337.420581414309 337.584204903175 337.740545944728 337.896886986282 338.042603733218 338.188320480154 338.310523251021 338.432726021888 338.544130098184 338.655534174480 338.765442959544 338.875351744608 338.966731194332 339.058110644055 339.128072537222 339.198034430388 339.268379466158 339.338724501929 339.424039477127 339.509354452325 339.599580850248 339.689807248171 339.786401591064 339.882995933956 339.997126144263 340.111256354571 340.235990611348 340.360724868124 340.473730387176 340.586735906228 340.689158181393 340.791580456557 340.900727360221 341.009874263885 341.154872396777 341.299870529668 341.483472761808 341.667074993948 341.868242153109 342.069409312271 342.264456214578 342.459503116885 342.650449487265 342.841395857645 343.021253454262 343.201111050878 343.354468577260 343.507826103641 343.637191865892 343.766557628142 343.886350911358 344.006144194574 344.115057152717 344.223970110859 344.323751742253 344.423533373646 344.526471089294 344.629408804942 344.735184447806 344.840960090669 344.940830075417 345.040700060165 345.142825602093 345.244951144021 345.353258414239 345.461565684457 345.561742598325 345.661919512192 345.747183300792 345.832447089391 345.918321519461 346.004195949531 346.100779004480 346.197362059429 346.299593561791 346.401825064153 346.512186102756 346.622547141358 346.744061652652 346.865576163946 346.988500888955 347.111425613964 347.232657355122 347.353889096279 347.483494160462 347.613099224645 347.764657761497 347.916216298350 348.111704887568 348.307193476786 348.500459990222 348.693726503657 348.890983069679 349.088239635700 349.291727004359 349.495214373018 349.702691602918
|
|
3
|
+
0.947488031764416 0.882390707609878 0.998157647188866 0.952418608947733 0.817270789273660 0.996885051409493 0.728085695102151 0.990851982305767 0.637884911652333 0.990072880709515 0.595762565698863 0.972711816092563 0.965856676712861 0.841499815568552 0.886523232980042 0.965146976116641 0.741668277822855 0.591837991311972 0.792193970670084 0.997641002602193 0.912710251366073 0.928807209498480 0.898237850232544 0.999266296612623 0.851886377313259 0.964215762668314 0.999965289268365 0.908641401475467 0.823260760767501 0.977611315847939 0.999974349583762 0.889567351166186 0.981488117858243 0.692616403426253 0.999799300554367 0.989320449671225 0.986825156618664 0.919534533638257 0.884225650135546 0.999760151806573 0.994414971601326 0.997951916161087 0.887100740140224 0.999562269234954 0.585920531538786 0.853141175940922 0.991037805096357 0.992152573567451 0.999332011617492 0.823545047162815 0.996443576352900 0.874206648004056 0.504419990571978 0.857988673199827 0.880110420364086 0.900473538757466 0.945823442689414 0.988574826041659 0.857348906652125 0.987735219949253 0.834849581142526 0.860603669548175 0.876631518580033 0.867738652432744 0.925070932274631 0.999566920332460 0.942733182253283 0.997254316870668 0.848459702866706 0.718212115789131 0.994416519154735 0.997261147415274 0.940630243662637 0.914072981316579 0.985177959596764 0.996989902230292 0.495901708774148 0.831373909694944 0.811585411032084 0.972422902327099 0.546506568978507 0.945466147984970 0.265167252870300 0.780073398100054 0.624172807601586 0.994702160557937 0.990471279156271 0.926205743368332 0.992289442371609 0.968178822977318 0.986639627303399 0.999350711154713 0.840996129164800 0.182032780051221 0.999420135756517 0.915674903525654 0.282483995465186 0.974733431424663 0.985662777337056 0.967728208260297 0.988995028275044 0.798827379975444 0.595406669550792 0.993177956795018 0.926237360222123 0.999240100281311 0.831292225438565 0.594676637874432 0.984614780462075 0.993177780480680 0.955256473058929 0.879741440104462 0.939025147345531 0.939963934618185 0.914610953158321 0.994504016810673 0.913917505780136 0.891820062027726 0.646946420306272 0.986077202293310 0.949554346292742 0.999015136479037 0.804778629910955 0.987191788882639 0.949946668384419 0.900620180118301 0.872239965517043 0.995495686190679 0.824647797345363 0.944362118100568 0.781994122827541 0.995249330421569 0.940627092330763 0.986553210751971 0.895606345391544 0.839413763311823 0.721710696426580 0.897831775256494 0.929067967094475 0.938208550669446 0.861961823764060 0.696886207592241 0.971019508447038 0.989448569084288 0.965598560947786 0.981922340009845 0.701411404420612 0.104600308707135 0.983164500996082 0.998017240668145 0.982302668213245 0.935061774348980 0.846669208371132 0.988816946876684 0.886186572312611 0.891383147444920 0.985965748104140 0.961108151167829 0.327083058925475 0.917596615233765 0.984060890087786 0.217730716921630 0.508169585980581 0.804407334152152 0.941039626571956 0.889349173772447 0.920662880643609 0.989464185236220 0.997605465838400 0.979297066753547 0.957172016172363 0.932132506038552 0.949575189698423 0.946393910886845 0.854301739303652 0.324800555892071 0.583488485913683 0.793620289936204 0.957643066851486 0.515647886050300 7.587988206706671E-002 0.996867521589717 0.997314974002783 0.989598664164471 0.999803995476453 0.791306381858096 0.999556931659123 0.812664089299613 0.999670698308111 0.900955932641992 0.992430774756580 0.976684805319769 0.975935960902469 0.839419982497541 0.541881541911152 0.650753760709872 0.920590438442358 0.885607859732216 0.993861723237700 0.983788524221900 0.863800866013731 0.980691997132225 0.995615856317590 0.998035099067158 0.999811714758427 0.964988924043206 0.941495533701261 0.985448935005877 0.989091806279382 0.984726628049543 0.983865840183017 0.978314689920389 0.989713063012797 0.748690231413388 0.850392645437634 0.926729817273933 0.964055523221987 0.858174179764269 0.971942532484567 0.975564067037686 0.973059960674606 0.982890313458071 0.974544248477077 0.880302225261198 0.823633420842541 0.939346535967996 0.952335555793629 0.972133528480918 0.882069362488227 0.957244691296426 0.931996655132657 0.998759581996794 0.730520403082927 0.987445088111139 0.943077947728455 0.830796266171339 0.995962582309898 0.976397067467406 0.999842741306551 0.994396191911193 0.980866711135284 0.974783154853931 0.958741105591158 0.887673550805615 0.821573436078837 0.981655775478357 0.994525598973073 0.902253799209502 0.820031864288388 0.927527981444656 0.969506633058071 0.993067476133660 0.762049414169121 0.976744402242804 0.321074631755019 0.996577589308699 0.999892202223984 0.867305592691854 0.999769688520615 0.999590503549087 0.998317218746979 0.990283611839627 0.984615565457175 0.873629200967010 0.925085621757581 0.649935802425659 0.752761503409992 0.915204684728773 0.822061098076671 0.999993666123134 0.794484496418216 0.679678268772147 0.512279128327064 0.907667460096379 0.998334986500622 0.999989606657796 0.932594465060915 0.824912191750586 0.668096805665421 0.994350905307304 0.999882062973145 0.932293050204642 0.953093150199542 0.980705662019545 0.892601959093755 0.945217594123708 0.991556051414152 0.930916373128400 0.850512922777617 0.994452949956122 0.554705372831289 0.997552734207715 0.889705207453044 0.867636913380303 0.733054431000601 0.189343181073062 0.826890203083248 0.999990059425219 0.873347393089944 0.937402781828196 0.991024201719805 0.960198320684717 0.931895807744739 0.836099661156904 0.994896965661127 0.999555521329668 0.925705142591297 0.999705976575668 0.337357688658303 0.927201512426989 0.997233493771151 0.996180780623039 0.965115513480998 0.945214504620896 0.400652424056755 0.978910246416990 0.988197843822846 0.983488981782400 0.944786229818726 0.898731539524064 0.958107295811193 0.997020194126272 0.975042358663448 0.949058616986090 0.966457179198796 0.971509076915235 0.866937202889720 0.990404613199697 0.992862121865960 0.985910200969520 0.978506742931891 0.892033009545154 0.536045358478198 0.994294000448681 0.991893801076137 0.999776834875325 0.943669866205371 0.887873440548549 0.835155315552832 0.865672193714240 0.957628043010181 0.958717440718235 0.765411016459182 0.998796714540848 0.872204491067018 0.752159100110003 0.982845533894251 0.980718422156039
|
|
4
|
+
2
|
|
5
|
+
-0.131254643360572 0.517069428110988 1.08182067584919 2.13326637358770 2.73375546874132 2.23997776835891 0.929699515256501 -1.11444707420409 -2.67573821245580 -2.87856330488736 -1.90209785853191 -0.960288396846844 -0.127309599506098 0.519123971628999 1.09845833595895 2.15377745092389 2.74178279831138 2.24390654627401 0.927034947716549 -1.11840574881049 -2.69421608102192 -2.89832584564760 -1.91246650303411 -0.960459432315943 -0.123323614572555 0.521218070743582 1.11513416618278 2.17432369186148 2.74984228497017 2.24786402897563 0.924395632660881 -1.12234432639115 -2.71267900802084 -2.91807776589179 -1.92282884807142 -0.960626468256293 -0.119335930046407 0.523313128694124 1.13181021448592 2.19487004452363 2.75790177699876 2.25182151436214 0.921756317605213 -1.12628287349056 -2.73114187405725 -2.93782932476804 -1.93319053133537 -0.960792260330621 -0.115346419561576 0.525409317286149 1.14848669811335 2.21541273077976 2.76595350089103 2.25576409079413 0.919094952776809 -1.13025072305144 -2.74964129524388 -2.95762143862511 -1.94359676941072 -0.961001887200727 -0.111663528494415 0.527041806941144 1.16591699949743 2.23572475782599 2.77558209490346 2.26100862164176 0.916135515605680 -1.13521851311631 -2.76853269649236 -2.97819389766126 -1.95447404199440 -0.961469922313078 -0.107860654017195 0.528789664565914 1.18345805341101 2.25613699562502 2.78530035789200 2.26632885561176 0.913237815703161 -1.14013859201476 -2.78739041267661 -2.99874290502888 -1.96533809630523 -0.961929689327584 -0.104054461617129 0.530539382987096 1.20099951099456 2.27654943992924 2.79501863022095 2.27164909425196 0.910340115800643 -1.14505866933551 -2.80624812570546 -3.01929207422557 -1.97620247742963 -0.962391060045031 -0.100251149809378 0.532282176944238 1.21853000024234 2.29694352384120 2.80471115010115 2.27693520990275 0.907399922368041 -1.15002684464570 -2.82515954118310 -3.03989522748933 -1.98712112423937 -0.962900372915200 -9.657142836469074E-002 0.534498365407295 1.23663732630217 2.31763175742248 2.81524537246355 2.28255222950951 0.904114205312311 -1.15622306722147 -2.84480538227602 -3.06114244925604 -1.99798435060416 -0.963119381070332 -9.277853340694056E-002 0.536818220007927 1.25483881112409 2.33840092971464 2.82584731348562 2.28822285505067 0.900867981465680 -1.16239113573849 -2.86443440846052 -3.08237920554798 -2.00884346092780 -0.963336241276855 -8.898545859310647E-002 0.539138117293329 1.27304020145945 2.35917005259441 2.83644925016947 2.29389347842271 0.897621757619049 -1.16855899932509 -2.88406302478419 -3.10361326189795 -2.01969758122834 -0.963539280350215 -8.516973153606598E-002 0.541497365991697 1.29129764237755 2.38001852337338 2.84715383206899 2.29969629210904 0.894537269185294 -1.17453457130906 -2.90346879331549 -3.12459976725524 -2.03027944733590 -0.963462846700087 -8.318137128004088E-002 0.543878489204072 1.30556339626621 2.40086196051427 2.86349649009386 2.30662073374580 0.887792197369512 -1.18505484118622 -2.92012438706208 -3.14478006706923 -2.03869306366509 -0.962605359946822 -8.131952325850622E-002 0.546127682709739 1.31969180297595 2.42157577340269 2.87971724679271 2.31343828536589 0.880955246846406 -1.19565053554695 -2.93683895106848 -3.16500456560729 -2.04713610718260 -0.961766931458097 -7.946636457773126E-002 0.548371989282220 1.33381912516008 2.44228903304709 2.89593798152913 2.32025582600477 0.874118296323300 -1.20624598994303 -2.95355303514557 -3.18522530648506 -2.05557211530886 -0.960907763364043 -7.757876207755750E-002 0.550679154240705 1.34803772029682 2.46313591470192 2.91233468733382 2.32730163521162 0.867561911867873 -1.21650906054413 -2.96988291770040 -3.20502042283162 -2.06354107589496 -0.959566493512336 -7.730219052307574E-002 0.550357935231823 1.35866963091892 2.48275195707344 2.93421989150165 2.33628089748756 0.856814006163955 -1.23047890110198 -2.98132949467161 -3.22331217432112 -2.07013703786751 -0.958710477392881 -7.712537096548630E-002 0.549926745420134 1.36918135193231 2.50225440094597 2.95599808828023 2.34516971366972 0.845992215661742 -1.24450607549849 -2.99281685452186 -3.24163247691575 -2.07674931917128 -0.957864064655987 -7.695143884179889E-002 0.549494111512313 1.37969307218733 2.52175684968569 2.97777629555155 2.35405853509826 0.835170425159529 -1.25853325141979 -3.00430421742168 -3.25995284087020 -2.08336172014513 -0.957018509640572 -7.677910249098863E-002 0.549056980809610 1.39019739462547 2.54124599510881 2.99953529400655 2.36292081968353 0.824314769787128 -1.27259985902399 -3.01583657881712 -3.27831926911595 -2.09002125120596 -0.956215270085997 -7.659034326116144E-002 0.548303305832745 1.40189817662095 2.56059404473287 3.02216449393333 2.37249493004545 0.812882095499132 -1.28748162695347 -3.02788580296873 -3.29720423678765 -2.09686072038062 -0.955524593494270 -7.628314748519592E-002 0.547660207647820 1.41370167565417 2.58003235634854 3.04487150080557 2.38213274135591 0.801499016162785 -1.30232615353693 -3.03991013937994 -3.31607228287846 -2.10369123413391 -0.954828406701153 -7.597388672782887E-002 0.547018261938728 1.42550541465766 2.59947079083601 3.06757851345117 2.39177055555305 0.790115936826438 -1.31717061821239 -3.05193435197517 -3.33493973126345 -2.11052067629153 -0.954130057562446 -7.566137287494026E-002 0.546380413810568 1.43731409572749 2.61891435943207 3.09029085224759 2.40141467672878 0.778740145296439 -1.33200480401653 -3.06394529463410 -3.35378911923031 -2.11732726754918 -0.953403895217910 -7.602451339366684E-002 0.546239431030685 1.44876402506477 2.63888506323516 3.11434437457380 2.41097132000011 0.766666985700136 -1.34824458585656 -3.07630256711981 -3.37247512859970 -2.12367610934446 -0.952323810467218 -7.635206537097822E-002 0.546127432756093 1.46023633487122 2.65887223673730 3.13840845582895 2.42053368947150 0.754594719575008 -1.36448616515093 -3.08866432798535 -3.39116569389385 -2.13002957460944 -0.951247011186959 -7.668156481946027E-002 0.546014308367871 1.47170833992158 2.67885925398630 3.16247252933391 2.43009605506780 0.742522453449881 -1.38072769771166 -3.10102599538363 -3.40985565045688 -2.13638191587943 -0.950167084869187 -7.700593418790123E-002 0.545910597459656 1.48319404185191 2.69886726639568 3.18656489627965 2.43969768454088 0.730500421637537 -1.39690705311059 -3.11331354277110 -3.42846020971871 -2.14263758255784 -0.948981962890397
|
|
6
|
+
315.551939725323 315.644481572339 315.737023419355 315.831019342046 315.925015264737 316.017871370006 316.110727475274 316.202019778800 316.293312082325 316.398957194298 316.504602306271 316.616296186514 316.727990066757 316.807925805088 316.887861543420 316.931141044518 316.974420545617 317.009107349720 317.043794153822 317.079929437763 317.116064721703 317.145706993755 317.175349265808 317.217323426682 317.259297587556 317.333560803835 317.407824020115 317.503177240388 317.598530460661 317.693263457235 317.787996453809 317.877822068571 317.967647683333 318.050567993525 318.133488303716 318.210189243932 318.286890184147 318.358269704290 318.429649224433 318.490571336492 318.551493448551 318.604564508711 318.657635568871 318.704397175047 318.751158781223 318.799682725764 318.848206670306 318.893484228360 318.938761786415 318.971518019974 319.004274253534 319.036061188978 319.067848124423 319.111286397250 319.154724670078 319.201719446166 319.248714222255 319.284772974631 319.320831727008 319.357168859719 319.393505992430 319.440776405357 319.488046818284 319.536228132110 319.584409445937 319.624345612089 319.664281778242 319.706243792973 319.748205807704 319.801823191447 319.855440575190 319.916163786446 319.976886997702 320.047187868720 320.117488739738 320.196045022062 320.274601304387 320.342956883518 320.411312462648 320.477813606563 320.544314750478 320.619464561199 320.694614371920 320.767508281146 320.840402190372 320.894032842931 320.947663495491 321.000026728028 321.052389960566 321.121797898971 321.191205837375 321.260278580623 321.329351323871 321.385910769388 321.442470214905 321.490940957730 321.539411700554 321.589024066796 321.638636433037 321.692886430859 321.747136428681 321.807410271409 321.867684114138 321.932598715622 321.997513317105 322.066533621897 322.135553926688 322.222083283889 322.308612641090 322.397868000686 322.487123360281 322.560303090986 322.633482821692 322.714287308585 322.795091795478 322.908260930562 323.021430065647 323.150533438384 323.279636811122 323.406616081944 323.533595352765 323.665482985268 323.797370617771 323.924810654343 324.052250690915 324.157171357403 324.262092023892 324.359777947874 324.457463871855 324.556286968541 324.655110065228 324.752623165413 324.850136265599 324.960022466940 325.069908668282 325.197416497441 325.324924326600 325.446490625083 325.568056923567 325.666280376633 325.764503829700 325.848864548032 325.933225266364 326.017071987071 326.100918707778 326.175668162728 326.250417617679 326.323939678717 326.397461739755 326.484110353133 326.570758966511 326.668917429334 326.767075892157 326.861808991717 326.956542091276 327.033335451432 327.110128811587 327.191444480573 327.272760149559 327.385077493485 327.497394837411 327.629431366801 327.761467896192 327.914194748308 328.066921600425 328.248121894823 328.429322189222 328.638637270929 328.847952352636 329.060521894823 329.273091437010 329.451251567944 329.629411698878 329.765631868200 329.901852037522 329.994517707214 330.087183376907 330.156606311731 330.226029246556 330.265626534325 330.305223822094 330.323366805201 330.341509788308 330.368464497079 330.395419205849 330.433097452678 330.470775699507 330.495675736625 330.520575773744 330.540523844711 330.560471915678 330.599888688797 330.639305461916 330.698128574083 330.756951686250 330.830811600114 330.904671513979 330.997038064963 331.089404615948 331.204868521107 331.320332426266 331.430668355679 331.541004285091 331.626799906690 331.712595528288 331.781599415280 331.850603302271 331.910793078141 331.970982854011 332.037661841837 332.104340829664 332.188600437080 332.272860044496 332.384631494404 332.496402944312 332.631585554655 332.766768164998 332.918152479767 333.069536794536 333.229342023138 333.389147251741 333.545866305444 333.702585359147 333.848586388998 333.994587418849 334.126660196019 334.258732973190 334.387671725906 334.516610478623 334.648009639758 334.779408800893 334.902725346625 335.026041892356 335.139417944176 335.252793995995 335.361458254357 335.470122512718 335.573790214812 335.677457916906 335.779511353252 335.881564789598 335.989359512658 336.097154235719 336.218065610003 336.338976984288 336.476357256642 336.613737528996 336.769402814510 336.925068100023 337.091013012734 337.256957925444 337.420581414309 337.584204903175 337.740545944728 337.896886986282 338.042603733218 338.188320480154 338.310523251021 338.432726021888 338.544130098184 338.655534174480 338.765442959544 338.875351744608 338.966731194332 339.058110644055 339.128072537222 339.198034430388 339.268379466158 339.338724501929 339.424039477127 339.509354452325 339.599580850248 339.689807248171 339.786401591064 339.882995933956 339.997126144263 340.111256354571 340.235990611348 340.360724868124 340.473730387176 340.586735906228 340.689158181393 340.791580456557 340.900727360221 341.009874263885 341.154872396777 341.299870529668 341.483472761808 341.667074993948 341.868242153109 342.069409312271 342.264456214578 342.459503116885 342.650449487265 342.841395857645 343.021253454262 343.201111050878 343.354468577260 343.507826103641 343.637191865892 343.766557628142 343.886350911358 344.006144194574 344.115057152717 344.223970110859 344.323751742253 344.423533373646 344.526471089294 344.629408804942 344.735184447806 344.840960090669 344.940830075417 345.040700060165 345.142825602093 345.244951144021 345.353258414239 345.461565684457 345.561742598325 345.661919512192 345.747183300792 345.832447089391 345.918321519461 346.004195949531 346.100779004480 346.197362059429 346.299593561791 346.401825064153 346.512186102756 346.622547141358 346.744061652652 346.865576163946 346.988500888955 347.111425613964 347.232657355122 347.353889096279 347.483494160462 347.613099224645 347.764657761497 347.916216298350 348.111704887568 348.307193476786 348.500459990222 348.693726503657 348.890983069679 349.088239635700 349.291727004359 349.495214373018 349.702691602918
|
|
7
|
+
0.947488031764416 0.882390707609878 0.998157647188866 0.952418608947733 0.817270789273660 0.996885051409493 0.728085695102151 0.990851982305767 0.637884911652333 0.990072880709515 0.595762565698863 0.972711816092563 0.965856676712861 0.841499815568552 0.886523232980042 0.965146976116641 0.741668277822855 0.591837991311972 0.792193970670084 0.997641002602193 0.912710251366073 0.928807209498480 0.898237850232544 0.999266296612623 0.851886377313259 0.964215762668314 0.999965289268365 0.908641401475467 0.823260760767501 0.977611315847939 0.999974349583762 0.889567351166186 0.981488117858243 0.692616403426253 0.999799300554367 0.989320449671225 0.986825156618664 0.919534533638257 0.884225650135546 0.999760151806573 0.994414971601326 0.997951916161087 0.887100740140224 0.999562269234954 0.585920531538786 0.853141175940922 0.991037805096357 0.992152573567451 0.999332011617492 0.823545047162815 0.996443576352900 0.874206648004056 0.504419990571978 0.857988673199827 0.880110420364086 0.900473538757466 0.945823442689414 0.988574826041659 0.857348906652125 0.987735219949253 0.834849581142526 0.860603669548175 0.876631518580033 0.867738652432744 0.925070932274631 0.999566920332460 0.942733182253283 0.997254316870668 0.848459702866706 0.718212115789131 0.994416519154735 0.997261147415274 0.940630243662637 0.914072981316579 0.985177959596764 0.996989902230292 0.495901708774148 0.831373909694944 0.811585411032084 0.972422902327099 0.546506568978507 0.945466147984970 0.265167252870300 0.780073398100054 0.624172807601586 0.994702160557937 0.990471279156271 0.926205743368332 0.992289442371609 0.968178822977318 0.986639627303399 0.999350711154713 0.840996129164800 0.182032780051221 0.999420135756517 0.915674903525654 0.282483995465186 0.974733431424663 0.985662777337056 0.967728208260297 0.988995028275044 0.798827379975444 0.595406669550792 0.993177956795018 0.926237360222123 0.999240100281311 0.831292225438565 0.594676637874432 0.984614780462075 0.993177780480680 0.955256473058929 0.879741440104462 0.939025147345531 0.939963934618185 0.914610953158321 0.994504016810673 0.913917505780136 0.891820062027726 0.646946420306272 0.986077202293310 0.949554346292742 0.999015136479037 0.804778629910955 0.987191788882639 0.949946668384419 0.900620180118301 0.872239965517043 0.995495686190679 0.824647797345363 0.944362118100568 0.781994122827541 0.995249330421569 0.940627092330763 0.986553210751971 0.895606345391544 0.839413763311823 0.721710696426580 0.897831775256494 0.929067967094475 0.938208550669446 0.861961823764060 0.696886207592241 0.971019508447038 0.989448569084288 0.965598560947786 0.981922340009845 0.701411404420612 0.104600308707135 0.983164500996082 0.998017240668145 0.982302668213245 0.935061774348980 0.846669208371132 0.988816946876684 0.886186572312611 0.891383147444920 0.985965748104140 0.961108151167829 0.327083058925475 0.917596615233765 0.984060890087786 0.217730716921630 0.508169585980581 0.804407334152152 0.941039626571956 0.889349173772447 0.920662880643609 0.989464185236220 0.997605465838400 0.979297066753547 0.957172016172363 0.932132506038552 0.949575189698423 0.946393910886845 0.854301739303652 0.324800555892071 0.583488485913683 0.793620289936204 0.957643066851486 0.515647886050300 7.587988206706671E-002 0.996867521589717 0.997314974002783 0.989598664164471 0.999803995476453 0.791306381858096 0.999556931659123 0.812664089299613 0.999670698308111 0.900955932641992 0.992430774756580 0.976684805319769 0.975935960902469 0.839419982497541 0.541881541911152 0.650753760709872 0.920590438442358 0.885607859732216 0.993861723237700 0.983788524221900 0.863800866013731 0.980691997132225 0.995615856317590 0.998035099067158 0.999811714758427 0.964988924043206 0.941495533701261 0.985448935005877 0.989091806279382 0.984726628049543 0.983865840183017 0.978314689920389 0.989713063012797 0.748690231413388 0.850392645437634 0.926729817273933 0.964055523221987 0.858174179764269 0.971942532484567 0.975564067037686 0.973059960674606 0.982890313458071 0.974544248477077 0.880302225261198 0.823633420842541 0.939346535967996 0.952335555793629 0.972133528480918 0.882069362488227 0.957244691296426 0.931996655132657 0.998759581996794 0.730520403082927 0.987445088111139 0.943077947728455 0.830796266171339 0.995962582309898 0.976397067467406 0.999842741306551 0.994396191911193 0.980866711135284 0.974783154853931 0.958741105591158 0.887673550805615 0.821573436078837 0.981655775478357 0.994525598973073 0.902253799209502 0.820031864288388 0.927527981444656 0.969506633058071 0.993067476133660 0.762049414169121 0.976744402242804 0.321074631755019 0.996577589308699 0.999892202223984 0.867305592691854 0.999769688520615 0.999590503549087 0.998317218746979 0.990283611839627 0.984615565457175 0.873629200967010 0.925085621757581 0.649935802425659 0.752761503409992 0.915204684728773 0.822061098076671 0.999993666123134 0.794484496418216 0.679678268772147 0.512279128327064 0.907667460096379 0.998334986500622 0.999989606657796 0.932594465060915 0.824912191750586 0.668096805665421 0.994350905307304 0.999882062973145 0.932293050204642 0.953093150199542 0.980705662019545 0.892601959093755 0.945217594123708 0.991556051414152 0.930916373128400 0.850512922777617 0.994452949956122 0.554705372831289 0.997552734207715 0.889705207453044 0.867636913380303 0.733054431000601 0.189343181073062 0.826890203083248 0.999990059425219 0.873347393089944 0.937402781828196 0.991024201719805 0.960198320684717 0.931895807744739 0.836099661156904 0.994896965661127 0.999555521329668 0.925705142591297 0.999705976575668 0.337357688658303 0.927201512426989 0.997233493771151 0.996180780623039 0.965115513480998 0.945214504620896 0.400652424056755 0.978910246416990 0.988197843822846 0.983488981782400 0.944786229818726 0.898731539524064 0.958107295811193 0.997020194126272 0.975042358663448 0.949058616986090 0.966457179198796 0.971509076915235 0.866937202889720 0.990404613199697 0.992862121865960 0.985910200969520 0.978506742931891 0.892033009545154 0.536045358478198 0.994294000448681 0.991893801076137 0.999776834875325 0.943669866205371 0.887873440548549 0.835155315552832 0.865672193714240 0.957628043010181 0.958717440718235 0.765411016459182 0.998796714540848 0.872204491067018 0.752159100110003 0.982845533894251 0.980718422156039
|
|
@@ -0,0 +1,7 @@
|
|
|
1
|
+
-0.1312545 0.5171010 1.081863 2.133283 2.733717 2.239956 0.9297554 -1.114420 -2.675813 -2.878533 -1.902167 -0.9602861 -0.1273115 0.5191532 1.098498 2.153796 2.741747 2.243891 0.9270887 -1.118383 -2.694288 -2.898298 -1.912531 -0.9604589 -0.1233273 0.5212451 1.115170 2.174344 2.749810 2.247854 0.9244474 -1.122326 -2.712748 -2.918053 -1.922890 -0.9606276 -0.1193415 0.5233380 1.131844 2.194892 2.757873 2.251818 0.9218062 -1.126269 -2.731208 -2.937807 -1.933248 -0.9607950 -0.1153538 0.5254319 1.148517 2.215437 2.765927 2.255767 0.9191428 -1.130242 -2.749705 -2.957601 -1.943650 -0.9610062 -0.1116722 0.5270622 1.165944 2.235751 2.775559 2.261016 0.9161811 -1.135213 -2.768593 -2.978176 -1.954524 -0.9614759 -0.1078705 0.5288079 1.183482 2.256165 2.785280 2.266341 0.9132813 -1.140137 -2.787447 -2.998728 -1.965384 -0.9619374 -0.1040655 0.5305554 1.201020 2.276579 2.795002 2.271666 0.9103814 -1.145061 -2.806301 -3.019279 -1.976244 -0.9624004 -0.1002634 0.5322961 1.218547 2.296975 2.804697 2.276956 0.9074391 -1.150033 -2.825208 -3.039885 -1.987159 -0.9629111 -9.6585110E-02 0.5345098 1.236650 2.317664 2.815235 2.282577 0.9041511 -1.156232 -2.844850 -3.061133 -1.998018 -0.9631324 -9.2793740E-02 0.5368271 1.254847 2.338434 2.825840 2.288252 0.9009027 -1.162404 -2.864474 -3.082371 -2.008873 -0.9633515 -8.9002162E-02 0.5391444 1.273044 2.359204 2.836445 2.293926 0.8976543 -1.168575 -2.884098 -3.103606 -2.019723 -0.9635568 -8.5187905E-02 0.5415012 1.291297 2.380054 2.847153 2.299733 0.8945679 -1.174553 -2.903498 -3.124593 -2.030300 -0.9634821 -8.3201811E-02 0.5438787 1.305556 2.400896 2.863496 2.306658 0.8878226 -1.185068 -2.920140 -3.144774 -2.038712 -0.9626333 -8.1342012E-02 0.5461245 1.319679 2.421609 2.879717 2.313476 0.8809856 -1.195659 -2.936841 -3.164999 -2.047152 -0.9618035 -7.9490878E-02 0.5483653 1.333800 2.442321 2.895938 2.320295 0.8741485 -1.206250 -2.953541 -3.185220 -2.055586 -0.9609529 -7.7605315E-02 0.5506690 1.348012 2.463167 2.912336 2.327341 0.8675922 -1.216508 -2.969857 -3.205015 -2.063553 -0.9596200 -7.7330299E-02 0.5503440 1.358636 2.482778 2.934216 2.336317 0.8568450 -1.230463 -2.981276 -3.223306 -2.070151 -0.9587819 -7.7154912E-02 0.5499091 1.369141 2.502276 2.955991 2.345202 0.8460240 -1.244474 -2.992736 -3.241627 -2.076766 -0.9579533 -7.6982312E-02 0.5494729 1.379645 2.521775 2.977765 2.354088 0.8352031 -1.258486 -3.004196 -3.259947 -2.083381 -0.9571256 -7.6811321E-02 0.5490321 1.390142 2.541260 2.999520 2.362947 0.8243483 -1.272537 -3.015701 -3.278312 -2.090042 -0.9563401 -7.6623864E-02 0.5482746 1.401835 2.560602 3.022145 2.372517 0.8129162 -1.287402 -3.027721 -3.297196 -2.096884 -0.9556683 -7.6317951E-02 0.5476279 1.413630 2.580036 3.044847 2.382150 0.8015336 -1.302230 -3.039716 -3.316063 -2.103717 -0.9549910 -7.6009966E-02 0.5469823 1.425426 2.599469 3.067550 2.391784 0.7901511 -1.317057 -3.051710 -3.334929 -2.110549 -0.9543115 -7.5698733E-02 0.5463407 1.437227 2.618907 3.090258 2.401424 0.7787760 -1.331874 -3.063692 -3.353776 -2.117358 -0.9536040 -7.6063856E-02 0.5461957 1.448669 2.638872 3.114306 2.410975 0.7667037 -1.348096 -3.076017 -3.372460 -2.123710 -0.9525445 -7.6393344E-02 0.5460796 1.460133 2.658854 3.138366 2.420532 0.7546321 -1.364318 -3.088346 -3.391149 -2.130067 -0.9514883 -7.6724812E-02 0.5459623 1.471596 2.678835 3.162424 2.430090 0.7425606 -1.380541 -3.100676 -3.409836 -2.136422 -0.9504288 -7.7051133E-02 0.5458546 1.483073 2.698837 3.186512 2.439687 0.7305394 -1.396702 -3.112931 -3.428438 -2.142681 -0.9492642
|
|
2
|
+
315.5519 315.6444 315.7370 315.8309 315.9249 316.0178 316.1107 316.2020 316.2933 316.3990 316.5046 316.6163 316.7279 316.8079 316.8878 316.9311 316.9743 317.0090 317.0437 317.0798 317.1160 317.1457 317.1753 317.2173 317.2592 317.3336 317.4079 317.5032 317.5985 317.6932 317.7880 317.8778 317.9677 318.0505 318.1334 318.2101 318.2868 318.3582 318.4297 318.4905 318.5515 318.6045 318.6576 318.7044 318.7512 318.7997 318.8482 318.8934 318.9387 318.9714 319.0042 319.0360 319.0678 319.1112 319.1547 319.2017 319.2487 319.2848 319.3208 319.3572 319.3935 319.4408 319.4881 319.5362 319.5844 319.6243 319.6643 319.7062 319.7481 319.8018 319.8554 319.9161 319.9768 320.0472 320.1175 320.1960 320.2744 320.3428 320.4112 320.4777 320.5442 320.6194 320.6946 320.7675 320.8404 320.8941 320.9478 321.0001 321.0525 321.1218 321.1912 321.2603 321.3293 321.3860 321.4426 321.4910 321.5395 321.5891 321.6387 321.6929 321.7472 321.8074 321.8677 321.9326 321.9974 322.0666 322.1357 322.2221 322.3085 322.3978 322.4871 322.5602 322.6334 322.7142 322.7951 322.9083 323.0215 323.1505 323.2796 323.4066 323.5336 323.6655 323.7974 323.9248 324.0522 324.1571 324.2621 324.3597 324.4574 324.5563 324.6551 324.7526 324.8501 324.9600 325.0699 325.1974 325.3250 325.4465 325.5681 325.6663 325.7645 325.8489 325.9333 326.0171 326.1010 326.1757 326.2505 326.3240 326.3975 326.4841 326.5707 326.6689 326.7671 326.8619 326.9567 327.0334 327.1101 327.1916 327.2731 327.3853 327.4976 327.6297 327.7618 327.9146 328.0674 328.2484 328.4293 328.6387 328.8480 329.0606 329.2731 329.4512 329.6294 329.7654 329.9015 329.9944 330.0872 330.1566 330.2260 330.2656 330.3052 330.3234 330.3416 330.3688 330.3960 330.4333 330.4705 330.4956 330.5207 330.5406 330.5605 330.5999 330.6394 330.6982 330.7570 330.8309 330.9047 330.9972 331.0897 331.2050 331.3204 331.4307 331.5411 331.6268 331.7126 331.7816 331.8506 331.9108 331.9710 332.0377 332.1044 332.1886 332.2728 332.3846 332.4964 332.6316 332.7668 332.9182 333.0696 333.2294 333.3892 333.5459 333.7026 333.8486 333.9946 334.1266 334.2587 334.3876 334.5166 334.6480 334.7794 334.9027 335.0260 335.1393 335.2527 335.3613 335.4700 335.5737 335.6774 335.7795 335.8816 335.9893 336.0971 336.2180 336.3389 336.4763 336.6138 336.7694 336.9250 337.0909 337.2568 337.4205 337.5841 337.7405 337.8969 338.0427 338.1884 338.3105 338.4326 338.5441 338.6555 338.7655 338.8755 338.9668 339.0580 339.1279 339.1978 339.2682 339.3386 339.4240 339.5094 339.5996 339.6898 339.7864 339.8830 339.9971 340.1111 340.2360 340.3608 340.4738 340.5868 340.6892 340.7917 340.9008 341.0099 341.1549 341.2999 341.4833 341.6667 341.8679 342.0690 342.2645 342.4600 342.6507 342.8414 343.0210 343.2006 343.3539 343.5072 343.6369 343.7666 343.8864 344.0062 344.1151 344.2241 344.3238 344.4236 344.5265 344.6294 344.7352 344.8410 344.9409 345.0409 345.1429 345.2449 345.3532 345.4615 345.5617 345.6619 345.7472 345.8325 345.9183 346.0041 346.1007 346.1973 346.2995 346.4017 346.5121 346.6226 346.7441 346.8656 346.9885 347.1115 347.2327 347.3539 347.4835 347.6131 347.7646 347.9162 348.1116 348.3071 348.4999 348.6928 348.8899 349.0870 349.2884 349.4897 349.7021
|
|
3
|
+
0.9473456 0.8820528 0.9981718 0.9523941 0.8169235 0.9968721 0.7275636 0.9908377 0.6374046 0.9900793 0.5952753 0.9727222 0.9657797 0.8412194 0.8863041 0.9650735 0.7412772 0.5911558 0.7918944 0.9976359 0.9126042 0.9287147 0.8981016 0.9992640 0.8516618 0.9641449 0.9999655 0.9084932 0.8229079 0.9775841 0.9999740 0.8893542 0.9814858 0.6921397 0.9997987 0.9893009 0.9868131 0.9194017 0.8840622 0.9997586 0.9944118 0.9979475 0.8869542 0.9995601 0.5853120 0.8529151 0.9910294 0.9921290 0.9993278 0.8233213 0.9964389 0.8739912 0.5036717 0.8577449 0.8799009 0.9003027 0.9457330 0.9885447 0.8570249 0.9877417 0.8345523 0.8603712 0.8764476 0.8674909 0.9248952 0.9995702 0.9426413 0.9972472 0.8482621 0.7178012 0.9944118 0.9972574 0.9405500 0.9139733 0.9851364 0.9969788 0.4954960 0.8312120 0.8113978 0.9723752 0.5457299 0.9453593 0.2643782 0.7797844 0.6239225 0.9947231 0.9904200 0.9260687 0.9922562 0.9681790 0.9866477 0.9993476 0.8408600 0.1812998 0.9994192 0.9156029 0.2818184 0.9746988 0.9856681 0.9676514 0.9889814 0.7985213 0.5947735 0.9931651 0.9261450 0.9992403 0.8309539 0.5941209 0.9845950 0.9931651 0.9551575 0.8795514 0.9389349 0.9398142 0.9144093 0.9944832 0.9138095 0.8916704 0.6464453 0.9860324 0.9495300 0.9990164 0.8044351 0.9871626 0.9498475 0.9004778 0.8719712 0.9955013 0.8243347 0.9442974 0.7816713 0.9952399 0.9405271 0.9865257 0.8954014 0.8391304 0.7213080 0.8975996 0.9288895 0.9381664 0.8617947 0.6964918 0.9709492 0.9894392 0.9655861 0.9818738 0.7010069 0.1039263 0.9831502 0.9979994 0.9821939 0.9348792 0.8464258 0.9887996 0.8859938 0.8913366 0.9858394 0.9611016 0.3264723 0.9174322 0.9840199 0.2168292 0.5075723 0.8040431 0.9408705 0.8893237 0.9207462 0.9894688 0.9976125 0.9792520 0.9570797 0.9320744 0.9494027 0.9463241 0.8543022 0.3240814 0.5827369 0.7932186 0.9575490 0.5153625 7.5613648E-02 0.9968399 0.9972875 0.9896256 0.9998081 0.7909694 0.9995499 0.8123233 0.9996711 0.9008568 0.9924324 0.9765879 0.9759859 0.8392387 0.5412561 0.6502172 0.9204307 0.8854656 0.9938530 0.9838016 0.8638172 0.9807110 0.9956102 0.9980421 0.9998107 0.9649314 0.9413726 0.9854034 0.9891018 0.9846898 0.9838746 0.9782504 0.9896938 0.7482746 0.8501588 0.9267029 0.9640011 0.8578818 0.9719143 0.9755108 0.9730045 0.9828760 0.9745318 0.8800598 0.8233213 0.9392827 0.9522703 0.9719781 0.8818638 0.9571581 0.9319032 0.9987569 0.7301404 0.9874538 0.9430911 0.8306218 0.9959595 0.9763831 0.9998436 0.9944333 0.9808433 0.9747442 0.9586537 0.8875102 0.8213251 0.9816285 0.9945189 0.9021925 0.8195850 0.9273343 0.9694217 0.9929974 0.7616304 0.9767481 0.3199174 0.9965921 0.9998900 0.8672252 0.9997818 0.9995781 0.9983168 0.9902315 0.9846542 0.8739912 0.9251511 0.6491545 0.7521552 0.9149528 0.8217025 0.9999934 0.7942597 0.6791515 0.5116771 0.9076493 0.9983443 0.9999844 0.9324158 0.8247468 0.6679452 0.9943252 0.9998839 0.9321719 0.9530309 0.9806179 0.8923972 0.9451830 0.9915662 0.9311679 0.8502639 0.9944261 0.5542958 0.9975463 0.8896603 0.8676236 0.7327968 0.1887221 0.8266525 0.9999914 0.8732761 0.9370639 0.9911562 0.9603313 0.9317075 0.8358660 0.9948754 0.9995601 0.9256614 0.9997001 0.3365756 0.9271325 0.9972321 0.9962590 0.9649491 0.9452492 0.4000554 0.9789069 0.9881831 0.9834098 0.9445416 0.8983963 0.9580159 0.9970051 0.9750311 0.9486145 0.9663744 0.9714338 0.8667263 0.9904200 0.9929007 0.9858394 0.9784330 0.8917311 0.5354862 0.9942672 0.9918803 0.9998053 0.9436064 0.8876644 0.8348811 0.8652251 0.9573147 0.9585381 0.7653055 0.9986812 0.8698055 0.7505626 0.9830258 0.9803512
|
|
4
|
+
2
|
|
5
|
+
-0.1312545 0.5171010 1.081863 2.133283 2.733717 2.239956 0.9297554 -1.114420 -2.675813 -2.878533 -1.902167 -0.9602861 -0.1273115 0.5191532 1.098498 2.153796 2.741747 2.243891 0.9270887 -1.118383 -2.694288 -2.898298 -1.912531 -0.9604589 -0.1233273 0.5212451 1.115170 2.174344 2.749810 2.247854 0.9244474 -1.122326 -2.712748 -2.918053 -1.922890 -0.9606276 -0.1193415 0.5233380 1.131844 2.194892 2.757873 2.251818 0.9218062 -1.126269 -2.731208 -2.937807 -1.933248 -0.9607950 -0.1153538 0.5254319 1.148517 2.215437 2.765927 2.255767 0.9191428 -1.130242 -2.749705 -2.957601 -1.943650 -0.9610062 -0.1116722 0.5270622 1.165944 2.235751 2.775559 2.261016 0.9161811 -1.135213 -2.768593 -2.978176 -1.954524 -0.9614759 -0.1078705 0.5288079 1.183482 2.256165 2.785280 2.266341 0.9132813 -1.140137 -2.787447 -2.998728 -1.965384 -0.9619374 -0.1040655 0.5305554 1.201020 2.276579 2.795002 2.271666 0.9103814 -1.145061 -2.806301 -3.019279 -1.976244 -0.9624004 -0.1002634 0.5322961 1.218547 2.296975 2.804697 2.276956 0.9074391 -1.150033 -2.825208 -3.039885 -1.987159 -0.9629111 -9.6585110E-02 0.5345098 1.236650 2.317664 2.815235 2.282577 0.9041511 -1.156232 -2.844850 -3.061133 -1.998018 -0.9631324 -9.2793740E-02 0.5368271 1.254847 2.338434 2.825840 2.288252 0.9009027 -1.162404 -2.864474 -3.082371 -2.008873 -0.9633515 -8.9002162E-02 0.5391444 1.273044 2.359204 2.836445 2.293926 0.8976543 -1.168575 -2.884098 -3.103606 -2.019723 -0.9635568 -8.5187905E-02 0.5415012 1.291297 2.380054 2.847153 2.299733 0.8945679 -1.174553 -2.903498 -3.124593 -2.030300 -0.9634821 -8.3201811E-02 0.5438787 1.305556 2.400896 2.863496 2.306658 0.8878226 -1.185068 -2.920140 -3.144774 -2.038712 -0.9626333 -8.1342012E-02 0.5461245 1.319679 2.421609 2.879717 2.313476 0.8809856 -1.195659 -2.936841 -3.164999 -2.047152 -0.9618035 -7.9490878E-02 0.5483653 1.333800 2.442321 2.895938 2.320295 0.8741485 -1.206250 -2.953541 -3.185220 -2.055586 -0.9609529 -7.7605315E-02 0.5506690 1.348012 2.463167 2.912336 2.327341 0.8675922 -1.216508 -2.969857 -3.205015 -2.063553 -0.9596200 -7.7330299E-02 0.5503440 1.358636 2.482778 2.934216 2.336317 0.8568450 -1.230463 -2.981276 -3.223306 -2.070151 -0.9587819 -7.7154912E-02 0.5499091 1.369141 2.502276 2.955991 2.345202 0.8460240 -1.244474 -2.992736 -3.241627 -2.076766 -0.9579533 -7.6982312E-02 0.5494729 1.379645 2.521775 2.977765 2.354088 0.8352031 -1.258486 -3.004196 -3.259947 -2.083381 -0.9571256 -7.6811321E-02 0.5490321 1.390142 2.541260 2.999520 2.362947 0.8243483 -1.272537 -3.015701 -3.278312 -2.090042 -0.9563401 -7.6623864E-02 0.5482746 1.401835 2.560602 3.022145 2.372517 0.8129162 -1.287402 -3.027721 -3.297196 -2.096884 -0.9556683 -7.6317951E-02 0.5476279 1.413630 2.580036 3.044847 2.382150 0.8015336 -1.302230 -3.039716 -3.316063 -2.103717 -0.9549910 -7.6009966E-02 0.5469823 1.425426 2.599469 3.067550 2.391784 0.7901511 -1.317057 -3.051710 -3.334929 -2.110549 -0.9543115 -7.5698733E-02 0.5463407 1.437227 2.618907 3.090258 2.401424 0.7787760 -1.331874 -3.063692 -3.353776 -2.117358 -0.9536040 -7.6063856E-02 0.5461957 1.448669 2.638872 3.114306 2.410975 0.7667037 -1.348096 -3.076017 -3.372460 -2.123710 -0.9525445 -7.6393344E-02 0.5460796 1.460133 2.658854 3.138366 2.420532 0.7546321 -1.364318 -3.088346 -3.391149 -2.130067 -0.9514883 -7.6724812E-02 0.5459623 1.471596 2.678835 3.162424 2.430090 0.7425606 -1.380541 -3.100676 -3.409836 -2.136422 -0.9504288 -7.7051133E-02 0.5458546 1.483073 2.698837 3.186512 2.439687 0.7305394 -1.396702 -3.112931 -3.428438 -2.142681 -0.9492642
|
|
6
|
+
315.5519 315.6444 315.7370 315.8309 315.9249 316.0178 316.1107 316.2020 316.2933 316.3990 316.5046 316.6163 316.7279 316.8079 316.8878 316.9311 316.9743 317.0090 317.0437 317.0798 317.1160 317.1457 317.1753 317.2173 317.2592 317.3336 317.4079 317.5032 317.5985 317.6932 317.7880 317.8778 317.9677 318.0505 318.1334 318.2101 318.2868 318.3582 318.4297 318.4905 318.5515 318.6045 318.6576 318.7044 318.7512 318.7997 318.8482 318.8934 318.9387 318.9714 319.0042 319.0360 319.0678 319.1112 319.1547 319.2017 319.2487 319.2848 319.3208 319.3572 319.3935 319.4408 319.4881 319.5362 319.5844 319.6243 319.6643 319.7062 319.7481 319.8018 319.8554 319.9161 319.9768 320.0472 320.1175 320.1960 320.2744 320.3428 320.4112 320.4777 320.5442 320.6194 320.6946 320.7675 320.8404 320.8941 320.9478 321.0001 321.0525 321.1218 321.1912 321.2603 321.3293 321.3860 321.4426 321.4910 321.5395 321.5891 321.6387 321.6929 321.7472 321.8074 321.8677 321.9326 321.9974 322.0666 322.1357 322.2221 322.3085 322.3978 322.4871 322.5602 322.6334 322.7142 322.7951 322.9083 323.0215 323.1505 323.2796 323.4066 323.5336 323.6655 323.7974 323.9248 324.0522 324.1571 324.2621 324.3597 324.4574 324.5563 324.6551 324.7526 324.8501 324.9600 325.0699 325.1974 325.3250 325.4465 325.5681 325.6663 325.7645 325.8489 325.9333 326.0171 326.1010 326.1757 326.2505 326.3240 326.3975 326.4841 326.5707 326.6689 326.7671 326.8619 326.9567 327.0334 327.1101 327.1916 327.2731 327.3853 327.4976 327.6297 327.7618 327.9146 328.0674 328.2484 328.4293 328.6387 328.8480 329.0606 329.2731 329.4512 329.6294 329.7654 329.9015 329.9944 330.0872 330.1566 330.2260 330.2656 330.3052 330.3234 330.3416 330.3688 330.3960 330.4333 330.4705 330.4956 330.5207 330.5406 330.5605 330.5999 330.6394 330.6982 330.7570 330.8309 330.9047 330.9972 331.0897 331.2050 331.3204 331.4307 331.5411 331.6268 331.7126 331.7816 331.8506 331.9108 331.9710 332.0377 332.1044 332.1886 332.2728 332.3846 332.4964 332.6316 332.7668 332.9182 333.0696 333.2294 333.3892 333.5459 333.7026 333.8486 333.9946 334.1266 334.2587 334.3876 334.5166 334.6480 334.7794 334.9027 335.0260 335.1393 335.2527 335.3613 335.4700 335.5737 335.6774 335.7795 335.8816 335.9893 336.0971 336.2180 336.3389 336.4763 336.6138 336.7694 336.9250 337.0909 337.2568 337.4205 337.5841 337.7405 337.8969 338.0427 338.1884 338.3105 338.4326 338.5441 338.6555 338.7655 338.8755 338.9668 339.0580 339.1279 339.1978 339.2682 339.3386 339.4240 339.5094 339.5996 339.6898 339.7864 339.8830 339.9971 340.1111 340.2360 340.3608 340.4738 340.5868 340.6892 340.7917 340.9008 341.0099 341.1549 341.2999 341.4833 341.6667 341.8679 342.0690 342.2645 342.4600 342.6507 342.8414 343.0210 343.2006 343.3539 343.5072 343.6369 343.7666 343.8864 344.0062 344.1151 344.2241 344.3238 344.4236 344.5265 344.6294 344.7352 344.8410 344.9409 345.0409 345.1429 345.2449 345.3532 345.4615 345.5617 345.6619 345.7472 345.8325 345.9183 346.0041 346.1007 346.1973 346.2995 346.4017 346.5121 346.6226 346.7441 346.8656 346.9885 347.1115 347.2327 347.3539 347.4835 347.6131 347.7646 347.9162 348.1116 348.3071 348.4999 348.6928 348.8899 349.0870 349.2884 349.4897 349.7021
|
|
7
|
+
0.9473456 0.8820528 0.9981718 0.9523941 0.8169235 0.9968721 0.7275636 0.9908377 0.6374046 0.9900793 0.5952753 0.9727222 0.9657797 0.8412194 0.8863041 0.9650735 0.7412772 0.5911558 0.7918944 0.9976359 0.9126042 0.9287147 0.8981016 0.9992640 0.8516618 0.9641449 0.9999655 0.9084932 0.8229079 0.9775841 0.9999740 0.8893542 0.9814858 0.6921397 0.9997987 0.9893009 0.9868131 0.9194017 0.8840622 0.9997586 0.9944118 0.9979475 0.8869542 0.9995601 0.5853120 0.8529151 0.9910294 0.9921290 0.9993278 0.8233213 0.9964389 0.8739912 0.5036717 0.8577449 0.8799009 0.9003027 0.9457330 0.9885447 0.8570249 0.9877417 0.8345523 0.8603712 0.8764476 0.8674909 0.9248952 0.9995702 0.9426413 0.9972472 0.8482621 0.7178012 0.9944118 0.9972574 0.9405500 0.9139733 0.9851364 0.9969788 0.4954960 0.8312120 0.8113978 0.9723752 0.5457299 0.9453593 0.2643782 0.7797844 0.6239225 0.9947231 0.9904200 0.9260687 0.9922562 0.9681790 0.9866477 0.9993476 0.8408600 0.1812998 0.9994192 0.9156029 0.2818184 0.9746988 0.9856681 0.9676514 0.9889814 0.7985213 0.5947735 0.9931651 0.9261450 0.9992403 0.8309539 0.5941209 0.9845950 0.9931651 0.9551575 0.8795514 0.9389349 0.9398142 0.9144093 0.9944832 0.9138095 0.8916704 0.6464453 0.9860324 0.9495300 0.9990164 0.8044351 0.9871626 0.9498475 0.9004778 0.8719712 0.9955013 0.8243347 0.9442974 0.7816713 0.9952399 0.9405271 0.9865257 0.8954014 0.8391304 0.7213080 0.8975996 0.9288895 0.9381664 0.8617947 0.6964918 0.9709492 0.9894392 0.9655861 0.9818738 0.7010069 0.1039263 0.9831502 0.9979994 0.9821939 0.9348792 0.8464258 0.9887996 0.8859938 0.8913366 0.9858394 0.9611016 0.3264723 0.9174322 0.9840199 0.2168292 0.5075723 0.8040431 0.9408705 0.8893237 0.9207462 0.9894688 0.9976125 0.9792520 0.9570797 0.9320744 0.9494027 0.9463241 0.8543022 0.3240814 0.5827369 0.7932186 0.9575490 0.5153625 7.5613648E-02 0.9968399 0.9972875 0.9896256 0.9998081 0.7909694 0.9995499 0.8123233 0.9996711 0.9008568 0.9924324 0.9765879 0.9759859 0.8392387 0.5412561 0.6502172 0.9204307 0.8854656 0.9938530 0.9838016 0.8638172 0.9807110 0.9956102 0.9980421 0.9998107 0.9649314 0.9413726 0.9854034 0.9891018 0.9846898 0.9838746 0.9782504 0.9896938 0.7482746 0.8501588 0.9267029 0.9640011 0.8578818 0.9719143 0.9755108 0.9730045 0.9828760 0.9745318 0.8800598 0.8233213 0.9392827 0.9522703 0.9719781 0.8818638 0.9571581 0.9319032 0.9987569 0.7301404 0.9874538 0.9430911 0.8306218 0.9959595 0.9763831 0.9998436 0.9944333 0.9808433 0.9747442 0.9586537 0.8875102 0.8213251 0.9816285 0.9945189 0.9021925 0.8195850 0.9273343 0.9694217 0.9929974 0.7616304 0.9767481 0.3199174 0.9965921 0.9998900 0.8672252 0.9997818 0.9995781 0.9983168 0.9902315 0.9846542 0.8739912 0.9251511 0.6491545 0.7521552 0.9149528 0.8217025 0.9999934 0.7942597 0.6791515 0.5116771 0.9076493 0.9983443 0.9999844 0.9324158 0.8247468 0.6679452 0.9943252 0.9998839 0.9321719 0.9530309 0.9806179 0.8923972 0.9451830 0.9915662 0.9311679 0.8502639 0.9944261 0.5542958 0.9975463 0.8896603 0.8676236 0.7327968 0.1887221 0.8266525 0.9999914 0.8732761 0.9370639 0.9911562 0.9603313 0.9317075 0.8358660 0.9948754 0.9995601 0.9256614 0.9997001 0.3365756 0.9271325 0.9972321 0.9962590 0.9649491 0.9452492 0.4000554 0.9789069 0.9881831 0.9834098 0.9445416 0.8983963 0.9580159 0.9970051 0.9750311 0.9486145 0.9663744 0.9714338 0.8667263 0.9904200 0.9929007 0.9858394 0.9784330 0.8917311 0.5354862 0.9942672 0.9918803 0.9998053 0.9436064 0.8876644 0.8348811 0.8652251 0.9573147 0.9585381 0.7653055 0.9986812 0.8698055 0.7505626 0.9830258 0.9803512
|
|
@@ -0,0 +1,3 @@
|
|
|
1
|
+
E (observations)
|
|
2
|
+
0.831000 1.045000 1.021000 0.970000 0.825000 0.891000 0.710000 0.801000 1.074000 1.148000 1.000000 0.928000 0.767000 0.701000 0.807000 0.902000 0.997000 1.224000 1.089000 0.973000 0.980000 0.665000 NOx (response)
|
|
3
|
+
4.818000 2.849000 3.275000 4.691000 4.255000 5.064000 2.118000 4.602000 2.286000 0.970000 3.965000 5.344000 3.834000 1.990000 5.199000 5.283000 3.752000 0.537000 1.640000 5.055000 4.937000 1.561000
|
|
@@ -0,0 +1,10 @@
|
|
|
1
|
+
loess(&gas):
|
|
2
|
+
4.883437 2.927639 3.572746 4.710510 4.815363 5.196193 2.512832 4.474939 2.149797 0.990044 4.089353 5.303627 3.864174 2.267121 4.575636 5.240293 4.154036 0.523682 1.853026 4.659184 4.521148 1.196414
|
|
3
|
+
loess(&gas_null):
|
|
4
|
+
4.892690 3.362733 3.728575 4.489023 4.831715 5.177112 2.535887 4.530370 2.873000 1.453502 4.059136 5.052650 3.943568 2.265232 4.614320 5.163913 4.107184 -0.276342 2.604373 4.445040 4.351615 1.051719
|
|
5
|
+
predict(gas_fit_E, m, &gas, &gas_pred, 0):
|
|
6
|
+
1.196414 5.068747 0.523682
|
|
7
|
+
pointwise(&gas_pred, m, coverage, &gas_ci):
|
|
8
|
+
0.407208 1.196414 1.985621 3.249187 3.679498 4.109808 4.631187 5.055708 5.480229 4.704010 5.135260 5.566510 2.759703 3.143656 3.527609 0.683247 1.196932 1.710617 -0.424684 0.523682 1.472049
|
|
9
|
+
anova(&gas_null, &gas, &gas_anova):
|
|
10
|
+
2.5531 15.663 10.1397 0.000860102
|
|
@@ -0,0 +1,4 @@
|
|
|
1
|
+
one_two (observations)
|
|
2
|
+
-0.957581 -2.809549 -0.696511 3.451000 0.509260 0.557854 0.052582 -2.050644 -1.115675 -1.183665 0.511959 0.334364 -2.057062 -0.121897 0.544238 0.600502 0.531074 0.495400 -1.608602 0.277371 0.290464 0.579894 -0.290441 1.306226 -0.482898 -0.716423 0.742413 -0.911613 1.279436 -0.189153 0.592293 0.952416 0.491436 -0.305681 -0.363871 -0.285424 -0.037209 -0.923529 1.138054 -1.331223 0.551234 -0.852726 1.196875 0.498782 0.320180 0.212447 1.009358 -0.900989 1.132164 0.018867 0.424170 -0.198621 0.955170 0.948321 0.473848 -0.699121 -0.612853 0.580431 1.277996 0.806797 -1.038559 1.008663 -0.578257 -0.323245 -0.756302 1.386352 0.722419 -1.216078 -0.498280 0.726247 -0.260119 -0.741135 -0.184111 0.307762 0.464568 -0.252531 -0.486504 0.426634 -1.303969 0.067149 1.771176 0.907249 0.432350 1.419897 -0.413389 2.442025 0.041138 0.509505 -0.282744 0.179882 -1.188083 0.982653 -1.042886 1.181365 -0.398340 -1.335565 -0.502790 0.484762 -0.806446 1.412077 -0.878874 -0.935197 -0.339255 0.164497 1.370018 -1.494684 1.380505 0.885084 0.835609 0.896235 -1.289541 0.233203 1.183198 -0.857934 -1.334234 -0.923351 0.769146 -0.377948 0.059114 -1.870615 -0.677868 0.038185 0.375301 0.964717 0.695051 -0.342140 -1.145463 -0.993246 -0.130573 1.213711 0.291241 1.106891 0.949571 0.463675 0.455723 0.398786 -0.015849 -1.397373 0.770062 0.083291 0.531798 0.049727 -0.734140 -0.963487 0.573561 -0.281942 -0.594503 0.770262 1.073983 -0.615706 -0.084794 -0.491630 -1.526969 -0.196881 0.165653 0.198357 0.974930 -0.954848 0.588474 -0.426883 0.177119 -0.916442 -1.885139 0.086894 0.456306 0.174285 -0.001308 -0.000585 0.284023 -0.365679 -0.548867 0.857848 0.699094 -0.140026 1.332454 1.601795 0.012415 0.243429 1.077369 1.859246 0.185910 0.033342 0.613008 1.068595 -0.683305 -0.128826 -1.655525 0.013086 0.062454 0.773042 0.127046 0.408652 1.195438 -0.185558 -1.299714 0.899675 -0.033648 -1.544602 0.655203 -0.713935
|
|
3
|
+
response
|
|
4
|
+
14.453553 6.622825 13.671414 14.197518 12.860530 12.522856 14.214638 7.924264 12.506938 13.734205 14.710855 13.596223 5.890019 13.558654 14.043167 13.931391 13.218920 17.090560 15.199322 13.261667 15.760636 12.083855 14.344907 12.609494 11.932959 13.408674 13.700765 13.013366 15.794999 14.600198 16.275751 11.564349 14.809023 12.982361 15.003502 14.737337 15.747677 11.674508 14.047278 14.666917 13.806240 13.611149 13.347149 14.225152 14.718846 14.217216 14.418058 14.719634 12.799716 13.933038 15.264603 14.660387 9.738691 14.443424 14.417284 15.184538 13.344938 15.372943 13.811554 15.103777 15.383834 14.368612 12.525202 14.325033 15.259658 13.004547 14.515988 15.176982 14.924187 13.872430 15.395366 13.428076 15.203430 14.186631 13.305833 14.074624 14.103092 13.499669 11.584675 14.264891 14.885616 13.967297 16.604680 10.367606 14.743473 16.308827 14.108673 13.590988 14.674546 15.294047 14.686723 13.611422 11.970270 13.884157 15.071776 12.589816 13.818745 14.245317 14.406530 14.347941
|
|
@@ -0,0 +1,16 @@
|
|
|
1
|
+
loess(&madeup):
|
|
2
|
+
13.898675 5.225583 14.265096 14.237045 13.366442 13.157273 14.189032 7.971348 13.215581 12.835208 13.967757 14.457845 6.039091 14.722007 13.758002 14.654508 14.243303 14.733353 13.696778 12.510036 15.026972 14.346972 14.547927 11.702517 14.222303 14.288345 13.755227 13.742749 14.250887 14.188543 14.388616 12.552611 14.100045 14.552719 14.514971 14.556744 14.496872 12.868609 12.871890 14.427522 14.410413 14.598345 13.415146 14.731999 14.657565 14.778805 13.988940 13.723749 11.932352 14.908485 14.457978 14.616147 11.100703 14.254687 14.411890 14.519007 13.734456 14.634786 13.080897 14.415189 14.623631 13.340049 12.584224 14.510466 14.293842 13.954341 14.281046 14.643596 14.492646 14.497077 14.497688 13.732467 14.573564 14.548255 13.560984 13.917026 14.496924 14.452028 11.687567 13.865395 13.908958 14.184016 14.514209 10.835782 14.203921 15.273614 13.702645 14.383224 14.506619 14.707063 14.595778 13.877506 12.326903 14.236236 13.962165 12.281695 14.477443 13.192765 13.985441 13.206705
|
|
3
|
+
loess(&madeup_new):
|
|
4
|
+
12.874079 10.889912 14.072472 14.688408 13.581655 13.741506 12.966195 11.778153 12.623840 12.391277 13.915575 14.549806 10.648311 13.714716 13.908051 14.338239 13.839956 14.661187 14.500423 12.892713 14.344422 14.557473 14.294094 14.258892 13.410545 14.058327 14.340850 12.675257 14.839750 12.706903 14.551439 14.002773 13.589104 14.094192 14.098484 14.246415 14.548803 11.925623 14.212327 14.501013 14.199258 14.517379 15.151176 14.202325 14.044125 14.549456 15.030205 12.957176 14.170346 14.166684 14.617600 14.160491 14.230993 14.821136 14.547624 14.489024 12.630700 14.292239 14.381246 14.834949 14.480650 14.747398 11.611371 14.519760 13.997474 14.591257 14.591759 14.478131 14.427821 14.785485 14.034793 12.821834 13.556879 14.606041 13.505266 12.561710 14.525273 14.555914 11.728344 13.296039 14.621341 14.570526 14.015254 14.451527 13.670446 15.113971 12.849321 14.570355 14.526773 13.616773 14.493652 14.496696 11.739043 14.886699 12.740301 12.265455 14.504240 13.584999 13.370259 15.378514
|
|
5
|
+
loess(&madeup_new) (family = symmetric):
|
|
6
|
+
13.853372 13.076392 14.410502 14.317606 13.530344 13.909561 13.483542 13.661465 13.738437 13.662424 14.078653 14.471725 13.362504 14.161261 14.056621 14.370400 13.827138 14.543779 14.563051 13.263340 14.405462 14.511807 14.391442 13.655949 13.919863 14.406599 14.307597 13.717301 14.561751 13.479815 14.423533 13.588048 13.644027 14.268385 14.283955 14.361030 14.528760 13.151054 13.826620 14.570806 14.120641 14.538117 14.640103 14.309962 14.069715 14.527437 14.606612 13.824681 13.607361 14.370144 14.535820 14.399329 14.122194 14.580712 14.478250 14.541334 13.593458 14.346338 14.016117 14.580121 14.566098 14.522580 12.578184 14.528995 14.305679 14.422916 14.523434 14.533691 14.490958 14.574831 14.352653 13.714577 13.911857 14.540189 13.522831 13.495538 14.528432 14.462617 13.481391 13.705501 14.399779 14.509204 14.010568 13.623165 14.211318 14.586042 13.392231 14.520972 14.525897 13.822231 14.577431 14.267078 13.423198 14.589490 13.573646 13.655627 14.531561 13.811793 13.979165 14.705119
|
|
7
|
+
loess(&madeup_new) (normalize = FALSE):
|
|
8
|
+
13.853372 13.076392 14.410502 14.317606 13.530344 13.909561 13.483542 13.661465 13.738437 13.662424 14.078653 14.471725 13.362504 14.161261 14.056621 14.370400 13.827138 14.543779 14.563051 13.263340 14.405462 14.511807 14.391442 13.655949 13.919863 14.406599 14.307597 13.717301 14.561751 13.479815 14.423533 13.588048 13.644027 14.268385 14.283955 14.361030 14.528760 13.151054 13.826620 14.570806 14.120641 14.538117 14.640103 14.309962 14.069715 14.527437 14.606612 13.824681 13.607361 14.370144 14.535820 14.399329 14.122194 14.580712 14.478250 14.541334 13.593458 14.346338 14.016117 14.580121 14.566098 14.522580 12.578184 14.528995 14.305679 14.422916 14.523434 14.533691 14.490958 14.574831 14.352653 13.714577 13.911857 14.540189 13.522831 13.495538 14.528432 14.462617 13.481391 13.705501 14.399779 14.509204 14.010568 13.623165 14.211318 14.586042 13.392231 14.520972 14.525897 13.822231 14.577431 14.267078 13.423198 14.589490 13.573646 13.655627 14.531561 13.811793 13.979165 14.705119
|
|
9
|
+
predict(newdata1, m, &madeup, &madeup_pred, 0):
|
|
10
|
+
8.156780 14.493587 14.854142
|
|
11
|
+
predict(newdata2, m, &madeup, &madeup_pred, 1):
|
|
12
|
+
14.491811 14.389726
|
|
13
|
+
pointwise(&madeup_pred, m, coverage, &madeup_ci):
|
|
14
|
+
13.761833 14.491811 15.221789 13.656419 14.389726 15.123034
|
|
15
|
+
anova(&madeup2, &madeup, &madeup_anova):
|
|
16
|
+
6.45034 81.2319 2.85933 0.0121449
|
|
@@ -0,0 +1,284 @@
|
|
|
1
|
+
import pickle
|
|
2
|
+
import tempfile
|
|
3
|
+
from pathlib import Path
|
|
4
|
+
|
|
5
|
+
import numpy as np
|
|
6
|
+
import numpy.testing as npt
|
|
7
|
+
import pytest
|
|
8
|
+
|
|
9
|
+
from skmisc.loess import loess, loess_anova
|
|
10
|
+
|
|
11
|
+
data_path = Path(__file__).parent
|
|
12
|
+
|
|
13
|
+
|
|
14
|
+
def madeup_data():
|
|
15
|
+
dfile = data_path / 'madeup_data'
|
|
16
|
+
rfile = data_path / 'madeup_result'
|
|
17
|
+
|
|
18
|
+
with open(dfile, 'r') as f:
|
|
19
|
+
f.readline()
|
|
20
|
+
x = np.fromiter(
|
|
21
|
+
(float(v) for v in f.readline().rstrip().split()),
|
|
22
|
+
np.float64).reshape(-1, 2)
|
|
23
|
+
f.readline()
|
|
24
|
+
y = np.fromiter(
|
|
25
|
+
(float(v) for v in f.readline().rstrip().split()),
|
|
26
|
+
np.float64)
|
|
27
|
+
|
|
28
|
+
results = []
|
|
29
|
+
with open(rfile, 'r') as f:
|
|
30
|
+
for i in range(8):
|
|
31
|
+
f.readline()
|
|
32
|
+
z = np.fromiter(
|
|
33
|
+
(float(v) for v in f.readline().rstrip().split()),
|
|
34
|
+
np.float64)
|
|
35
|
+
results.append(z)
|
|
36
|
+
|
|
37
|
+
newdata1 = np.array([[-2.5, 0.], [2.5, 0.], [0., 0.]])
|
|
38
|
+
newdata2 = np.array([[-0.5, 0.5], [0., 0.]])
|
|
39
|
+
return (x, y, results, newdata1, newdata2)
|
|
40
|
+
|
|
41
|
+
|
|
42
|
+
def gas_data():
|
|
43
|
+
NOx = np.array([4.818, 2.849, 3.275, 4.691, 4.255, 5.064, 2.118, 4.602,
|
|
44
|
+
2.286, 0.970, 3.965, 5.344, 3.834, 1.990, 5.199, 5.283,
|
|
45
|
+
3.752, 0.537, 1.640, 5.055, 4.937, 1.561])
|
|
46
|
+
E = np.array([0.831, 1.045, 1.021, 0.970, 0.825, 0.891, 0.71, 0.801,
|
|
47
|
+
1.074, 1.148, 1.000, 0.928, 0.767, 0.701, 0.807, 0.902,
|
|
48
|
+
0.997, 1.224, 1.089, 0.973, 0.980, 0.665])
|
|
49
|
+
gas_fit_E = np.array([0.665, 0.949, 1.224])
|
|
50
|
+
newdata = np.array([0.6650000, 0.7581667, 0.8513333, 0.9445000,
|
|
51
|
+
1.0376667, 1.1308333, 1.2240000])
|
|
52
|
+
alpha = 0.01
|
|
53
|
+
|
|
54
|
+
rfile = data_path / 'gas_result'
|
|
55
|
+
results = []
|
|
56
|
+
with open(rfile, 'r') as f:
|
|
57
|
+
for i in range(8):
|
|
58
|
+
f.readline()
|
|
59
|
+
z = np.fromiter(
|
|
60
|
+
(float(v) for v in f.readline().rstrip().split()),
|
|
61
|
+
np.float64)
|
|
62
|
+
results.append(z)
|
|
63
|
+
return (E, NOx, gas_fit_E, newdata, alpha, results)
|
|
64
|
+
|
|
65
|
+
|
|
66
|
+
class TestLoess2d(object):
|
|
67
|
+
"Test class for lowess."
|
|
68
|
+
d = madeup_data()
|
|
69
|
+
|
|
70
|
+
def test_2dbasic(self):
|
|
71
|
+
"2D standard"
|
|
72
|
+
(x, y, results, _, _) = self.d
|
|
73
|
+
madeup = loess(x, y)
|
|
74
|
+
madeup.model.span = 0.5
|
|
75
|
+
madeup.model.normalize = True
|
|
76
|
+
madeup.fit()
|
|
77
|
+
npt.assert_almost_equal(madeup.outputs.fitted_values, results[0], 5)
|
|
78
|
+
npt.assert_almost_equal(madeup.outputs.enp, 14.9, 1)
|
|
79
|
+
npt.assert_almost_equal(madeup.outputs.residual_scale, 0.9693, 4)
|
|
80
|
+
|
|
81
|
+
def test_2d_modflags(self):
|
|
82
|
+
"2D - modification of model flags"
|
|
83
|
+
(x, y, results, _, _) = self.d
|
|
84
|
+
madeup = loess(x, y)
|
|
85
|
+
madeup.model.span = 0.8
|
|
86
|
+
madeup.model.drop_square = [True, False]
|
|
87
|
+
madeup.model.parametric = [True, False]
|
|
88
|
+
npt.assert_equal(madeup.model.parametric[:2], [1, 0])
|
|
89
|
+
madeup.fit()
|
|
90
|
+
npt.assert_almost_equal(madeup.outputs.fitted_values, results[1], 5)
|
|
91
|
+
npt.assert_almost_equal(madeup.outputs.enp, 6.9, 1)
|
|
92
|
+
npt.assert_almost_equal(madeup.outputs.residual_scale, 1.4804, 4)
|
|
93
|
+
|
|
94
|
+
def test_2d_modfamily(self):
|
|
95
|
+
"2D - family modification"
|
|
96
|
+
(x, y, results, _, _) = self.d
|
|
97
|
+
madeup = loess(x, y)
|
|
98
|
+
madeup.model.span = 0.8
|
|
99
|
+
madeup.model.drop_square = [True, False]
|
|
100
|
+
madeup.model.parametric = [True, False]
|
|
101
|
+
madeup.model.family = "symmetric"
|
|
102
|
+
madeup.fit()
|
|
103
|
+
npt.assert_almost_equal(madeup.outputs.fitted_values, results[2], 5)
|
|
104
|
+
npt.assert_almost_equal(madeup.outputs.enp, 6.9, 1)
|
|
105
|
+
npt.assert_almost_equal(madeup.outputs.residual_scale, 1.0868, 4)
|
|
106
|
+
|
|
107
|
+
def test_2d_modnormalize(self):
|
|
108
|
+
"2D - normalization modification"
|
|
109
|
+
(x, y, results, _, _) = self.d
|
|
110
|
+
madeup = loess(x, y)
|
|
111
|
+
madeup.model.span = 0.8
|
|
112
|
+
madeup.model.drop_square = [True, False]
|
|
113
|
+
madeup.model.parametric = [True, False]
|
|
114
|
+
madeup.model.family = "symmetric"
|
|
115
|
+
madeup.model.normalize = False
|
|
116
|
+
madeup.fit()
|
|
117
|
+
npt.assert_almost_equal(madeup.outputs.fitted_values, results[3], 5)
|
|
118
|
+
npt.assert_almost_equal(madeup.outputs.enp, 6.9, 1)
|
|
119
|
+
npt.assert_almost_equal(madeup.outputs.residual_scale, 1.0868, 4)
|
|
120
|
+
|
|
121
|
+
def test_2d_pred_nostderr(self):
|
|
122
|
+
"2D prediction - no stderr"
|
|
123
|
+
(x, y, results, newdata1, _) = self.d
|
|
124
|
+
madeup = loess(x, y)
|
|
125
|
+
madeup.model.span = 0.5
|
|
126
|
+
madeup.model.normalize = True
|
|
127
|
+
prediction = madeup.predict(newdata1, stderror=False)
|
|
128
|
+
npt.assert_almost_equal(prediction.values, results[4], 5)
|
|
129
|
+
#
|
|
130
|
+
prediction = madeup.predict(newdata1, stderror=False)
|
|
131
|
+
npt.assert_almost_equal(prediction.values, results[4], 5)
|
|
132
|
+
|
|
133
|
+
def test_2d_pred_nodata(self):
|
|
134
|
+
"2D prediction - nodata"
|
|
135
|
+
(x, y, _, _, _) = self.d
|
|
136
|
+
madeup = loess(x, y)
|
|
137
|
+
try:
|
|
138
|
+
madeup.predict(None)
|
|
139
|
+
except ValueError:
|
|
140
|
+
pass
|
|
141
|
+
else:
|
|
142
|
+
raise AssertionError("The test should have failed")
|
|
143
|
+
|
|
144
|
+
def test_2d_pred_stderr(self):
|
|
145
|
+
"2D prediction - w/ stderr"
|
|
146
|
+
(x, y, results, _, newdata2) = self.d
|
|
147
|
+
madeup = loess(x, y)
|
|
148
|
+
madeup.model.span = 0.5
|
|
149
|
+
madeup.model.normalize = True
|
|
150
|
+
prediction = madeup.predict(newdata2, stderror=True)
|
|
151
|
+
npt.assert_almost_equal(prediction.values, results[5], 5)
|
|
152
|
+
npt.assert_almost_equal(prediction.stderr, [0.276746, 0.278009], 5)
|
|
153
|
+
npt.assert_almost_equal(prediction.residual_scale, 0.969302, 6)
|
|
154
|
+
npt.assert_almost_equal(prediction.df, 81.2319, 4)
|
|
155
|
+
|
|
156
|
+
# Direct access
|
|
157
|
+
prediction = madeup.predict(newdata2, stderror=True)
|
|
158
|
+
npt.assert_almost_equal(prediction.values, results[5], 5)
|
|
159
|
+
npt.assert_almost_equal(prediction.stderr, [0.276746, 0.278009], 5)
|
|
160
|
+
npt.assert_almost_equal(prediction.residual_scale, 0.969302, 6)
|
|
161
|
+
npt.assert_almost_equal(prediction.df, 81.2319, 4)
|
|
162
|
+
|
|
163
|
+
def test_2d_pred_confinv(self):
|
|
164
|
+
"2D prediction - confidence"
|
|
165
|
+
(x, y, results, _, newdata2) = self.d
|
|
166
|
+
madeup = loess(x, y)
|
|
167
|
+
madeup.model.span = 0.5
|
|
168
|
+
madeup.model.normalize = True
|
|
169
|
+
prediction = madeup.predict(newdata2, stderror=True)
|
|
170
|
+
ci = prediction.confidence(alpha=0.01)
|
|
171
|
+
npt.assert_almost_equal(ci.lower, results[6][::3], 5)
|
|
172
|
+
npt.assert_almost_equal(ci.fit, results[6][1::3], 5)
|
|
173
|
+
npt.assert_almost_equal(ci.upper, results[6][2::3], 5)
|
|
174
|
+
|
|
175
|
+
|
|
176
|
+
class TestLoessGas(object):
|
|
177
|
+
"Test class for lowess."
|
|
178
|
+
|
|
179
|
+
d = gas_data()
|
|
180
|
+
|
|
181
|
+
def test_1dbasic(self):
|
|
182
|
+
"Basic test 1d"
|
|
183
|
+
(E, NOx, _, _, _, results) = self.d
|
|
184
|
+
gas = loess(E, NOx)
|
|
185
|
+
gas.model.span = 2./3.
|
|
186
|
+
gas.fit()
|
|
187
|
+
npt.assert_almost_equal(gas.outputs.fitted_values, results[0], 6)
|
|
188
|
+
npt.assert_almost_equal(gas.outputs.enp, 5.5, 1)
|
|
189
|
+
npt.assert_almost_equal(gas.outputs.residual_scale, 0.3404, 4)
|
|
190
|
+
|
|
191
|
+
def test_1dbasic_alt(self):
|
|
192
|
+
"Basic test 1d - part #2"
|
|
193
|
+
(E, NOx, _, _, _, results) = self.d
|
|
194
|
+
gas_null = loess(E, NOx)
|
|
195
|
+
gas_null.model.span = 1.0
|
|
196
|
+
gas_null.fit()
|
|
197
|
+
npt.assert_almost_equal(gas_null.outputs.fitted_values, results[1], 6)
|
|
198
|
+
npt.assert_almost_equal(gas_null.outputs.enp, 3.5, 1)
|
|
199
|
+
npt.assert_almost_equal(gas_null.outputs.residual_scale, 0.5197, 4)
|
|
200
|
+
|
|
201
|
+
def test_1dpredict(self):
|
|
202
|
+
"Basic test 1d - prediction"
|
|
203
|
+
(E, NOx, gas_fit_E, _, _, results) = self.d
|
|
204
|
+
gas = loess(E, NOx, span=2./3.)
|
|
205
|
+
gas.fit()
|
|
206
|
+
predicted = gas.predict(gas_fit_E, stderror=False).values
|
|
207
|
+
npt.assert_almost_equal(predicted, results[2], 6)
|
|
208
|
+
|
|
209
|
+
def test_1dpredict_2(self):
|
|
210
|
+
"Basic test 1d - new predictions"
|
|
211
|
+
(E, NOx, _, newdata, _, results) = self.d
|
|
212
|
+
# gas = loess(E, NOx, span=2./3.)
|
|
213
|
+
gas = loess(E, NOx)
|
|
214
|
+
gas.model.span = 2./3.
|
|
215
|
+
prediction = gas.predict(newdata, stderror=True)
|
|
216
|
+
ci = prediction.confidence(alpha=0.01)
|
|
217
|
+
npt.assert_almost_equal(ci.lower, results[3][0::3], 6)
|
|
218
|
+
npt.assert_almost_equal(ci.fit, results[3][1::3], 6)
|
|
219
|
+
npt.assert_almost_equal(ci.upper, results[3][2::3], 6)
|
|
220
|
+
|
|
221
|
+
def test_anova(self):
|
|
222
|
+
"Tests anova"
|
|
223
|
+
(E, NOx, _, _, _, results) = self.d
|
|
224
|
+
gas = loess(E, NOx, span=2./3.)
|
|
225
|
+
gas.fit()
|
|
226
|
+
gas_null = loess(E, NOx, span=1.0)
|
|
227
|
+
gas_null.fit()
|
|
228
|
+
gas_anova = loess_anova(gas, gas_null)
|
|
229
|
+
gas_anova_theo = results[4]
|
|
230
|
+
npt.assert_almost_equal(gas_anova.dfn, gas_anova_theo[0], 5)
|
|
231
|
+
npt.assert_almost_equal(gas_anova.dfd, gas_anova_theo[1], 5)
|
|
232
|
+
npt.assert_almost_equal(gas_anova.F_value, gas_anova_theo[2], 5)
|
|
233
|
+
npt.assert_almost_equal(gas_anova.Pr_F, gas_anova_theo[3], 5)
|
|
234
|
+
|
|
235
|
+
def test_failures(self):
|
|
236
|
+
"Tests failures"
|
|
237
|
+
(E, NOx, gas_fit_E, _, _, _) = self.d
|
|
238
|
+
gas = loess(E, NOx, span=2./3.)
|
|
239
|
+
# This one should fail (all parametric)
|
|
240
|
+
gas.model.parametric = True
|
|
241
|
+
with pytest.raises(ValueError):
|
|
242
|
+
gas.fit()
|
|
243
|
+
|
|
244
|
+
# This one also (all drop_square)
|
|
245
|
+
gas.model.drop_square = True
|
|
246
|
+
with pytest.raises(ValueError):
|
|
247
|
+
gas.fit()
|
|
248
|
+
|
|
249
|
+
gas.model.degree = 1
|
|
250
|
+
with pytest.raises(ValueError):
|
|
251
|
+
gas.fit()
|
|
252
|
+
|
|
253
|
+
# This one should not (revert to std)
|
|
254
|
+
gas.model.parametric = False
|
|
255
|
+
gas.model.drop_square = False
|
|
256
|
+
gas.model.degree = 2
|
|
257
|
+
gas.fit()
|
|
258
|
+
|
|
259
|
+
# Now, for predict .................
|
|
260
|
+
prediction = gas.predict(gas_fit_E, stderror=False)
|
|
261
|
+
# This one should fail (extrapolation & blending)
|
|
262
|
+
with pytest.raises(ValueError):
|
|
263
|
+
gas.predict(prediction.values, stderror=False)
|
|
264
|
+
|
|
265
|
+
# But this one should not ..........
|
|
266
|
+
gas.predict(gas_fit_E, stderror=False)
|
|
267
|
+
|
|
268
|
+
def test_pickling(self):
|
|
269
|
+
rs = np.random.RandomState(seed=123)
|
|
270
|
+
n = 500
|
|
271
|
+
np.arange(n)
|
|
272
|
+
x = np.linspace(0, 500, n)
|
|
273
|
+
y = rs.randint(1, 1000, n)
|
|
274
|
+
lo = loess(y, x, span=0.3)
|
|
275
|
+
lo.fit()
|
|
276
|
+
|
|
277
|
+
# pickle lo, unpickle it as lo2 and compare
|
|
278
|
+
# predictions from both
|
|
279
|
+
with tempfile.TemporaryFile() as f:
|
|
280
|
+
pickle.dump(lo, f)
|
|
281
|
+
f.seek(0)
|
|
282
|
+
lo2 = pickle.load(f)
|
|
283
|
+
|
|
284
|
+
npt.assert_array_equal(lo.predict(x).values, lo2.predict(x).values)
|