Mplus VERSION 7.2
MUTHEN & MUTHEN
05/07/2014   2:03 PM

INPUT INSTRUCTIONS

  TITLE:		this is an example of a Monte Carlo
  			simulation study this is an example of a path analysis
  			with a combination of censored, categorical,
  			and unordered categorical dependent variables
  MONTECARLO:			
  			NAMES = y1 u1 u2 x1-x3;
  			NOBSERVATIONS = 500;
  			NREPS = 1;
  			SEED = 53487;
  			GENERATE = y1(ca 1) u1(1) u2(n 2);
  			CENSORED = y1(a);
  			CATEGORICAL = u1;
  			NOMINAL = u2;
  			SAVE = ex3.15.dat;	
  MODEL POPULATION:
  			[x1-x3@0];
  			x1-x3@1;
  			y1 ON x1*-1 x2*.5 x3*1;
  			y1*1;
  			[y1*.5];
  			u1 ON x1*1 x2*-.5 x3*-1;
  			u2#1 ON y1*.5 u1*1 x2*-1;
  			u2#2 ON y1*.-5 u1*-1 x2*1;
  MODEL:
  			y1 ON x1*-1 x2*.5 x3*1;
  			y1*1;
  			[y1*.5];
  			u1 ON x1*1 x2*-.5 x3*-1;
  			u2#1 ON y1*.5 u1*1 x2*-1;
  			u2#2 ON y1*.-5 u1*-1 x2*1;
  OUTPUT:		TECH8 TECH9;



INPUT READING TERMINATED NORMALLY



this is an example of a Monte Carlo
simulation study this is an example of a path analysis
with a combination of censored, categorical,
and unordered categorical dependent variables

SUMMARY OF ANALYSIS

Number of groups                                                 1
Number of observations                                         500

Number of replications
    Requested                                                    1
    Completed                                                    1
Value of seed                                                53487

Number of dependent variables                                    3
Number of independent variables                                  3
Number of continuous latent variables                            0

Observed dependent variables

  Censored
   Y1

  Binary and ordered categorical (ordinal)
   U1

  Unordered categorical (nominal)
   U2

Observed independent variables
   X1          X2          X3


Estimator                                                      MLR
Information matrix                                        OBSERVED
Optimization Specifications for the Quasi-Newton Algorithm for
Continuous Outcomes
  Maximum number of iterations                                 100
  Convergence criterion                                  0.100D-05
Optimization Specifications for the EM Algorithm
  Maximum number of iterations                                 500
  Convergence criteria
    Loglikelihood change                                 0.100D-02
    Relative loglikelihood change                        0.100D-05
    Derivative                                           0.100D-02
Optimization Specifications for the M step of the EM Algorithm for
Categorical Latent variables
  Number of M step iterations                                    1
  M step convergence criterion                           0.100D-02
  Basis for M step termination                           ITERATION
Optimization Specifications for the M step of the EM Algorithm for
Censored, Binary or Ordered Categorical (Ordinal), Unordered
Categorical (Nominal) and Count Outcomes
  Number of M step iterations                                    1
  M step convergence criterion                           0.100D-02
  Basis for M step termination                           ITERATION
  Maximum value for logit thresholds                            15
  Minimum value for logit thresholds                           -15
  Minimum expected cell size for chi-square              0.100D-01
Optimization algorithm                                         EMA
Integration Specifications
  Type                                                    STANDARD
  Number of integration points                                  15
  Dimensions of numerical integration                            0
  Adaptive quadrature                                           ON
Link                                                         LOGIT
Cholesky                                                       OFF


SUMMARY OF CENSORED LIMITS

      Y1                 1.000


SAMPLE STATISTICS FOR THE FIRST REPLICATION


     SAMPLE STATISTICS


           Means
              X1            X2            X3
              ________      ________      ________
 1             -0.009        -0.014        -0.029


           Covariances
              X1            X2            X3
              ________      ________      ________
 X1             1.026
 X2             0.013         0.972
 X3            -0.086         0.020         1.047


           Correlations
              X1            X2            X3
              ________      ________      ________
 X1             1.000
 X2             0.013         1.000
 X3            -0.083         0.020         1.000




MODEL FIT INFORMATION

Number of Free Parameters                       17

