Mplus VERSION 7.4
MUTHEN & MUTHEN
06/02/2016 5:16 PM
INPUT INSTRUCTIONS
Title:
Clinical Trials data
M (intent) original 4-category mediator version
data:
file = smoking.txt;
variable:
names = intent tx ciguse;
usev = tx ciguse intent;
categorical = ciguse intent;
Analysis:
estimator = ml;
bootstrap = 1000;
link = probit;
model:
ciguse on intent tx;
intent on tx;
model indirect: ciguse ind intent tx;
output:
tech1 tech8 sampstat patterns cinterval(bootstrap);
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) original 4-category mediator version
SUMMARY OF ANALYSIS
Number of groups 1
Number of observations 864
Number of dependent variables 2
Number of independent variables 1
Number of continuous latent variables 0
Observed dependent variables
Binary and ordered categorical (ordinal)
CIGUSE INTENT
Observed independent variables
TX
Estimator ML
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
Number of bootstrap draws
Requested 1000
Completed 1000
Optimization algorithm EMA
Integration Specifications
Type STANDARD
Number of integration points 15
Dimensions of numerical integration 0
Adaptive quadrature ON
Link PROBIT
Cholesky OFF
Input data file(s)
smoking.txt
Input data format FREE
UNIVARIATE PROPORTIONS AND COUNTS FOR CATEGORICAL VARIABLES
CIGUSE
Category 1 0.819 708.000
Category 2 0.181 156.000
INTENT
Category 1 0.745 644.000
Category 2 0.119 103.000
Category 3 0.069 60.000
Category 4 0.066 57.000
SAMPLE STATISTICS
SAMPLE STATISTICS
Means
TX
________
1 0.571
Covariances
TX
________
TX 0.245
Correlations
TX
________
TX 1.000
UNIVARIATE SAMPLE STATISTICS
UNIVARIATE HIGHER-ORDER MOMENT DESCRIPTIVE STATISTICS
Variable/ Mean/ Skewness/ Minimum/ % with Percentiles
Sample Size Variance Kurtosis Maximum Min/Max 20%/60% 40%/80% Median
TX 0.571 -0.285 0.000 42.94% 0.000 0.000 1.000
864.000 0.245 -1.919 1.000 57.06% 1.000 1.000
THE MODEL ESTIMATION TERMINATED NORMALLY
MODEL FIT INFORMATION
Number of Free Parameters 7
Loglikelihood
H0 Value -1055.162
Information Criteria
Akaike (AIC) 2124.325
Bayesian (BIC) 2157.656
Sample-Size Adjusted BIC 2135.426
(n* = (n + 2) / 24)
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
CIGUSE ON
INTENT 0.608 0.053 11.571 0.000
TX -0.203 0.106 -1.913 0.056
INTENT ON
TX -0.246 0.088 -2.794 0.005
Thresholds
CIGUSE$1 1.189 0.089 13.365 0.000
INTENT$1 0.525 0.066 7.914 0.000
INTENT$2 0.970 0.072 13.425 0.000
INTENT$3 1.378 0.081 17.037 0.000
TOTAL, INDIRECT, AND DIRECT EFFECTS BASED ON COUNTERFACTUALS (CAUSALLY-DEFINED EFFECTS)
Two-Tailed
Estimate S.E. Est./S.E. P-Value
Effects from TX to CIGUSE
Tot natural IE -0.030 0.011 -2.710 0.007
Pure natural DE -0.046 0.024 -1.903 0.057
Total effect -0.076 0.026 -2.911 0.004
Odds ratios for binary Y
Tot natural IE 0.803 0.064 12.494 0.000
Pure natural DE 0.751 0.116 6.463 0.000
Total effect 0.603 0.108 5.610 0.000
Other effects
Pure natural IE -0.033 0.012 -2.728 0.006
Tot natural DE -0.043 0.023 -1.898 0.058
Total effect -0.076 0.026 -2.911 0.004
Odds ratios for other effects for binary Y
Pure natural IE 0.819 0.060 13.711 0.000
Tot natural DE 0.736 0.122 6.029 0.000
Total effect 0.603 0.108 5.610 0.000
CONFIDENCE INTERVALS OF MODEL RESULTS
Lower .5% Lower 2.5% Lower 5% Estimate Upper 5% Upper 2.5% Upper .5%
CIGUSE ON
INTENT 0.468 0.509 0.527 0.608 0.700 0.716 0.743
TX -0.481 -0.410 -0.368 -0.203 -0.016 0.018 0.068
INTENT ON
TX -0.477 -0.418 -0.394 -0.246 -0.097 -0.074 -0.026
Thresholds
CIGUSE$1 0.988 1.018 1.049 1.189 1.349 1.374 1.433
INTENT$1 0.348 0.395 0.418 0.525 0.638 0.659 0.708
INTENT$2 0.794 0.836 0.861 0.970 1.095 1.114 1.159
INTENT$3 1.184 1.235 1.256 1.378 1.524 1.554 1.603
CONFIDENCE INTERVALS OF TOTAL, INDIRECT, AND DIRECT EFFECTS BASED ON COUNTERFACTUALS (CAUSALLY-DEFINED EFFECTS)
Lower .5% Lower 2.5% Lower 5% Estimate Upper 5% Upper 2.5% Upper .5%
Effects from TX to CIGUSE
Tot natural IE -0.059 -0.052 -0.048 -0.030 -0.012 -0.009 -0.003
Pure natural DE -0.106 -0.093 -0.085 -0.046 -0.004 0.004 0.015
Total effect -0.146 -0.125 -0.118 -0.076 -0.030 -0.023 -0.006
Odds ratios for binary Y
Tot natural IE 0.646 0.685 0.705 0.803 0.916 0.938 0.976
Pure natural DE 0.511 0.559 0.592 0.751 0.976 1.027 1.104
Total effect 0.385 0.436 0.455 0.603 0.815 0.858 0.958
Other effects
Pure natural IE -0.066 -0.057 -0.053 -0.033 -0.013 -0.010 -0.004
Tot natural DE -0.099 -0.088 -0.079 -0.043 -0.003 0.004 0.015
Total effect -0.146 -0.125 -0.118 -0.076 -0.030 -0.023 -0.006
Odds ratios for other effects for binary Y
Pure natural IE 0.681 0.707 0.726 0.819 0.924 0.943 0.979
Tot natural DE 0.488 0.540 0.570 0.736 0.974 1.028 1.108
Total effect 0.385 0.436 0.455 0.603 0.815 0.858 0.958
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
TAU
CIGUSE$1 INTENT$1 INTENT$2 INTENT$3
________ ________ ________ ________
1 4 5 6 7
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 0 0
BETA
CIGUSE INTENT TX
________ ________ ________
CIGUSE 0 1 2
INTENT 0 0 3
TX 0 0 0
PSI
CIGUSE INTENT TX
________ ________ ________
CIGUSE 0
INTENT 0 0
TX 0 0 0
STARTING VALUES
TAU
CIGUSE$1 INTENT$1 INTENT$2 INTENT$3
________ ________ ________ ________
1 0.840 0.597 1.030 1.472
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 1.000
TX 0.000 0.000 0.123
TECHNICAL 8 OUTPUT
E STEP ITER LOGLIKELIHOOD ABS CHANGE REL CHANGE ALGORITHM
1 -0.11336167D+04 0.0000000 0.0000000 EM
2 -0.10557277D+04 77.8889909 0.0687084 EM
3 -0.10551628D+04 0.5649177 0.0005351 EM
4 -0.10551624D+04 0.0003701 0.0000004 EM
PLOT INFORMATION
The following plots are available:
Histograms (sample values)
Scatterplots (sample values)
Sample proportions
Bootstrap distributions
DIAGRAM INFORMATION
Use View Diagram under the Diagram menu in the Mplus Editor to view the diagram.
If running Mplus from the Mplus Diagrammer, the diagram opens automatically.
Diagram output
c:\users\gryphon\desktop\chapter8\ex8.26.dgm
Beginning Time: 17:16:58
Ending Time: 17:17:01
Elapsed Time: 00:00:03
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-2015 Muthen & Muthen
Back to examples