Skip to content

Commit

Permalink
keep table values for nonsig observations in localknox hotspot (#136)
Browse files Browse the repository at this point in the history
  • Loading branch information
knaaptime committed Jun 14, 2024
1 parent 8ea6cbd commit 261768a
Showing 1 changed file with 9 additions and 1 deletion.
10 changes: 9 additions & 1 deletion pointpats/spacetime.py
Original file line number Diff line number Diff line change
Expand Up @@ -28,6 +28,7 @@
import libpysal as lps
import numpy as np
import pandas
import pandas as pd
import scipy.stats as stats
from libpysal import cg
from libpysal.graph import Graph
Expand Down Expand Up @@ -1352,6 +1353,13 @@ def hotspots(self, crit=0.05, inference="permutation", keep_neighbors=True):
(self.adjlist.focal.isin(pdf_sig.index.values))
| self.adjlist.neighbor.isin(pdf_sig.index.values)
]
pdf_sig = pd.concat([pdf_sig,
self._gdf[self._gdf.index.isin(temp_neighbors.neighbor.values)

][[col, "time"]].rename(
columns={col: "pvalue", "time": "focal_time"}
)
])

pdf_sig = pdf_sig.merge(
temp_neighbors, how='outer', left_index=True, right_on="focal"
Expand All @@ -1361,7 +1369,7 @@ def hotspots(self, crit=0.05, inference="permutation", keep_neighbors=True):
graph = Graph.from_adjacency(pdf_sig.assign(weight=1))
pdf_sig["cluster"] = graph.component_labels.values
if not keep_neighbors :
pdf_sig = pdf_sig.dropna(subset=['pvalue'])
pdf_sig = pdf_sig[pdf_sig.pvalue<=crit]

return self._gdf[["geometry"]].merge(
pdf_sig.copy(), left_index=True, right_on="focal"
Expand Down

0 comments on commit 261768a

Please sign in to comment.