全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
16072 10
2015-01-21
各位,我做logit模型(有5个显著),回归后用mfx求边际效应,发现没有一个显著的;但是我又用margins求边际效应却发现出现了5个变量显著,请问这个该怎么办?是继续用mfx还是用margins?谢谢!
二维码

扫码加我 拉你入群

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

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

全部回复
2015-1-21 19:54:43
没有什么不一样的啊

. webuse lbw
(Hosmer & Lemeshow data)

. logit low age lwt  smoke ptl ht ui

Iteration 0:   log likelihood =   -117.336  
Iteration 1:   log likelihood =  -104.7002  
Iteration 2:   log likelihood = -104.39685  
Iteration 3:   log likelihood = -104.39591  
Iteration 4:   log likelihood = -104.39591  

Logistic regression                               Number of obs   =        189
                                                  LR chi2(6)      =      25.88
                                                  Prob > chi2     =     0.0002
Log likelihood = -104.39591                       Pseudo R2       =     0.1103

------------------------------------------------------------------------------
         low |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         age |  -.0422544   .0345836    -1.22   0.222     -.110037    .0255281
         lwt |  -.0142885   .0066523    -2.15   0.032    -.0273267   -.0012502
       smoke |   .5506312   .3436293     1.60   0.109    -.1228697    1.224132
         ptl |   .5932548   .3484191     1.70   0.089     -.089634    1.276144
          ht |   1.862491   .6862291     2.71   0.007     .5175064    3.207475
          ui |   .7367904   .4564882     1.61   0.107      -.15791    1.631491
       _cons |    1.37896   1.088893     1.27   0.205    -.7552312    3.513151
------------------------------------------------------------------------------

. mfx

Marginal effects after logit
      y  = Pr(low) (predict)
         =  .28930847
------------------------------------------------------------------------------
variable |      dy/dx    Std. Err.     z    P>|z|  [    95% C.I.   ]      X
---------+--------------------------------------------------------------------
     age |  -.0086879      .00708   -1.23   0.220  -.022558  .005183   23.2381
     lwt |  -.0029378      .00135   -2.18   0.029  -.005579 -.000297    129.82
   smoke*|   .1156273      .07302    1.58   0.113  -.027487  .258742   .391534
     ptl |   .1219786      .07191    1.70   0.090  -.018965  .262922   .195767
      ht*|   .4340082      .14329    3.03   0.002   .153163  .714854   .063492
      ui*|   .1652437      .10863    1.52   0.128  -.047676  .378164   .148148
------------------------------------------------------------------------------
(*) dy/dx is for discrete change of dummy variable from 0 to 1


. margins,dydx(*) atmean

Conditional marginal effects                      Number of obs   =        189
Model VCE    : OIM

Expression   : Pr(low), predict()
dy/dx w.r.t. : age lwt smoke ptl ht ui
at           : age             =     23.2381 (mean)
               lwt             =    129.8201 (mean)
               smoke           =    .3915344 (mean)
               ptl             =    .1957672 (mean)
               ht              =    .0634921 (mean)
               ui              =    .1481481 (mean)

------------------------------------------------------------------------------
             |            Delta-method
             |      dy/dx   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         age |  -.0086879    .007077    -1.23   0.220    -.0225585    .0051827
         lwt |  -.0029378   .0013475    -2.18   0.029    -.0055789   -.0002968
       smoke |   .1132148   .0703391     1.61   0.107    -.0246472    .2510768
         ptl |   .1219786   .0719113     1.70   0.090     -.018965    .2629222
          ht |    .382945     .14083     2.72   0.007     .1069232    .6589668
          ui |   .1514908    .093994     1.61   0.107    -.0327341    .3357157
------------------------------------------------------------------------------



二维码

扫码加我 拉你入群

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

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

2015-1-22 09:48:34
耶,那我为什么两者不一样呢?用mfx都没有显著的,而用margins却有好几个显著,是怎么回事,没有搞明白
二维码

扫码加我 拉你入群

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

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

2015-1-22 13:25:14
应该一样才对
二维码

扫码加我 拉你入群

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

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

2016-3-6 09:42:45
蓝色 发表于 2015-1-22 13:25
应该一样才对
我想请教一下,当因变量是四分类定序变量(ologit)的时候,用margins dydx(*)求边际效应时,是不是跟mlogit的求法是一样的呀,就是说还是得一个一个的求,是吗?
二维码

扫码加我 拉你入群

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

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

