Correlation

comrbidities

Make TBI Status (In comorbidity section we have 188, while 209 reported TBI )

Make (QOLIE10)

Make (VR-12)

Making scales

## 
##   1   2   4 
## 250 250  11

Join Scales to LEEP

# step1: lm(Y~X)

	Model1-QOLIE10			Model2-QOLIE10			Model3-QOLIE10			Model4-QOLIE10			Model5-QOLIE10			Model6-QOLIE10
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	26.51	24.69 – 28.32	<0.001	31.17	30.04 – 32.29	<0.001	25.00	14.03 – 35.97	<0.001	17.75	15.40 – 20.10	<0.001	19.10	16.65 – 21.55	<0.001	16.61	14.14 – 19.08	<0.001
Psychiatric Comorbidity	4.73	2.59 – 6.86	<0.001				6.23	-4.79 – 17.26	0.267
Somatic Comorbidity				-4.62	-6.78 – -2.46	<0.001	1.55	-9.57 – 12.67	0.784
PHQ										2.64	2.16 – 3.12	<0.001				1.97	1.28 – 2.66	<0.001
GAD2													2.17	1.71 – 2.63	<0.001	0.85	0.21 – 1.48	0.009
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.068 / 0.065			0.064 / 0.061			0.069 / 0.061			0.314 / 0.312			0.250 / 0.247			0.332 / 0.327

	Model1-VR12			Model2-VR12			Model3-VR12			Model4-VR12			Model5-VR12			Model6-VR12
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	43.95	42.21 – 45.68	<0.001	41.38	40.30 – 42.47	<0.001	54.00	43.54 – 64.46	<0.001	55.01	52.97 – 57.06	<0.001	53.56	51.38 – 55.74	<0.001	56.23	54.09 – 58.36	<0.001
Psychiatric Comorbidity	-2.70	-4.75 – -0.65	0.010				-12.75	-23.26 – -2.24	0.018
Somatic Comorbidity				2.28	0.20 – 4.35	0.032	-10.34	-20.94 – 0.26	0.056
PHQ										-2.82	-3.24 – -2.41	<0.001				-2.11	-2.70 – -1.51	<0.001
GAD2													-2.32	-2.73 – -1.91	<0.001	-0.90	-1.45 – -0.36	0.001
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.025 / 0.021			0.018 / 0.014			0.039 / 0.032			0.408 / 0.406			0.325 / 0.322			0.432 / 0.427

	Model1-MCS			Model2-MCS			Model3-MCS			Model4-MCS			Model5-MCS			Model6-MCS
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	17.96	16.91 – 19.01	<0.001	16.73	16.08 – 17.38	<0.001	23.00	16.69 – 29.31	<0.001	23.82	22.49 – 25.14	<0.001	22.79	21.39 – 24.19	<0.001	24.30	22.89 – 25.70	<0.001
Psychiatric Comorbidity	-1.29	-2.53 – -0.06	0.040				-6.34	-12.68 – 0.01	0.050
Somatic Comorbidity				1.09	-0.16 – 2.33	0.088	-5.18	-11.59 – 1.22	0.112
PHQ										-1.47	-1.74 – -1.20	<0.001				-1.19	-1.58 – -0.80	<0.001
GAD2													-1.16	-1.42 – -0.89	<0.001	-0.36	-0.72 – 0.00	0.051
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.016 / 0.012			0.011 / 0.007			0.026 / 0.018			0.309 / 0.307			0.225 / 0.222			0.319 / 0.314

	Model1-PCS			Model2-PCS			Model3-PCS			Model4-PCS			Model5-PCS			Model6-PCS
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	25.99	25.10 – 26.87	<0.001	24.65	24.10 – 25.20	<0.001	31.00	25.70 – 36.30	<0.001	31.20	30.12 – 32.27	<0.001	30.77	29.66 – 31.88	<0.001	31.93	30.81 – 33.05	<0.001
Psychiatric Comorbidity	-1.40	-2.44 – -0.36	0.008				-6.41	-11.74 – -1.09	0.018
Somatic Comorbidity				1.19	0.14 – 2.24	0.026	-5.15	-10.53 – 0.22	0.060
PHQ										-1.35	-1.57 – -1.13	<0.001				-0.92	-1.23 – -0.61	<0.001
GAD2													-1.16	-1.37 – -0.95	<0.001	-0.55	-0.83 – -0.26	<0.001
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.027 / 0.023			0.019 / 0.015			0.040 / 0.032			0.363 / 0.361			0.317 / 0.314			0.396 / 0.391

