Results: Simple Mediation Model - Vale and Maurelli (1983) - Skewness = 2, Kurtosis = 7 - Complete Data - Parametric Bootstrap Confidence Intervals Assuming Multivariate Normal Distribution

results_vm_mod_pb.mvn_ci

Format

A data frame with the following variables

taskid

Simulation task identification number.

n

Sample size.

simreps

Monte Carlo replications.

taudot

Population slope of path from x to y \(\left( \dot{\tau} \right)\)

beta

Population slope of path from m to y \(\left( \beta \right)\)

alpha

Population slope of path from x to m \(\left( \alpha \right)\)

alphabeta

Population indirect effect of x on y through m \(\left( \alpha \beta \right)\)

sigma2x

Population variance of x \(\left( \sigma_{x}^{2} \right)\)

sigma2epsilonm

Population error variance of m \(\left( \sigma_{\varepsilon_{m}}^{2} \right)\)

sigma2epsilony

Population error variance of y \(\left( \sigma_{\varepsilon_{y}}^{2} \right)\)

mux

Population mean of x \(\left( \mu_x \right)\).

deltam

Population intercept of m \(\left( \delta_m \right)\).

deltay

Population intercept of y \(\left( \delta_y \right)\).

est

Mean of the estimate of the indirect effect \(\left( \hat{\alpha} \hat{\beta} \right)\).

se

Mean of the estimate of standard error of the indirect effect \(\left( \hat{\alpha} \hat{\beta} \right)\).

reps

Monte Carlo method of bootstrap replications.

ci_0.05

Mean of the lower limit confidence interval for the 99.9% confidence interval.

ci_0.5

Mean of the lower limit confidence interval for the 99% confidence interval.

ci_2.5

Mean of the lower limit confidence interval for the 95% confidence interval.

ci_97.5

Mean of the upper limit confidence interval for the 95% confidence interval.

ci_99.5

Mean of the upper limit confidence interval for the 99% confidence interval.

ci_99.95

Mean of the upper limit confidence interval for the 99.9% confidence interval.

zero_hit_99.9

Mean zero hit for the 99.9% confidence interval.

zero_hit_99

Mean zero hit for the 99% confidence interval.

zero_hit_95

Mean zero hit for the 95% confidence interval.

len_99.9

Mean confidence interval length for the 99.9% confidence interval.

len_99

Mean confidence interval length for the 99% confidence interval.

len_95

Mean confidence interval length for the 95% confidence interval.

shape_99.9

Mean confidence interval shape for the 99.9% confidence interval.

shape_99

Mean confidence interval shape for the 99% confidence interval.

shape_95

Mean confidence interval shape for the 95% confidence interval.

theta_hit_99.9

Mean theta hit for the 99.9% confidence interval.

theta_hit_99

Mean theta hit for the 99% confidence interval.

theta_hit_95

Mean theta hit for the 95% confidence interval.

theta_miss_99.9

Mean theta miss for the 99.9% confidence interval.

theta_miss_99

Mean theta miss for the 99% confidence interval.

theta_miss_95

Mean theta miss for the 95% confidence interval.

theta

Population parameter \(\alpha \beta\).

power_99.9

Mean power for the 99.9% confidence interval.

power_99

Mean power for the 99% confidence interval.

power_95

Mean power for the 95% confidence interval.

liberal_ll_99.9

Lower limit of the liberal criteria for the 99.9% confidence interval.

liberal_ul_99.9

Upper limit of the liberal criteria for the 99.9% confidence interval.

moderate_ll_99.9

Lower limit of the moderate criteria for the 99.9% confidence interval.

moderate_ul_99.9

Upper limit of the moderate criteria for the 99.9% confidence interval.

strict_ll_99.9

Lower limit of the strict criteria for the 99.9% confidence interval.

strict_ul_99.9

Upper limit of the strict criteria for the 99.9% confidence interval.

liberal_ll_99

Lower limit of the liberal criteria for the 99% confidence interval.

liberal_ul_99

Upper limit of the liberal criteria for the 99% confidence interval.

moderate_ll_99

Lower limit of the moderate criteria for the 99% confidence interval.

moderate_ul_99

Upper limit of the moderate criteria for the 99% confidence interval.

strict_ll_99

Lower limit of the strict criteria for the 99% confidence interval.

strict_ul_99

Upper limit of the strict criteria for the 99% confidence interval.

liberal_ll_95

