如何通过probit模型计算或显示出的逆米尔斯比率

蓝色 发表于 2012-3-9 16:58
1、如果predict 命令 不了解，还是应该先找本书，把基础的知识补充补充吧。这是最基本的命令之一的。否则， ...

我刚才又翻了一下书， predict命令里 只写明xb是模型拟合值，确实没有gw的含义，我看的书是人大版王天夫的《stata实用教程》

sungmoo 发表于 2012-3-9 17:06
*想偷懒的话，可以用以下直接得到变量imr：
heckman y x*, sel(g=w*) two m(imr)
*其中，g是选择方程中的 ...

版主，我刚才看help手册时候明白了偷懒的方法，也试过了这种方法，有米尔斯比的显著性等信息，我不太想偷懒，想看看按照两步走怎么弄

alabala 发表于 2012-3-9 17:12 我刚才又翻了一下书， predict命令里 只写明xb是模型拟合值，确实没有gw的含义，我看的书是人大版王天夫的 ...

help probit_postestimation

上面的gw，就是prob回归中自变量与系数估计值的线性组合。

补充一句：gw换成其他符号也可以。

gw 是自己定义的变量

没有笑话的意思。学软件应该先把基础的有了解。否则，全是问题，别人根本解释不完的

本问题的关键，与使用的软件无关。只是一个计量理论的问题。

help predict                                                                                     dialog:  predict
-----------------------------------------------------------------------------------------------------------------
Title
    [R] predict -- Obtain predictions, residuals, etc., after estimation

Syntax
    After single-equation (SE) models
        predict [type] newvar [if] [in] [, single_options]

    After multiple-equation (ME) models
        predict [type] newvar [if] [in] [, multiple_options]
        predict [type] {stub*|newvar1 ... newvarq} [if] [in] , scores

sungmoo 发表于 2012-3-9 17:23
本问题的关键，与使用的软件无关。只是一个计量理论的问题。

恩，版主说的对，我确实该好好补一下基础理论，谢谢各位老师，蓝色和版主，你们太让人感动了。是要学好计量理论，不然别人的论文都看不懂

蓝色 发表于 2012-3-9 17:24
help predict                                                                                     d ...

结合两位老师的讲解，确实感到理论不行，刚才又翻了一下IMR的计算，修改了一下命令如下，不知对否？
另外请两位老师给推荐本书，我已经看过人大王天夫版的《stata实用教程》
      
我用两种命令得出的结果为什么不一样呢？一个是偷懒的heckman方法，一个用的两步法
  1.      sysuse auto, clear
          probit foreign weight mpg
          predict gw, xb
            g lambda=normalden(gw)/normal(gw) if foreign==1
            reg price weight headroom rep78 lambda if foreign== 1
2. heckman price weight headroom rep78, two sel(foreign= weight mpg)


另外，是否可以在第一步中直接求出针对选择方程中 foreign=0的IMR？
命令为：replace lambda=-normalden(gw)/1-normal(gw) if foreign==0
不知道这样写对吗？

两种结果：
   

Variable m1	m2
_
weight 5.0226671***
headroom -77.976666
rep78 -692.51437
lambda 818.32881
_cons -2651.2977
price
weight	4.6592634***
headroom	-96.762281
rep78	-712.16727
_cons	-1891.3727
foreign
weight	-.00244374***
mpg	-.09083704*
_cons	8.1625526***
mills
lambda	1095.4711
legend: * p<.1; **	p<.05; *** p<.01

蓝色 发表于 2012-3-9 19:07
所以你需要看计量的书，wooldridge的书P582 中文
理论不懂软件的结果就是错误的

看了您给的书上的介绍，我马上找这本书，刚接触计量，真是博大精深，很感兴趣，谢谢您的指点

replace lambda=-normalden(gw)/1-normal(gw) if foreign==0

1-normal(gw)=normal(-gw)