There are a lot of errors with the local function

[vangjunxiu]

---

model = Gravity(flows, Oi, Dj, Dij, 'exp') local= model.local(loc_index=Origin, locs=np.unique(Origin))

Exception Traceback (most recent call last) Input In [19], in 1 model = Gravity(flows, Oi, Dj, Dij, 'exp') ----> 2 local= model.local(loc_index=Origin, locs=np.unique(Origin))

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\gravity.py:478, in Gravity.local(self, loc_index, locs) 476 d_vars = self.dv[subset.reshape(self.dv.shape[0]),:] 477 dij = self.reshape(self.c[subset]) --> 478 model = Gravity(f, o_vars, d_vars, dij, self.cf, 479 constant=False) 480 results['AIC'].append(model.AIC) 481 results['deviance'].append(model.deviance)

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\gravity.py:429, in Gravity.init(self, flows, o_vars, d_vars, cost, cost_func, constant, framework, SF, CD, Lag, Quasi) 426 self.c = np.reshape(cost, (-1,1)) 427 #User.check_arrays(self.f, self.ov, self.dv, self.c) --> 429 BaseGravity.init(self, self.f, self.c, 430 cost_func=cost_func, o_vars=self.ov, d_vars=self.dv, constant=constant, 431 framework=framework, SF=SF, CD=CD, Lag=Lag, Quasi=Quasi)

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\gravity.py:245, in BaseGravity.init(self, flows, cost, cost_func, o_vars, d_vars, origins, destinations, constant, framework, SF, CD, Lag, Quasi) 243 if (framework.lower() == 'glm'): 244 if not Quasi: --> 245 results = self.fit(framework='glm') 246 else: 247 results = self.fit(framework='glm', Quasi=True)

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\count_model.py:87, in CountModel.fit(self, framework, Quasi) 85 if (framework.lower() == 'glm'): 86 if not Quasi: ---> 87 results = GLM(self.y, self.X, family = Poisson(), constant=self.constant).fit() 88 else: 89 results = GLM(self.y, self.X, family = QuasiPoisson(), constant=self.constant).fit()

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\glm\glm.py:81, in GLM.init(self, y, X, family, offset, y_fix, constant) 76 def init(self, y, X, family=family.Gaussian(), offset=None, y_fix = None, 77 constant=True): 78 """ 79 Initialize class 80 """ ---> 81 self.n = USER.check_arrays(y, X) 82 USER.check_y(y, self.n) 83 self.y = y

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\spreg\user_output.py:358, in check_arrays(*arrays) 356 raise Exception("all input arrays must have exactly two dimensions") 357 if shape[0] < shape[1]: --> 358 raise Exception("one or more input arrays have more columns than rows") 359 if not spu.spisfinite(i): 360 raise Exception("one or more input arrays have missing/NaN values")

Exception: one or more input arrays have more columns than rows

---

ZeroDivisionError Traceback (most recent call last) Input In [5], in 1 model = Gravity(flows, Oi, Dj, Dij, 'exp') ----> 2 local= model.local(loc_index=Origin, locs=np.unique(Origin))

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\gravity.py:478, in Gravity.local(self, loc_index, locs) 476 d_vars = self.dv[subset.reshape(self.dv.shape[0]),:] 477 dij = self.reshape(self.c[subset]) --> 478 model = Gravity(f, o_vars, d_vars, dij, self.cf, 479 constant=False) 480 results['AIC'].append(model.AIC) 481 results['deviance'].append(model.deviance)

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\gravity.py:429, in Gravity.init(self, flows, o_vars, d_vars, cost, cost_func, constant, framework, SF, CD, Lag, Quasi) 426 self.c = np.reshape(cost, (-1,1)) 427 #User.check_arrays(self.f, self.ov, self.dv, self.c) --> 429 BaseGravity.init(self, self.f, self.c, 430 cost_func=cost_func, o_vars=self.ov, d_vars=self.dv, constant=constant, 431 framework=framework, SF=SF, CD=CD, Lag=Lag, Quasi=Quasi)

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\gravity.py:245, in BaseGravity.init(self, flows, cost, cost_func, o_vars, d_vars, origins, destinations, constant, framework, SF, CD, Lag, Quasi) 243 if (framework.lower() == 'glm'): 244 if not Quasi: --> 245 results = self.fit(framework='glm') 246 else: 247 results = self.fit(framework='glm', Quasi=True)

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\count_model.py:90, in CountModel.fit(self, framework, Quasi) 88 else: 89 results = GLM(self.y, self.X, family = QuasiPoisson(), constant=self.constant).fit() ---> 90 return CountModelResults(results) 92 else: 93 raise NotImplemented('Poisson GLM is the only count model currently implemented')

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\spint\count_model.py:186, in CountModelResults.init(self, results) 184 self.tvalues = results.tvalues 185 self.D2 = results.D2 --> 186 self.adj_D2 = results.adj_D2 187 self.pseudoR2 = results.pseudoR2 188 self.adj_pseudoR2 = results.adj_pseudoR2

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\glm\utils.py:370, in CachedAttribute.get(self, obj, type) 367 # print("[_cachedval=%s]" % _cachedval) 368 if _cachedval is None: 369 # Call the "fget" function --> 370 _cachedval = self.fget(obj) 371 # Set the attribute in obj 372 # print("Setting %s in cache to %s" % (name, _cachedval)) 373 try:

File D:\Anaconda3\envs\tensorflow\lib\site-packages\pysal\contrib\glm\glm.py:370, in GLMResults.adj_D2(self) 368 @cache_readonly 369 def adj_D2(self): --> 370 return 1.0 - (float(self.n) - 1.0)/(float(self.n) - float(self.k)) * (1.0-self.D2)

ZeroDivisionError: float division by zero

[ljwolf]

Can you post your data? These errors seem to be related to the structure of your input data.

[vangjunxiu]

I can give you my data, hope you can point me to solve these problems, thank you.

[vangjunxiu]

I have sent the data to your mailbox，thank you.

[ljwolf]

I uploaded the data here. I will look into this as best I can in the coming week.

Workday_OD.csv