[R] wrong signs using MNP.

Kosuke Imai kimai at Princeton.Edu
Sun Apr 3 14:36:56 CEST 2005


Dear Bill,
  You might want to check whether the parameterization of MNP differs from
the book you mention, which I don't have. In particular, I would check how
the choice specific covariate is calculated. More importantly, I would
also make sure that your chain has converged. Based on my experiences of
the multinomial probit model, 1500 draws is not sufficient. Of course,
this depends on one's model, data, and starting values. On this issue, you
might consult the accompanying paper as well as my Journal of Econometrics
article (both available at
http://www.princeton.edu/~kimai/research/MNP.html). The former will walk
you through how one might use MNP and coda packages to conduct convergence
diagnostics using multiple chains via some detailed examples.
Hope this helps,
Kosuke

---------------------------------------------------------
Kosuke Imai               Office: Corwin Hall 041
Assistant Professor       Phone: 609-258-6601 
Department of Politics    eFax:  973-556-1929
Princeton University      Email: kimai at Princeton.Edu
Princeton, NJ 08544-1012  http://www.princeton.edu/~kimai
---------------------------------------------------------

> Hi,
>   I recently found the MNP package.  Out of curiosity, I tried to reproduce
> results from Greene (Econometric Analysis, fourth edition) on page 874.
> 
> The signs of the estimates are all opposite those of Greene's table.  Might
> anyone be able to tell me what I am doing wrong?
> 
> I have attached the function call, the coefficients, and a few rows of the
> data.  The dataset has 210 observations.
> 
> Thanks,
> 
> --Bill
> 
> 
> res1<-mnp(mc ~ 1,
>   choiceX = list(a=cbind(ttme.1,gc.1,hinc),
>                  b=cbind(ttme.2,gc.2,0),
>                  c=cbind(ttme.3,gc.3,0),
>                  d=cbind(ttme.4,gc.4,0)
>             ),
>   cXnames = c("ttme","gc","hinc"),
>   data = wide,
>   n.draws = 1500,
>   burnin = 50,
>   base = "d",
>   verbose = TRUE
> )
> 
> Coefficients:
>                    mean  std.dev.      2.5%  97.5%
> (Intercept):a -1.003078  0.508752 -2.292813 -0.207
> (Intercept):b -1.511963  0.582518 -2.875069 -0.705
> (Intercept):c -1.241368  0.411084 -2.231260 -0.646
> ttme           0.022973  0.007515  0.011386  0.042
> gc             0.002175  0.001981 -0.001518  0.006
> hinc          -0.030616  0.006030 -0.042062 -0.018
> 
> Covariances:
>        mean std.dev.    2.5%  97.5%
> a:a  1.0000   0.0000  1.0000  1.000
> a:b -0.8026   0.4655 -1.8569 -0.252
> a:c -0.7246   0.3721 -1.8966 -0.305
> b:b  1.4315   1.9552  0.1265  7.973
> b:c  1.1526   1.7011  0.1585  7.276
> c:c  1.1281   1.6327  0.2184  7.682
> 
> Number of alternatives: 4
> Number of observations: 210
> Number of Gibbs draws: 1450
> 
> 
> The data:
> 
> mc ttme.1 ttme.2 ttme.3 ttme.4 gc.1 gc.2 gc.3 gc.4 hinc
>  d     69     34     35      0   70   71   70   30   35
>  d     64     44     53      0   68   84   85   50   30
>  d     69     34     35      0  129  195  149  101   40
>  d     64     44     53      0   59   79   81   32   70
>  d     64     44     53      0   82   93   94   99   45
>  b     69     40     35      0   70   57   58   43   20
>  a     45     34     35      0  160  213  167  125   45
>  d     69     34     35      0  137  149  146  135   12
>  d     69     34     35      0   70   71   70   40   40
>  d     69     34     35      0   65   69   68   30   70
>  d     64     44     53      0   68   70   73   36   15
>  d     64     44     53      0   79   90   91   44   35
>  d     64     44     53      0   63   81   83   41   50
>  d     64     44     53      0   72   85   86   58   40
>  d     64     44     53      0  109  214  189  199   26
>  b     69     20     35      0   73   55   72   52   26
> 
> .
> .
> .
> 
> 
> 
> 
>




More information about the R-help mailing list