def test_knox(self): result = interaction.knox(self.events.space, self.events.t, delta=20, tau=5, permutations=1) self.assertEqual(result['stat'], 13.0)
import pysal import pysal.spatial_dynamics.interaction as interaction np.random.seed(100) # In[14]: ###knox running path = "C:/Users/titog/Desktop/city lab/knox test/files/time" events = SpaceTimeEvents(path,'time') # In[30]: result = interaction.knox(events.space, events.t ,delta=600,tau=1,permutations=99) # In[31]: print(result['stat']) # In[27]: print("%2.8f"%result['pvalue']) # In[113]: """
def test_knox(self): result = interaction.knox(self.events.space, self.events.t, delta=20, tau=5, permutations=1) self.assertEquals(result["stat"], 13.0)
provider = layer.dataProvider() ### Input for attribute table fields fields = provider.fields() fields.append(QgsField('STAT', QVariant.Float)) fields.append(QgsField('PVALUE', QVariant.Float)) fields.append(QgsField('COUNTS', QVariant.Int)) writer = self.getOutputFromName(self.OUTPUT).getVectorWriter( fields, provider.geometryType(), layer.crs() ) events = interaction.SpaceTimeEvents(self.INPUT,'T') f = pysal.open(filename.replace('.shp','.dbf')) y=np.array(f.by_col[str(field)]) kt = interaction.knox(events.space, events.t, delta=20, tau=5, permutations=99) self.setOutputValue(self.P_SIM, str(kt.p_sim)) outFeat = QgsFeature() i = 0 for inFeat in processing.features(layer): inGeom = inFeat.geometry() outFeat.setGeometry(inGeom) attrs = inFeat.attributes() attrs.append(float(kt.p_sim[i])) outFeat.setAttributes(attrs) writer.addFeature(outFeat) i+=1 del writer
def test_knox(self): result = interaction.knox( self.events, delta=20, tau=5, permutations=1) self.assertEquals(result['stat'], 13.0)