Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • MIMIC (SEM) Models with panel (longitudinal) data

    IS it possible to estimate a MIMIC (SEM) model for a panel of countries and years? Could I replicate using STATA theTable 6 of this paper: Dybka, P., Kowalczuk, M., Olesiński, B., Torój, A., & Rozkrut, M. (2019). Currency demand and MIMIC models: towards a structured hybrid method of measuring the shadow economy. International Tax and Public Finance, 26(1), 4-40.
    Thanks

    I

  • #2
    Hallo Antonio,

    Have you found a solution on applying MIMIC model for a panel of countries and years? I would like to run a similar model. It would be grate if you could share the solution.

    Thank you,
    Rattiya

    Comment


    • #3
      It wouldn't be something like the following?

      .ÿ
      .ÿversionÿ17.0

      .ÿ
      .ÿclearÿ*

      .ÿ
      .ÿsetÿseedÿ`=strreverse("1619609")'

      .ÿ
      .ÿquietlyÿsetÿobsÿ500

      .ÿ
      .ÿgenerateÿintÿpidÿ=ÿ_n

      .ÿgenerateÿdoubleÿpid_uÿ=ÿrnormal()

      .ÿgenerateÿdoubleÿpreÿ=ÿruniform(-0.5,ÿ0.5)

      .ÿ
      .ÿquietlyÿexpandÿ10

      .ÿbysortÿpid:ÿgenerateÿbyteÿtimÿ=ÿ_n

      .ÿ
      .ÿquietlyÿdrawnormÿe1ÿe2ÿe3,ÿdoubleÿcorr(1ÿ0.5ÿ0.5ÿ\ÿ0.5ÿ1ÿ0.5ÿ\ÿ0.5ÿ0.5ÿ1)

      .ÿ
      .ÿforvaluesÿiÿ=ÿ1/3ÿ{
      ÿÿ2.ÿÿÿÿÿÿÿÿÿgenerateÿdoubleÿout`i'ÿ=ÿ1ÿ+ÿ1ÿ*ÿpreÿ+ÿ(timÿ-ÿ5.5)ÿ/ÿ9ÿ+ÿ///
      >ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿpid_uÿ+ÿ///
      >ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿe`i'
      ÿÿ3.ÿ}

      .ÿ
      .ÿ*
      .ÿ*ÿBeginÿhere
      .ÿ*
      .ÿgsemÿ(out1ÿout2ÿout3ÿ<-ÿF)ÿ(Fÿ<-ÿc.preÿi.tim)ÿ///
      >ÿÿÿÿÿÿÿÿÿ(out1@1ÿ<-ÿM[pid])ÿ(out2@1ÿ<-ÿM[pid])ÿ(out3@1ÿ<-ÿM[pid]),ÿ///
      >ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿnocnsreportÿnodvheaderÿnolog

      GeneralizedÿstructuralÿequationÿmodelÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿobsÿ=ÿ5,000
      Logÿlikelihoodÿ=ÿ-20074.355

      ------------------------------------------------------------------------------
      ÿÿÿÿÿÿÿÿÿÿÿÿÿ|ÿCoefficientÿÿStd.ÿerr.ÿÿÿÿÿÿzÿÿÿÿP>|z|ÿÿÿÿÿ[95%ÿconf.ÿinterval]
      -------------+----------------------------------------------------------------
      out1ÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿM[pid]ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿÿFÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ.5690399ÿÿÿ.0571988ÿÿÿÿÿ9.95ÿÿÿ0.000ÿÿÿÿÿ.4569323ÿÿÿÿ.6811474
      -------------+----------------------------------------------------------------
      out2ÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿM[pid]ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿÿFÿ|ÿÿÿ1.031858ÿÿÿ.0226704ÿÿÿÿ45.52ÿÿÿ0.000ÿÿÿÿÿ.9874251ÿÿÿÿ1.076292
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ.5620843ÿÿÿ.0577623ÿÿÿÿÿ9.73ÿÿÿ0.000ÿÿÿÿÿ.4488722ÿÿÿÿ.6752964
      -------------+----------------------------------------------------------------
      out3ÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿM[pid]ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿÿFÿ|ÿÿÿ1.003519ÿÿÿ.0221415ÿÿÿÿ45.32ÿÿÿ0.000ÿÿÿÿÿ.9601219ÿÿÿÿ1.046915
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿ.573707ÿÿÿ.0572028ÿÿÿÿ10.03ÿÿÿ0.000ÿÿÿÿÿ.4615916ÿÿÿÿ.6858225
      -------------+----------------------------------------------------------------
      Fÿÿÿÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿpreÿ|ÿÿÿ.9899723ÿÿÿ.1547052ÿÿÿÿÿ6.40ÿÿÿ0.000ÿÿÿÿÿ.6867556ÿÿÿÿ1.293189
      ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿtimÿ|
      ÿÿÿÿÿÿÿÿÿÿ2ÿÿ|ÿÿÿÿ.126801ÿÿÿ.0501352ÿÿÿÿÿ2.53ÿÿÿ0.011ÿÿÿÿÿ.0285378ÿÿÿÿ.2250642
      ÿÿÿÿÿÿÿÿÿÿ3ÿÿ|ÿÿÿ.2507597ÿÿÿ.0502802ÿÿÿÿÿ4.99ÿÿÿ0.000ÿÿÿÿÿ.1522123ÿÿÿÿÿ.349307
      ÿÿÿÿÿÿÿÿÿÿ4ÿÿ|ÿÿÿ.2923716ÿÿÿ.0502811ÿÿÿÿÿ5.81ÿÿÿ0.000ÿÿÿÿÿ.1938225ÿÿÿÿ.3909207
      ÿÿÿÿÿÿÿÿÿÿ5ÿÿ|ÿÿÿ.3975033ÿÿÿ.0504394ÿÿÿÿÿ7.88ÿÿÿ0.000ÿÿÿÿÿ.2986439ÿÿÿÿ.4963626
      ÿÿÿÿÿÿÿÿÿÿ6ÿÿ|ÿÿÿ.5029159ÿÿÿ.0505696ÿÿÿÿÿ9.95ÿÿÿ0.000ÿÿÿÿÿ.4038013ÿÿÿÿ.6020306
      ÿÿÿÿÿÿÿÿÿÿ7ÿÿ|ÿÿÿ.6766229ÿÿÿ.0510235ÿÿÿÿ13.26ÿÿÿ0.000ÿÿÿÿÿ.5766186ÿÿÿÿ.7766272
      ÿÿÿÿÿÿÿÿÿÿ8ÿÿ|ÿÿÿ.7972289ÿÿÿ.0512257ÿÿÿÿ15.56ÿÿÿ0.000ÿÿÿÿÿ.6968284ÿÿÿÿ.8976294
      ÿÿÿÿÿÿÿÿÿÿ9ÿÿ|ÿÿÿ.9301085ÿÿÿ.0516572ÿÿÿÿ18.01ÿÿÿ0.000ÿÿÿÿÿ.8288623ÿÿÿÿ1.031355
      ÿÿÿÿÿÿÿÿÿ10ÿÿ|ÿÿÿ.9792638ÿÿÿ.0517185ÿÿÿÿ18.93ÿÿÿ0.000ÿÿÿÿÿ.8778973ÿÿÿÿÿ1.08063
      -------------+----------------------------------------------------------------
      ÿÿvar(M[pid])|ÿÿÿ.9508162ÿÿÿ.0642172ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.8329276ÿÿÿÿÿ1.08539
      ÿÿÿÿÿvar(e.F)|ÿÿÿÿ.467868ÿÿÿ.0180779ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4337443ÿÿÿÿ.5046762
      -------------+----------------------------------------------------------------
      ÿÿvar(e.out1)|ÿÿÿ.5052615ÿÿÿ.0147036ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4772495ÿÿÿÿ.5349177
      ÿÿvar(e.out2)|ÿÿÿ.4718743ÿÿÿ.0146545ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4440086ÿÿÿÿ.5014888
      ÿÿvar(e.out3)|ÿÿÿÿÿ.49997ÿÿÿ.0146032ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4721521ÿÿÿÿ.5294268
      ------------------------------------------------------------------------------

      .ÿ
      .ÿexit

      endÿofÿdo-file


      .


      I suppose that if you believed that the regression coefficient for the random effect differs for each manifest variable, then you could let that float.

      Comment

      Working...
      X