Mplus VERSION 7.3
MUTHEN & MUTHEN
09/22/2014   5:18 PM

INPUT INSTRUCTIONS

  title:	this is an example of a discrete-time
  	survival analysis
  montecarlo:
  	names = u1-u4 x;
  	generate = u1-u4(1);
  	categorical = u1-u4;
      missing = u1-u4;
  	nobs = 500;
  	nreps = 1;
  	save = ex6.19.dat;

  model population:
      [x@0]; x@1;
      [u1$1*2 u2$1*1.5 u3$1*1 u4$1*1];
  	f by u1-u4@1;
      f on x*.5;
      f@0;

  model missing:
      [u2-u4@-15]; ! Low logit intercepts give P(missing) = 0
      u2-u4 on u1@30; ! if u1=1, later u's are missing with certainty
      u3-u4 on u2@30; ! if u2=1, ...
      u4 on u3@30;  ! if u3=1, ...

  analysis:
      estimator = mlr;

  model:
      [u1$1*2 u2$1*1.5 u3$1*1 u4$1*1];
  	f by u1-u4@1;
      f on x*.5;
      f@0;

  output:
  	tech8 tech9;




INPUT READING TERMINATED NORMALLY



this is an example of a discrete-time
survival analysis

SUMMARY OF ANALYSIS

Number of groups                                                 1
Number of observations                                         500

Number of replications
    Requested                                                    1
    Completed                                                    1
Value of seed                                                    0

Number of dependent variables                                    4
Number of independent variables                                  1
Number of continuous latent variables                            1

Observed dependent variables

  Binary and ordered categorical (ordinal)
   U1          U2          U3          U4

Observed independent variables
   X

Continuous latent variables
   F


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
Maximum number of iterations for H1                           2000
Convergence criterion for H1                             0.100D-03
Optimization algorithm                                         EMA
Integration Specifications
  Type                                                    STANDARD
  Number of integration points                                  15
  Dimensions of numerical integration                            0
  Adaptive quadrature                                           ON
Link                                                         LOGIT
Cholesky                                                        ON


SUMMARY OF DATA FOR THE FIRST REPLICATION

     Number of missing data patterns             7


SUMMARY OF MISSING DATA PATTERNS FOR THE FIRST REPLICATION


     MISSING DATA PATTERNS FOR U (x = not missing)

           1  2  3  4  5  6  7
 U1              x  x  x  x
 U2        x  x  x  x  x     x
 U3        x     x  x        x
 U4        x        x


     MISSING DATA PATTERN FREQUENCIES FOR U

    Pattern   Frequency     Pattern   Frequency     Pattern   Frequency
          1         148           4         112           7          51
          2          40           5          59
          3          55           6          35


COVARIANCE COVERAGE OF DATA FOR THE FIRST REPLICATION

Minimum covariance coverage value   0.100


     PROPORTION OF DATA PRESENT FOR U


           Covariance Coverage
              U1            U2            U3            U4
              ________      ________      ________      ________
 U1             0.522
 U2             0.452         0.930
 U3             0.334         0.732         0.732
 U4             0.224         0.520         0.520         0.520


SAMPLE STATISTICS FOR THE FIRST REPLICATION


     SAMPLE STATISTICS


           Means
              X
              ________
 1              0.033


           Covariances
              X
              ________
 X              0.978


           Correlations
              X
              ________
 X              1.000




MODEL FIT INFORMATION

Number of Free Parameters                        5

Loglikelihood

    H0 Value

        Mean                              -678.342
        Std Dev                              0.000
        Number of successful computations        1

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

Information Criteria

    Akaike (AIC)

        Mean                              1366.683
        Std Dev                              0.000
        Number of successful computations        1

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

    Bayesian (BIC)

        Mean                              1387.756
        Std Dev                              0.000
        Number of successful computations        1

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

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

        Mean                              1371.886
        Std Dev                              0.000
        Number of successful computations        1

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



MODEL RESULTS

                              ESTIMATES              S. E.     M. S. E.  95%  % Sig
                 Population   Average   Std. Dev.   Average             Cover Coeff
 F        BY
  U1                  1.000     1.0000     0.0000     0.0000     0.0000 1.000 0.000
  U2                  1.000     1.0000     0.0000     0.0000     0.0000 1.000 0.000
  U3                  1.000     1.0000     0.0000     0.0000     0.0000 1.000 0.000
  U4                  1.000     1.0000     0.0000     0.0000     0.0000 1.000 0.000

 F          ON
  X                   0.500     0.5348     0.0000     0.0689     0.0012 1.000 1.000

 Thresholds
  U1$1                2.000     1.9970     0.0000     0.1869     0.0000 1.000 1.000
  U2$1                1.500     1.3828     0.0000     0.1189     0.0137 1.000 1.000
  U3$1                1.000     0.9016     0.0000     0.1195     0.0097 1.000 1.000
  U4$1                1.000     1.0834     0.0000     0.1451     0.0070 1.000 1.000

 Residual Variances
  F                   0.000     0.0000     0.0000     0.0000     0.0000 1.000 0.000