2016-3-6 09:46:34
sunshinewp 发表于 2016-3-6 09:42
我想请教一下,当因变量是四分类定序变量(ologit)的时候,用margins dydx(*)求边际效应时,是不是跟mlogi ...
Title

    [R] ologit postestimation -- Postestimation tools for ologit


Description

    The following postestimation commands are available after ologit:

    Command              Description
    ----------------------------------------------------------------------------------------------------------------------
        contrast         contrasts and ANOVA-style joint tests of estimates
        estat ic         Akaike's and Schwarz's Bayesian information criteria (AIC and BIC)
        estat summarize  summary statistics for the estimation sample
        estat vce        variance-covariance matrix of the estimators (VCE)
        estat (svy)      postestimation statistics for survey data
        estimates        cataloging estimation results
    (1) forecast         dynamic forecasts and simulations
        lincom           point estimates, standard errors, testing, and inference for linear combinations of coefficients
        linktest         link test for model specification
    (2) lrtest           likelihood-ratio test
        margins          marginal means, predictive margins, marginal effects, and average marginal effects
        marginsplot      graph the results from margins (profile plots, interaction plots, etc.)
        nlcom            point estimates, standard errors, testing, and inference for nonlinear combinations of
                           coefficients
        predict          predictions, residuals, influence statistics, and other diagnostic measures
        predictnl        point estimates, standard errors, testing, and inference for generalized predictions
        pwcompare        pairwise comparisons of estimates
        suest            seemingly unrelated estimation
        test             Wald tests of simple and composite linear hypotheses
        testnl           Wald tests of nonlinear hypotheses
    ----------------------------------------------------------------------------------------------------------------------
    (1) forecast is not appropriate with mi or svy estimation results.
    (2) lrtest is not appropriate with svy estimation results.


Syntax for predict

        predict [type] {stub* | newvar | newvarlist} [if] [in] [, statistic outcome(outcome) nooffset]

        predict [type] {stub* | newvarlist} [if] [in] , scores

    statistic          Description
    ----------------------------------------------------------------------------------------------------------------------
    Main
      pr               predicted probabilities; the default
      xb               linear prediction
      stdp             standard error of the linear prediction
    ----------------------------------------------------------------------------------------------------------------------
    If you do not specify outcome(), pr (with one new variable specified) assumes outcome(#1).
    You specify one or k new variables with pr, where k is the number of outcomes.
    You specify one new variable with xb and stdp.
    These statistics are available both in and out of sample; type predict ... if e(sample) ... if wanted only for the
      estimation sample.


Menu for predict

    Statistics > Postestimation > Predictions, residuals, etc.


Options for predict

        +------+
    ----+ Main +----------------------------------------------------------------------------------------------------------

    pr, the default, calculates the predicted probabilities.  If you do not also specify the outcome() option, you specify
        k new variables, where k is the number of categories of the dependent variable.  Say that you fit a model by
        typing ologit result x1 x2, and result takes on three values.  Then you could type predict p1 p2 p3 to obtain all
        three predicted probabilities.  If you specify the outcome() option, you must specify one new variable.  Say that
        result takes on the values 1, 2, and 3.  Typing predict p1, outcome(1) would produce the same p1.

    xb calculates the linear prediction.  You specify one new variable, for example, predict linear, xb.  The linear
        prediction is defined, ignoring the contribution of the estimated cutpoints.

    stdp calculates the standard error of the linear prediction.  You specify one new variable, for example, predict se,
        stdp.

    outcome(outcome) specifies for which outcome the predicted probabilities are to be calculated.  outcome() should
        contain either one value of the dependent variable or one of #1, #2, ..., with #1 meaning the first category of
        the dependent variable, #2 meaning the second category, etc.

    nooffset is relevant only if you specified offset(varname) for ologit.  It modifies the calculations made by predict
        so that they ignore the offset variable; the linear prediction is treated as xb rather than as xb + offset.

    scores calculates equation-level score variables.  The number of score variables created will equal the number of
        outcomes in the model.  If the number of outcomes in the model was k, then

        The first new variable will contain the derivative of the log likelihood with respect to the regression equation.

        The other new variables will contain the derivative of the log likelihood with respect to the cutpoints.


Examples

    Setup
        . webuse fullauto
        . ologit rep77 i.foreign length mpg

    Predicted probabilities for each of the five outcomes
        . predict poor fair avg good exc

    Average marginal effects on the probability of an excellent repair record
        . margins, dydx(*) predict(outcome(5))


    Report information criteria
        . estat ic
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群