Announcement

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

  • #16
    Originally posted by Diep Nguyen View Post
    Here is my data file. Could you check the SEM model with your code? I've tried but it failed....
    With Newton-Raphson, sometimes you have to coax the likelihood function into the range where the algorithm can converge. With maximum-likeihood estimators like SEM, if you're having trouble with convergence sometimes you can fit a somewhat more constrained model, and then use those estimates as starting values for the less-constrained final model.

    That tactic works here. The outpout is shown below.

    The do-file and log-file are attached for your convenience. I followed the diagram in your initial post; I hope that I got everything correct. (The model is too Rube Goldberg for me and so nothing that I'd ever contemplate.)

    .ÿ
    .ÿversionÿ18.0

    .ÿ
    .ÿlogÿcloseÿ_all

    .ÿlogÿusingÿsem.smcl,ÿnomsgÿname(lo)ÿreplace

    .ÿ
    .ÿclearÿ*

    .ÿ
    .ÿquietlyÿimportÿdelimitedÿ"Dataÿ20230922.csv",ÿasdoubleÿ

    .ÿ
    .ÿquietlyÿsemÿ///
    >ÿÿÿÿÿÿÿÿÿ(c1x?ÿ<-ÿSE1)ÿ(c2x?ÿ<-ÿSE2)ÿ(c3x?ÿ<-ÿSE3)ÿ(SE1ÿSE2ÿSE3ÿ<-ÿSE)ÿ///
    >ÿÿÿÿÿÿÿÿÿ(b1x1ÿb1x2ÿ<-ÿGC1)ÿ(b1x3-b2x2ÿ<-ÿGC2)ÿ(b2x3-ÿb3x5ÿ<-ÿGC3)ÿ(GC1ÿGC2ÿGC3ÿ<-ÿGC)ÿ///
    >ÿÿÿÿÿÿÿÿÿ(SE@1ÿGC@1ÿe1x?ÿ<-ÿCL)

    .ÿassertÿe(converged)

    .ÿ
    .ÿtempnameÿB

    .ÿmatrixÿdefineÿ`B'ÿ=ÿe(b)

    .ÿ
    .ÿsemÿ///
    >ÿÿÿÿÿÿÿÿÿ(c1x?ÿ<-ÿSE1)ÿ(c2x?ÿ<-ÿSE2)ÿ(c3x?ÿ<-ÿSE3)ÿ(SE1ÿSE2ÿSE3ÿ<-ÿSE)ÿ///
    >ÿÿÿÿÿÿÿÿÿ(b1x1ÿb1x2ÿ<-ÿGC1)ÿ(b1x3-b2x2ÿ<-ÿGC2)ÿ(b2x3-ÿb3x5ÿ<-ÿGC3)ÿ(GC1ÿGC2ÿGC3ÿ<-ÿGC)ÿ///
    >ÿÿÿÿÿÿÿÿÿ(SEÿGCÿe1x?ÿ<-ÿCL),ÿfrom(`B')ÿnocnsreportÿnodescribeÿnolog

    StructuralÿequationÿmodelÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿobsÿ=ÿ494
    Estimationÿmethod:ÿml

    Logÿlikelihoodÿ=ÿ-11691.205

    ------------------------------------------------------------------------------
    ÿÿÿÿÿÿÿÿÿÿÿÿÿ|ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿOIM
    ÿÿÿÿÿÿÿÿÿÿÿÿÿ|ÿCoefficientÿÿstd.ÿerr.ÿÿÿÿÿÿzÿÿÿÿP>|z|ÿÿÿÿÿ[95%ÿconf.ÿinterval]
    -------------+----------------------------------------------------------------
    Structuralÿÿÿ|
    ÿÿSE1ÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿSEÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿSE2ÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿSEÿ|ÿÿÿ1.051997ÿÿÿ.0367389ÿÿÿÿ28.63ÿÿÿ0.000ÿÿÿÿÿ.9799895ÿÿÿÿ1.124004
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿSE3ÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿSEÿ|ÿÿÿ1.062761ÿÿÿ.0331458ÿÿÿÿ32.06ÿÿÿ0.000ÿÿÿÿÿ.9977968ÿÿÿÿ1.127726
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿSEÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿ.9161776ÿÿÿ.0368929ÿÿÿÿ24.83ÿÿÿ0.000ÿÿÿÿÿ.8438689ÿÿÿÿ.9884863
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿGC1ÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿGCÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿGC2ÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿGCÿ|ÿÿÿ1.278963ÿÿÿ.0930113ÿÿÿÿ13.75ÿÿÿ0.000ÿÿÿÿÿ1.096664ÿÿÿÿ1.461261
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿGC3ÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿGCÿ|ÿÿÿ1.203971ÿÿÿ.0915942ÿÿÿÿ13.14ÿÿÿ0.000ÿÿÿÿÿÿ1.02445ÿÿÿÿ1.383493
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿGCÿÿÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿ.5432622ÿÿÿ.0415703ÿÿÿÿ13.07ÿÿÿ0.000ÿÿÿÿÿÿ.461786ÿÿÿÿ.6247384
    -------------+----------------------------------------------------------------
    Measurementÿÿ|
    ÿÿc1x1ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE1ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.708502ÿÿÿ.0521633ÿÿÿÿ71.09ÿÿÿ0.000ÿÿÿÿÿ3.606264ÿÿÿÿÿ3.81074
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿc1x2ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE1ÿ|ÿÿÿ1.031135ÿÿÿ.0248369ÿÿÿÿ41.52ÿÿÿ0.000ÿÿÿÿÿ.9824552ÿÿÿÿ1.079814
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.718623ÿÿÿ.0527051ÿÿÿÿ70.56ÿÿÿ0.000ÿÿÿÿÿ3.615323ÿÿÿÿ3.821923
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿc1x3ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE1ÿ|ÿÿÿ1.038696ÿÿÿÿÿ.02402ÿÿÿÿ43.24ÿÿÿ0.000ÿÿÿÿÿ.9916175ÿÿÿÿ1.085774
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.773279ÿÿÿ.0525773ÿÿÿÿ71.77ÿÿÿ0.000ÿÿÿÿÿÿ3.67023ÿÿÿÿ3.876329
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿc2x1ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE2ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.591093ÿÿÿ.0537639ÿÿÿÿ66.79ÿÿÿ0.000ÿÿÿÿÿ3.485718ÿÿÿÿ3.696468
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿc2x2ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE2ÿ|ÿÿÿ.9363871ÿÿÿ.0259911ÿÿÿÿ36.03ÿÿÿ0.000ÿÿÿÿÿ.8854455ÿÿÿÿ.9873286
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿÿ3.6417ÿÿÿ.0519096ÿÿÿÿ70.15ÿÿÿ0.000ÿÿÿÿÿ3.539959ÿÿÿÿ3.743441
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿc2x3ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE2ÿ|ÿÿÿ.9859987ÿÿÿ.0266454ÿÿÿÿ37.00ÿÿÿ0.000ÿÿÿÿÿ.9337746ÿÿÿÿ1.038223
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.633603ÿÿÿ.0543432ÿÿÿÿ66.86ÿÿÿ0.000ÿÿÿÿÿ3.527093ÿÿÿÿ3.740114
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿc3x1ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE3ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.696356ÿÿÿ.0513888ÿÿÿÿ71.93ÿÿÿ0.000ÿÿÿÿÿ3.595636ÿÿÿÿ3.797077
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿc3x2ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿSE3ÿ|ÿÿÿ.9850818ÿÿÿ.0204865ÿÿÿÿ48.08ÿÿÿ0.000ÿÿÿÿÿÿ.944929ÿÿÿÿ1.025235
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.633603ÿÿÿ.0517951ÿÿÿÿ70.15ÿÿÿ0.000ÿÿÿÿÿ3.532087ÿÿÿÿÿ3.73512
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb1x1ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC1ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿ3.01417ÿÿÿ.0506033ÿÿÿÿ59.56ÿÿÿ0.000ÿÿÿÿÿ2.914989ÿÿÿÿ3.113351
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb1x2ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC1ÿ|ÿÿÿ1.061108ÿÿÿÿ.069233ÿÿÿÿ15.33ÿÿÿ0.000ÿÿÿÿÿ.9254141ÿÿÿÿ1.196802
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.277328ÿÿÿ.0499556ÿÿÿÿ65.60ÿÿÿ0.000ÿÿÿÿÿ3.179417ÿÿÿÿ3.375239
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb1x3ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.777328ÿÿÿ.0464047ÿÿÿÿ81.40ÿÿÿ0.000ÿÿÿÿÿ3.686376ÿÿÿÿ3.868279
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb1x4ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿ1.054854ÿÿÿ.0380935ÿÿÿÿ27.69ÿÿÿ0.000ÿÿÿÿÿ.9801918ÿÿÿÿ1.129516
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.759109ÿÿÿ.0475745ÿÿÿÿ79.02ÿÿÿ0.000ÿÿÿÿÿ3.665865ÿÿÿÿ3.852354
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb1x5ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿÿ1.03348ÿÿÿ.0429075ÿÿÿÿ24.09ÿÿÿ0.000ÿÿÿÿÿÿ.949383ÿÿÿÿ1.117578
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.615385ÿÿÿ.0496229ÿÿÿÿ72.86ÿÿÿ0.000ÿÿÿÿÿ3.518125ÿÿÿÿ3.712644
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb2x1ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿ.9341177ÿÿÿÿ.040127ÿÿÿÿ23.28ÿÿÿ0.000ÿÿÿÿÿ.8554702ÿÿÿÿ1.012765
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.668016ÿÿÿ.0454251ÿÿÿÿ80.75ÿÿÿ0.000ÿÿÿÿÿ3.578985ÿÿÿÿ3.757048
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb2x2ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿÿ.875959ÿÿÿ.0388336ÿÿÿÿ22.56ÿÿÿ0.000ÿÿÿÿÿ.7998465ÿÿÿÿ.9520714
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.791498ÿÿÿ.0433938ÿÿÿÿ87.37ÿÿÿ0.000ÿÿÿÿÿ3.706448ÿÿÿÿ3.876548
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb2x3ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.421053ÿÿÿÿ.048685ÿÿÿÿ70.27ÿÿÿ0.000ÿÿÿÿÿ3.325632ÿÿÿÿ3.516473
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb2x4ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ.8743251ÿÿÿÿ.042508ÿÿÿÿ20.57ÿÿÿ0.000ÿÿÿÿÿ.7910109ÿÿÿÿ.9576392
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.346154ÿÿÿ.0450744ÿÿÿÿ74.24ÿÿÿ0.000ÿÿÿÿÿÿ3.25781ÿÿÿÿ3.434498
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb3x1ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ.9838105ÿÿÿ.0424863ÿÿÿÿ23.16ÿÿÿ0.000ÿÿÿÿÿÿ.900539ÿÿÿÿ1.067082
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.267206ÿÿÿ.0466347ÿÿÿÿ70.06ÿÿÿ0.000ÿÿÿÿÿ3.175804ÿÿÿÿ3.358609
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb3x2ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿÿ1.07211ÿÿÿ.0432244ÿÿÿÿ24.80ÿÿÿ0.000ÿÿÿÿÿ.9873915ÿÿÿÿ1.156828
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.370445ÿÿÿ.0487181ÿÿÿÿ69.18ÿÿÿ0.000ÿÿÿÿÿÿ3.27496ÿÿÿÿ3.465931
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb3x3ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿÿ.986871ÿÿÿ.0430915ÿÿÿÿ22.90ÿÿÿ0.000ÿÿÿÿÿ.9024132ÿÿÿÿ1.071329
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.481781ÿÿÿ.0470114ÿÿÿÿ74.06ÿÿÿ0.000ÿÿÿÿÿ3.389641ÿÿÿÿ3.573922
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb3x4ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ1.067303ÿÿÿ.0416187ÿÿÿÿ25.64ÿÿÿ0.000ÿÿÿÿÿ.9857322ÿÿÿÿ1.148875
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.370445ÿÿÿ.0474397ÿÿÿÿ71.05ÿÿÿ0.000ÿÿÿÿÿ3.277465ÿÿÿÿ3.463425
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿb3x5ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ1.026561ÿÿÿ.0428075ÿÿÿÿ23.98ÿÿÿ0.000ÿÿÿÿÿ.9426597ÿÿÿÿ1.110462
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.532389ÿÿÿ.0476024ÿÿÿÿ74.21ÿÿÿ0.000ÿÿÿÿÿÿ3.43909ÿÿÿÿ3.625688
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿe1x1ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.548583ÿÿÿ.0495986ÿÿÿÿ71.55ÿÿÿ0.000ÿÿÿÿÿ3.451372ÿÿÿÿ3.645794
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿe1x2ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿ1.017133ÿÿÿ.0290793ÿÿÿÿ34.98ÿÿÿ0.000ÿÿÿÿÿ.9601387ÿÿÿÿ1.074127
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿ3.58502ÿÿÿ.0492502ÿÿÿÿ72.79ÿÿÿ0.000ÿÿÿÿÿ3.488492ÿÿÿÿ3.681549
    ÿÿ-----------+----------------------------------------------------------------
    ÿÿe1x3ÿÿÿÿÿÿÿ|
    ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿ.9974863ÿÿÿ.0279925ÿÿÿÿ35.63ÿÿÿ0.000ÿÿÿÿÿÿ.942622ÿÿÿÿ1.052351
    ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿ3.61336ÿÿÿ.0479527ÿÿÿÿ75.35ÿÿÿ0.000ÿÿÿÿÿ3.519375ÿÿÿÿ3.707346
    -------------+----------------------------------------------------------------
    ÿÿvar(e.c1x1)|ÿÿÿ.1829607ÿÿÿ.0147854ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1561603ÿÿÿÿ.2143607
    ÿÿvar(e.c1x2)|ÿÿÿ.1375924ÿÿÿ.0126478ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1149079ÿÿÿÿ.1647551
    ÿÿvar(e.c1x3)|ÿÿÿ.1127763ÿÿÿ.0115301ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.092298ÿÿÿÿ.1377983
    ÿÿvar(e.c2x1)|ÿÿÿ.1721754ÿÿÿ.0165614ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1425919ÿÿÿÿ.2078966
    ÿÿvar(e.c2x2)|ÿÿÿ.2300579ÿÿÿÿÿ.01881ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1959933ÿÿÿÿ.2700432
    ÿÿvar(e.c2x3)|ÿÿÿ.2380278ÿÿÿ.0197887ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2022377ÿÿÿÿ.2801517
    ÿÿvar(e.c3x1)|ÿÿÿ.0838395ÿÿÿ.0114939ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0640847ÿÿÿÿ.1096839
    ÿÿvar(e.c3x2)|ÿÿÿ.1406958ÿÿÿ.0133321ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1168483ÿÿÿÿ.1694103
    ÿÿvar(e.b1x1)|ÿÿÿ.5055991ÿÿÿ.0504007ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4158668ÿÿÿÿ.6146931
    ÿÿvar(e.b1x2)|ÿÿÿ.3777787ÿÿÿ.0498613ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2916696ÿÿÿÿ.4893096
    ÿÿvar(e.b1x3)|ÿÿÿ.2702206ÿÿÿÿÿ.02173ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2308173ÿÿÿÿ.3163504
    ÿÿvar(e.b1x4)|ÿÿÿ.2350812ÿÿÿ.0203308ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1984279ÿÿÿÿ.2785052
    ÿÿvar(e.b1x5)|ÿÿÿÿÿ.36886ÿÿÿ.0275195ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.3186809ÿÿÿÿ.4269403
    ÿÿvar(e.b2x1)|ÿÿÿ.3269028ÿÿÿ.0242152ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2827263ÿÿÿÿ.3779819
    ÿÿvar(e.b2x2)|ÿÿÿ.3213116ÿÿÿ.0235175ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2783718ÿÿÿÿ.3708751
    ÿÿvar(e.b2x3)|ÿÿÿ.3756743ÿÿÿ.0265939ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.3270056ÿÿÿÿ.4315865
    ÿÿvar(e.b2x4)|ÿÿÿÿ.395759ÿÿÿ.0271897ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.3459003ÿÿÿÿ.4528044
    ÿÿvar(e.b3x1)|ÿÿÿ.3046711ÿÿÿÿ.021963ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2645271ÿÿÿÿ.3509072
    ÿÿvar(e.b3x2)|ÿÿÿ.2584475ÿÿÿ.0196891ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2226003ÿÿÿÿ.3000674
    ÿÿvar(e.b3x3)|ÿÿÿ.3172986ÿÿÿ.0227518ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2756976ÿÿÿÿ.3651769
    ÿÿvar(e.b3x4)|ÿÿÿ.2058956ÿÿÿ.0165136ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1759452ÿÿÿÿ.2409444
    ÿÿvar(e.b3x5)|ÿÿÿ.2813735ÿÿÿ.0207878ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2434426ÿÿÿÿ.3252143
    ÿÿvar(e.e1x1)|ÿÿÿ.2167468ÿÿÿÿ.017411ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1851726ÿÿÿÿ.2537048
    ÿÿvar(e.e1x2)|ÿÿÿÿ.165225ÿÿÿ.0146709ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1388336ÿÿÿÿ.1966332
    ÿÿvar(e.e1x3)|ÿÿÿ.1424442ÿÿÿ.0132339ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1187307ÿÿÿÿ.1708939
    ÿÿÿvar(e.SE1)|ÿÿÿ.1437629ÿÿÿ.0166918ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.114503ÿÿÿÿ.1804998
    ÿÿÿvar(e.SE2)|ÿÿÿ.1297473ÿÿÿ.0185665ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0980152ÿÿÿÿ.1717526
    ÿÿÿvar(e.SE3)|ÿÿÿ.0715456ÿÿÿ.0147052ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0478223ÿÿÿÿ.1070375
    ÿÿÿÿvar(e.SE)|ÿÿÿ.1793293ÿÿÿÿ.020207ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1437928ÿÿÿÿ.2236482
    ÿÿÿvar(e.GC1)|ÿÿÿ.3171735ÿÿÿ.0434836ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2424373ÿÿÿÿ.4149485
    ÿÿÿvar(e.GC2)|ÿÿÿ.0702159ÿÿÿ.0192738ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0410002ÿÿÿÿ.1202499
    ÿÿÿvar(e.GC3)|ÿÿÿ.1542162ÿÿÿ.0214327ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1174441ÿÿÿÿ.2025016
    ÿÿÿÿvar(e.GC)|ÿÿÿ.1475165ÿÿÿ.0228982ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1088211ÿÿÿÿ.1999714
    ÿÿÿÿÿÿvar(CL)|ÿÿÿ.9985042ÿÿÿ.0768234ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.8587361ÿÿÿÿ1.161021
    ------------------------------------------------------------------------------
    LRÿtestÿofÿmodelÿvs.ÿsaturated:ÿchi2(267)ÿ=ÿ1129.29ÿÿÿÿÿÿÿProbÿ>ÿchi2ÿ=ÿ0.0000

    .ÿassertÿe(converged)

    .ÿ
    .ÿmata:ÿrowsum(st_matrix("e(ilog)")ÿ:!=ÿ0)ÿ//ÿConvergedÿinÿsixÿiterations
    ÿÿ6

    .ÿ
    .ÿquietlyÿlogÿcloseÿlo

    .ÿ
    .ÿexit

    endÿofÿdo-file


    .


    The data seem artificial.
    Attached Files
    Last edited by Joseph Coveney; 26 Sep 2023, 20:49.

    Comment


    • #17
      Sorry, I got the arrows pointing the wrong way between CL and SE and GC. This might be closer to the original diagram. Do-file is again attached if I've still got it wrong and you need to tweak it.

      .ÿ
      .ÿversionÿ18.0

      .ÿ
      .ÿlogÿcloseÿ_all

      .ÿlogÿusingÿsem2.smcl,ÿnomsgÿname(lo)ÿreplace

      .ÿ
      .ÿclearÿ*

      .ÿ
      .ÿquietlyÿimportÿdelimitedÿ"Dataÿ20230922.csv",ÿasdoubleÿ

      .ÿ
      .ÿquietlyÿsemÿ///
      >ÿÿÿÿÿ(e1x?ÿ<-ÿCL)ÿ///
      >ÿÿÿÿÿ(c1x?ÿ<-ÿSE1)ÿ(c2x?ÿ<-ÿSE2)ÿ(c3x?ÿ<-ÿSE3)ÿ(SE1@1ÿSE2@1ÿSE3@1ÿCL@1ÿ<-ÿSE)ÿ///
      >ÿÿÿÿÿ(b1x1ÿb1x2ÿ<-ÿGC1)ÿ(b1x3-b2x2ÿ<-ÿGC2)ÿ(b2x3-b3x5ÿ<-ÿGC3)ÿ(GC1@1ÿGC2@1ÿGC3@1ÿCL@1ÿ<-ÿGC),ÿ///
      >ÿÿÿÿÿÿÿÿÿcovariance(SE*GC@0)

      .ÿ
      .ÿassertÿe(converged)

      .ÿ
      .ÿtempnameÿB

      .ÿmatrixÿdefineÿ`B'ÿ=ÿe(b)

      .ÿ
      .ÿsemÿ///
      >ÿÿÿÿÿ(e1x?ÿ<-ÿCL)ÿ///
      >ÿÿÿÿÿ(c1x?ÿ<-ÿSE1)ÿ(c2x?ÿ<-ÿSE2)ÿ(c3x?ÿ<-ÿSE3)ÿ(SE1ÿSE2ÿSE3ÿCLÿ<-ÿSE)ÿ///
      >ÿÿÿÿÿ(b1x1ÿb1x2ÿ<-ÿGC1)ÿ(b1x3-b2x2ÿ<-ÿGC2)ÿ(b2x3-b3x5ÿ<-ÿGC3)ÿ(GC1ÿGC2ÿGC3ÿCLÿ<-ÿGC),ÿ///
      >ÿÿÿÿÿÿÿÿÿcovariance(SE*GC@0)ÿ///
      >ÿÿÿÿÿÿÿÿÿÿÿÿÿfrom(`B')ÿnocnsreportÿnodescribeÿnolog

      StructuralÿequationÿmodelÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿobsÿ=ÿ494
      Estimationÿmethod:ÿml

      Logÿlikelihoodÿ=ÿ-11894.907

      ------------------------------------------------------------------------------
      ÿÿÿÿÿÿÿÿÿÿÿÿÿ|ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿOIM
      ÿÿÿÿÿÿÿÿÿÿÿÿÿ|ÿCoefficientÿÿstd.ÿerr.ÿÿÿÿÿÿzÿÿÿÿP>|z|ÿÿÿÿÿ[95%ÿconf.ÿinterval]
      -------------+----------------------------------------------------------------
      Structuralÿÿÿ|
      ÿÿCLÿÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿSEÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿÿÿÿGCÿ|ÿÿÿ.2207291ÿÿÿ.0669043ÿÿÿÿÿ3.30ÿÿÿ0.001ÿÿÿÿÿ.0895991ÿÿÿÿÿ.351859
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿSE1ÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿSEÿ|ÿÿÿ1.297417ÿÿÿ.0808953ÿÿÿÿ16.04ÿÿÿ0.000ÿÿÿÿÿ1.138865ÿÿÿÿ1.455969
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿSE2ÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿSEÿ|ÿÿÿ1.368674ÿÿÿ.0843209ÿÿÿÿ16.23ÿÿÿ0.000ÿÿÿÿÿ1.203408ÿÿÿÿÿ1.53394
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿSE3ÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿSEÿ|ÿÿÿ1.377896ÿÿÿ.0811375ÿÿÿÿ16.98ÿÿÿ0.000ÿÿÿÿÿÿ1.21887ÿÿÿÿ1.536923
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿGC1ÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿGCÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿGC2ÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿGCÿ|ÿÿÿ1.290642ÿÿÿÿ.099343ÿÿÿÿ12.99ÿÿÿ0.000ÿÿÿÿÿ1.095933ÿÿÿÿ1.485351
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿGC3ÿÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿGCÿ|ÿÿÿ1.190267ÿÿÿ.0916615ÿÿÿÿ12.99ÿÿÿ0.000ÿÿÿÿÿ1.010614ÿÿÿÿÿ1.36992
      -------------+----------------------------------------------------------------
      Measurementÿÿ|
      ÿÿe1x1ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.548583ÿÿÿ.0459973ÿÿÿÿ77.15ÿÿÿ0.000ÿÿÿÿÿÿ3.45843ÿÿÿÿ3.638736
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿe1x2ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿ1.018518ÿÿÿ.0285954ÿÿÿÿ35.62ÿÿÿ0.000ÿÿÿÿÿ.9624722ÿÿÿÿ1.074564
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿ3.58502ÿÿÿ.0454802ÿÿÿÿ78.83ÿÿÿ0.000ÿÿÿÿÿ3.495881ÿÿÿÿÿ3.67416
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿe1x3ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿÿCLÿ|ÿÿÿ.9991152ÿÿÿ.0275719ÿÿÿÿ36.24ÿÿÿ0.000ÿÿÿÿÿ.9450752ÿÿÿÿ1.053155
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿ3.61336ÿÿÿ.0442245ÿÿÿÿ81.70ÿÿÿ0.000ÿÿÿÿÿ3.526682ÿÿÿÿ3.700039
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc1x1ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE1ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.708502ÿÿÿ.0521633ÿÿÿÿ71.09ÿÿÿ0.000ÿÿÿÿÿ3.606264ÿÿÿÿÿ3.81074
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc1x2ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE1ÿ|ÿÿÿ1.031661ÿÿÿ.0248958ÿÿÿÿ41.44ÿÿÿ0.000ÿÿÿÿÿ.9828658ÿÿÿÿ1.080456
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.718623ÿÿÿ.0527051ÿÿÿÿ70.56ÿÿÿ0.000ÿÿÿÿÿ3.615323ÿÿÿÿ3.821923
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc1x3ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE1ÿ|ÿÿÿ1.039584ÿÿÿ.0240609ÿÿÿÿ43.21ÿÿÿ0.000ÿÿÿÿÿ.9924255ÿÿÿÿ1.086743
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.773279ÿÿÿ.0525773ÿÿÿÿ71.77ÿÿÿ0.000ÿÿÿÿÿÿ3.67023ÿÿÿÿ3.876329
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc2x1ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE2ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.591093ÿÿÿ.0537639ÿÿÿÿ66.79ÿÿÿ0.000ÿÿÿÿÿ3.485718ÿÿÿÿ3.696468
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc2x2ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE2ÿ|ÿÿÿ.9358522ÿÿÿÿ.025967ÿÿÿÿ36.04ÿÿÿ0.000ÿÿÿÿÿ.8849577ÿÿÿÿ.9867466
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿÿ3.6417ÿÿÿ.0519096ÿÿÿÿ70.15ÿÿÿ0.000ÿÿÿÿÿ3.539959ÿÿÿÿ3.743441
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc2x3ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE2ÿ|ÿÿÿ.9853883ÿÿÿ.0266218ÿÿÿÿ37.01ÿÿÿ0.000ÿÿÿÿÿ.9332105ÿÿÿÿ1.037566
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.633603ÿÿÿ.0543432ÿÿÿÿ66.86ÿÿÿ0.000ÿÿÿÿÿ3.527093ÿÿÿÿ3.740114
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc3x1ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE3ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.696356ÿÿÿ.0513888ÿÿÿÿ71.93ÿÿÿ0.000ÿÿÿÿÿ3.595636ÿÿÿÿ3.797077
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿc3x2ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿSE3ÿ|ÿÿÿ.9801784ÿÿÿ.0205158ÿÿÿÿ47.78ÿÿÿ0.000ÿÿÿÿÿ.9399683ÿÿÿÿ1.020389
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.633603ÿÿÿ.0517951ÿÿÿÿ70.15ÿÿÿ0.000ÿÿÿÿÿ3.532087ÿÿÿÿÿ3.73512
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb1x1ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC1ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿÿ3.01417ÿÿÿ.0506033ÿÿÿÿ59.56ÿÿÿ0.000ÿÿÿÿÿ2.914989ÿÿÿÿ3.113351
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb1x2ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC1ÿ|ÿÿÿ1.072046ÿÿÿ.0707843ÿÿÿÿ15.15ÿÿÿ0.000ÿÿÿÿÿ.9333116ÿÿÿÿ1.210781
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.277328ÿÿÿ.0499556ÿÿÿÿ65.60ÿÿÿ0.000ÿÿÿÿÿ3.179417ÿÿÿÿ3.375239
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb1x3ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.777328ÿÿÿ.0464047ÿÿÿÿ81.40ÿÿÿ0.000ÿÿÿÿÿ3.686376ÿÿÿÿ3.868279
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb1x4ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿ1.055056ÿÿÿÿÿ.03767ÿÿÿÿ28.01ÿÿÿ0.000ÿÿÿÿÿ.9812239ÿÿÿÿ1.128888
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.759109ÿÿÿ.0475745ÿÿÿÿ79.02ÿÿÿ0.000ÿÿÿÿÿ3.665865ÿÿÿÿ3.852354
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb1x5ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿ1.035841ÿÿÿ.0425775ÿÿÿÿ24.33ÿÿÿ0.000ÿÿÿÿÿ.9523905ÿÿÿÿ1.119291
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.615385ÿÿÿ.0496229ÿÿÿÿ72.86ÿÿÿ0.000ÿÿÿÿÿ3.518125ÿÿÿÿ3.712644
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb2x1ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿ.9290195ÿÿÿ.0400142ÿÿÿÿ23.22ÿÿÿ0.000ÿÿÿÿÿ.8505931ÿÿÿÿ1.007446
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.668016ÿÿÿ.0454251ÿÿÿÿ80.75ÿÿÿ0.000ÿÿÿÿÿ3.578985ÿÿÿÿ3.757048
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb2x2ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC2ÿ|ÿÿÿ.8619442ÿÿÿ.0388385ÿÿÿÿ22.19ÿÿÿ0.000ÿÿÿÿÿ.7858222ÿÿÿÿ.9380663
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.791498ÿÿÿ.0433938ÿÿÿÿ87.37ÿÿÿ0.000ÿÿÿÿÿ3.706448ÿÿÿÿ3.876548
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb2x3ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿÿÿÿÿÿÿÿ1ÿÿ(constrained)
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.421053ÿÿÿÿ.048685ÿÿÿÿ70.27ÿÿÿ0.000ÿÿÿÿÿ3.325632ÿÿÿÿ3.516473
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb2x4ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ.8718875ÿÿÿ.0424867ÿÿÿÿ20.52ÿÿÿ0.000ÿÿÿÿÿ.7886151ÿÿÿÿ.9551598
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.346154ÿÿÿ.0450744ÿÿÿÿ74.24ÿÿÿ0.000ÿÿÿÿÿÿ3.25781ÿÿÿÿ3.434498
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb3x1ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ.9843613ÿÿÿ.0424042ÿÿÿÿ23.21ÿÿÿ0.000ÿÿÿÿÿ.9012505ÿÿÿÿ1.067472
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.267206ÿÿÿ.0466347ÿÿÿÿ70.06ÿÿÿ0.000ÿÿÿÿÿ3.175804ÿÿÿÿ3.358609
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb3x2ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ1.072927ÿÿÿÿ.043126ÿÿÿÿ24.88ÿÿÿ0.000ÿÿÿÿÿ.9884011ÿÿÿÿ1.157452
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.370445ÿÿÿ.0487181ÿÿÿÿ69.18ÿÿÿ0.000ÿÿÿÿÿÿ3.27496ÿÿÿÿ3.465931
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb3x3ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿ.9849571ÿÿÿ.0430638ÿÿÿÿ22.87ÿÿÿ0.000ÿÿÿÿÿ.9005537ÿÿÿÿ1.069361
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.481781ÿÿÿ.0470114ÿÿÿÿ74.06ÿÿÿ0.000ÿÿÿÿÿ3.389641ÿÿÿÿ3.573922
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb3x4ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿÿ1.06821ÿÿÿ.0415305ÿÿÿÿ25.72ÿÿÿ0.000ÿÿÿÿÿ.9868116ÿÿÿÿ1.149608
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.370445ÿÿÿ.0474397ÿÿÿÿ71.05ÿÿÿ0.000ÿÿÿÿÿ3.277465ÿÿÿÿ3.463425
      ÿÿ-----------+----------------------------------------------------------------
      ÿÿb3x5ÿÿÿÿÿÿÿ|
      ÿÿÿÿÿÿÿÿÿGC3ÿ|ÿÿÿÿ1.02276ÿÿÿ.0428086ÿÿÿÿ23.89ÿÿÿ0.000ÿÿÿÿÿ.9388568ÿÿÿÿ1.106663
      ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ3.532389ÿÿÿ.0476024ÿÿÿÿ74.21ÿÿÿ0.000ÿÿÿÿÿÿ3.43909ÿÿÿÿ3.625688
      -------------+----------------------------------------------------------------
      ÿÿvar(e.e1x1)|ÿÿÿ.2101613ÿÿÿ.0172828ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1788766ÿÿÿÿ.2469176
      ÿÿvar(e.e1x2)|ÿÿÿÿ.155578ÿÿÿ.0147006ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.129276ÿÿÿÿ.1872312
      ÿÿvar(e.e1x3)|ÿÿÿ.1326231ÿÿÿ.0133488ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1088789ÿÿÿÿ.1615453
      ÿÿvar(e.c1x1)|ÿÿÿ.1841062ÿÿÿ.0148514ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1571825ÿÿÿÿ.2156417
      ÿÿvar(e.c1x2)|ÿÿÿ.1375514ÿÿÿ.0126749ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1148232ÿÿÿÿ.1647784
      ÿÿvar(e.c1x3)|ÿÿÿ.1118709ÿÿÿ.0115265ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0914144ÿÿÿÿ.1369051
      ÿÿvar(e.c2x1)|ÿÿÿ.1712331ÿÿÿ.0165219ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1417285ÿÿÿÿ.2068799
      ÿÿvar(e.c2x2)|ÿÿÿ.2304902ÿÿÿ.0188342ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1963802ÿÿÿÿÿ.270525
      ÿÿvar(e.c2x3)|ÿÿÿ.2386239ÿÿÿÿ.019817ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2027796ÿÿÿÿ.2808042
      ÿÿvar(e.c3x1)|ÿÿÿ.0777328ÿÿÿ.0116565ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0579378ÿÿÿÿÿ.104291
      ÿÿvar(e.c3x2)|ÿÿÿ.1465922ÿÿÿ.0137781ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1219289ÿÿÿÿ.1762442
      ÿÿvar(e.b1x1)|ÿÿÿÿ.513347ÿÿÿ.0509308ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4226302ÿÿÿÿ.6235359
      ÿÿvar(e.b1x2)|ÿÿÿ.3689651ÿÿÿÿ.050677ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2818861ÿÿÿÿ.4829441
      ÿÿvar(e.b1x3)|ÿÿÿ.2650962ÿÿÿ.0215926ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2259808ÿÿÿÿ.3109821
      ÿÿvar(e.b1x4)|ÿÿÿ.2290387ÿÿÿ.0201902ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1926965ÿÿÿÿÿ.272235
      ÿÿvar(e.b1x5)|ÿÿÿ.3594852ÿÿÿ.0272056ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.3099294ÿÿÿÿ.4169648
      ÿÿvar(e.b2x1)|ÿÿÿ.3300177ÿÿÿ.0246332ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2851031ÿÿÿÿ.3820081
      ÿÿvar(e.b2x2)|ÿÿÿ.3368326ÿÿÿ.0244336ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2921921ÿÿÿÿÿ.388293
      ÿÿvar(e.b2x3)|ÿÿÿ.3745263ÿÿÿ.0265833ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.3258858ÿÿÿÿ.4304268
      ÿÿvar(e.b2x4)|ÿÿÿ.3982713ÿÿÿ.0273655ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.3480906ÿÿÿÿÿ.455686
      ÿÿvar(e.b3x1)|ÿÿÿ.3026967ÿÿÿ.0218739ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2627224ÿÿÿÿ.3487533
      ÿÿvar(e.b3x2)|ÿÿÿ.2557326ÿÿÿ.0195659ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2201209ÿÿÿÿ.2971056
      ÿÿvar(e.b3x3)|ÿÿÿ.3191859ÿÿÿ.0228988ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2773176ÿÿÿÿ.3673754
      ÿÿvar(e.b3x4)|ÿÿÿ.2030464ÿÿÿ.0164061ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1733079ÿÿÿÿ.2378877
      ÿÿvar(e.b3x5)|ÿÿÿ.2863667ÿÿÿ.0211159ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2478319ÿÿÿÿ.3308932
      ÿÿÿÿvar(e.CL)|ÿÿÿ.2087328ÿÿÿ.0206698ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1719096ÿÿÿÿ.2534435
      ÿÿÿvar(e.SE1)|ÿÿÿ.1421682ÿÿÿ.0168967ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1126256ÿÿÿÿ.1794602
      ÿÿÿvar(e.SE2)|ÿÿÿ.1239188ÿÿÿ.0188404ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0919861ÿÿÿÿ.1669369
      ÿÿÿvar(e.SE3)|ÿÿÿ.0787265ÿÿÿ.0159291ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0529533ÿÿÿÿ.1170438
      ÿÿÿvar(e.GC1)|ÿÿÿ.3087343ÿÿÿÿ.043016ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.2349562ÿÿÿÿ.4056792
      ÿÿÿvar(e.GC2)|ÿÿÿ.0609172ÿÿÿ.0277761ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.0249242ÿÿÿÿ.1488874
      ÿÿÿvar(e.GC3)|ÿÿÿ.1688942ÿÿÿ.0280936ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.1219066ÿÿÿÿ.2339927
      ÿÿÿÿÿÿvar(SE)|ÿÿÿ.6047105ÿÿÿ.0783173ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4691444ÿÿÿÿ.7794503
      ÿÿÿÿÿÿvar(GC)|ÿÿÿ.4429002ÿÿÿ.0646439ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.3327116ÿÿÿÿ.5895814
      ------------------------------------------------------------------------------
      LRÿtestÿofÿmodelÿvs.ÿsaturated:ÿchi2(267)ÿ=ÿ1536.70ÿÿÿÿÿÿÿProbÿ>ÿchi2ÿ=ÿ0.0000

      .ÿassertÿe(converged)

      .ÿ
      .ÿmata:ÿrowsum(st_matrix("e(ilog)")ÿ:!=ÿ0)
      ÿÿ7

      .ÿ
      .ÿquietlyÿlogÿcloseÿlo

      .ÿ
      .ÿexit

      endÿofÿdo-file


      .
      Attached Files

      Comment


      • #18
        Big thanks to Erik for your advice. I will recheck everything!
        Diep

        Comment


        • #19
          Thank you very much for your patience. I will follow your advice!!!

          Comment


          • #20
            Joseph,

            Thank you very much for your support!!!

            Amazingly you made it work!

            Diep

            Comment

            Working...
            X