Mplus DEVELOPMENT (Dev 10/23/2011) MUTHEN & MUTHEN 10/23/2011 3:18 PM INPUT INSTRUCTIONS Title: Clinical Trials data M (intent) treated as continuous data: file = 4cat m.dat; variable: names = intent tx ciguse w; usev = tx ciguse intent; categorical = ciguse; freqweight = w; define: intent = (intent-1.456)/0.8854; Analysis: estimator = mlr; link = probit; model: [ciguse$1] (mbeta0); ciguse on tx (beta2) intent (beta1); [intent] (gamma0); intent on tx (gamma1); intent (sig2); model constraint: new(ind dir arg11 arg10 arg00 v1 v0 probit11 probit10 probit00 indirect direct orind ordir arg probit oddrat); dir=beta2; ind=beta1*gamma1; arg11=-mbeta0+beta2+beta1*(gamma0+gamma1); arg10=-mbeta0+beta2+beta1*gamma0; arg00=-mbeta0+beta1*gamma0; v1=beta1^2*sig2+1; v0=beta1^2*sig2+1; probit11=arg11/sqrt(v1); probit10=arg10/sqrt(v1); probit00=arg00/sqrt(v0); ! Phi function needed below: indirect=phi(probit11)-phi(probit10); direct=phi(probit10)-phi(probit00); orind=(phi(probit11)/(1-phi(probit11)))/(phi(probit10)/(1-phi(probit10))); ordir=(phi(probit10)/(1-phi(probit10)))/(phi(probit00)/(1-phi(probit00))); ! for naive indirect effect: arg=-mbeta0+beta1*gamma0+beta1*gamma1; probit=arg/sqrt(v0); oddrat=(phi(probit)/(1-phi(probit)))/(phi(probit00)/(1-phi(probit00))); output: tech1 tech8 sampstat patterns standardized; plot: type = plot3; *** WARNING in OUTPUT command PATTERNS option is available only when there are missing data. Request for PATTERNS is ignored. 1 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS Clinical Trials data M (intent) treated as continuous SUMMARY OF ANALYSIS Number of groups 1 Number of observations 864 Number of patterns 16 Number of dependent variables 2 Number of independent variables 1 Number of continuous latent variables 0 Observed dependent variables Continuous INTENT Binary and ordered categorical (ordinal) CIGUSE Observed independent variables TX Variables with special functions Weight variable W 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 10 Minimum value for logit thresholds -10 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 PROBIT Cholesky ON Input data file(s) 4cat m.dat Input data format FREE UNIVARIATE PROPORTIONS AND COUNTS FOR CATEGORICAL VARIABLES CIGUSE Category 1 0.819 708.000 Category 2 0.181 156.000 SAMPLE STATISTICS SAMPLE STATISTICS Means INTENT TX ________ ________ 1 0.000 0.571 Covariances INTENT TX ________ ________ INTENT 1.000 TX -0.046 0.245 Correlations INTENT TX ________ ________ INTENT 1.000 TX -0.092 1.000 THE MODEL ESTIMATION TERMINATED NORMALLY MODEL FIT INFORMATION Number of Free Parameters 6 Loglikelihood H0 Value -1557.334 H0 Scaling Correction Factor 1.195 for MLR Information Criteria Akaike (AIC) 3126.668 Bayesian (BIC) 3155.237 Sample-Size Adjusted BIC 3136.183 (n* = (n + 2) / 24) MODEL RESULTS Two-Tailed Estimate S.E. Est./S.E. P-Value CIGUSE ON TX -0.203 0.109 -1.867 0.062 INTENT 0.538 0.048 11.227 0.000 INTENT ON TX -0.186 0.070 -2.664 0.008 Intercepts INTENT 0.106 0.056 1.906 0.057 Thresholds CIGUSE$1 0.912 0.080 11.432 0.000 Residual Variances INTENT 0.990 0.069 14.291 0.000 New/Additional Parameters IND -0.100 0.038 -2.602 0.009 DIR -0.203 0.109 -1.867 0.062 ARG11 -1.158 0.079 -14.579 0.000 ARG10 -1.058 0.081 -13.072 0.000 ARG00 -0.855 0.085 -10.105 0.000 V1 1.287 0.055 23.545 0.000 V0 1.287 0.055 23.545 0.000 PROBIT11 -1.021 0.072 -14.240 0.000 PROBIT10 -0.933 0.075 -12.514 0.000 PROBIT00 -0.754 0.076 -9.947 0.000 INDIRECT -0.022 0.009 -2.548 0.011 DIRECT -0.050 0.027 -1.853 0.064 ORIND 0.853 0.051 16.587 0.000 ORDIR 0.731 0.123 5.941 0.000 ARG -0.955 0.084 -11.419 0.000 PROBIT -0.842 0.074 -11.447 0.000 ODDRAT 0.858 0.050 17.226 0.000 STANDARDIZED MODEL RESULTS STDYX Standardization Two-Tailed Estimate S.E. Est./S.E. P-Value CIGUSE ON TX -0.088 0.047 -1.882 0.060 INTENT 0.470 0.035 13.458 0.000 INTENT ON TX -0.092 0.034 -2.679 0.007 Intercepts INTENT 0.106 0.053 1.988 0.047 Thresholds CIGUSE$1 0.797 0.073 10.985 0.000 Residual Variances INTENT 0.992 0.006 156.871 0.000 STDY Standardization Two-Tailed Estimate S.E. Est./S.E. P-Value CIGUSE ON TX -0.177 0.094 -1.884 0.060 INTENT 0.470 0.035 13.458 0.000 INTENT ON TX -0.186 0.069 -2.685 0.007 Intercepts INTENT 0.106 0.053 1.988 0.047 Thresholds CIGUSE$1 0.797 0.073 10.985 0.000 Residual Variances INTENT 0.992 0.006 156.871 0.000 STD Standardization Two-Tailed Estimate S.E. Est./S.E. P-Value CIGUSE ON TX -0.203 0.109 -1.867 0.062 INTENT 0.538 0.048 11.227 0.000 INTENT ON TX -0.186 0.070 -2.664 0.008 Intercepts INTENT 0.106 0.056 1.906 0.057 Thresholds CIGUSE$1 0.912 0.080 11.432 0.000 Residual Variances INTENT 0.990 0.069 14.291 0.000 R-SQUARE Observed Two-Tailed Variable Estimate S.E. Est./S.E. P-Value CIGUSE 0.236 0.034 7.020 0.000 INTENT 0.008 0.006 1.340 0.180 QUALITY OF NUMERICAL RESULTS Condition Number for the Information Matrix 0.215E-03 (ratio of smallest to largest eigenvalue) TECHNICAL 1 OUTPUT PARAMETER SPECIFICATION TAU CIGUSE$1 ________ 1 6 NU CIGUSE INTENT TX ________ ________ ________ 1 0 0 0 LAMBDA CIGUSE INTENT TX ________ ________ ________ CIGUSE 0 0 0 INTENT 0 0 0 TX 0 0 0 THETA CIGUSE INTENT TX ________ ________ ________ CIGUSE 0 INTENT 0 0 TX 0 0 0 ALPHA CIGUSE INTENT TX ________ ________ ________ 1 0 1 0 BETA CIGUSE INTENT TX ________ ________ ________ CIGUSE 0 2 3 INTENT 0 0 4 TX 0 0 0 PSI CIGUSE INTENT TX ________ ________ ________ CIGUSE 0 INTENT 0 5 TX 0 0 0 PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS NEW/ADDITIONAL PARAMETERS IND DIR ARG11 ARG10 ARG00 ________ ________ ________ ________ ________ 1 7 8 9 10 11 NEW/ADDITIONAL PARAMETERS V1 V0 PROBIT11 PROBIT10 PROBIT00 ________ ________ ________ ________ ________ 1 12 13 14 15 16 NEW/ADDITIONAL PARAMETERS INDIRECT DIRECT ORIND ORDIR ARG ________ ________ ________ ________ ________ 1 17 18 19 20 21 NEW/ADDITIONAL PARAMETERS PROBIT ODDRAT ________ ________ 1 22 23 STARTING VALUES TAU CIGUSE$1 ________ 1 0.840 NU CIGUSE INTENT TX ________ ________ ________ 1 0.000 0.000 0.000 LAMBDA CIGUSE INTENT TX ________ ________ ________ CIGUSE 1.000 0.000 0.000 INTENT 0.000 1.000 0.000 TX 0.000 0.000 1.000 THETA CIGUSE INTENT TX ________ ________ ________ CIGUSE 0.000 INTENT 0.000 0.000 TX 0.000 0.000 0.000 ALPHA CIGUSE INTENT TX ________ ________ ________ 1 0.000 0.000 0.000 BETA CIGUSE INTENT TX ________ ________ ________ CIGUSE 0.000 0.000 0.000 INTENT 0.000 0.000 0.000 TX 0.000 0.000 0.000 PSI CIGUSE INTENT TX ________ ________ ________ CIGUSE 1.000 INTENT 0.000 0.500 TX 0.000 0.000 0.123 STARTING VALUES FOR THE ADDITIONAL PARAMETERS NEW/ADDITIONAL PARAMETERS IND DIR ARG11 ARG10 ARG00 ________ ________ ________ ________ ________ 1 0.500 0.500 0.500 0.500 0.500 NEW/ADDITIONAL PARAMETERS V1 V0 PROBIT11 PROBIT10 PROBIT00 ________ ________ ________ ________ ________ 1 0.500 0.500 0.500 0.500 0.500 NEW/ADDITIONAL PARAMETERS INDIRECT DIRECT ORIND ORDIR ARG ________ ________ ________ ________ ________ 1 0.500 0.500 0.500 0.500 0.500 NEW/ADDITIONAL PARAMETERS PROBIT ODDRAT ________ ________ 1 0.500 0.500 TECHNICAL 8 OUTPUT E STEP ITER LOGLIKELIHOOD ABS CHANGE REL CHANGE ALGORITHM 1 -0.17664884D+04 0.0000000 0.0000000 EM 2 -0.17128090D+04 53.6793732 0.0303876 FS 3 -0.15573624D+04 155.4465808 0.0907554 EM 4 -0.15573339D+04 0.0284823 0.0000183 FS 5 -0.15573338D+04 0.0001253 0.0000001 EM PLOT INFORMATION The following plots are available: Histograms (sample values) Scatterplots (sample values) Sample proportions Beginning Time: 15:18:07 Ending Time: 15:18:07 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-2011 Muthen & Muthen