Lower limit of the liberal criteria for the 95% confidence interval.

liberal_ul_95

Upper limit of the liberal criteria for the 95% confidence interval.

moderate_ll_95

Lower limit of the moderate criteria for the 95% confidence interval.

moderate_ul_95

Upper limit of the moderate criteria for the 95% confidence interval.

strict_ll_95

Lower limit of the strict criteria for the 95% confidence interval.

strict_ul_95

Upper limit of the strict criteria for the 95% confidence interval.

serlin_ll_95

Lower limit of the Serlin criteria for the 95% confidence interval.

serlin_ul_95

Upper limit of the Serlin criteria for the 95% confidence interval.

liberal_99.9

Logical. 1 if miss rate is inside the liberal robustness criteria for 99.9% confidence interval.

liberal_99

Logical. 1 if miss rate is inside the liberal robustness criteria for 99% confidence interval.

liberal_95

Logical. 1 if miss rate is inside the liberal robustness criteria for 95% confidence interval.

moderate_99.9

Logical. 1 if miss rate is inside the moderate robustness criteria for 99.9% confidence interval.

moderate_99

Logical. 1 if miss rate is inside the moderate robustness criteria for 99% confidence interval.

moderate_95

Logical. 1 if miss rate is inside the moderate robustness criteria for 95% confidence interval.

strict_99.9

Logical. 1 if miss rate is inside the strict robustness criteria for 99.9% confidence interval.

strict_99

Logical. 1 if miss rate is inside the strict robustness criteria for 99% confidence interval.

strict_95

Logical. 1 if miss rate is inside the strict robustness criteria for 95% confidence interval.

serlin_95

Logical. 1 if miss rate is inside the Serlin robustness criteria for 95% confidence interval.

missing

Type of missingness.

std

Standardized vs. unstandardize indirect effect.

Method

Method used. Fit in this case.

n_label

Sample size labels.

alpha_label

\(\alpha\) labels.

beta_label

\(\beta\) labels.

taudot_label

\(\dot{\tau}\) labels.

theta_label

\(\theta\) labels.

Details

The simple mediation model is given by $$ y_i = \delta_y + \dot{\tau} x_i + \beta m_i + \varepsilon_{y_{i}} $$

$$ m_i = \delta_m + \alpha x_i + \varepsilon_{m_{i}} $$

The parameters for the mean structure are $$ \boldsymbol{\theta}_{\text{mean structure}} = \left\{ \mu_x, \delta_m, \delta_y \right\} . $$

The parameters for the covariance structure are $$ \boldsymbol{\theta}_{\text{covariance structure}} = \left\{ \dot{\tau}, \beta, \alpha, \sigma_{x}^{2}, \sigma_{\varepsilon_{m}}^{2}, \sigma_{\varepsilon_{y}}^{2} \right\} . $$

See also

Examples

