Multi Resolution Fit.
Project description
Multiple Resolution Goodness of Fit.
The thesis of [cuevs2013]_ is that:
"... there is no one 'proper' resolution, but rather a range of
resolutions is necessary to adequately describe the fit of
models with reality."
Usage
#!/usr/bin/env python3
# -*- Coding: UTF-8 -*-
"""
Example from [costanza89]_.
References
----------
.. [costanza89] COSTANZA, Robert. Model goodness of fit: a multiple resolution
procedure. Ecological modelling, v. 47, n. 3-4,
p. 199-215, 1989.
"""
from multiresolutionfit import Multiresoutionfit
from numpy import arange, array
from numpy.random import randint
import matplotlib.pyplot as plt
scene1 = array([[1, 1, 1, 1, 2, 2, 2, 3, 3, 3],
[1, 1, 1, 2, 2, 2, 3, 3, 3, 3],
[1, 1, 2, 2, 2, 3, 3, 3, 3, 3],
[3, 3, 2, 2, 3, 3, 3, 3, 3, 3],
[1, 3, 3, 3, 3, 3, 3, 3, 3, 3],
[1, 1, 1, 3, 3, 3, 3, 3, 3, 3],
[2, 2, 2, 2, 2, 2, 2, 2, 3, 3],
[3, 3, 3, 3, 3, 3, 3, 3, 3, 3],
[3, 3, 3, 3, 2, 2, 3, 3, 3, 3],
[3, 3, 3, 3, 2, 2, 2, 2, 3, 3]
])
scene2 = array([[1, 1, 2, 2, 2, 2, 2, 2, 3, 3],
[1, 1, 1, 1, 2, 3, 3, 3, 3, 3],
[1, 1, 1, 2, 3, 3, 3, 3, 3, 3],
[3, 1, 2, 2, 3, 3, 3, 4, 4, 4],
[3, 3, 3, 3, 3, 3, 3, 3, 3, 3],
[1, 1, 1, 3, 3, 3, 3, 3, 3, 3],
[1, 1, 2, 2, 2, 2, 2, 2, 3, 3],
[1, 2, 2, 3, 3, 2, 2, 3, 3, 3],
[3, 3, 3, 3, 2, 2, 2, 3, 3, 3],
[3, 3, 3, 3, 2, 2, 2, 2, 3, 3]
])
obj = Multiresoutionfit(scene1, scene2, verbose=True)
MAXW = 10
k = 0.1
wins = arange(1, 11, 1, dtype=int)
ftot, fw, wins = obj.ft(k=k, wins=wins)
print(f"\nWeighted fit: {ftot:.2f}\n")
z, fit = obj.zvalue(k=k, wins=wins, permutations=30)
print(f"z value {z:.2f}.")
plt.plot(wins, fw, marker='D')
plt.xticks(wins)
plt.ylim(ymax=0.95, ymin=0.75)
plt.xlim(xmax=MAXW, xmin=1)
plt.grid(True)
plt.show()
References
COSTANZA, Robert. Model goodness of fit: a multiple resolution procedure. Ecological modelling, v. 47, n. 3-4, p. 199-215, 1989.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Close
Hashes for multiresolutionfit-1.0.0-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 5b4c3b53c0d8791c7021c9204a64ee155a81d2675755e9b6bba6ac62f15f4a92 |
|
MD5 | beb4f28e746a116fd26a6b8a60173a8d |
|
BLAKE2b-256 | 15d784ba45152f71f21d1bb964822d0bc6fa9e5b28ebaeaa765dce59cbdabc78 |