我正在研究推动欺诈行为的个人特征。 我正在使用具有随机效果的面板概率估算器。
library(pglm)
re.probit <- pglm(fraudulent ~ gender + earn + svo_angle + gender*earn +
risk_av , model=("random"), effect=("twoways"),
index=c("indiv", "Period"), family=binomial(link="probit"),
data=FINAL)
summary(re.probit)
这是输出：
Estimates:
Estimate Std. error t value Pr(> t)
(Intercept) 0.380148 1.008398 0.377 0.70619
gender -0.005578 0.464420 -0.012 0.99042 *
earn -0.335599 1.090525 -0.308 0.75828
svo_angle -0.008280 0.010688 -0.775 0.43855
risk_av -0.209750 0.096782 -2.167 0.03022 *
gender:earn -0.392419 0.676833 -0.580 0.56206 **
sigma 0.841900 0.314763 2.675 0.00748 **
我想测试模型的适合度。函数PseudoR2()执行此操作，但在pglm()函数后它不起作用。你知道一个与pglm()一起使用的类似命令吗？或者如何评估该模型的拟合优度？
然后我有第二个问题。如何将我的模型导出到Latex或excel（我尝试使用Stargazer包，但对pglm()模型也不起作用...）
谢谢你的帮助，
马克