我在SAS和R中运行了相同的概率回归,而我的系数估计值(基本上)是等效的,报告的测试统计数据是不同的.具体而言,SAS将测试统计报告为t统计,而R将测试统计报告为z统计.
我检查了我的计量经济学文本,发现(很少详细说明)它以t统计数据报告了probit结果.
哪个统计数据合适？为什么R与SAS不同？
这是我的SAS代码:
proc qlim data=DavesData;
model y = x1 x2 x3/ discrete(d=probit);
run;
quit;
这是我的R代码:
> model.1 <- glm(y ~ x1 + x2 + x3, family=binomial(link="probit"))
> summary(model.1)