Loglikelihood

    H0 Value

        Mean                             -1181.376
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       0.000        -1181.376      -1181.376
           0.980       0.000        -1181.376      -1181.376
           0.950       0.000        -1181.376      -1181.376
           0.900       0.000        -1181.376      -1181.376
           0.800       0.000        -1181.376      -1181.376
           0.700       0.000        -1181.376      -1181.376
           0.500       0.000        -1181.376      -1181.376
           0.300       0.000        -1181.376      -1181.376
           0.200       0.000        -1181.376      -1181.376
           0.100       0.000        -1181.376      -1181.376
           0.050       0.000        -1181.376      -1181.376
           0.020       0.000        -1181.376      -1181.376
           0.010       0.000        -1181.376      -1181.376

Information Criteria

    Akaike (AIC)

        Mean                              2396.751
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       0.000         2396.751       2396.751
           0.980       0.000         2396.751       2396.751
           0.950       0.000         2396.751       2396.751
           0.900       0.000         2396.751       2396.751
           0.800       0.000         2396.751       2396.751
           0.700       0.000         2396.751       2396.751
           0.500       0.000         2396.751       2396.751
           0.300       0.000         2396.751       2396.751
           0.200       0.000         2396.751       2396.751
           0.100       0.000         2396.751       2396.751
           0.050       0.000         2396.751       2396.751
           0.020       0.000         2396.751       2396.751
           0.010       0.000         2396.751       2396.751

    Bayesian (BIC)

        Mean                              2468.399
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       0.000         2468.399       2468.399
           0.980       0.000         2468.399       2468.399
           0.950       0.000         2468.399       2468.399
           0.900       0.000         2468.399       2468.399
           0.800       0.000         2468.399       2468.399
           0.700       0.000         2468.399       2468.399
           0.500       0.000         2468.399       2468.399
           0.300       0.000         2468.399       2468.399
           0.200       0.000         2468.399       2468.399
           0.100       0.000         2468.399       2468.399
           0.050       0.000         2468.399       2468.399
           0.020       0.000         2468.399       2468.399
           0.010       0.000         2468.399       2468.399

    Sample-Size Adjusted BIC (n* = (n + 2) / 24)

        Mean                              2414.440
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       0.000         2414.440       2414.440
           0.980       0.000         2414.440       2414.440
           0.950       0.000         2414.440       2414.440
           0.900       0.000         2414.440       2414.440
           0.800       0.000         2414.440       2414.440
           0.700       0.000         2414.440       2414.440
           0.500       0.000         2414.440       2414.440
           0.300       0.000         2414.440       2414.440
           0.200       0.000         2414.440       2414.440
           0.100       0.000         2414.440       2414.440
           0.050       0.000         2414.440       2414.440
           0.020       0.000         2414.440       2414.440
           0.010       0.000         2414.440       2414.440



MODEL RESULTS

                              ESTIMATES              S. E.     M. S. E.  95%  % Sig
                 Population   Average   Std. Dev.   Average             Cover Coeff
 Y1         ON
  X1                 -1.000    -0.9916     0.0000     0.0542     0.0001 1.000 1.000
  X2                  0.500     0.4883     0.0000     0.0483     0.0001 1.000 1.000
  X3                  1.000     0.9200     0.0000     0.0489     0.0064 1.000 1.000

 U1         ON
  X1                  1.000     1.0306     0.0000     0.1210     0.0009 1.000 1.000
  X2                 -0.500    -0.5988     0.0000     0.1241     0.0098 1.000 1.000
  X3                 -1.000    -0.9251     0.0000     0.1165     0.0056 1.000 1.000

 U2#1       ON
  Y1                  0.500     0.4504     0.0000     0.0965     0.0025 1.000 1.000
  U1                  1.000     1.5030     0.0000     0.2577     0.2530 1.000 1.000
  X2                 -1.000    -0.6650     0.0000     0.1645     0.1122 0.000 1.000

 U2#2       ON
  Y1                  0.000    -0.0203     0.0000     0.1535     0.0004 1.000 0.000
  U1                 -1.000    -0.8865     0.0000     0.3595     0.0129 1.000 1.000
  X2                  1.000     1.7306     0.0000     0.2189     0.5337 0.000 1.000

 Intercepts
  Y1                  0.500     0.4168     0.0000     0.0535     0.0069 1.000 1.000
  U2#1                0.000    -0.4963     0.0000     0.1757     0.2463 0.000 1.000
  U2#2                0.000    -0.5441     0.0000     0.2005     0.2961 0.000 1.000

 Thresholds
  U1$1                0.000     0.0941     0.0000     0.1086     0.0089 1.000 0.000

 Residual Variances
  Y1                  1.000     0.9042     0.0000     0.0726     0.0092 1.000 1.000


QUALITY OF NUMERICAL RESULTS

     Average Condition Number for the Information Matrix      0.100E-01
       (ratio of smallest to largest eigenvalue)


