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)
