请教下stata16构建空间权重矩阵不同维错误的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如附件图所示，在构建31*31的空间权重矩阵后，和原面板数据求莫兰指数，但出现了不同维的错误，请教一下怎么解决？
命令代码：spatgsa upgrade,weights(W) moran
错误提醒：Matrix W is 31x31, the dataset in use has 496 obs.
To run -spatgsa- weights matrix dimension must equal N. of obs

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

moran只能做截面

一年一年的可以

一知半解统计人 发表于 2020-3-18 20:42
一年一年的可以

好的，现在能做出来了，谢谢

遇到这个问题是因为空间权重矩阵（W）的维度与你的数据集观测值数量不匹配。在你的情况下，W是一个31x31的矩阵，而数据集有496个观测值。这通常意味着W矩阵是根据一个包含较少单元格或观察点的数据集构建的。

解决这个问题的方法主要有以下几点：

### 1. 确认空间权重矩阵（W）是否正确构建

- 首先检查你的W矩阵构建过程，确保它是基于与你当前使用的数据集相同的空间单位。例如，如果数据是31个省份的时间序列数据，则W矩阵应该是这31个省份之间的空间权重。

### 2. 数据集的观测值数量问题

- 如果W矩阵仅包含空间上的信息（如省份间的关系），而你的面板数据包含了时间和空间的信息，那么每个时间点上你都有31个观测值。这种情况下，可能需要在运行`spatgsa`命令之前对数据进行处理。

### 3. 使用面板数据构建W矩阵

- 如果你打算基于包含时间信息的完整面板数据集来构建W矩阵，则需要先聚合数据到一个只包括空间单位（如省份）的横截面数据集中，然后在此基础上构建W。完成后，再将W与原始的面板数据结合。

### 4. 数据重组或汇总

- 另一种方法是考虑是否可以按时间平均你的数据或者选择特定的时间点来创建一个符合W矩阵维度的数据集。

在具体操作上，你可以尝试以下步骤：

```stata
* 假设省份ID是"prov_id"
* 时间变量为"time"

// 将面板数据汇总到横截面级别以构建W
egen prov_data = sum(your_variable), by(prov_id)
bysort prov_id: keep if _n == 1 // 或者选择特定的时间点

// 构建W矩阵
spwcreate W, idvar(prov_id) wtype(q)

// 然后使用这个W矩阵与你的原始数据计算莫兰指数，确保在执行`spatgsa upgrade`命令时数据集是按时间展开的面板数据。
```

关键在于确保你的空间权重矩阵和你打算进行空间分析的数据集在单位上一致。希望这能帮助到你！

此文本由CAIE学术大模型生成，添加下方二维码，优先体验功能试用