QUALITY OF NUMERICAL RESULTS

     Average Condition Number for the Information Matrix      0.136E+00
       (ratio of smallest to largest eigenvalue)


TECHNICAL 1 OUTPUT


     PARAMETER SPECIFICATION


           TAU
              U1$1          U2$1          U3$1          U4$1
              ________      ________      ________      ________
 1                  2             3             4             5


           NU
              U1            U2            U3            U4            X
              ________      ________      ________      ________      ________
 1                  0             0             0             0             0


           LAMBDA
              F             X
              ________      ________
 U1                 0             0
 U2                 0             0
 U3                 0             0
 U4                 0             0
 X                  0             0


           THETA
              U1            U2            U3            U4            X
              ________      ________      ________      ________      ________
 U1                 0
 U2                 0             0
 U3                 0             0             0
 U4                 0             0             0             0
 X                  0             0             0             0             0


           ALPHA
              F             X
              ________      ________
 1                  0             0


           BETA
              F             X
              ________      ________
 F                  0             1
 X                  0             0


           PSI
              F             X
              ________      ________
 F                  0
 X                  0             0


     STARTING VALUES


           TAU
              U1$1          U2$1          U3$1          U4$1
              ________      ________      ________      ________
 1              2.000         1.500         1.000         1.000


           NU
              U1            U2            U3            U4            X
              ________      ________      ________      ________      ________
 1              0.000         0.000         0.000         0.000         0.000


           LAMBDA
              F             X
              ________      ________
 U1             1.000         0.000
 U2             1.000         0.000
 U3             1.000         0.000
 U4             1.000         0.000
 X              0.000         1.000


           THETA
              U1            U2            U3            U4            X
              ________      ________      ________      ________      ________
 U1             1.000
 U2             0.000         1.000
 U3             0.000         0.000         1.000
 U4             0.000         0.000         0.000         1.000
 X              0.000         0.000         0.000         0.000         0.000


           ALPHA
              F             X
              ________      ________
 1              0.000         0.000


           BETA
              F             X
              ________      ________
 F              0.000         0.500
 X              0.000         0.000


           PSI
              F             X
              ________      ________
 F              0.000
 X              0.000         0.500


     POPULATION VALUES


           TAU
              U1$1          U2$1          U3$1          U4$1
              ________      ________      ________      ________
 1              2.000         1.500         1.000         1.000


           NU
              U1            U2            U3            U4            X
              ________      ________      ________      ________      ________
 1              0.000         0.000         0.000         0.000         0.000


           LAMBDA
              F             X
              ________      ________
 U1             1.000         0.000
 U2             1.000         0.000
 U3             1.000         0.000
 U4             1.000         0.000
 X              0.000         1.000


           THETA
              U1            U2            U3            U4            X
              ________      ________      ________      ________      ________
 U1             0.000
 U2             0.000         0.000
 U3             0.000         0.000         0.000
 U4             0.000         0.000         0.000         0.000
 X              0.000         0.000         0.000         0.000         0.000


           ALPHA
              F             X
              ________      ________
 1              0.000         0.000


           BETA
              F             X
              ________      ________
 F              0.000         0.500
 X              0.000         0.000


           PSI
              F             X
              ________      ________
 F              0.000
 X              0.000         1.000


TECHNICAL 8 OUTPUT


  TECHNICAL 8 OUTPUT FOR REPLICATION 1


   E STEP  ITER  LOGLIKELIHOOD    ABS CHANGE   REL CHANGE  ALGORITHM
              1 -0.67957786D+03    0.0000000    0.0000000  EM
              2 -0.67834260D+03    1.2352647    0.0018177  EM
              3 -0.67834172D+03    0.0008733    0.0000013  EM
              4 -0.67834172D+03    0.0000000    0.0000000  EM


TECHNICAL 9 OUTPUT

  Error messages for each replication (if any)



SAVEDATA INFORMATION

  Order of variables

    U1
    U2
    U3
    U4
    X

  Save file
    ex6.19.dat

  Save file format           Free
  Save file record length    10000
  Missing designated by 999


     Beginning Time:  17:18:27
        Ending Time:  17:18:27
       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