[R] 2SLS for panel data, re

Millo Giovanni Giovanni_Millo at Generali.com
Fri Nov 15 17:31:28 CET 2013


Dear Chanita,

impossible to tell without a reproducible example. You do not even
include your Stata call.

Assuming you meant 'ivreg2' w/o "t", there are four rows of possible
arguments to it in the help page, but I don't seem to find any switch
for random effects. Are you sure you are not comparing a RE model with a
pooled one?

Best wishes,
Giovanni

Giovanni Millo, PhD
Research Dept.,
Assicurazioni Generali SpA
Via Machiavelli 3,
34132 Trieste (Italy)
tel. +39 040 671184
fax  +39 040 671160


-------------- original message ------------------

Date: Thu, 14 Nov 2013 08:12:33 -0800
From: Chanita Holmes <chanita.holmes at gmail.com>
To: r-help at r-project.org
Subject: [R] 2SLS for panel data, re
Message-ID:
	
<CAEaD0=6SkDOTODbPZZ19LE2ca6YxByXWp+q-GiL1g+8eFB22DA at mail.gmail.com>
Content-Type: text/plain

Hi,

I am trying to estimate a 2sls using panel data (random effect model). I
tried the same estimation in STATA using the ivtreg2 command. However
STATA
and R are giving me two different results. I figure there is something
with
my R code:

iv=plm(formula=wecon~fdistockgdp +trade + polrightsreversed +lnrgdpch +
execleft + muslim2+c100rat +c111rat +yeardum|
polrightsreversed+lnrgdpch+
execleft+muslim2+c100rat+c111rat+yeardum
+lnpop+lnarea+devcountrycomlanguage+bitcum,
data = women, index = c("country", "year"), random.method = c("swar"),
inst.method = c("bvk"), model="random")
summary(iv)

Coefficients :
                    Estimate Std. Error t-value  Pr(>|t|)
(Intercept)       -0.2258528  0.2951301 -0.7653 0.4441954
fdistockgdp       -0.0067207  0.0077315 -0.8693 0.3847993
trade              0.0068462  0.0023687  2.8903 0.0038863 **
polrightsreversed  0.0092366  0.0106174  0.8699 0.3844229
lnrgdpch           0.1246679  0.0389043  3.2045 0.0013724 **
execleft           0.1118046  0.0340817  3.2805 0.0010524 **
muslim2           -0.0044742  0.0012433 -3.5986 0.0003270 ***
c100rat            0.0226208  0.0595134  0.3801 0.7039114
c111rat            0.0165951  0.0618339  0.2684 0.7884310
yeardum1982        0.1479947  0.0588824  2.5134 0.0120282 *
yeardum1983        0.1783255  0.0606153  2.9419 0.0032958 **
yeardum1984        0.0344572  0.0597167  0.5770 0.5639907
yeardum1985        0.2206961  0.0610344  3.6159 0.0003060 ***
yeardum1986        0.2428015  0.0649779  3.7367 0.0001912 ***
yeardum1987        0.0489043  0.0615708  0.7943 0.4271186
yeardum1988        0.2243599  0.0605343  3.7063 0.0002155 ***
yeardum1989        0.2215060  0.0624042  3.5495 0.0003940 ***
yeardum1990        0.0688333  0.0607056  1.1339 0.2569648
yeardum1991        0.1370871  0.0638830  2.1459 0.0319892 *
yeardum1992        0.1851857  0.0630868  2.9354 0.0033655 **
yeardum1993        0.0904620  0.0698526  1.2950 0.1954420
yeardum1994        0.1003735  0.0737431  1.3611 0.1736137
yeardum1995        0.1164818  0.0721240  1.6150 0.1064494
yeardum1996        0.0482520  0.0787232  0.6129 0.5399837
yeardum1997        0.1049161  0.0895001  1.1722 0.2412247
yeardum1998        0.2191887  0.1109757  1.9751 0.0483807 *
yeardum1999        0.1573342  0.1397150  1.1261 0.2602422
yeardum2000        0.1532796  0.1627206  0.9420 0.3463059
---
 However STATA gives me
------------------------------------------------------------
-----------------------
            wecon |      Coef.   Std. Err.      z    P>|z|     [95%
Conf.
Interval]
------------------+-----------------------------------------
-----------------------
            trade |   .0093915   .0027483     3.42   0.001      .004005
.014778
      fdistockgdp |  -.0169171   .0092405    -1.83   0.067    -.0350281
.0011938
polrightsreversed |   .0165855   .0119176     1.39   0.164    -.0067726
.0399436
         lnrgdpch |   .1045675   .0431179     2.43   0.015     .0200579
.189077
         execleft |   .1373652   .0384442     3.57   0.000     .0620159
.2127145
          muslim2 |  -.0043645   .0013551    -3.22   0.001    -.0070205
-.0017085
          c100rat |   .0480539   .0657304     0.73   0.465    -.0807752
.1768831
          c111rat |   .0170048   .0676272     0.25   0.801    -.1155421
.1495516

Really would appreciate any help explaining why the results are so
different

	[[alternative HTML version deleted]]



--------------------- end original message ---------

 
Ai sensi del D.Lgs. 196/2003 si precisa che le informazi...{{dropped:12}}



More information about the R-help mailing list