__author__ = "Xun Li <lixun910@gmail.com>"
from ..libgeoda import gda_localg, gda_localgstar, VecBool
from .lisa import lisa
import multiprocessing
[docs]def local_g(w, data, **kwargs):
"""Local Getis-Ord's G Statistics
The function to apply Getis-Ord's local G statistics
Args:
w (Weight): An instance of Weight class.
data (tuple/list/pandas.Series): A list of numeric values of selected variable
undefs (list, optional): A list of boolean values to indicate which value is undefined or null
permutations (int, optional): The number of permutations for the LISA computation
permutation_method (str, optional): The permutation method used for the LISA computation. Options are {'complete', 'lookup-table'}. Default is 'complete'.
significance_cutoff (float, optional): A cutoff value for significance p-values to filter not-significant clusters
cpu_threads (int, optional): The number of cpu threads used for parallel LISA computation
seed (int, optional): The seed for random number generator
Returns:
lisa: An instance of lisa class
"""
if w == None:
raise ValueError("Weights is None.")
undefs = VecBool() if 'undefs' not in kwargs else kwargs['undefs']
significance_cutoff = 0.05 if 'significance_cutoff' not in kwargs else kwargs['significance_cutoff']
permutations = 999 if 'permutations' not in kwargs else kwargs['permutations']
permutation_method = 'complete' if 'permutation_method' not in kwargs else kwargs['permutation_method']
cpu_threads = 6 if 'cpu_threads' not in kwargs else kwargs['cpu_threads']
seed = 123456789 if 'seed' not in kwargs else kwargs['seed']
lisa_obj = gda_localg(w.gda_w, list(data), list(undefs), significance_cutoff, cpu_threads, permutations, permutation_method, seed)
return lisa(lisa_obj)
[docs]def local_gstar(w, data, **kwargs):
"""Local Getis-Ord's G* Statistics
The function to apply Getis-Ord's local G* statistics
Args:
w (Weight): An instance of Weight class.
data (tuple/list/pandas.Series): A list of numeric values of selected variable
undefs (list, optional): A list of boolean values to indicate which value is undefined or null
permutations (int, optional): The number of permutations for the LISA computation
permutation_method (str, optional): The permutation method used for the LISA computation. Options are {'complete', 'lookup-table'}. Default is 'complete'.
significance_cutoff (float, optional): A cutoff value for significance p-values to filter not-significant clusters
cpu_threads (int, optional): The number of cpu threads used for parallel LISA computation
seed (int, optional): The seed for random number generator
Returns:
lisa: An instance of lisa class
"""
if w == None:
raise ValueError("Weights is None.")
undefs = VecBool() if 'undefs' not in kwargs else kwargs['undefs']
significance_cutoff = 0.05 if 'significance_cutoff' not in kwargs else kwargs['significance_cutoff']
permutations = 999 if 'permutations' not in kwargs else kwargs['permutations']
permutation_method = 'complete' if 'permutation_method' not in kwargs else kwargs['permutation_method']
cpu_threads = 6 if 'cpu_threads' not in kwargs else kwargs['cpu_threads']
seed = 123456789 if 'seed' not in kwargs else kwargs['seed']
lisa_obj = gda_localgstar(w.gda_w, list(data), list(undefs), significance_cutoff, cpu_threads, permutations, permutation_method, seed)
return lisa(lisa_obj)