예제 #1
0
def semiGWR(semiGWRMod):
    """
    get diagnostics for semiGWR model
    """
    
     
    if semiGWRMod.mType == 0: # Gaussian model
        semiGWRMod.aic = Diagnostics.get_AIC_GWR(semiGWRMod)
        semiGWRMod.aicc = Diagnostics.get_AICc_GWR(semiGWRMod) 
        semiGWRMod.bic = Diagnostics.get_BIC_GWR(semiGWRMod) 
        semiGWRMod.cv = Diagnostics.get_CV_GWR(semiGWRMod)
        semiGWRMod.R2 = Diagnostics.r2_GWR(semiGWRMod)
        semiGWRMod.R2_adj = Diagnostics.ar2_GWR(semiGWRMod)
    else:
        semiGWRMod.aic = Diagnostics.get_AIC_GWGLM(semiGWRMod)
        semiGWRMod.aicc = Diagnostics.get_AICc_GWGLM(semiGWRMod) 
        semiGWRMod.bic = Diagnostics.get_BIC_GWGLM(semiGWRMod) 
        semiGWRMod.dev_null = Diagnostics.dev_mod_GLM(semiGWRMod)
        semiGWRMod.pdev = 1.0 - semiGWRMod.dev_res/semiGWRMod.dev_null 

    summary_semiGWR(semiGWRMod)
예제 #2
0
def GWGLM(GWRMod):
    """
    get diagnostics for GWGLM model
    """
    if GWRMod.mType == 0:
        GWRMod.aic = Diagnostics.get_AIC_GWR(GWRMod)
        GWRMod.aicc = Diagnostics.get_AICc_GWR(GWRMod)
        GWRMod.bic = Diagnostics.get_BIC_GWR(GWRMod) 
        GWRMod.cv = Diagnostics.get_CV_GWR(GWRMod)
        GWRMod.R2 = Diagnostics.r2_GWR(GWRMod)
        GWRMod.R2_adj = Diagnostics.ar2_GWR(GWRMod)
    else:
        GWRMod.aic = Diagnostics.get_AIC_GWGLM(GWRMod)
        GWRMod.aicc = Diagnostics.get_AICc_GWGLM(GWRMod) 
        GWRMod.bic = Diagnostics.get_BIC_GWGLM(GWRMod) 
        #GWGLMMod.dev_res = Diagnostics.dev_res_GWGLM(GWGLMMod)
        GWRMod.dev_null = Diagnostics.dev_mod_GLM(GWRMod)
        GWRMod.pdev = 1.0 - GWRMod.dev_res/GWRMod.dev_null    
    

    summary_GWGLM(GWRMod)