コード例 #1
0
ファイル: M_semiGWR.py プロジェクト: TaylorOshan/GIS596_GWR
 def __init__(self, y, x_global, x_local, kernel, mType=0, offset=None, sigma2_v1=False, tol=1.0e-6, maxIter=200, y_name="", y_off_name="", x_name_glob=[], x_name_loc=[], fle_name="", summaryGLM=False):
     """
     Initialize class
     """
     semiGWR_Base.__init__(self, y, x_global, x_local, kernel, mType, offset, sigma2_v1, tol, maxIter)     
     self.y_name = y_name
     self.x_name_glob = x_name_glob
     self.x_name_loc = x_name_loc
     if y_off_name is None:
         y_off_name = " "
     self.y_off_name = y_off_name
     
     # add x names-global
     n_xname = len(self.x_name_glob)
     if n_xname < self.nVars_glob:
         for i in range(n_xname,self.nVars_glob):
             self.x_name_glob.append("name" + str(i))
     # add x names-local        
     n_xname = len(self.x_name_loc)
     if n_xname < self.nVars_loc:
         for i in range(n_xname,self.nVars_loc-1):
             self.x_name_loc.append("name" + str(i))          
     if len(self.x_name_loc) < self.nVars_loc:
         self.x_name_loc.insert(0,'Intercept')
     self.fle_name = fle_name    
     self.mName = semiGWR_Names[self.mType]
     
     if summaryGLM:
         x = np.hstack((x_global,x_local)) 
         x_name = [ elem for elem in x_name_glob]
         x_name.extend(x_name_loc)
         #if mType == 0:
             #self.OLS = OLS(y,x,sigma2_v1,y_name,x_name,fle_name)
         #else:
         self.GLM = GLM(y, x, mType, offset, sigma2_v1, tol, maxIter, y_name, y_off_name, x_name, fle_name)        
     
     Summary.semiGWR(self)