If we had previously saved out CBPS weights, we could read them in instead of re-creating them.
@@ -895,7 +908,7 @@ Next, we look inside the output to summarize the weighting process (drawing on the summary method from weightIt, for any given imputed dataset.
This summary also provides the effective sample size, or the sample size that results from applying the weights to the original sample size, for each time point. Weighting can often result is an effective or weighted sample size that is smaller than the orignal sample size and is something to keep in mind when evaluating weighting methods. In this example, we see that our original 1,292 sample is reduced to 605 upon weighting with the CBPS method.
We then view a distribution plot of the weights for any one imputed dataset. The user has the option to supply save.out
to save plots to the home directory.
@@ -1191,22 +1204,149 @@
-# method <- "gbm"
-#
-# weights.gbm <- createWeights(obj = obj,
-# data = data,
-# method = method,
-# formulas = formulas,
-# save.out = save.out)
-#
-# print(weights.gbm,
-# i = 1)
-#
-# summary(weights.gbm[[1]])
-#
-# plot(weights.gbm,
-# i = 1,
-# save.out = save.out)
+
method <- "gbm"
+
+weights.gbm <- createWeights(obj = obj,
+ data = data,
+ method = method,
+ formulas = formulas,
+ save.out = save.out)
+
+print(weights.gbm,
+ i = 1)
+#>
+#> For imputation 1 and the `gbm` weighting method, the median weight value is 1.26 (SD = 8.18; range = 0-216).
+
+summary(weights.gbm[[1]])
+#> Time 1
+#> Summary of weights
+#>
+#> - Weight ranges:
+#>
+#> Min Max
+#> all 0 |---------------------------| 216.1534
+#>
+#> - Units with the 5 most extreme weights:
+#>
+#> 342 30 475 1145 36
+#> all 43.0076 43.1447 55.5631 150.016 216.1534
+#>
+#> - Weight statistics:
+#>
+#> Coef of Var MAD Entropy # Zeros
+#> all 3.383 0.916 1.016 0
+#>
+#> - Effective Sample Sizes:
+#>
+#> Total
+#> Unweighted 1292.
+#> Weighted 103.87
+#>
+#> Time 2
+#> Summary of weights
+#>
+#> - Weight ranges:
+#>
+#> Min Max
+#> all 0 |---------------------------| 216.1534
+#>
+#> - Units with the 5 most extreme weights:
+#>
+#> 342 30 475 1145 36
+#> all 43.0076 43.1447 55.5631 150.016 216.1534
+#>
+#> - Weight statistics:
+#>
+#> Coef of Var MAD Entropy # Zeros
+#> all 3.383 0.916 1.016 0
+#>
+#> - Effective Sample Sizes:
+#>
+#> Total
+#> Unweighted 1292.
+#> Weighted 103.87
+#>
+#> Time 3
+#> Summary of weights
+#>
+#> - Weight ranges:
+#>
+#> Min Max
+#> all 0 |---------------------------| 216.1534
+#>
+#> - Units with the 5 most extreme weights:
+#>
+#> 342 30 475 1145 36
+#> all 43.0076 43.1447 55.5631 150.016 216.1534
+#>
+#> - Weight statistics:
+#>
+#> Coef of Var MAD Entropy # Zeros
+#> all 3.383 0.916 1.016 0
+#>
+#> - Effective Sample Sizes:
+#>
+#> Total
+#> Unweighted 1292.
+#> Weighted 103.87
+#>
+#> Time 4
+#> Summary of weights
+#>
+#> - Weight ranges:
+#>
+#> Min Max
+#> all 0 |---------------------------| 216.1534
+#>
+#> - Units with the 5 most extreme weights:
+#>
+#> 342 30 475 1145 36
+#> all 43.0076 43.1447 55.5631 150.016 216.1534
+#>
+#> - Weight statistics:
+#>
+#> Coef of Var MAD Entropy # Zeros
+#> all 3.383 0.916 1.016 0
+#>
+#> - Effective Sample Sizes:
+#>
+#> Total
+#> Unweighted 1292.
+#> Weighted 103.87
+#>
+#> Time 5
+#> Summary of weights
+#>
+#> - Weight ranges:
+#>
+#> Min Max
+#> all 0 |---------------------------| 216.1534
+#>
+#> - Units with the 5 most extreme weights:
+#>
+#> 342 30 475 1145 36
+#> all 43.0076 43.1447 55.5631 150.016 216.1534
+#>
+#> - Weight statistics:
+#>
+#> Coef of Var MAD Entropy # Zeros
+#> all 3.383 0.916 1.016 0
+#>
+#> - Effective Sample Sizes:
+#>
+#> Total
+#> Unweighted 1292.
+#> Weighted 103.87
+
+plot(weights.gbm,
+ i = 1,
+ save.out = save.out)
+The GBM method produces a similar mean as GLM and a similarly large range (0-216).
summary(balance_stats.cbps,
save.out = save.out)
-#> Averaging across imputed datasets using cbps weighting method.
-#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 0 out of 241 (0.0%) covariates across time points remain imbalanced:
+#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 2 out of 241 (0.8%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.07):
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
@@ -1741,7 +1884,7 @@ Workflow: Continuous Exposure
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 0 |
+#> | ESETA1.35 | 59 | 2 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
@@ -1751,7 +1894,8 @@ Workflow: Continuous Exposure
summary(balance_stats.cbps,
i = 1,
save.out = save.out)
-#> USER ALERT: For imputation 1 using `cbps` weighting method: As shown below, 4 out of 241 (1.7%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.10):
+#> No covariates remain imbalaned for imputation 1 using `cbps` weighting method.
+#> USER ALERT: For imputation 1 using `cbps` weighting method: As shown below, 0 out of 241 (0.0%) covariates across time points remain imbalanced:
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
@@ -1761,9 +1905,9 @@ Workflow: Continuous Exposure
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 2 |
+#> | ESETA1.35 | 59 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.58 | 69 | 2 |
+#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
#>
#> Table: Imbalanced Covariates for imputation 1 using `cbps`
@@ -1805,123 +1949,123 @@ Bart
@@ -1994,8 +2138,7 @@ Workflow: Continuous Exposure
summary(balance_stats.cbps,
save.out = save.out)
-#> Averaging across imputed datasets using cbps weighting method.
-#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 0 out of 241 (0.0%) covariates across time points remain imbalanced:
+#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 2 out of 241 (0.8%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.07):
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
@@ -2005,7 +2148,7 @@ Workflow: Continuous Exposure
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 0 |
+#> | ESETA1.35 | 59 | 2 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
@@ -2031,24 +2174,40 @@ Workflow: Continuous Exposure
#>
#> Table: Imbalanced Covariates Averaging Across Imputed Datasetrs using `glm`
-# summary(balance_stats.gbm,
-# save.out = save.out)
+summary(balance_stats.gbm,
+ save.out = save.out)
+#> USER ALERT: Averaging across imputed datasets using `gbm` weighting method: As shown below, 27 out of 241 (11.2%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.12 (max: 0.22):
+#> +-----------+-----------------------+----------------------------+
+#> | Exposure | Total # of covariates | # of imbalanced covariates |
+#> +===========+=======================+============================+
+#> | ESETA1.6 | 28 | 7 |
+#> +-----------+-----------------------+----------------------------+
+#> | ESETA1.15 | 38 | 3 |
+#> +-----------+-----------------------+----------------------------+
+#> | ESETA1.24 | 47 | 10 |
+#> +-----------+-----------------------+----------------------------+
+#> | ESETA1.35 | 59 | 5 |
+#> +-----------+-----------------------+----------------------------+
+#> | ESETA1.58 | 69 | 2 |
+#> +-----------+-----------------------+----------------------------+
+#>
+#> Table: Imbalanced Covariates Averaging Across Imputed Datasetrs using `gbm`
summary(balance_stats.bart,
save.out = save.out)
-#> USER ALERT: Averaging across imputed datasets using `bart` weighting method: As shown below, 16 out of 241 (6.6%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.13 (max: 0.19):
+#> USER ALERT: Averaging across imputed datasets using `bart` weighting method: As shown below, 28 out of 241 (11.6%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.13 (max: 0.24):
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
-#> | ESETA1.6 | 28 | 3 |
+#> | ESETA1.6 | 28 | 7 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.15 | 38 | 9 |
+#> | ESETA1.15 | 38 | 13 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.24 | 47 | 4 |
+#> | ESETA1.24 | 47 | 7 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 0 |
+#> | ESETA1.35 | 59 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.58 | 69 | 0 |
+#> | ESETA1.58 | 69 | 1 |
#> +-----------+-----------------------+----------------------------+
#>
#> Table: Imbalanced Covariates Averaging Across Imputed Datasetrs using `bart`
@@ -2080,7 +2239,7 @@ Workflow: Continuous Exposure
#> $effective.sample.size
#> Total
#> Unweighted 1292.0000
-#> Weighted 567.0924
+#> Weighted 605.0426
summary(weights.glm[[1]])[[1]][6]
#> $effective.sample.size
@@ -2088,19 +2247,23 @@ Workflow: Continuous Exposure
#> Unweighted 1292.00000
#> Weighted 82.74434
-# summary(weights.gbm[[1]])[[1]][6]
+summary(weights.gbm[[1]])[[1]][6]
+#> $effective.sample.size
+#> Total
+#> Unweighted 1292.0000
+#> Weighted 103.8726
summary(weights.bart[[1]])[[1]][6]
#> $effective.sample.size
#> Total
#> Unweighted 1292.00000
-#> Weighted 36.66269
+#> Weighted 27.65344
summary(weights.super[[1]])[[1]][6]
#> $effective.sample.size
-#> Total
-#> Unweighted 1292.0000
-#> Weighted 83.0753
+#> Total
+#> Unweighted 1292.00000
+#> Weighted 83.05617
From this, we also find that the CBPS method yields the highest effective sample size of 605.
From these inspections, we identify the best performing weighting method as CBPS.
@@ -2119,487 +2282,487 @@ Workflow: Continuous Exposure
#> +-----------+----------------+---------------+------------+----------+
#> | exposure | covariate | std_bal_stats | bal_thresh | balanced |
#> +===========+================+===============+============+==========+
-#> | ESETA1.6 | state_1 | -0.015841 | 0.1 | 1 |
+#> | ESETA1.6 | state_1 | 0.0056014 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | BioDadInHH2_1 | 0.049787 | 0.1 | 1 |
+#> | ESETA1.6 | BioDadInHH2_1 | -0.0010128 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | PmAge2 | 0.021011 | 0.1 | 1 |
+#> | ESETA1.6 | PmAge2 | -0.0014006 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | PmBlac2_1 | -0.033318 | 0.1 | 1 |
+#> | ESETA1.6 | PmBlac2_1 | 0.0056993 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | TcBlac2_1 | -0.017018 | 0.1 | 1 |
+#> | ESETA1.6 | TcBlac2_1 | 0.0057264 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | PmMrSt2_1 | 0.026148 | 0.1 | 1 |
+#> | ESETA1.6 | PmMrSt2_1 | -0.003174 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | PmEd2 | 0.001702 | 0.05 | 1 |
+#> | ESETA1.6 | PmEd2 | -0.0126755 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | KFASTScr | -0.016403 | 0.1 | 1 |
+#> | ESETA1.6 | KFASTScr | -0.0123196 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | RMomAgeU | 0.072699 | 0.1 | 1 |
+#> | ESETA1.6 | RMomAgeU | -0.0050029 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | RHealth_1 | -0.013994 | 0.1 | 1 |
+#> | ESETA1.6 | RHealth_1 | -0.0058639 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | RHealth_2 | -0.020431 | 0.1 | 1 |
+#> | ESETA1.6 | RHealth_2 | -0.000575 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | RHealth_3 | 0.02928 | 0.1 | 1 |
+#> | ESETA1.6 | RHealth_3 | 0.00236 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | RHealth_4 | 0.01127 | 0.1 | 1 |
+#> | ESETA1.6 | RHealth_4 | 0.0062365 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | RHealth_5 | 0.012702 | 0.1 | 1 |
+#> | ESETA1.6 | RHealth_5 | 0.0165904 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | HomeOwnd_1 | 0.003785 | 0.1 | 1 |
+#> | ESETA1.6 | HomeOwnd_1 | -0.0015327 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | HomeOwnd_2 | 0.00374 | 0.1 | 1 |
+#> | ESETA1.6 | HomeOwnd_2 | -0.0004367 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | HomeOwnd_3 | -0.005058 | 0.1 | 1 |
+#> | ESETA1.6 | HomeOwnd_3 | 0.0036865 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | HomeOwnd_4 | -0.002294 | 0.1 | 1 |
+#> | ESETA1.6 | HomeOwnd_4 | 0.0058464 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | SWghtLB | -0.030219 | 0.1 | 1 |
+#> | ESETA1.6 | SWghtLB | -0.0076889 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | SurpPreg_1 | -0.004567 | 0.1 | 1 |
+#> | ESETA1.6 | SurpPreg_1 | 0.0027974 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | SmokTotl_1 | 0.02189 | 0.1 | 1 |
+#> | ESETA1.6 | SmokTotl_1 | 0.0041922 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | SmokTotl_2 | -0.027643 | 0.1 | 1 |
+#> | ESETA1.6 | SmokTotl_2 | -0.0021733 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | SmokTotl_3 | 0.012046 | 0.1 | 1 |
+#> | ESETA1.6 | SmokTotl_3 | 0.0127002 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | DrnkFreq_1 | -0.026237 | 0.1 | 1 |
+#> | ESETA1.6 | DrnkFreq_1 | -0.0056479 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | DrnkFreq_2 | 0.013558 | 0.1 | 1 |
+#> | ESETA1.6 | DrnkFreq_2 | 0.0062915 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | peri_health | -0.035748 | 0.1 | 1 |
+#> | ESETA1.6 | peri_health | 0.0018227 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | caregiv_health | -0.003855 | 0.1 | 1 |
+#> | ESETA1.6 | caregiv_health | 0.0063088 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.6 | gov_assist | -0.022816 | 0.1 | 1 |
+#> | ESETA1.6 | gov_assist | 0.0142948 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | state_1 | 0.036174 | 0.1 | 1 |
+#> | ESETA1.15 | state_1 | 0.0093982 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | BioDadInHH2_1 | -0.012204 | 0.1 | 1 |
+#> | ESETA1.15 | BioDadInHH2_1 | -0.0060919 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmAge2 | 0.004854 | 0.1 | 1 |
+#> | ESETA1.15 | PmAge2 | -0.00208 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmBlac2_1 | 0.042683 | 0.1 | 1 |
+#> | ESETA1.15 | PmBlac2_1 | 0.0078345 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | TcBlac2_1 | 0.025346 | 0.1 | 1 |
+#> | ESETA1.15 | TcBlac2_1 | 0.004728 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmMrSt2_1 | 0.004899 | 0.1 | 1 |
+#> | ESETA1.15 | PmMrSt2_1 | -0.003418 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmEd2 | -0.029365 | 0.05 | 1 |
+#> | ESETA1.15 | PmEd2 | -0.0143257 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | KFASTScr | -0.037699 | 0.1 | 1 |
+#> | ESETA1.15 | KFASTScr | -0.0059526 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RMomAgeU | -0.041403 | 0.1 | 1 |
+#> | ESETA1.15 | RMomAgeU | -0.0017348 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_1 | -0.035278 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_1 | -0.0033083 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_2 | 0.001438 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_2 | -0.0021593 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_3 | -0.011183 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_3 | 0.0017789 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_4 | 0.010747 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_4 | 0.0109611 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_5 | -0.049295 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_5 | -0.0110053 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_1 | -0.033142 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_1 | -0.007286 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_2 | 0.044428 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_2 | -0.0015007 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_3 | -0.016749 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_3 | 0.0027981 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_4 | 0.000614 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_4 | 0.00847 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SWghtLB | 0.020957 | 0.1 | 1 |
+#> | ESETA1.15 | SWghtLB | -0.0013944 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SurpPreg_1 | 0.027651 | 0.1 | 1 |
+#> | ESETA1.15 | SurpPreg_1 | 0.0011906 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_1 | -0.050611 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_1 | 0.0014261 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_2 | 0.069412 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_2 | 0.0056452 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_3 | -0.026007 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_3 | -0.0058635 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | DrnkFreq_1 | -0.011565 | 0.1 | 1 |
+#> | ESETA1.15 | DrnkFreq_1 | 0.0061577 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | DrnkFreq_2 | 0.003636 | 0.1 | 1 |
+#> | ESETA1.15 | DrnkFreq_2 | -0.0064732 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | peri_health | 0.011326 | 0.1 | 1 |
+#> | ESETA1.15 | peri_health | -0.001696 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | caregiv_health | 0.007437 | 0.1 | 1 |
+#> | ESETA1.15 | caregiv_health | 0.0067301 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | gov_assist | -0.015918 | 0.1 | 1 |
+#> | ESETA1.15 | gov_assist | 0.0086264 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SAAmylase.6 | -0.038498 | 0.1 | 1 |
+#> | ESETA1.15 | SAAmylase.6 | 0.0031339 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | MDI.6 | 0.045319 | 0.1 | 1 |
+#> | ESETA1.15 | MDI.6 | -0.0008854 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHasSO.6_1 | -0.02082 | 0.1 | 1 |
+#> | ESETA1.15 | RHasSO.6_1 | 0.0072932 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | WndNbrhood.6 | 0.004196 | 0.1 | 1 |
+#> | ESETA1.15 | WndNbrhood.6 | -0.0053101 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | IBRAttn.6 | -0.019396 | 0.1 | 1 |
+#> | ESETA1.15 | IBRAttn.6 | -0.0030064 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | B18Raw.6 | 0.033228 | 0.1 | 1 |
+#> | ESETA1.15 | B18Raw.6 | 0.020191 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HOMEETA1.6 | -0.008945 | 0.1 | 1 |
+#> | ESETA1.15 | HOMEETA1.6 | -0.0133514 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | InRatioCor.6 | -0.045347 | 0.05 | 1 |
+#> | ESETA1.15 | InRatioCor.6 | -0.0212558 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | CORTB.6 | 0.010691 | 0.1 | 1 |
+#> | ESETA1.15 | CORTB.6 | -0.0004624 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | ESETA1.6 | 0.043475 | 0.1 | 1 |
+#> | ESETA1.15 | ESETA1.6 | 0.0190618 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | state_1 | 0.006241 | 0.1 | 1 |
+#> | ESETA1.24 | state_1 | 0.0178116 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | BioDadInHH2_1 | -0.035814 | 0.1 | 1 |
+#> | ESETA1.24 | BioDadInHH2_1 | -0.0113743 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | PmAge2 | -0.017532 | 0.1 | 1 |
+#> | ESETA1.24 | PmAge2 | -0.0025439 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | PmBlac2_1 | 0.015892 | 0.1 | 1 |
+#> | ESETA1.24 | PmBlac2_1 | 0.0137147 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | TcBlac2_1 | 0.030768 | 0.1 | 1 |
+#> | ESETA1.24 | TcBlac2_1 | 0.0194831 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | PmMrSt2_1 | 0.062369 | 0.1 | 1 |
+#> | ESETA1.24 | PmMrSt2_1 | 0.0056732 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | PmEd2 | 0.00653 | 0.05 | 1 |
+#> | ESETA1.24 | PmEd2 | -0.0088749 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | KFASTScr | 0.032766 | 0.1 | 1 |
+#> | ESETA1.24 | KFASTScr | 0.0027367 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RMomAgeU | 0.032742 | 0.1 | 1 |
+#> | ESETA1.24 | RMomAgeU | -0.0066146 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RHealth_1 | 0.00826 | 0.1 | 1 |
+#> | ESETA1.24 | RHealth_1 | -0.0024297 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RHealth_2 | 0.025263 | 0.1 | 1 |
+#> | ESETA1.24 | RHealth_2 | 0.0132047 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RHealth_3 | -0.029311 | 0.1 | 1 |
+#> | ESETA1.24 | RHealth_3 | -0.0131262 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RHealth_4 | -0.020001 | 0.1 | 1 |
+#> | ESETA1.24 | RHealth_4 | -0.0077858 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RHealth_5 | 0.057299 | 0.1 | 1 |
+#> | ESETA1.24 | RHealth_5 | 0.0232206 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | HomeOwnd_1 | -0.032929 | 0.1 | 1 |
+#> | ESETA1.24 | HomeOwnd_1 | -0.0102612 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | HomeOwnd_2 | 0.008234 | 0.1 | 1 |
+#> | ESETA1.24 | HomeOwnd_2 | -0.0164388 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | HomeOwnd_3 | 0.000591 | 0.1 | 1 |
+#> | ESETA1.24 | HomeOwnd_3 | 0.0056505 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | HomeOwnd_4 | 0.022929 | 0.1 | 1 |
+#> | ESETA1.24 | HomeOwnd_4 | 0.0031887 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | SWghtLB | -0.008564 | 0.1 | 1 |
+#> | ESETA1.24 | SWghtLB | 0.0047081 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | SurpPreg_1 | -0.037616 | 0.1 | 1 |
+#> | ESETA1.24 | SurpPreg_1 | 0.0027554 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | SmokTotl_1 | 0.011311 | 0.1 | 1 |
+#> | ESETA1.24 | SmokTotl_1 | 0.0039094 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | SmokTotl_2 | -0.000889 | 0.1 | 1 |
+#> | ESETA1.24 | SmokTotl_2 | -0.0184712 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | SmokTotl_3 | 0.002071 | 0.1 | 1 |
+#> | ESETA1.24 | SmokTotl_3 | 0.0039583 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | DrnkFreq_1 | -0.011027 | 0.1 | 1 |
+#> | ESETA1.24 | DrnkFreq_1 | 0.0133345 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | DrnkFreq_2 | 0.012054 | 0.1 | 1 |
+#> | ESETA1.24 | DrnkFreq_2 | -0.0110411 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | peri_health | -0.040005 | 0.1 | 1 |
+#> | ESETA1.24 | peri_health | -0.0071191 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | caregiv_health | 0.002807 | 0.1 | 1 |
+#> | ESETA1.24 | caregiv_health | -0.0052552 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | gov_assist | 0.017188 | 0.1 | 1 |
+#> | ESETA1.24 | gov_assist | -0.0017485 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | SAAmylase.6 | 0.015602 | 0.1 | 1 |
+#> | ESETA1.24 | SAAmylase.6 | 0.0159673 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | SAAmylase.15 | 0.043926 | 0.1 | 1 |
+#> | ESETA1.24 | SAAmylase.15 | 0.0082825 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | MDI.6 | 0.02066 | 0.1 | 1 |
+#> | ESETA1.24 | MDI.6 | 0.0094737 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | MDI.15 | 0.002145 | 0.1 | 1 |
+#> | ESETA1.24 | MDI.15 | -0.0034069 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RHasSO.6_1 | -0.027332 | 0.1 | 1 |
+#> | ESETA1.24 | RHasSO.6_1 | -0.0140331 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | RHasSO.15_1 | 0.000214 | 0.1 | 1 |
+#> | ESETA1.24 | RHasSO.15_1 | -0.001231 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | WndNbrhood.6 | -0.000298 | 0.1 | 1 |
+#> | ESETA1.24 | WndNbrhood.6 | 0.0017775 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | IBRAttn.6 | 0.015692 | 0.1 | 1 |
+#> | ESETA1.24 | IBRAttn.6 | -0.0016465 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | IBRAttn.15 | 0.004147 | 0.1 | 1 |
+#> | ESETA1.24 | IBRAttn.15 | -0.000343 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | B18Raw.6 | 0.093979 | 0.1 | 1 |
+#> | ESETA1.24 | B18Raw.6 | 0.0822361 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | B18Raw.15 | -0.003861 | 0.1 | 1 |
+#> | ESETA1.24 | B18Raw.15 | 0.0005663 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | HOMEETA1.6 | 0.006479 | 0.1 | 1 |
+#> | ESETA1.24 | HOMEETA1.6 | 0.0113979 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | HOMEETA1.15 | 0.016107 | 0.1 | 1 |
+#> | ESETA1.24 | HOMEETA1.15 | 0.0067273 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | InRatioCor.6 | -0.043574 | 0.05 | 1 |
+#> | ESETA1.24 | InRatioCor.6 | -0.0374541 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | InRatioCor.15 | -0.002573 | 0.05 | 1 |
+#> | ESETA1.24 | InRatioCor.15 | 0.0098439 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | CORTB.6 | 0.040583 | 0.1 | 1 |
+#> | ESETA1.24 | CORTB.6 | 0.0563049 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | CORTB.15 | -0.007733 | 0.1 | 1 |
+#> | ESETA1.24 | CORTB.15 | -0.005924 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | ESETA1.6 | 0.038786 | 0.1 | 1 |
+#> | ESETA1.24 | ESETA1.6 | 0.0273448 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.24 | ESETA1.15 | 0.016269 | 0.1 | 1 |
+#> | ESETA1.24 | ESETA1.15 | 0.0061281 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | state_1 | 0.012297 | 0.1 | 1 |
+#> | ESETA1.35 | state_1 | 0.0159812 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | BioDadInHH2_1 | -0.007364 | 0.1 | 1 |
+#> | ESETA1.35 | BioDadInHH2_1 | -0.0216829 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | PmAge2 | -0.011717 | 0.1 | 1 |
+#> | ESETA1.35 | PmAge2 | -0.0074947 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | PmBlac2_1 | -0.001426 | 0.1 | 1 |
+#> | ESETA1.35 | PmBlac2_1 | 0.0103925 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | TcBlac2_1 | -0.009506 | 0.1 | 1 |
+#> | ESETA1.35 | TcBlac2_1 | 0.0163783 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | PmMrSt2_1 | 0.004603 | 0.1 | 1 |
+#> | ESETA1.35 | PmMrSt2_1 | -0.0010987 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | PmEd2 | -0.041457 | 0.05 | 1 |
+#> | ESETA1.35 | PmEd2 | -0.0389852 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | KFASTScr | -0.036504 | 0.1 | 1 |
+#> | ESETA1.35 | KFASTScr | -0.0250932 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RMomAgeU | 0.014532 | 0.1 | 1 |
+#> | ESETA1.35 | RMomAgeU | -0.0138737 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHealth_1 | -0.014423 | 0.1 | 1 |
+#> | ESETA1.35 | RHealth_1 | 0.0182202 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHealth_2 | -0.057253 | 0.1 | 1 |
+#> | ESETA1.35 | RHealth_2 | -0.0231654 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHealth_3 | 0.065055 | 0.1 | 1 |
+#> | ESETA1.35 | RHealth_3 | -0.0001511 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHealth_4 | 0.019599 | 0.1 | 1 |
+#> | ESETA1.35 | RHealth_4 | 0.0240331 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHealth_5 | 0.009704 | 0.1 | 1 |
+#> | ESETA1.35 | RHealth_5 | 0.005529 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | HomeOwnd_1 | 0.011584 | 0.1 | 1 |
+#> | ESETA1.35 | HomeOwnd_1 | -0.0251607 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | HomeOwnd_2 | 0.020508 | 0.1 | 1 |
+#> | ESETA1.35 | HomeOwnd_2 | 0.0322371 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | HomeOwnd_3 | -0.030714 | 0.1 | 1 |
+#> | ESETA1.35 | HomeOwnd_3 | -0.0198839 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | HomeOwnd_4 | 0.005852 | 0.1 | 1 |
+#> | ESETA1.35 | HomeOwnd_4 | 0.0124531 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SWghtLB | -0.019471 | 0.1 | 1 |
+#> | ESETA1.35 | SWghtLB | -0.0030648 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SurpPreg_1 | 0.059723 | 0.1 | 1 |
+#> | ESETA1.35 | SurpPreg_1 | 0.0266638 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SmokTotl_1 | 0.017933 | 0.1 | 1 |
+#> | ESETA1.35 | SmokTotl_1 | -0.0129774 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SmokTotl_2 | -0.036241 | 0.1 | 1 |
+#> | ESETA1.35 | SmokTotl_2 | -0.013428 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SmokTotl_3 | 0.0324 | 0.1 | 1 |
+#> | ESETA1.35 | SmokTotl_3 | 0.009222 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | DrnkFreq_1 | -0.007962 | 0.1 | 1 |
+#> | ESETA1.35 | DrnkFreq_1 | -0.0225268 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | DrnkFreq_2 | 0.006528 | 0.1 | 1 |
+#> | ESETA1.35 | DrnkFreq_2 | 0.0090972 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | peri_health | 0.024272 | 0.1 | 1 |
+#> | ESETA1.35 | peri_health | 0.0014589 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | caregiv_health | 0.001126 | 0.1 | 1 |
+#> | ESETA1.35 | caregiv_health | -0.0051781 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | gov_assist | 0.001664 | 0.1 | 1 |
+#> | ESETA1.35 | gov_assist | 0.0318938 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SAAmylase.6 | -0.02835 | 0.1 | 1 |
+#> | ESETA1.35 | SAAmylase.6 | -0.0262552 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SAAmylase.15 | -0.016089 | 0.1 | 1 |
+#> | ESETA1.35 | SAAmylase.15 | -0.0163234 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | SAAmylase.24 | 0.002003 | 0.1 | 1 |
+#> | ESETA1.35 | SAAmylase.24 | 0.0100267 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | MDI.6 | -0.003452 | 0.1 | 1 |
+#> | ESETA1.35 | MDI.6 | -0.0185937 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | MDI.15 | 0.016688 | 0.1 | 1 |
+#> | ESETA1.35 | MDI.15 | 0.0252672 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHasSO.6_1 | 0.027168 | 0.1 | 1 |
+#> | ESETA1.35 | RHasSO.6_1 | 0.0185998 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHasSO.15_1 | 0.066367 | 0.1 | 1 |
+#> | ESETA1.35 | RHasSO.15_1 | 0.0371793 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | RHasSO.24_1 | 0.059345 | 0.1 | 1 |
+#> | ESETA1.35 | RHasSO.24_1 | 0.0131871 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | WndNbrhood.6 | -0.045594 | 0.1 | 1 |
+#> | ESETA1.35 | WndNbrhood.6 | -0.0321717 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | WndNbrhood.24 | -0.023004 | 0.1 | 1 |
+#> | ESETA1.35 | WndNbrhood.24 | -0.022859 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | IBRAttn.6 | -0.051726 | 0.1 | 1 |
+#> | ESETA1.35 | IBRAttn.6 | -0.0604992 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | IBRAttn.15 | -0.039767 | 0.1 | 1 |
+#> | ESETA1.35 | IBRAttn.15 | -0.0174987 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | IBRAttn.24 | -0.001796 | 0.1 | 1 |
+#> | ESETA1.35 | IBRAttn.24 | 0.0059855 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | B18Raw.6 | 0.053306 | 0.1 | 1 |
+#> | ESETA1.35 | B18Raw.6 | 0.0402665 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | B18Raw.15 | 0.04418 | 0.1 | 1 |
+#> | ESETA1.35 | B18Raw.15 | 0.0266371 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | B18Raw.24 | 0.030505 | 0.1 | 1 |
+#> | ESETA1.35 | B18Raw.24 | 0.0130311 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | HOMEETA1.6 | -0.034064 | 0.1 | 1 |
+#> | ESETA1.35 | HOMEETA1.6 | -0.0573771 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | HOMEETA1.15 | -0.037334 | 0.1 | 1 |
+#> | ESETA1.35 | HOMEETA1.15 | -0.020653 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | HOMEETA1.24 | -0.016869 | 0.1 | 1 |
+#> | ESETA1.35 | HOMEETA1.24 | -0.0066419 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | InRatioCor.6 | -0.03603 | 0.05 | 1 |
+#> | ESETA1.35 | InRatioCor.6 | -0.0650277 | 0.05 | 0 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | InRatioCor.15 | -0.022606 | 0.05 | 1 |
+#> | ESETA1.35 | InRatioCor.15 | -0.0609491 | 0.05 | 0 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | InRatioCor.24 | -0.014519 | 0.05 | 1 |
+#> | ESETA1.35 | InRatioCor.24 | -0.036796 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | CORTB.6 | 0.034283 | 0.1 | 1 |
+#> | ESETA1.35 | CORTB.6 | 0.0382809 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | CORTB.15 | 0.00186 | 0.1 | 1 |
+#> | ESETA1.35 | CORTB.15 | 0.0095674 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | CORTB.24 | 0.021881 | 0.1 | 1 |
+#> | ESETA1.35 | CORTB.24 | 0.0266067 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | EARS_TJo.24 | -0.009914 | 0.1 | 1 |
+#> | ESETA1.35 | EARS_TJo.24 | -0.0122612 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | LESMnPos.24 | -0.015144 | 0.1 | 1 |
+#> | ESETA1.35 | LESMnPos.24 | -0.0106439 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | LESMnNeg.24 | 0.032104 | 0.1 | 1 |
+#> | ESETA1.35 | LESMnNeg.24 | 0.0081802 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | ESETA1.6 | 0.084487 | 0.1 | 1 |
+#> | ESETA1.35 | ESETA1.6 | 0.0844947 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | ESETA1.15 | 0.08033 | 0.1 | 1 |
+#> | ESETA1.35 | ESETA1.15 | 0.0775117 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.35 | ESETA1.24 | 0.003504 | 0.1 | 1 |
+#> | ESETA1.35 | ESETA1.24 | 0.0204742 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | state_1 | 0.046707 | 0.1 | 1 |
+#> | ESETA1.58 | state_1 | 0.0173188 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | BioDadInHH2_1 | -0.032556 | 0.1 | 1 |
+#> | ESETA1.58 | BioDadInHH2_1 | -0.029331 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | PmAge2 | -0.011158 | 0.1 | 1 |
+#> | ESETA1.58 | PmAge2 | -0.0022438 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | PmBlac2_1 | 0.045969 | 0.1 | 1 |
+#> | ESETA1.58 | PmBlac2_1 | 0.0167733 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | TcBlac2_1 | 0.062286 | 0.1 | 1 |
+#> | ESETA1.58 | TcBlac2_1 | 0.0285232 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | PmMrSt2_1 | -0.02112 | 0.1 | 1 |
+#> | ESETA1.58 | PmMrSt2_1 | -0.0308851 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | PmEd2 | 0.003704 | 0.05 | 1 |
+#> | ESETA1.58 | PmEd2 | -0.005835 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | KFASTScr | 0.001911 | 0.1 | 1 |
+#> | ESETA1.58 | KFASTScr | 0.0001801 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RMomAgeU | 0.042917 | 0.1 | 1 |
+#> | ESETA1.58 | RMomAgeU | 0.0200869 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHealth_1 | 0.001836 | 0.1 | 1 |
+#> | ESETA1.58 | RHealth_1 | -0.0046526 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHealth_2 | 0.011966 | 0.1 | 1 |
+#> | ESETA1.58 | RHealth_2 | 0.0058821 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHealth_3 | -0.000588 | 0.1 | 1 |
+#> | ESETA1.58 | RHealth_3 | 0.0057173 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHealth_4 | -0.031635 | 0.1 | 1 |
+#> | ESETA1.58 | RHealth_4 | -0.0099921 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHealth_5 | 0.03408 | 0.1 | 1 |
+#> | ESETA1.58 | RHealth_5 | -0.0018151 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HomeOwnd_1 | -0.020218 | 0.1 | 1 |
+#> | ESETA1.58 | HomeOwnd_1 | -0.0129367 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HomeOwnd_2 | 0.024941 | 0.1 | 1 |
+#> | ESETA1.58 | HomeOwnd_2 | 0.0088235 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HomeOwnd_3 | 0.002759 | 0.1 | 1 |
+#> | ESETA1.58 | HomeOwnd_3 | 0.029316 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HomeOwnd_4 | 0.002021 | 0.1 | 1 |
+#> | ESETA1.58 | HomeOwnd_4 | -0.0281944 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SWghtLB | -0.022619 | 0.1 | 1 |
+#> | ESETA1.58 | SWghtLB | -0.0300973 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SurpPreg_1 | -0.018775 | 0.1 | 1 |
+#> | ESETA1.58 | SurpPreg_1 | -0.0039986 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SmokTotl_1 | -0.044081 | 0.1 | 1 |
+#> | ESETA1.58 | SmokTotl_1 | 0.002585 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SmokTotl_2 | 0.003514 | 0.1 | 1 |
+#> | ESETA1.58 | SmokTotl_2 | -0.0052944 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SmokTotl_3 | -0.008938 | 0.1 | 1 |
+#> | ESETA1.58 | SmokTotl_3 | -0.0020945 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | DrnkFreq_1 | 0.010825 | 0.1 | 1 |
+#> | ESETA1.58 | DrnkFreq_1 | -0.0068799 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | DrnkFreq_2 | -0.012394 | 0.1 | 1 |
+#> | ESETA1.58 | DrnkFreq_2 | 0.0052751 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | peri_health | -0.019048 | 0.1 | 1 |
+#> | ESETA1.58 | peri_health | 0.0042015 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | caregiv_health | -0.001676 | 0.1 | 1 |
+#> | ESETA1.58 | caregiv_health | 0.0116077 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | gov_assist | 0.01349 | 0.1 | 1 |
+#> | ESETA1.58 | gov_assist | 0.007267 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SAAmylase.6 | 0.028656 | 0.1 | 1 |
+#> | ESETA1.58 | SAAmylase.6 | 0.0316109 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SAAmylase.15 | -0.02358 | 0.1 | 1 |
+#> | ESETA1.58 | SAAmylase.15 | -0.030204 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | SAAmylase.24 | -0.017473 | 0.1 | 1 |
+#> | ESETA1.58 | SAAmylase.24 | -0.0039453 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | MDI.6 | 0.035384 | 0.1 | 1 |
+#> | ESETA1.58 | MDI.6 | 0.0458769 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | MDI.15 | -0.069623 | 0.1 | 1 |
+#> | ESETA1.58 | MDI.15 | -0.046456 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHasSO.6_1 | 0.066159 | 0.1 | 1 |
+#> | ESETA1.58 | RHasSO.6_1 | 0.0649412 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHasSO.15_1 | 0.005755 | 0.1 | 1 |
+#> | ESETA1.58 | RHasSO.15_1 | -0.0122701 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHasSO.24_1 | -0.025656 | 0.1 | 1 |
+#> | ESETA1.58 | RHasSO.24_1 | -0.0169821 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | RHasSO.35_1 | -0.013387 | 0.1 | 1 |
+#> | ESETA1.58 | RHasSO.35_1 | 0.0116231 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | WndNbrhood.6 | -0.042805 | 0.1 | 1 |
+#> | ESETA1.58 | WndNbrhood.6 | -0.0373402 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | WndNbrhood.24 | -0.006412 | 0.1 | 1 |
+#> | ESETA1.58 | WndNbrhood.24 | -0.001917 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | WndNbrhood.35 | 0.003724 | 0.1 | 1 |
+#> | ESETA1.58 | WndNbrhood.35 | 0.0023001 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | IBRAttn.6 | -0.04382 | 0.1 | 1 |
+#> | ESETA1.58 | IBRAttn.6 | -0.0320238 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | IBRAttn.15 | -0.01477 | 0.1 | 1 |
+#> | ESETA1.58 | IBRAttn.15 | -0.0253196 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | IBRAttn.24 | -0.048202 | 0.1 | 1 |
+#> | ESETA1.58 | IBRAttn.24 | -0.0374467 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | B18Raw.6 | 0.03759 | 0.1 | 1 |
+#> | ESETA1.58 | B18Raw.6 | 0.016658 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | B18Raw.15 | 0.078278 | 0.1 | 1 |
+#> | ESETA1.58 | B18Raw.15 | 0.0519722 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | B18Raw.24 | 0.051997 | 0.1 | 1 |
+#> | ESETA1.58 | B18Raw.24 | 0.0765176 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HOMEETA1.6 | -0.030969 | 0.1 | 1 |
+#> | ESETA1.58 | HOMEETA1.6 | -0.0038499 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HOMEETA1.15 | 0.004609 | 0.1 | 1 |
+#> | ESETA1.58 | HOMEETA1.15 | 0.0141684 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HOMEETA1.24 | -0.047911 | 0.1 | 1 |
+#> | ESETA1.58 | HOMEETA1.24 | -0.0396149 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | HOMEETA1.35 | -0.07415 | 0.1 | 1 |
+#> | ESETA1.58 | HOMEETA1.35 | -0.0181722 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | InRatioCor.6 | -0.036626 | 0.05 | 1 |
+#> | ESETA1.58 | InRatioCor.6 | -0.021147 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | InRatioCor.15 | -0.046776 | 0.05 | 1 |
+#> | ESETA1.58 | InRatioCor.15 | -0.0288403 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | InRatioCor.24 | -0.044794 | 0.05 | 1 |
+#> | ESETA1.58 | InRatioCor.24 | -0.0306187 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | InRatioCor.35 | -0.021949 | 0.05 | 1 |
+#> | ESETA1.58 | InRatioCor.35 | -0.0000517 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | CORTB.6 | 0.024018 | 0.1 | 1 |
+#> | ESETA1.58 | CORTB.6 | 0.034015 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | CORTB.15 | -0.008257 | 0.1 | 1 |
+#> | ESETA1.58 | CORTB.15 | 0.0126726 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | CORTB.24 | 0.004373 | 0.1 | 1 |
+#> | ESETA1.58 | CORTB.24 | -0.0234976 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | EARS_TJo.24 | -0.018775 | 0.1 | 1 |
+#> | ESETA1.58 | EARS_TJo.24 | -0.0205811 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | EARS_TJo.35 | 0.023288 | 0.1 | 1 |
+#> | ESETA1.58 | EARS_TJo.35 | -0.0102582 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | LESMnPos.24 | -0.019433 | 0.1 | 1 |
+#> | ESETA1.58 | LESMnPos.24 | -0.0325432 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | LESMnPos.35 | -0.016828 | 0.1 | 1 |
+#> | ESETA1.58 | LESMnPos.35 | -0.0198356 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | LESMnNeg.24 | 0.027503 | 0.1 | 1 |
+#> | ESETA1.58 | LESMnNeg.24 | 0.0397038 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | LESMnNeg.35 | 0.054757 | 0.1 | 1 |
+#> | ESETA1.58 | LESMnNeg.35 | 0.0159532 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | StrDif_Tot.35 | 0.022732 | 0.1 | 1 |
+#> | ESETA1.58 | StrDif_Tot.35 | 0.025157 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | fscore.35 | -0.059882 | 0.1 | 1 |
+#> | ESETA1.58 | fscore.35 | -0.0372574 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | ESETA1.6 | 0.039731 | 0.1 | 1 |
+#> | ESETA1.58 | ESETA1.6 | 0.0391978 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | ESETA1.15 | 0.031636 | 0.1 | 1 |
+#> | ESETA1.58 | ESETA1.15 | 0.0257108 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | ESETA1.24 | 0.014185 | 0.1 | 1 |
+#> | ESETA1.58 | ESETA1.24 | 0.0109497 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.58 | ESETA1.35 | 0.019732 | 0.1 | 1 |
+#> | ESETA1.58 | ESETA1.35 | 0.0073936 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
#>
#> Table: Balance Stats for All Exposure Time Points Averaging Across Imputed Datasets
@@ -2649,7 +2812,7 @@ Workflow: Continuous Exposure
#> ESETA1.24 ~ state + BioDadInHH2 + PmAge2 + PmBlac2 + TcBlac2 + PmMrSt2 + PmEd2 + KFASTScr + RMomAgeU + RHealth + HomeOwnd + SWghtLB + SurpPreg + SmokTotl + DrnkFreq + peri_health + caregiv_health + gov_assist + SAAmylase.15 + MDI.15 + RHasSO.15 + IBRAttn.15 + B18Raw.15 + HOMEETA1.15 + InRatioCor.15 + CORTB.15 + ESETA1.15
#> USER ALERT: Please manually inspect the updated balancing formula below that includes time-varying confounders at t-1 and those greater at further lags that remained imbalanced:
#> At time point 35, the update formula for ESETA1.35 is:
-#> ESETA1.35 ~ state + BioDadInHH2 + PmAge2 + PmBlac2 + TcBlac2 + PmMrSt2 + PmEd2 + KFASTScr + RMomAgeU + RHealth + HomeOwnd + SWghtLB + SurpPreg + SmokTotl + DrnkFreq + peri_health + caregiv_health + gov_assist + SAAmylase.24 + RHasSO.24 + WndNbrhood.24 + IBRAttn.24 + B18Raw.24 + HOMEETA1.24 + InRatioCor.24 + CORTB.24 + EARS_TJo.24 + LESMnPos.24 + LESMnNeg.24 + ESETA1.24
+#> ESETA1.35 ~ state + BioDadInHH2 + PmAge2 + PmBlac2 + TcBlac2 + PmMrSt2 + PmEd2 + KFASTScr + RMomAgeU + RHealth + HomeOwnd + SWghtLB + SurpPreg + SmokTotl + DrnkFreq + peri_health + caregiv_health + gov_assist + SAAmylase.24 + RHasSO.24 + WndNbrhood.24 + IBRAttn.24 + B18Raw.24 + HOMEETA1.24 + InRatioCor.24 + CORTB.24 + EARS_TJo.24 + LESMnPos.24 + LESMnNeg.24 + ESETA1.24 + InRatioCor.6 + InRatioCor.15
#> USER ALERT: Please manually inspect the updated balancing formula below that includes time-varying confounders at t-1 and those greater at further lags that remained imbalanced:
#> At time point 58, the update formula for ESETA1.58 is:
#> ESETA1.58 ~ state + BioDadInHH2 + PmAge2 + PmBlac2 + TcBlac2 + PmMrSt2 + PmEd2 + KFASTScr + RMomAgeU + RHealth + HomeOwnd + SWghtLB + SurpPreg + SmokTotl + DrnkFreq + peri_health + caregiv_health + gov_assist + RHasSO.35 + WndNbrhood.35 + HOMEETA1.35 + InRatioCor.35 + EARS_TJo.35 + LESMnPos.35 + LESMnNeg.35 + StrDif_Tot.35 + fscore.35 + ESETA1.35
@@ -2673,7 +2836,7 @@ Workflow: Continuous Exposure
obj = obj,
method = method,
formulas = formulas,
- max.it = 1, # testing only
+ # max.it = 1, # testing only
save.out = save.out)
We then inspect these final weights by imputed dataset. The user has the option to supply save.out
to save printed or plot output to the home directory.
@@ -2682,7 +2845,7 @@ Workflow: Continuous Exposure
print(final_weights,
i = 1)
#>
-#> For imputation 1 and the `cbps` weighting method, the median weight value is 0.77 (SD = 1.18; range = 0-9).
+#> For imputation 1 and the `cbps` weighting method, the median weight value is 0.77 (SD = 1.17; range = 0-9).
plot(final_weights,
i = 1,
@@ -2776,7 +2939,7 @@ Workflow: Continuous Exposure
print(trim_weights.s2,
i = 1)
#>
-#> For imputation 1 and the `cbps` weighting method, after trimming at 98th quantile, the median weight value is 0.77 (SD = 1.06; range = 0-5).
+#> For imputation 1 and the `cbps` weighting method, after trimming at 98th quantile, the median weight value is 0.77 (SD = 1.05; range = 0-5).
plot(trim_weights.s2,
i = 1,
@@ -2816,17 +2979,17 @@ Workflow: Continuous Exposure
summary(final_balance_stats_trim,
save.out = save.out)
-#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 2 out of 241 (0.8%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.06):
+#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 1 out of 241 (0.4%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.05 (max: 0.05):
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
#> | ESETA1.6 | 28 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.15 | 38 | 0 |
+#> | ESETA1.15 | 38 | 1 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 2 |
+#> | ESETA1.35 | 59 | 0 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
@@ -2840,81 +3003,81 @@ Workflow: Continuous Exposure
#> +-----------+----------------+---------------+------------+----------+
#> | exposure | covariate | std_bal_stats | bal_thresh | balanced |
#> +===========+================+===============+============+==========+
-#> | ESETA1.15 | state_1 | 0.018002 | 0.1 | 1 |
+#> | ESETA1.15 | state_1 | 0.017188 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | BioDadInHH2_1 | -0.01641 | 0.1 | 1 |
+#> | ESETA1.15 | BioDadInHH2_1 | -0.0163 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmAge2 | -0.009705 | 0.1 | 1 |
+#> | ESETA1.15 | PmAge2 | -0.008023 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmBlac2_1 | 0.017384 | 0.1 | 1 |
+#> | ESETA1.15 | PmBlac2_1 | 0.016619 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | TcBlac2_1 | 0.017685 | 0.1 | 1 |
+#> | ESETA1.15 | TcBlac2_1 | 0.018256 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmMrSt2_1 | -0.00618 | 0.1 | 1 |
+#> | ESETA1.15 | PmMrSt2_1 | -0.00511 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmEd2 | -0.036545 | 0.05 | 1 |
+#> | ESETA1.15 | PmEd2 | -0.037445 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | KFASTScr | -0.020649 | 0.1 | 1 |
+#> | ESETA1.15 | KFASTScr | -0.019385 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RMomAgeU | -0.005992 | 0.1 | 1 |
+#> | ESETA1.15 | RMomAgeU | -0.004175 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_1 | -0.001056 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_1 | 0.002172 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_2 | -0.0065 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_2 | -0.008696 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_3 | 0.004525 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_3 | 0.004039 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_4 | 0.006969 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_4 | 0.007655 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_5 | -0.011918 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_5 | -0.011764 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_1 | -0.003678 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_1 | -0.004551 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_2 | 0.003725 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_2 | 0.003749 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_3 | 0.006049 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_3 | 0.004181 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_4 | 0.008359 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_4 | 0.010203 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SWghtLB | -0.000996 | 0.1 | 1 |
+#> | ESETA1.15 | SWghtLB | -0.000544 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SurpPreg_1 | 0.008394 | 0.1 | 1 |
+#> | ESETA1.15 | SurpPreg_1 | 0.00931 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_1 | 0.010939 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_1 | 0.010631 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_2 | 0.009158 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_2 | 0.008245 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_3 | -0.005918 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_3 | -0.005169 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | DrnkFreq_1 | 0.007173 | 0.1 | 1 |
+#> | ESETA1.15 | DrnkFreq_1 | 0.009757 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | DrnkFreq_2 | -0.017437 | 0.1 | 1 |
+#> | ESETA1.15 | DrnkFreq_2 | -0.01846 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | peri_health | -0.000413 | 0.1 | 1 |
+#> | ESETA1.15 | peri_health | -0.001012 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | caregiv_health | 0.008989 | 0.1 | 1 |
+#> | ESETA1.15 | caregiv_health | 0.009529 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | gov_assist | 0.015211 | 0.1 | 1 |
+#> | ESETA1.15 | gov_assist | 0.014943 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SAAmylase.6 | 0.010964 | 0.1 | 1 |
+#> | ESETA1.15 | SAAmylase.6 | 0.009242 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | MDI.6 | 0.01114 | 0.1 | 1 |
+#> | ESETA1.15 | MDI.6 | 0.010527 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHasSO.6_1 | 0.018068 | 0.1 | 1 |
+#> | ESETA1.15 | RHasSO.6_1 | 0.02075 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | WndNbrhood.6 | -0.017711 | 0.1 | 1 |
+#> | ESETA1.15 | WndNbrhood.6 | -0.016214 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | IBRAttn.6 | -0.002143 | 0.1 | 1 |
+#> | ESETA1.15 | IBRAttn.6 | -0.003919 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | B18Raw.6 | 0.025159 | 0.1 | 1 |
+#> | ESETA1.15 | B18Raw.6 | 0.027468 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HOMEETA1.6 | -0.024603 | 0.1 | 1 |
+#> | ESETA1.15 | HOMEETA1.6 | -0.023461 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | InRatioCor.6 | -0.047661 | 0.05 | 1 |
+#> | ESETA1.15 | InRatioCor.6 | -0.050616 | 0.05 | 0 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | CORTB.6 | -0.005651 | 0.1 | 1 |
+#> | ESETA1.15 | CORTB.6 | -0.005384 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | ESETA1.6 | 0.03469 | 0.1 | 1 |
+#> | ESETA1.15 | ESETA1.6 | 0.034089 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
#>
#> Table: Balance Stats for Exposure Time 2 Averaging Across Imputed Datasets
@@ -2925,120 +3088,120 @@ Workflow: Continuous Exposure
#>
#> - Weight ranges:
#>
-#> Min Max
-#> all 0 |---------------------------| 3.528
+#> Min Max
+#> all 0 |---------------------------| 3.5109
#>
#> - Units with the 5 most extreme weights:
-#>
-#> 89 79 30 27 11
-#> all 3.528 3.528 3.528 3.528 3.528
+#>
+#> 89 79 30 27 11
+#> all 3.5109 3.5109 3.5109 3.5109 3.5109
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 0.896 0.687 0.372 0
+#> all 0.893 0.685 0.371 0
#>
#> - Effective Sample Sizes:
#>
-#> Total
-#> Unweighted 1292.
-#> Weighted 717.08
+#> Total
+#> Unweighted 1292.
+#> Weighted 718.9
#>
#> Time 2
#> Summary of weights
#>
#> - Weight ranges:
#>
-#> Min Max
-#> all 0 |---------------------------| 3.528
+#> Min Max
+#> all 0 |---------------------------| 3.5109
#>
#> - Units with the 5 most extreme weights:
-#>
-#> 89 79 30 27 11
-#> all 3.528 3.528 3.528 3.528 3.528
+#>
+#> 89 79 30 27 11
+#> all 3.5109 3.5109 3.5109 3.5109 3.5109
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 0.896 0.687 0.372 0
+#> all 0.893 0.685 0.371 0
#>
#> - Effective Sample Sizes:
#>
-#> Total
-#> Unweighted 1292.
-#> Weighted 717.08
+#> Total
+#> Unweighted 1292.
+#> Weighted 718.9
#>
#> Time 3
#> Summary of weights
#>
#> - Weight ranges:
#>
-#> Min Max
-#> all 0 |---------------------------| 3.528
+#> Min Max
+#> all 0 |---------------------------| 3.5109
#>
#> - Units with the 5 most extreme weights:
-#>
-#> 89 79 30 27 11
-#> all 3.528 3.528 3.528 3.528 3.528
+#>
+#> 89 79 30 27 11
+#> all 3.5109 3.5109 3.5109 3.5109 3.5109
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 0.896 0.687 0.372 0
+#> all 0.893 0.685 0.371 0
#>
#> - Effective Sample Sizes:
#>
-#> Total
-#> Unweighted 1292.
-#> Weighted 717.08
+#> Total
+#> Unweighted 1292.
+#> Weighted 718.9
#>
#> Time 4
#> Summary of weights
#>
#> - Weight ranges:
#>
-#> Min Max
-#> all 0 |---------------------------| 3.528
+#> Min Max
+#> all 0 |---------------------------| 3.5109
#>
#> - Units with the 5 most extreme weights:
-#>
-#> 89 79 30 27 11
-#> all 3.528 3.528 3.528 3.528 3.528
+#>
+#> 89 79 30 27 11
+#> all 3.5109 3.5109 3.5109 3.5109 3.5109
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 0.896 0.687 0.372 0
+#> all 0.893 0.685 0.371 0
#>
#> - Effective Sample Sizes:
#>
-#> Total
-#> Unweighted 1292.
-#> Weighted 717.08
+#> Total
+#> Unweighted 1292.
+#> Weighted 718.9
#>
#> Time 5
#> Summary of weights
#>
#> - Weight ranges:
#>
-#> Min Max
-#> all 0 |---------------------------| 3.528
+#> Min Max
+#> all 0 |---------------------------| 3.5109
#>
#> - Units with the 5 most extreme weights:
-#>
-#> 89 79 30 27 11
-#> all 3.528 3.528 3.528 3.528 3.528
+#>
+#> 89 79 30 27 11
+#> all 3.5109 3.5109 3.5109 3.5109 3.5109
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 0.896 0.687 0.372 0
+#> all 0.893 0.685 0.371 0
#>
#> - Effective Sample Sizes:
#>
-#> Total
-#> Unweighted 1292.
-#> Weighted 717.08
+#> Total
+#> Unweighted 1292.
+#> Weighted 718.9
As shown above, the trimmed weights result in one imbalanced confounder in relation to exposure at 15 months and an effective sample size of 719.
Untrimmed
@@ -3060,7 +3223,8 @@ Workflow: Continuous Exposure
summary(final_balance_stats_untrim,
save.out = save.out)
-#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 2 out of 241 (0.8%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.07):
+#> Averaging across imputed datasets using cbps weighting method.
+#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 0 out of 241 (0.0%) covariates across time points remain imbalanced:
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
@@ -3070,7 +3234,7 @@ Workflow: Continuous Exposure
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 2 |
+#> | ESETA1.35 | 59 | 0 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
@@ -3084,23 +3248,23 @@ Workflow: Continuous Exposure
#> - Weight ranges:
#>
#> Min Max
-#> all 0 |---------------------------| 8.5138
+#> all 0 |---------------------------| 8.5407
#>
#> - Units with the 5 most extreme weights:
#>
-#> 1246 1145 294 697 919
-#> all 7.4581 7.6788 7.9902 8.2032 8.5138
+#> 877 1145 697 294 919
+#> all 7.1323 7.3944 7.7172 8.5237 8.5407
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 1.066 0.725 0.443 0
+#> all 1.056 0.722 0.439 0
#>
#> - Effective Sample Sizes:
#>
#> Total
#> Unweighted 1292.
-#> Weighted 605.04
+#> Weighted 611.17
#>
#> Time 2
#> Summary of weights
@@ -3108,23 +3272,23 @@ Workflow: Continuous Exposure
#> - Weight ranges:
#>
#> Min Max
-#> all 0 |---------------------------| 8.5138
+#> all 0 |---------------------------| 8.5407
#>
#> - Units with the 5 most extreme weights:
#>
-#> 1246 1145 294 697 919
-#> all 7.4581 7.6788 7.9902 8.2032 8.5138
+#> 877 1145 697 294 919
+#> all 7.1323 7.3944 7.7172 8.5237 8.5407
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 1.066 0.725 0.443 0
+#> all 1.056 0.722 0.439 0
#>
#> - Effective Sample Sizes:
#>
#> Total
#> Unweighted 1292.
-#> Weighted 605.04
+#> Weighted 611.17
#>
#> Time 3
#> Summary of weights
@@ -3132,23 +3296,23 @@ Workflow: Continuous Exposure
#> - Weight ranges:
#>
#> Min Max
-#> all 0 |---------------------------| 8.5138
+#> all 0 |---------------------------| 8.5407
#>
#> - Units with the 5 most extreme weights:
#>
-#> 1246 1145 294 697 919
-#> all 7.4581 7.6788 7.9902 8.2032 8.5138
+#> 877 1145 697 294 919
+#> all 7.1323 7.3944 7.7172 8.5237 8.5407
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 1.066 0.725 0.443 0
+#> all 1.056 0.722 0.439 0
#>
#> - Effective Sample Sizes:
#>
#> Total
#> Unweighted 1292.
-#> Weighted 605.04
+#> Weighted 611.17
#>
#> Time 4
#> Summary of weights
@@ -3156,23 +3320,23 @@ Workflow: Continuous Exposure
#> - Weight ranges:
#>
#> Min Max
-#> all 0 |---------------------------| 8.5138
+#> all 0 |---------------------------| 8.5407
#>
#> - Units with the 5 most extreme weights:
#>
-#> 1246 1145 294 697 919
-#> all 7.4581 7.6788 7.9902 8.2032 8.5138
+#> 877 1145 697 294 919
+#> all 7.1323 7.3944 7.7172 8.5237 8.5407
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 1.066 0.725 0.443 0
+#> all 1.056 0.722 0.439 0
#>
#> - Effective Sample Sizes:
#>
#> Total
#> Unweighted 1292.
-#> Weighted 605.04
+#> Weighted 611.17
#>
#> Time 5
#> Summary of weights
@@ -3180,23 +3344,23 @@ Workflow: Continuous Exposure
#> - Weight ranges:
#>
#> Min Max
-#> all 0 |---------------------------| 8.5138
+#> all 0 |---------------------------| 8.5407
#>
#> - Units with the 5 most extreme weights:
#>
-#> 1246 1145 294 697 919
-#> all 7.4581 7.6788 7.9902 8.2032 8.5138
+#> 877 1145 697 294 919
+#> all 7.1323 7.3944 7.7172 8.5237 8.5407
#>
#> - Weight statistics:
#>
#> Coef of Var MAD Entropy # Zeros
-#> all 1.066 0.725 0.443 0
+#> all 1.056 0.722 0.439 0
#>
#> - Effective Sample Sizes:
#>
#> Total
#> Unweighted 1292.
-#> Weighted 605.04
+#> Weighted 611.17
We see that the non-trimmed weights result in no imbalanced confounders and an effective sample size of 611.
In this scenario, given that, for the trimmed weights, the imbalanced confounder (income at 6 months) is only related to exposure at -0.051 and the trimmed weights yield a higher effective sample size, we proceed with the trimmed weights.
@@ -3205,17 +3369,17 @@ Workflow: Continuous Exposure
summary(final_balance_stats_trim,
save.out = save.out)
-#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 2 out of 241 (0.8%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.06):
+#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 1 out of 241 (0.4%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.05 (max: 0.05):
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
#> | ESETA1.6 | 28 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.15 | 38 | 0 |
+#> | ESETA1.15 | 38 | 1 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 2 |
+#> | ESETA1.35 | 59 | 0 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
@@ -3255,7 +3419,7 @@ Workflow: Continuous Exposure
summary(final_balance_stats.s1,
save.out = save.out)
-#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 3 out of 241 (1.2%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.06):
+#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 1 out of 241 (0.4%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.06):
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
@@ -3265,7 +3429,7 @@ Workflow: Continuous Exposure
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 2 |
+#> | ESETA1.35 | 59 | 0 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
@@ -3279,81 +3443,81 @@ Workflow: Continuous Exposure
#> +-----------+----------------+---------------+------------+----------+
#> | exposure | covariate | std_bal_stats | bal_thresh | balanced |
#> +===========+================+===============+============+==========+
-#> | ESETA1.15 | state_1 | 0.02248 | 0.1 | 1 |
+#> | ESETA1.15 | state_1 | 0.02175 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | BioDadInHH2_1 | -0.02024 | 0.1 | 1 |
+#> | ESETA1.15 | BioDadInHH2_1 | -0.01857 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmAge2 | -0.01219 | 0.1 | 1 |
+#> | ESETA1.15 | PmAge2 | -0.00978 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmBlac2_1 | 0.0244 | 0.1 | 1 |
+#> | ESETA1.15 | PmBlac2_1 | 0.02279 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | TcBlac2_1 | 0.02265 | 0.1 | 1 |
+#> | ESETA1.15 | TcBlac2_1 | 0.02269 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmMrSt2_1 | -0.00906 | 0.1 | 1 |
+#> | ESETA1.15 | PmMrSt2_1 | -0.00654 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | PmEd2 | -0.04175 | 0.05 | 1 |
+#> | ESETA1.15 | PmEd2 | -0.04287 | 0.05 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | KFASTScr | -0.02587 | 0.1 | 1 |
+#> | ESETA1.15 | KFASTScr | -0.02611 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RMomAgeU | -0.01202 | 0.1 | 1 |
+#> | ESETA1.15 | RMomAgeU | -0.00973 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_1 | -0.00346 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_1 | -0.00117 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_2 | -0.00661 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_2 | -0.00813 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_3 | 0.00912 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_3 | 0.00942 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_4 | 0.00351 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_4 | 0.00349 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHealth_5 | -0.01245 | 0.1 | 1 |
+#> | ESETA1.15 | RHealth_5 | -0.01243 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_1 | -0.00579 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_1 | -0.00602 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_2 | 0.00296 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_2 | 0.00153 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_3 | 0.00995 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_3 | 0.00942 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HomeOwnd_4 | 0.00817 | 0.1 | 1 |
+#> | ESETA1.15 | HomeOwnd_4 | 0.00997 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SWghtLB | -0.00333 | 0.1 | 1 |
+#> | ESETA1.15 | SWghtLB | -0.00312 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SurpPreg_1 | 0.01431 | 0.1 | 1 |
+#> | ESETA1.15 | SurpPreg_1 | 0.01335 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_1 | 0.01323 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_1 | 0.01352 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_2 | 0.0114 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_2 | 0.01077 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SmokTotl_3 | -0.00775 | 0.1 | 1 |
+#> | ESETA1.15 | SmokTotl_3 | -0.00753 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | DrnkFreq_1 | 0.0125 | 0.1 | 1 |
+#> | ESETA1.15 | DrnkFreq_1 | 0.01492 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | DrnkFreq_2 | -0.02183 | 0.1 | 1 |
+#> | ESETA1.15 | DrnkFreq_2 | -0.02308 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | peri_health | -0.00145 | 0.1 | 1 |
+#> | ESETA1.15 | peri_health | -0.00114 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | caregiv_health | 0.01119 | 0.1 | 1 |
+#> | ESETA1.15 | caregiv_health | 0.01331 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | gov_assist | 0.01836 | 0.1 | 1 |
+#> | ESETA1.15 | gov_assist | 0.01817 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | SAAmylase.6 | 0.01106 | 0.1 | 1 |
+#> | ESETA1.15 | SAAmylase.6 | 0.01058 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | MDI.6 | 0.016 | 0.1 | 1 |
+#> | ESETA1.15 | MDI.6 | 0.01466 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | RHasSO.6_1 | 0.0208 | 0.1 | 1 |
+#> | ESETA1.15 | RHasSO.6_1 | 0.02011 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | WndNbrhood.6 | -0.02262 | 0.1 | 1 |
+#> | ESETA1.15 | WndNbrhood.6 | -0.02011 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | IBRAttn.6 | -0.0034 | 0.1 | 1 |
+#> | ESETA1.15 | IBRAttn.6 | -0.00386 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | B18Raw.6 | 0.02608 | 0.1 | 1 |
+#> | ESETA1.15 | B18Raw.6 | 0.0289 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | HOMEETA1.6 | -0.02883 | 0.1 | 1 |
+#> | ESETA1.15 | HOMEETA1.6 | -0.02621 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | InRatioCor.6 | -0.06024 | 0.05 | 0 |
+#> | ESETA1.15 | InRatioCor.6 | -0.06194 | 0.05 | 0 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | CORTB.6 | -0.00713 | 0.1 | 1 |
+#> | ESETA1.15 | CORTB.6 | -0.00842 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
-#> | ESETA1.15 | ESETA1.6 | 0.0427 | 0.1 | 1 |
+#> | ESETA1.15 | ESETA1.6 | 0.04248 | 0.1 | 1 |
#> +-----------+----------------+---------------+------------+----------+
#>
#> Table: Balance Stats for Exposure Time 2 Averaging Across Imputed Datasets
@@ -3374,7 +3538,8 @@ Workflow: Continuous Exposure
summary(final_balance_stats.s2,
save.out = save.out)
-#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 2 out of 241 (0.8%) covariates across time points remain imbalanced with a remaining median absolute correlation of 0.06 (max: 0.06):
+#> Averaging across imputed datasets using cbps weighting method.
+#> USER ALERT: Averaging across imputed datasets using `cbps` weighting method: As shown below, 0 out of 241 (0.0%) covariates across time points remain imbalanced:
#> +-----------+-----------------------+----------------------------+
#> | Exposure | Total # of covariates | # of imbalanced covariates |
#> +===========+=======================+============================+
@@ -3384,7 +3549,7 @@ Workflow: Continuous Exposure
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.24 | 47 | 0 |
#> +-----------+-----------------------+----------------------------+
-#> | ESETA1.35 | 59 | 2 |
+#> | ESETA1.35 | 59 | 0 |
#> +-----------+-----------------------+----------------------------+
#> | ESETA1.58 | 69 | 0 |
#> +-----------+-----------------------+----------------------------+
@@ -3457,7 +3622,7 @@ Workflow: Continuous Exposure
#> Please inspect the Wald test to determine if the exposures collectively predict significant variation in the outcome compared to a model without exposure terms.
#> We strongly suggest only conducting history comparisons if the test is significant.
#> test statistic df1 df2 dfcom p.value riv
-#> 1 ~~ 2 -1.181783 3 0.4300163 NA 1 2.021014
+#> 1 ~~ 2 -1.295578 3 0.4041257 NA 1 2.226662
#>
#>
#> The marginal model, m0, pooled across imputed datasets is summarized below:
@@ -3466,13 +3631,13 @@ Workflow: Continuous Exposure
#> +--------------------+-------+-----------------+--------+
#> | | Est. | CI | p |
#> +====================+=======+=================+========+
-#> | (Intercept) | 0.483 | [0.452, 0.513] | <0.001 |
+#> | (Intercept) | 0.482 | [0.448, 0.516] | <0.001 |
#> +--------------------+-------+-----------------+--------+
-#> | ESETA1.Childhood | 0.028 | [-0.112, 0.169] | 0.442 |
+#> | ESETA1.Childhood | 0.029 | [-0.104, 0.163] | 0.418 |
#> +--------------------+-------+-----------------+--------+
-#> | ESETA1.Infancy | 0.003 | [-0.042, 0.048] | 0.891 |
+#> | ESETA1.Infancy | 0.004 | [-0.036, 0.044] | 0.845 |
#> +--------------------+-------+-----------------+--------+
-#> | ESETA1.Toddlerhood | 0.012 | [-0.054, 0.079] | 0.650 |
+#> | ESETA1.Toddlerhood | 0.011 | [-0.047, 0.069] | 0.670 |
#> +--------------------+-------+-----------------+--------+
#> | Num.Obs. | 1292 | | |
#> +--------------------+-------+-----------------+--------+
@@ -3501,8 +3666,8 @@ Workflow: Continuous Exposure
save.out = save.out)
#> Please inspect the Wald test to determine if the exposures collectively predict significant variation in the outcome compared to a model without exposure terms.
#> We strongly suggest only conducting history comparisons if the test is significant.
-#> test statistic df1 df2 dfcom p.value riv
-#> 1 ~~ 2 -0.42805 3 0.6223825 NA 1 1.252613
+#> test statistic df1 df2 dfcom p.value riv
+#> 1 ~~ 2 -0.3020476 3 0.6594817 NA 1 1.174877
#>
#>
#> The marginal model, m0, pooled across imputed datasets is summarized below:
@@ -3511,13 +3676,13 @@ Workflow: Continuous Exposure
#> +--------------------+-------+-----------------+--------+
#> | | Est. | CI | p |
#> +====================+=======+=================+========+
-#> | (Intercept) | 0.481 | [0.450, 0.513] | <0.001 |
+#> | (Intercept) | 0.481 | [0.448, 0.514] | <0.001 |
#> +--------------------+-------+-----------------+--------+
-#> | ESETA1.Childhood | 0.028 | [-0.085, 0.141] | 0.400 |
+#> | ESETA1.Childhood | 0.029 | [-0.072, 0.130] | 0.357 |
#> +--------------------+-------+-----------------+--------+
-#> | ESETA1.Infancy | 0.006 | [-0.040, 0.051] | 0.784 |
+#> | ESETA1.Infancy | 0.006 | [-0.039, 0.051] | 0.782 |
#> +--------------------+-------+-----------------+--------+
-#> | ESETA1.Toddlerhood | 0.015 | [-0.061, 0.091] | 0.601 |
+#> | ESETA1.Toddlerhood | 0.014 | [-0.064, 0.092] | 0.628 |
#> +--------------------+-------+-----------------+--------+
#> | Num.Obs. | 1292 | | |
#> +--------------------+-------+-----------------+--------+
@@ -3543,26 +3708,26 @@ Workflow: Continuous Exposure
save.out = save.out)
#> Please inspect the Wald test to determine if the exposures collectively predict significant variation in the outcome compared to a model without exposure terms.
#> We strongly suggest only conducting history comparisons if the test is significant.
-#> test statistic df1 df2 dfcom p.value riv
-#> 1 ~~ 2 -1.653469 3 0.3504546 NA 1 2.86164
+#> test statistic df1 df2 dfcom p.value riv
+#> 1 ~~ 2 -1.758818 3 0.3310847 NA 1 3.208959
#>
#>
#> The marginal model, m0, pooled across imputed datasets is summarized below:
-#> +--------------------+--------+-----------------+--------+
-#> | | (1) |
-#> +--------------------+--------+-----------------+--------+
-#> | | Est. | CI | p |
-#> +====================+========+=================+========+
-#> | (Intercept) | 0.485 | [0.454, 0.515] | <0.001 |
-#> +--------------------+--------+-----------------+--------+
-#> | ESETA1.Childhood | 0.029 | [-0.150, 0.209] | 0.487 |
-#> +--------------------+--------+-----------------+--------+
-#> | ESETA1.Infancy | -0.002 | [-0.055, 0.052] | 0.934 |
-#> +--------------------+--------+-----------------+--------+
-#> | ESETA1.Toddlerhood | 0.009 | [-0.061, 0.080] | 0.748 |
-#> +--------------------+--------+-----------------+--------+
-#> | Num.Obs. | 1292 | | |
-#> +--------------------+--------+-----------------+--------+
+#> +--------------------+-------+-----------------+--------+
+#> | | (1) |
+#> +--------------------+-------+-----------------+--------+
+#> | | Est. | CI | p |
+#> +====================+=======+=================+========+
+#> | (Intercept) | 0.483 | [0.448, 0.518] | <0.001 |
+#> +--------------------+-------+-----------------+--------+
+#> | ESETA1.Childhood | 0.030 | [-0.142, 0.202] | 0.466 |
+#> +--------------------+-------+-----------------+--------+
+#> | ESETA1.Infancy | 0.000 | [-0.045, 0.044] | 0.982 |
+#> +--------------------+-------+-----------------+--------+
+#> | ESETA1.Toddlerhood | 0.007 | [-0.053, 0.067] | 0.778 |
+#> +--------------------+-------+-----------------+--------+
+#> | Num.Obs. | 1292 | | |
+#> +--------------------+-------+-----------------+--------+
With a comparable result.
@@ -3643,7 +3808,7 @@ Workflow: Continuous Exposure
#> +=======+==========+===========+==========+===========+
#> | l-l-l | 0.47 | 0.014 | 0.44 | 0.49 |
#> +-------+----------+-----------+----------+-----------+
-#> | h-h-h | 0.5 | 0.019 | 0.43 | 0.57 |
+#> | h-h-h | 0.5 | 0.02 | 0.42 | 0.57 |
#> +-------+----------+-----------+----------+-----------+
#>
#> Conducting multiple comparison correction for all pairings between comparison histories and each refernece history using the BH method.
@@ -3652,7 +3817,7 @@ Workflow: Continuous Exposure
#> +-------------------+------------+------------+-----------+-------------+------------+--------------+
#> | term | estimate | std.error | p.value | conf.low | conf.high | p.value_corr |
#> +===================+============+============+===========+=============+============+==============+
-#> | (h-h-h) - (l-l-l) | 0.02991272 | 0.01979075 | 0.2172254 | -0.02903923 | 0.08886468 | 0.2172254 |
+#> | (h-h-h) - (l-l-l) | 0.03035383 | 0.01948627 | 0.1991802 | -0.02526992 | 0.08597759 | 0.1991802 |
#> +-------------------+------------+------------+-----------+-------------+------------+--------------+
We first confirm the distribution of our sample across our exposure histories is reasonable.
@@ -3661,10 +3826,10 @@ Workflow: Continuous Exposure
summary(results,
save.out = save.out)
-#> term estimate std.error statistic df p.value
-#> 1 (h-h-h) - (l-l-l) 0.02991272 0.01979075 1.511449 3.40476 0.2172254
+#> term estimate std.error statistic df p.value
+#> 1 (h-h-h) - (l-l-l) 0.03035383 0.01948627 1.557704 3.739199 0.1991802
#> conf.low conf.high dose p.value_corr
-#> 1 -0.02903923 0.08886468 3 - 0 0.2172254
+#> 1 -0.02526992 0.08597759 3 - 0 0.1991802
We then inspect the history comparison and conclude that there is no evidence that children who experienced different histories of exposure to economic strain over infancy, toddlerhood, and early childhood differ in their behavioral problems in early childhood.
Lastly, we plot the results. The user has the option to supply save.out
to save plot output to the home directory.
@@ -3697,10 +3862,10 @@ Workflow: Continuous Exposure
save.out = save.out)
summary(results.s1,
save.out = save.out)
-#> term estimate std.error statistic df p.value
-#> 1 (h-h-h) - (l-l-l) 0.03275407 0.01735569 1.887224 5.210738 0.1154298
-#> conf.low conf.high dose p.value_corr
-#> 1 -0.01132286 0.07683101 3 - 0 0.1154298
+#> term estimate std.error statistic df p.value
+#> 1 (h-h-h) - (l-l-l) 0.03327812 0.01659957 2.004757 6.875701 0.08575804
+#> conf.low conf.high dose p.value_corr
+#> 1 -0.006117941 0.07267418 3 - 0 0.08575804
As shown above, results indicate a marginal but non-significant contrast between “l-l-l” and “h-h-h” histories of economic strain exposure in relation to behavior problems in early childhood.
@@ -3718,10 +3883,10 @@ Workflow: Continuous Exposure
save.out = save.out)
summary(results.s2,
save.out = save.out)
-#> term estimate std.error statistic df p.value
-#> 1 (h-h-h) - (l-l-l) 0.02588619 0.0234072 1.105907 2.462813 0.3649947
+#> term estimate std.error statistic df p.value
+#> 1 (h-h-h) - (l-l-l) 0.0264079 0.02274075 1.161259 2.655186 0.339228
#> conf.low conf.high dose p.value_corr
-#> 1 -0.05869734 0.1104697 3 - 0 0.3649947
+#> 1 -0.05158198 0.1043978 3 - 0 0.339228
Similarly, we find no evidence for differences in behavioral problems as function of history of exposure to economic strain.
@@ -3729,13 +3894,13 @@ Workflow: Continuous Exposure