data(results_vm_mod_pb.mvn_ci, package = "jeksterslabRmedsimple") head(results_vm_mod_pb.mvn_ci)
#> taskid n simreps taudot beta alpha alphabeta sigma2x #> 1 1 1000 5000 0.1414214 0.7140742 0.7140742 0.509902 225 #> 2 2 500 5000 0.1414214 0.7140742 0.7140742 0.509902 225 #> 3 3 250 5000 0.1414214 0.7140742 0.7140742 0.509902 225 #> 4 4 200 5000 0.1414214 0.7140742 0.7140742 0.509902 225 #> 5 5 150 5000 0.1414214 0.7140742 0.7140742 0.509902 225 #> 6 6 100 5000 0.1414214 0.7140742 0.7140742 0.509902 225 #> sigma2epsilonm sigma2epsilony mux deltam deltay est se reps #> 1 110.2721 73.3221 100 28.59258 14.45045 0.5100761 0.02434670 5000 #> 2 110.2721 73.3221 100 28.59258 14.45045 0.5118519 0.03465958 5000 #> 3 110.2721 73.3221 100 28.59258 14.45045 0.5120330 0.04924292 5000 #> 4 110.2721 73.3221 100 28.59258 14.45045 0.5138135 0.05550023 5000 #> 5 110.2721 73.3221 100 28.59258 14.45045 0.5142247 0.06427220 5000 #> 6 110.2721 73.3221 100 28.59258 14.45045 0.5175581 0.08005077 5000 #> ci_0.05 ci_0.5 ci_2.5 ci_97.5 ci_99.5 ci_99.95 zero_hit_99.9 #> 1 0.4336865 0.4490725 0.4631972 0.5585365 0.5741798 0.5916634 0e+00 #> 2 0.4047542 0.4259028 0.4455765 0.5813219 0.6040027 0.6295337 0e+00 #> 3 0.3626318 0.3915594 0.4186994 0.6115778 0.6449828 0.6827734 0e+00 #> 4 0.3463307 0.3786731 0.4090317 0.6264247 0.6644933 0.7076204 0e+00 #> 5 0.3219603 0.3588668 0.3935560 0.6452264 0.6900724 0.7409502 0e+00 #> 6 0.2808694 0.3259985 0.3684624 0.6820302 0.7393363 0.8054666 6e-04 #> zero_hit_99 zero_hit_95 len_99.9 len_99 len_95 shape_99.9 shape_99 #> 1 0 0 0.1579769 0.1251073 0.09533931 1.070486 1.051596 #> 2 0 0 0.2247795 0.1780999 0.13574540 1.101275 1.073007 #> 3 0 0 0.3201417 0.2534234 0.19287841 1.145559 1.104569 #> 4 0 0 0.3612896 0.2858202 0.21739300 1.160701 1.116392 #> 5 0 0 0.4189899 0.3312056 0.25167040 1.183254 1.133677 #> 6 0 0 0.5245973 0.4133377 0.31356774 1.222197 1.160656 #> shape_95 theta_hit_99.9 theta_hit_99 theta_hit_95 theta_miss_99.9 #> 1 1.034156 0.9502 0.8778 0.7630 0.0498 #> 2 1.048573 0.9498 0.8812 0.7702 0.0502 #> 3 1.067186 0.9522 0.8842 0.7718 0.0478 #> 4 1.075464 0.9526 0.8876 0.7770 0.0474 #> 5 1.086647 0.9546 0.8980 0.7850 0.0454 #> 6 1.104857 0.9588 0.8994 0.7864 0.0412 #> theta_miss_99 theta_miss_95 theta power_99.9 power_99 power_95 #> 1 0.1222 0.2370 0.509902 1.0000 1 1 #> 2 0.1188 0.2298 0.509902 1.0000 1 1 #> 3 0.1158 0.2282 0.509902 1.0000 1 1 #> 4 0.1124 0.2230 0.509902 1.0000 1 1 #> 5 0.1020 0.2150 0.509902 1.0000 1 1 #> 6 0.1006 0.2136 0.509902 0.9994 1 1 #> liberal_ll_99.9 liberal_ul_99.9 moderate_ll_99.9 moderate_ul_99.9 #> 1 5e-04 0.0015 8e-04 0.0012 #> 2 5e-04 0.0015 8e-04 0.0012 #> 3 5e-04 0.0015 8e-04 0.0012 #> 4 5e-04 0.0015 8e-04 0.0012 #> 5 5e-04 0.0015 8e-04 0.0012 #> 6 5e-04 0.0015 8e-04 0.0012 #> strict_ll_99.9 strict_ul_99.9 liberal_ll_99 liberal_ul_99 moderate_ll_99 #> 1 9e-04 0.0011 0.005 0.015 0.008 #> 2 9e-04 0.0011 0.005 0.015 0.008 #> 3 9e-04 0.0011 0.005 0.015 0.008 #> 4 9e-04 0.0011 0.005 0.015 0.008 #> 5 9e-04 0.0011 0.005 0.015 0.008 #> 6 9e-04 0.0011 0.005 0.015 0.008 #> moderate_ul_99 strict_ll_99 strict_ul_99 liberal_ll_95 liberal_ul_95 #> 1 0.012 0.009 0.011 0.025 0.075 #> 2 0.012 0.009 0.011 0.025 0.075 #> 3 0.012 0.009 0.011 0.025 0.075 #> 4 0.012 0.009 0.011 0.025 0.075 #> 5 0.012 0.009 0.011 0.025 0.075 #> 6 0.012 0.009 0.011 0.025 0.075 #> moderate_ll_95 moderate_ul_95 strict_ll_95 strict_ul_95 serlin_ll_95 #> 1 0.04 0.06 0.045 0.055 0.035 #> 2 0.04 0.06 0.045 0.055 0.035 #> 3 0.04 0.06 0.045 0.055 0.035 #> 4 0.04 0.06 0.045 0.055 0.035 #> 5 0.04 0.06 0.045 0.055 0.035 #> 6 0.04 0.06 0.045 0.055 0.035 #> serlin_ul_95 liberal_99.9 liberal_99 liberal_95 moderate_99.9 moderate_99 #> 1 0.065 0 0 0 0 0 #> 2 0.065 0 0 0 0 0 #> 3 0.065 0 0 0 0 0 #> 4 0.065 0 0 0 0 0 #> 5 0.065 0 0 0 0 0 #> 6 0.065 0 0 0 0 0 #> moderate_95 strict_99.9 strict_99 strict_95 serlin_95 missing std #> 1 0 0 0 0 0 Complete Unstandardized #> 2 0 0 0 0 0 Complete Unstandardized #> 3 0 0 0 0 0 Complete Unstandardized #> 4 0 0 0 0 0 Complete Unstandardized #> 5 0 0 0 0 0 Complete Unstandardized #> 6 0 0 0 0 0 Complete Unstandardized #> Method n_label alpha_label beta_label taudot_label theta_label #> 1 PBPC n: 1000 α: 0.71 β: 0.71 τ̇: 0.14 0.51(0.71,0.71) #> 2 PBPC n: 500 α: 0.71 β: 0.71 τ̇: 0.14 0.51(0.71,0.71) #> 3 PBPC n: 250 α: 0.71 β: 0.71 τ̇: 0.14 0.51(0.71,0.71) #> 4 PBPC n: 200 α: 0.71 β: 0.71 τ̇: 0.14 0.51(0.71,0.71) #> 5 PBPC n: 150 α: 0.71 β: 0.71 τ̇: 0.14 0.51(0.71,0.71) #> 6 PBPC n: 100 α: 0.71 β: 0.71 τ̇: 0.14 0.51(0.71,0.71)
str(results_vm_mod_pb.mvn_ci)
#> 'data.frame': 1062 obs. of 79 variables: #> $ taskid : num 1 2 3 4 5 6 7 8 9 10 ... #> $ n : num 1000 500 250 200 150 100 75 50 20 1000 ... #> $ simreps : num 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 ... #> $ taudot : num 0.141 0.141 0.141 0.141 0.141 ... #> $ beta : num 0.714 0.714 0.714 0.714 0.714 ... #> $ alpha : num 0.714 0.714 0.714 0.714 0.714 ... #> $ alphabeta : num 0.51 0.51 0.51 0.51 0.51 ... #> $ sigma2x : num 225 225 225 225 225 225 225 225 225 225 ... #> $ sigma2epsilonm : num 110 110 110 110 110 ... #> $ sigma2epsilony : num 73.3 73.3 73.3 73.3 73.3 ... #> $ mux : num 100 100 100 100 100 100 100 100 100 100 ... #> $ deltam : num 28.6 28.6 28.6 28.6 28.6 ... #> $ deltay : num 14.5 14.5 14.5 14.5 14.5 ... #> $ est : num 0.51 0.512 0.512 0.514 0.514 ... #> $ se : num 0.0243 0.0347 0.0492 0.0555 0.0643 ... #> $ reps : num 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 ... #> $ ci_0.05 : num 0.434 0.405 0.363 0.346 0.322 ... #> $ ci_0.5 : num 0.449 0.426 0.392 0.379 0.359 ... #> $ ci_2.5 : num 0.463 0.446 0.419 0.409 0.394 ... #> $ ci_97.5 : num 0.559 0.581 0.612 0.626 0.645 ... #> $ ci_99.5 : num 0.574 0.604 0.645 0.664 0.69 ... #> $ ci_99.95 : num 0.592 0.63 0.683 0.708 0.741 ... #> $ zero_hit_99.9 : num 0 0 0 0 0 ... #> $ zero_hit_99 : num 0 0 0 0 0 ... #> $ zero_hit_95 : num 0 0 0 0 0 0 0.0004 0.0056 0.179 0 ... #> $ len_99.9 : num 0.158 0.225 0.32 0.361 0.419 ... #> $ len_99 : num 0.125 0.178 0.253 0.286 0.331 ... #> $ len_95 : num 0.0953 0.1357 0.1929 0.2174 0.2517 ... #> $ shape_99.9 : num 1.07 1.1 1.15 1.16 1.18 ... #> $ shape_99 : num 1.05 1.07 1.1 1.12 1.13 ... #> $ shape_95 : num 1.03 1.05 1.07 1.08 1.09 ... #> $ theta_hit_99.9 : num 0.95 0.95 0.952 0.953 0.955 ... #> $ theta_hit_99 : num 0.878 0.881 0.884 0.888 0.898 ... #> $ theta_hit_95 : num 0.763 0.77 0.772 0.777 0.785 ... #> $ theta_miss_99.9 : num 0.0498 0.0502 0.0478 0.0474 0.0454 0.0412 0.0354 0.024 0.0156 0.0464 ... #> $ theta_miss_99 : num 0.122 0.119 0.116 0.112 0.102 ... #> $ theta_miss_95 : num 0.237 0.23 0.228 0.223 0.215 ... #> $ theta : num 0.51 0.51 0.51 0.51 0.51 ... #> $ power_99.9 : num 1 1 1 1 1 ... #> $ power_99 : num 1 1 1 1 1 ... #> $ power_95 : num 1 1 1 1 1 ... #> $ liberal_ll_99.9 : num 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 ... #> $ liberal_ul_99.9 : num 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 ... #> $ moderate_ll_99.9: num 8e-04 8e-04 8e-04 8e-04 8e-04 8e-04 8e-04 8e-04 8e-04 8e-04 ... #> $ moderate_ul_99.9: num 0.0012 0.0012 0.0012 0.0012 0.0012 0.0012 0.0012 0.0012 0.0012 0.0012 ... #> $ strict_ll_99.9 : num 9e-04 9e-04 9e-04 9e-04 9e-04 9e-04 9e-04 9e-04 9e-04 9e-04 ... #> $ strict_ul_99.9 : num 0.0011 0.0011 0.0011 0.0011 0.0011 0.0011 0.0011 0.0011 0.0011 0.0011 ... #> $ liberal_ll_99 : num 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 ... #> $ liberal_ul_99 : num 0.015 0.015 0.015 0.015 0.015 0.015 0.015 0.015 0.015 0.015 ... #> $ moderate_ll_99 : num 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 ... #> $ moderate_ul_99 : num 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 ... #> $ strict_ll_99 : num 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 ... #> $ strict_ul_99 : num 0.011 0.011 0.011 0.011 0.011 0.011 0.011 0.011 0.011 0.011 ... #> $ liberal_ll_95 : num 0.025 0.025 0.025 0.025 0.025 0.025 0.025 0.025 0.025 0.025 ... #> $ liberal_ul_95 : num 0.075 0.075 0.075 0.075 0.075 0.075 0.075 0.075 0.075 0.075 ... #> $ moderate_ll_95 : num 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 ... #> $ moderate_ul_95 : num 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 ... #> $ strict_ll_95 : num 0.045 0.045 0.045 0.045 0.045 0.045 0.045 0.045 0.045 0.045 ... #> $ strict_ul_95 : num 0.055 0.055 0.055 0.055 0.055 0.055 0.055 0.055 0.055 0.055 ... #> $ serlin_ll_95 : num 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 ... #> $ serlin_ul_95 : num 0.065 0.065 0.065 0.065 0.065 0.065 0.065 0.065 0.065 0.065 ... #> $ liberal_99.9 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ liberal_99 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ liberal_95 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ moderate_99.9 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ moderate_99 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ moderate_95 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ strict_99.9 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ strict_99 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ strict_95 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ serlin_95 : num 0 0 0 0 0 0 0 0 0 0 ... #> $ missing : chr "Complete" "Complete" "Complete" "Complete" ... #> $ std : chr "Unstandardized" "Unstandardized" "Unstandardized" "Unstandardized" ... #> $ Method : chr "PBPC" "PBPC" "PBPC" "PBPC" ... #> $ n_label : Factor w/ 9 levels "n: 20","n: 50",..: 9 8 7 6 5 4 3 2 1 9 ... #> $ alpha_label : Factor w/ 4 levels "α: 0.00","α: 0.38",..: 4 4 4 4 4 4 4 4 4 4 ... #> $ beta_label : Factor w/ 4 levels "β: 0.00","β: 0.38",..: 4 4 4 4 4 4 4 4 4 4 ... #> $ taudot_label : Factor w/ 4 levels "τ̇: 0.00","τ̇: 0.14",..: 2 2 2 2 2 2 2 2 2 1 ... #> $ theta_label : chr "0.51(0.71,0.71)" "0.51(0.71,0.71)" "0.51(0.71,0.71)" "0.51(0.71,0.71)" ...