Some descriptive of Illness Perception

## 
## 1\nNot at all severe   2\nA little severe   3\nSomewhat severe 
##                   39                   65                   52 
## 4\nModerately severe      5\nQuite severe       6\nVery severe 
##                   54                   29                   16 
##  7\nExtremely severe 
##                    6

## 
##  1  2  3  4  5  6  7 
## 39 65 52 54 29 16  6

step2: lm(M~X)

	Model1-Severity			Model2-Severity			Model3-Severity			Model4-severity			Model5-Severity			Model6-Severity
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	2.93	2.58 – 3.29	<0.001	3.23	3.01 – 3.45	<0.001	2.00	-0.15 – 4.15	0.069	1.62	1.12 – 2.13	<0.001	1.86	1.35 – 2.37	<0.001	1.52	0.98 – 2.05	<0.001
Psychiatric Comorbidity	0.31	-0.11 – 0.73	0.143				1.24	-0.92 – 3.41	0.259
Somatic Comorbidity				-0.27	-0.70 – 0.15	0.205	0.96	-1.23 – 3.14	0.389
PHQ										0.33	0.23 – 0.43	<0.001				0.27	0.12 – 0.42	<0.001
GAD2													0.26	0.16 – 0.36	<0.001	0.08	-0.06 – 0.22	0.256
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.008 / 0.004			0.006 / 0.002			0.011 / 0.003			0.137 / 0.134			0.099 / 0.096			0.142 / 0.135

step3: lm(Y ~X+M)

fit1 <- lm(SUM_QOLIE10~psy_com+severe, data=LEEP)
#summary(fit1) # significant
fit2 <- lm(SUM_QOLIE10~somatic_com+severe, data=LEEP)
#summary(fit1) # significant
fit3 <- lm(SUM_QOLIE10~psy_com+somatic_com+severe, data=LEEP)
#summary(fit1) # significant
fit4 <- lm(SUM_QOLIE10~PHQ+severe, data=LEEP)
#summary(fit1) # significant
fit5 <- lm(SUM_QOLIE10~GAD2+severe, data=LEEP)
#summary(fit1) # significant
fit6<- lm(SUM_QOLIE10~PHQ+GAD2+severe, data=LEEP)
#summary(fit1) # significant

pl <- c(
  `(Intercept)` = "Intercept",psy_com = "Psychiatric Comorbidity ",somatic_com  ="Somatic Comorbidity")