TECHNICAL 1 OUTPUT


     PARAMETER SPECIFICATION


           TAU
              U1$1
              ________
 1                 17


           NU
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 1                  0             0             0             0             0


           NU
              X1            X2            X3
              ________      ________      ________
 1                  0             0             0


           LAMBDA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1               0             0             0             0             0
 Y1                 0             0             0             0             0
 U2#1               0             0             0             0             0
 U2#2               0             0             0             0             0
 U1                 0             0             0             0             0
 X1                 0             0             0             0             0
 X2                 0             0             0             0             0
 X3                 0             0             0             0             0


           LAMBDA
              X1            X2            X3
              ________      ________      ________
 Y1#1               0             0             0
 Y1                 0             0             0
 U2#1               0             0             0
 U2#2               0             0             0
 U1                 0             0             0
 X1                 0             0             0
 X2                 0             0             0
 X3                 0             0             0


           THETA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1               0
 Y1                 0             0
 U2#1               0             0             0
 U2#2               0             0             0             0
 U1                 0             0             0             0             0
 X1                 0             0             0             0             0
 X2                 0             0             0             0             0
 X3                 0             0             0             0             0


           THETA
              X1            X2            X3
              ________      ________      ________
 X1                 0
 X2                 0             0
 X3                 0             0             0


           ALPHA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 1                  0             1             2             3             0


           ALPHA
              X1            X2            X3
              ________      ________      ________
 1                  0             0             0


           BETA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1               0             0             0             0             0
 Y1                 0             0             0             0             0
 U2#1               0             7             0             0             8
 U2#2               0            10             0             0            11
 U1                 0             0             0             0             0
 X1                 0             0             0             0             0
 X2                 0             0             0             0             0
 X3                 0             0             0             0             0


           BETA
              X1            X2            X3
              ________      ________      ________
 Y1#1               0             0             0
 Y1                 4             5             6
 U2#1               0             9             0
 U2#2               0            12             0
 U1                13            14            15
 X1                 0             0             0
 X2                 0             0             0
 X3                 0             0             0


           PSI
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1               0
 Y1                 0            16
 U2#1               0             0             0
 U2#2               0             0             0             0
 U1                 0             0             0             0             0
 X1                 0             0             0             0             0
 X2                 0             0             0             0             0
 X3                 0             0             0             0             0


           PSI
              X1            X2            X3
              ________      ________      ________
 X1                 0
 X2                 0             0
 X3                 0             0             0


     STARTING VALUES


           TAU
              U1$1
              ________
 1              0.000


           NU
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 1              0.000         0.000         0.000         0.000         0.000


           NU
              X1            X2            X3
              ________      ________      ________
 1              0.000         0.000         0.000


           LAMBDA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           1.000         0.000         0.000         0.000         0.000
 Y1             0.000         1.000         0.000         0.000         0.000
 U2#1           0.000         0.000         1.000         0.000         0.000
 U2#2           0.000         0.000         0.000         1.000         0.000
 U1             0.000         0.000         0.000         0.000         1.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           LAMBDA
              X1            X2            X3
              ________      ________      ________
 Y1#1           0.000         0.000         0.000
 Y1             0.000         0.000         0.000
 U2#1           0.000         0.000         0.000
 U2#2           0.000         0.000         0.000
 U1             0.000         0.000         0.000
 X1             1.000         0.000         0.000
 X2             0.000         1.000         0.000
 X3             0.000         0.000         1.000


           THETA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           0.000
 Y1             0.000         0.000
 U2#1           0.000         0.000         0.000
 U2#2           0.000         0.000         0.000         0.000
 U1             0.000         0.000         0.000         0.000         0.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           THETA
              X1            X2            X3
              ________      ________      ________
 X1             0.000
 X2             0.000         0.000
 X3             0.000         0.000         0.000


           ALPHA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 1            -20.000         0.500         0.000         0.000         0.000


           ALPHA
              X1            X2            X3
              ________      ________      ________
 1              0.000         0.000         0.000


           BETA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           0.000         0.000         0.000         0.000         0.000
 Y1             0.000         0.000         0.000         0.000         0.000
 U2#1           0.000         0.500         0.000         0.000         1.000
 U2#2           0.000         0.000         0.000         0.000        -1.000
 U1             0.000         0.000         0.000         0.000         0.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           BETA
              X1            X2            X3
              ________      ________      ________
 Y1#1           0.000         0.000         0.000
 Y1            -1.000         0.500         1.000
 U2#1           0.000        -1.000         0.000
 U2#2           0.000         1.000         0.000
 U1             1.000        -0.500        -1.000
 X1             0.000         0.000         0.000
 X2             0.000         0.000         0.000
 X3             0.000         0.000         0.000


           PSI
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           0.000
 Y1             0.000         1.000
 U2#1           0.000         0.000         0.000
 U2#2           0.000         0.000         0.000         0.000
 U1             0.000         0.000         0.000         0.000         1.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           PSI
              X1            X2            X3
              ________      ________      ________
 X1             0.500
 X2             0.000         0.500
 X3             0.000         0.000         0.500


     POPULATION VALUES


           TAU
              U1$1
              ________
 1              0.000


           NU
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 1              0.000         0.000         0.000         0.000         0.000


           NU
              X1            X2            X3
              ________      ________      ________
 1              0.000         0.000         0.000


           LAMBDA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           1.000         0.000         0.000         0.000         0.000
 Y1             0.000         1.000         0.000         0.000         0.000
 U2#1           0.000         0.000         1.000         0.000         0.000
 U2#2           0.000         0.000         0.000         1.000         0.000
 U1             0.000         0.000         0.000         0.000         1.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           LAMBDA
              X1            X2            X3
              ________      ________      ________
 Y1#1           0.000         0.000         0.000
 Y1             0.000         0.000         0.000
 U2#1           0.000         0.000         0.000
 U2#2           0.000         0.000         0.000
 U1             0.000         0.000         0.000
 X1             1.000         0.000         0.000
 X2             0.000         1.000         0.000
 X3             0.000         0.000         1.000


           THETA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           0.000
 Y1             0.000         0.000
 U2#1           0.000         0.000         0.000
 U2#2           0.000         0.000         0.000         0.000
 U1             0.000         0.000         0.000         0.000         0.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           THETA
              X1            X2            X3
              ________      ________      ________
 X1             0.000
 X2             0.000         0.000
 X3             0.000         0.000         0.000


           ALPHA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 1            -20.000         0.500         0.000         0.000         0.000


           ALPHA
              X1            X2            X3
              ________      ________      ________
 1              0.000         0.000         0.000


           BETA
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           0.000         0.000         0.000         0.000         0.000
 Y1             0.000         0.000         0.000         0.000         0.000
 U2#1           0.000         0.500         0.000         0.000         1.000
 U2#2           0.000         0.000         0.000         0.000        -1.000
 U1             0.000         0.000         0.000         0.000         0.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           BETA
              X1            X2            X3
              ________      ________      ________
 Y1#1           0.000         0.000         0.000
 Y1            -1.000         0.500         1.000
 U2#1           0.000        -1.000         0.000
 U2#2           0.000         1.000         0.000
 U1             1.000        -0.500        -1.000
 X1             0.000         0.000         0.000
 X2             0.000         0.000         0.000
 X3             0.000         0.000         0.000


           PSI
              Y1#1          Y1            U2#1          U2#2          U1
              ________      ________      ________      ________      ________
 Y1#1           0.000
 Y1             0.000         1.000
 U2#1           0.000         0.000         0.000
 U2#2           0.000         0.000         0.000         0.000
 U1             0.000         0.000         0.000         0.000         0.000
 X1             0.000         0.000         0.000         0.000         0.000
 X2             0.000         0.000         0.000         0.000         0.000
 X3             0.000         0.000         0.000         0.000         0.000


           PSI
              X1            X2            X3
              ________      ________      ________
 X1             1.000
 X2             0.000         1.000
 X3             0.000         0.000         1.000


TECHNICAL 8 OUTPUT


  TECHNICAL 8 OUTPUT FOR REPLICATION 1


   E STEP  ITER  LOGLIKELIHOOD    ABS CHANGE   REL CHANGE  ALGORITHM
              1 -0.11978265D+04    0.0000000    0.0000000  EM
              2 -0.11818891D+04   15.9373501    0.0133052  EM
              3 -0.11813763D+04    0.5128394    0.0004339  EM
              4 -0.11813756D+04    0.0007050    0.0000006  EM


TECHNICAL 9 OUTPUT

  Error messages for each replication (if any)



SAVEDATA INFORMATION

  Order of variables

    Y1
    U2
    U1
    X1
    X2
    X3

  Save file
    ex3.15.dat

  Save file format           Free
  Save file record length    10000


     Beginning Time:  14:03:44
        Ending Time:  14:03:44
       Elapsed Time:  00:00:00



MUTHEN & MUTHEN
3463 Stoner Ave.
Los Angeles, CA  90066

Tel: (310) 391-9971
Fax: (310) 391-8971
Web: www.StatModel.com
Support: Support@StatModel.com

Copyright (c) 1998-2014 Muthen & Muthen

Back to examples