def test_jacquez(self):
result = interaction.jacquez(self.events.space,
self.events.t, k=3, permutations=1)
if scp_version > 11:
self.assertEquals(result['stat'], 12)
else:
self.assertEquals(result['stat'], 13)
def test_jacquez(self):
result = interaction.jacquez(self.events.space,
self.events.t, k=3, permutations=1)
if scp_version > 11:
self.assertEquals(result['stat'], 12)
else:
self.assertEquals(result['stat'], 13)
filename = self.getParameterValue(self.INPUT)
layer = dataobjects.getObjectFromUri(filename)
filename = dataobjects.exportVectorLayer(layer)
provider = layer.dataProvider()
fields = provider.fields()
fields.append(QgsField('STAT', QVariant.Float))
fields.append(QgsField('PVALUE', QVariant.Float))
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)])
jt = interaction.jacquez(events.space, events.t ,k=3,permutations=99)
self.setOutputValue(self.P_SIM, str(jt.p_sim))
outFeat = QgsFeature()
i = 0
for inFeat in processing.features(layer):
inGeom = inFeat.geometry()
outFeat.setGeometry(inGeom)
attrs = inFeat.attributes()
attrs.append(float(jt.p_sim[i]))
outFeat.setAttributes(attrs)
writer.addFeature(outFeat)
i+=1
del writer
def test_jacquez(self):
result = interaction.jacquez(self.events.space,
self.events.t, k=3, permutations=1)
self.assertEquals(result['stat'], 13)
