非线性似不相关回归nlsur中如何得到对数行列式log-det值？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

最近在用nlsur估计参数，发现有些文章中给出了log-det值，这个值在stata中要如何得到呢？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：NLSUR 非线性 SUR nls Log

楼主知道怎么得到了吗

同问，请问楼主知道了吗

求助！！！

capt mat drop A  
mat A=e(V)   
dis log(det(A))

阿蛋在呢 发表于 2022-12-13 15:19
capt mat drop A  
mat A=e(V)   
dis log(det(A))

请问一下log_det值有啥评判标准吗？越大越好吗？如果为负数咋办

在Stata中使用`nlsur`命令进行非线性似不相关回归（Nonlinear Seemingly Unrelated Regression）后获取对数行列式（log-determinant）的值，并不是一个直接输出的结果。这是因为`nlsur`主要关注于参数估计和标准误等，而不是直接提供协方差矩阵或其对数行列式的值。

但是，你可以通过以下步骤计算得到该值：

1. 首先运行你的`nlsur`模型。
2. 使用命令`estat vce`来显示模型的估计误差方差-协方差矩阵（VCE）。
3. 然后使用矩阵操作命令如`mat`和`ln`来计算该矩阵的行列式并取对数。

具体来说，你可以按照以下步骤进行：

```
nlsur (eq1: y1 = {b0} + {b1}*x1, nograd) (eq2: y2 = {c0} + {c1}*x2), ///
       tech(bhhh nr) iterate(50)

estat vce

mat V = e(V)
scalar log_det = ln(det(V))
di "Log-Determinant is: `log_det'"
```

但是需要注意的是，`nlsur`所使用的VCE通常是估计的，并且可能经过某些调整（如White标准误），因此计算行列式时使用原始矩阵可能不完全反映最终模型的实际精度。此外，在处理大规模数据或复杂模型时，协方差矩阵可能会变得不稳定或奇异，这可能会影响对数行列式的准确性。

在学术论文中报告log-det值通常是为了说明模型的某些特定性质（例如，信息量、模型的稳定性等），但这一数值并不是`nlsur`模型估计过程中的直接输出。因此，在实际应用时，应根据具体研究目的来判断是否需要计算并报告此值。
  
请注意，上述操作在Stata 13以上的版本中应该可以正常运行。如果你使用的Stata版本较低，命令和函数的可用性可能会有所不同。

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