Python实现地理加权回归的代码

以下是Python实现地理加权回归的代码：

import numpy as np
import pysal as ps
import statsmodels.api as sm

# 加载数据
db = ps.open('data.dbf', 'r')
y_var = 'crime'
x_vars = ['ind', 'pop']

# 构建空间权重矩阵
w = ps.weights.Queen.from_shapefile('data.shp')
w.transform = 'R'

# 构建地理加权回归模型
model = ps.spreg.GM_Lag(y=db.by_col(y_var), x=db.by_col(x_vars), w=w)
print(model.summary)