def GLM(GLMMod):
"""
get diagnostics for GLM model
"""
if GLMMod.mType == 0:
GLMMod.t_stat = Diagnostics.tstat_GLM(GLMMod)
GLMMod.f_stat = Diagnostics.fstat_OLS(GLMMod)
GLMMod.r2 = Diagnostics.r2_OLS(GLMMod)
GLMMod.ar2 = Diagnostics.ar2_OLS(GLMMod)
GLMMod.cv = Diagnostics.get_CV_OLS(GLMMod)
GLMMod.sig2ML = GLMMod.sigma2_n
GLMMod.mulColli = Diagnostics.ci_OLS(GLMMod)
#summary_OLS(GLMMod)
else:
GLMMod.dev_null = Diagnostics.dev_mod_GLM(GLMMod)
GLMMod.pdev = 1.0 - GLMMod.dev_res/GLMMod.dev_null
GLMMod.t_stat = Diagnostics.tstat_GLM(GLMMod, True)
#GLMMod.std_err = Diagnostics.se_betas_GLM(GLMMod)
GLMMod.aic = Diagnostics.get_AIC_GLM(GLMMod)
GLMMod.aicc = Diagnostics.get_AICc_GLM(GLMMod)
GLMMod.bic = Diagnostics.get_BIC_GLM(GLMMod)
summary_GLM(GLMMod)
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)
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)