tab_model(
  fit1,fit2,fit3,fit4, fit5,fit6,pred.labels = pl,
    dv.labels = c("Model1-QOLIE10Severity", "Model2-QOLIE10Severity","Model3-QOLIE10Severity", "Model4-QOLIE10Severity", "Model5-QOLIE10Severity","Model6-QOLIE10Severity" ),
  string.pred = "Coeffcient",
  string.ci = "Conf. Int (95%)",
  string.p = "P-Value"
)

	Model1-QOLIE10Severity			Model2-QOLIE10Severity			Model3-QOLIE10Severity			Model4-QOLIE10Severity			Model5-QOLIE10Severity			Model6-QOLIE10Severity
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	18.32	16.17 – 20.48	<0.001	22.10	20.17 – 24.03	<0.001	19.41	10.16 – 28.66	<0.001	14.35	12.09 – 16.62	<0.001	14.83	12.48 – 17.19	<0.001	13.50	11.15 – 15.86	<0.001
Psychiatric Comorbidity	3.85	2.06 – 5.65	<0.001				2.76	-6.50 – 12.02	0.558
severe	2.79	2.27 – 3.31	<0.001	2.81	2.28 – 3.33	<0.001	2.79	2.27 – 3.32	<0.001	2.09	1.58 – 2.60	<0.001	2.29	1.78 – 2.81	<0.001	2.05	1.54 – 2.56	<0.001
Somatic Comorbidity				-3.85	-5.66 – -2.04	<0.001	-1.13	-10.46 – 8.20	0.812
PHQ										1.94	1.48 – 2.40	<0.001				1.42	0.78 – 2.05	<0.001
GAD2													1.57	1.15 – 2.00	<0.001	0.69	0.12 – 1.26	0.019
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.349 / 0.344			0.348 / 0.343			0.349 / 0.342			0.451 / 0.447			0.422 / 0.418			0.463 / 0.457

	Model1-MCS-Severity			Model2-MCS-Severity			Model3-MCS-Severity			Model4-MCS-Severity			Model5-MCS-Severity			Model6-MCS-Severity
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	22.21	20.91 – 23.50	<0.001	21.44	20.27 – 22.60	<0.001	25.88	20.33 – 31.42	<0.001	25.44	24.12 – 26.76	<0.001	24.92	23.53 – 26.32	<0.001	25.79	24.41 – 27.17	<0.001
Psychiatric Comorbidity	-0.84	-1.92 – 0.24	0.127				-4.55	-10.10 – 1.01	0.108
severe	-1.45	-1.76 – -1.14	<0.001	-1.46	-1.77 – -1.14	<0.001	-1.44	-1.75 – -1.12	<0.001	-1.00	-1.30 – -0.70	<0.001	-1.14	-1.45 – -0.84	<0.001	-0.98	-1.28 – -0.68	<0.001
Somatic Comorbidity				0.69	-0.40 – 1.78	0.216	-3.81	-9.40 – 1.79	0.182
PHQ										-1.14	-1.41 – -0.87	<0.001				-0.92	-1.30 – -0.55	<0.001
GAD2													-0.86	-1.11 – -0.61	<0.001	-0.28	-0.62 – 0.05	0.101
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.255 / 0.249			0.253 / 0.247			0.260 / 0.252			0.408 / 0.404			0.360 / 0.355			0.415 / 0.408

	Model1-PCS-Severity			Model2-PCS-Severity			Model3-PCS-Severity			Model4-PCS-Severity			Model5-PCS-Severity			Model6-PCS-Severity
Coeffcient	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value	Estimates	Conf. Int (95%)	P-Value
Intercept	28.55	27.38 – 29.72	<0.001	27.51	26.45 – 28.56	<0.001	32.72	27.72 – 37.73	<0.001	31.85	30.71 – 32.99	<0.001	31.72	30.54 – 32.91	<0.001	32.49	31.32 – 33.66	<0.001
Psychiatric Comorbidity	-1.13	-2.10 – -0.15	0.024				-5.34	-10.36 – -0.33	0.037
severe	-0.88	-1.16 – -0.59	<0.001	-0.88	-1.17 – -0.60	<0.001	-0.86	-1.15 – -0.58	<0.001	-0.40	-0.66 – -0.14	0.002	-0.51	-0.77 – -0.25	<0.001	-0.37	-0.62 – -0.12	0.004
Somatic Comorbidity				0.95	-0.04 – 1.94	0.059	-4.33	-9.38 – 0.72	0.093
PHQ										-1.22	-1.45 – -0.98	<0.001				-0.82	-1.13 – -0.50	<0.001
GAD2													-1.03	-1.24 – -0.81	<0.001	-0.52	-0.80 – -0.23	<0.001
Observations	261			261			261			261			261			261
R2 / R2 adjusted	0.148 / 0.142			0.143 / 0.137			0.158 / 0.148			0.385 / 0.381			0.355 / 0.350			0.415 / 0.408

CI for a*b:

## Warning: package 'RMediation' was built under R version 4.0.4

## Loading required package: MASS

## 
## Attaching package: 'MASS'

## The following object is masked from 'package:dplyr':
## 
##     select

## Loading required package: lavaan

## Warning: package 'lavaan' was built under R version 4.0.4

## This is lavaan 0.6-7

## lavaan is BETA software! Please report any bugs.

## Loading required package: e1071

## 
## Attaching package: 'e1071'

## The following object is masked from 'package:gtools':
## 
##     permutations

## Warning in arrows(MedCI[[1]][1], yci, MedCI[[1]][2], yci, length = smidge, :
## 'length', 'angle', or 'code' greater than length 1; values after the first are
## ignored

## $`97.5% CI`
## [1] 0.4129502 0.9126848
## 
## $Estimate
## [1] 0.6463027
## 
## $SE
## [1] 0.1277206

As shown in the resulted figure, the mean of ab is 0.64, and standard error 0.128. The corresponding 95% CI is (0.41,0.93).