Main model
Depression
h7_phq_results <- list()
h7_phq_models <- list()
for (h7_diff_var in h7_diff_vars) {
formula <- as.formula(paste0("phq_total ~ t1_", h7_diff_var, "_z * t1_age_z * time + (1 | ID)"))
model <- lmer(formula, data = full_data_long)
tidy_model <- tidy(model, conf.int = TRUE) %>%
filter(!str_detect(term, "Intercept|sd")) %>% # removes intercept terms
rename(
Variable = term,
B = estimate,
SE = std.error,
t = statistic,
p = p.value,
CI_Lower = conf.low,
CI_Upper = conf.high
) %>%
dplyr::select(Variable, B, SE, t, p, CI_Lower, CI_Upper) %>%
mutate(predictor = h7_diff_var)
h7_phq_results[[h7_diff_var]] <- tidy_model
h7_phq_models[[h7_diff_var]] <- model
h7_phq_results_df <- bind_rows(h7_phq_results) %>%
mutate(p.sig = ifelse(p < 0.05, "*", ""))
}
print(h7_phq_results_df)
## # A tibble: 66 × 9
## Variable B SE t p
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z 0.0440 0.184 0.240 8.11e- 1
## 2 t1_age_z -1.59 0.194 -8.20 6.42e-16
## 3 time2 -0.541 0.241 -2.24 2.53e- 2
## 4 time3 -0.655 0.262 -2.50 1.26e- 2
## 5 t1_f2f_freq_diff_family_z:t1_age_z -0.103 0.192 -0.535 5.93e- 1
## 6 t1_f2f_freq_diff_family_z:time2 0.948 0.262 3.62 3.21e- 4
## 7 t1_f2f_freq_diff_family_z:time3 0.630 0.291 2.16 3.10e- 2
## 8 t1_age_z:time2 -0.158 0.237 -0.666 5.06e- 1
## 9 t1_age_z:time3 -0.323 0.269 -1.20 2.31e- 1
## 10 t1_f2f_freq_diff_family_z:t1_age_z:time2 -0.123 0.238 -0.516 6.06e- 1
## 11 t1_f2f_freq_diff_family_z:t1_age_z:time3 -0.202 0.304 -0.664 5.07e- 1
## 12 t1_phone_freq_diff_family_z 0.601 0.185 3.25 1.18e- 3
## 13 t1_age_z -1.59 0.194 -8.20 6.27e-16
## 14 time2 -0.543 0.241 -2.25 2.46e- 2
## 15 time3 -0.625 0.262 -2.39 1.71e- 2
## 16 t1_phone_freq_diff_family_z:t1_age_z -0.136 0.190 -0.714 4.76e- 1
## 17 t1_phone_freq_diff_family_z:time2 -0.160 0.242 -0.664 5.07e- 1
## 18 t1_phone_freq_diff_family_z:time3 0.722 0.280 2.58 1.02e- 2
## 19 t1_age_z:time2 -0.242 0.237 -1.02 3.07e- 1
## 20 t1_age_z:time3 -0.322 0.262 -1.23 2.21e- 1
## 21 t1_phone_freq_diff_family_z:t1_age_z:time2 0.0201 0.230 0.0873 9.30e- 1
## 22 t1_phone_freq_diff_family_z:t1_age_z:time3 -0.801 0.311 -2.58 1.02e- 2
## 23 t1_elec_freq_diff_family_z 0.388 0.183 2.12 3.46e- 2
## 24 t1_age_z -1.58 0.194 -8.15 9.10e-16
## 25 time2 -0.558 0.243 -2.30 2.18e- 2
## 26 time3 -0.650 0.263 -2.47 1.37e- 2
## 27 t1_elec_freq_diff_family_z:t1_age_z -0.286 0.207 -1.38 1.67e- 1
## 28 t1_elec_freq_diff_family_z:time2 0.290 0.244 1.19 2.33e- 1
## 29 t1_elec_freq_diff_family_z:time3 -0.0580 0.323 -0.180 8.58e- 1
## 30 t1_age_z:time2 -0.206 0.238 -0.864 3.88e- 1
## 31 t1_age_z:time3 -0.263 0.265 -0.992 3.22e- 1
## 32 t1_elec_freq_diff_family_z:t1_age_z:time2 0.175 0.228 0.769 4.42e- 1
## 33 t1_elec_freq_diff_family_z:t1_age_z:time3 -0.425 0.319 -1.33 1.83e- 1
## 34 t1_f2f_freq_diff_friend_z 0.374 0.190 1.97 4.87e- 2
## 35 t1_age_z -1.82 0.188 -9.71 1.61e-21
## 36 time2 -0.656 0.241 -2.73 6.56e- 3
## 37 time3 -1.01 0.266 -3.81 1.52e- 4
## 38 t1_f2f_freq_diff_friend_z:t1_age_z -0.416 0.170 -2.44 1.48e- 2
## 39 t1_f2f_freq_diff_friend_z:time2 0.0473 0.244 0.194 8.47e- 1
## 40 t1_f2f_freq_diff_friend_z:time3 0.129 0.270 0.479 6.32e- 1
## 41 t1_age_z:time2 -0.0964 0.233 -0.413 6.80e- 1
## 42 t1_age_z:time3 -0.224 0.256 -0.872 3.83e- 1
## 43 t1_f2f_freq_diff_friend_z:t1_age_z:time2 -0.00479 0.213 -0.0225 9.82e- 1
## 44 t1_f2f_freq_diff_friend_z:t1_age_z:time3 -0.157 0.229 -0.684 4.94e- 1
## 45 t1_phone_freq_diff_friend_z 0.443 0.183 2.41 1.59e- 2
## 46 t1_age_z -1.84 0.184 -10.0 8.04e-23
## 47 time2 -0.576 0.236 -2.44 1.51e- 2
## 48 time3 -0.914 0.257 -3.55 4.14e- 4
## 49 t1_phone_freq_diff_friend_z:t1_age_z -0.368 0.185 -1.99 4.72e- 2
## 50 t1_phone_freq_diff_friend_z:time2 -0.0795 0.241 -0.330 7.42e- 1
## 51 t1_phone_freq_diff_friend_z:time3 -0.263 0.249 -1.06 2.91e- 1
## 52 t1_age_z:time2 -0.0708 0.230 -0.308 7.58e- 1
## 53 t1_age_z:time3 -0.290 0.248 -1.17 2.42e- 1
## 54 t1_phone_freq_diff_friend_z:t1_age_z:time2 0.382 0.225 1.70 8.95e- 2
## 55 t1_phone_freq_diff_friend_z:t1_age_z:time3 0.192 0.227 0.849 3.96e- 1
## 56 t1_elec_freq_diff_friend_z 0.507 0.183 2.77 5.72e- 3
## 57 t1_age_z -1.84 0.182 -10.1 4.85e-23
## 58 time2 -0.638 0.233 -2.74 6.37e- 3
## 59 time3 -0.921 0.257 -3.58 3.64e- 4
## 60 t1_elec_freq_diff_friend_z:t1_age_z -0.214 0.183 -1.17 2.42e- 1
## 61 t1_elec_freq_diff_friend_z:time2 0.105 0.228 0.458 6.47e- 1
## 62 t1_elec_freq_diff_friend_z:time3 -0.00760 0.266 -0.0285 9.77e- 1
## 63 t1_age_z:time2 -0.130 0.224 -0.578 5.63e- 1
## 64 t1_age_z:time3 -0.284 0.246 -1.15 2.49e- 1
## 65 t1_elec_freq_diff_friend_z:t1_age_z:time2 0.0499 0.250 0.199 8.42e- 1
## 66 t1_elec_freq_diff_friend_z:t1_age_z:time3 0.0977 0.290 0.337 7.36e- 1
## CI_Lower CI_Upper predictor p.sig
## <dbl> <dbl> <chr> <chr>
## 1 -0.316 0.404 f2f_freq_diff_family ""
## 2 -1.97 -1.21 f2f_freq_diff_family "*"
## 3 -1.01 -0.0671 f2f_freq_diff_family "*"
## 4 -1.17 -0.141 f2f_freq_diff_family "*"
## 5 -0.479 0.274 f2f_freq_diff_family ""
## 6 0.433 1.46 f2f_freq_diff_family "*"
## 7 0.0577 1.20 f2f_freq_diff_family "*"
## 8 -0.623 0.308 f2f_freq_diff_family ""
## 9 -0.852 0.206 f2f_freq_diff_family ""
## 10 -0.590 0.345 f2f_freq_diff_family ""
## 11 -0.799 0.395 f2f_freq_diff_family ""
## 12 0.239 0.964 phone_freq_diff_family "*"
## 13 -1.97 -1.21 phone_freq_diff_family "*"
## 14 -1.02 -0.0699 phone_freq_diff_family "*"
## 15 -1.14 -0.112 phone_freq_diff_family "*"
## 16 -0.509 0.237 phone_freq_diff_family ""
## 17 -0.635 0.314 phone_freq_diff_family ""
## 18 0.171 1.27 phone_freq_diff_family "*"
## 19 -0.707 0.223 phone_freq_diff_family ""
## 20 -0.837 0.194 phone_freq_diff_family ""
## 21 -0.432 0.473 phone_freq_diff_family ""
## 22 -1.41 -0.191 phone_freq_diff_family "*"
## 23 0.0281 0.747 elec_freq_diff_family "*"
## 24 -1.96 -1.20 elec_freq_diff_family "*"
## 25 -1.03 -0.0814 elec_freq_diff_family "*"
## 26 -1.17 -0.134 elec_freq_diff_family "*"
## 27 -0.692 0.120 elec_freq_diff_family ""
## 28 -0.188 0.769 elec_freq_diff_family ""
## 29 -0.692 0.576 elec_freq_diff_family ""
## 30 -0.673 0.262 elec_freq_diff_family ""
## 31 -0.785 0.258 elec_freq_diff_family ""
## 32 -0.273 0.624 elec_freq_diff_family ""
## 33 -1.05 0.201 elec_freq_diff_family ""
## 34 0.00222 0.746 f2f_freq_diff_friend "*"
## 35 -2.19 -1.45 f2f_freq_diff_friend "*"
## 36 -1.13 -0.184 f2f_freq_diff_friend "*"
## 37 -1.53 -0.491 f2f_freq_diff_friend "*"
## 38 -0.750 -0.0816 f2f_freq_diff_friend "*"
## 39 -0.433 0.527 f2f_freq_diff_friend ""
## 40 -0.400 0.659 f2f_freq_diff_friend ""
## 41 -0.555 0.362 f2f_freq_diff_friend ""
## 42 -0.727 0.280 f2f_freq_diff_friend ""
## 43 -0.422 0.413 f2f_freq_diff_friend ""
## 44 -0.607 0.293 f2f_freq_diff_friend ""
## 45 0.0829 0.802 phone_freq_diff_friend "*"
## 46 -2.21 -1.48 phone_freq_diff_friend "*"
## 47 -1.04 -0.112 phone_freq_diff_friend "*"
## 48 -1.42 -0.408 phone_freq_diff_friend "*"
## 49 -0.731 -0.00462 phone_freq_diff_friend "*"
## 50 -0.553 0.394 phone_freq_diff_friend ""
## 51 -0.751 0.225 phone_freq_diff_friend ""
## 52 -0.522 0.380 phone_freq_diff_friend ""
## 53 -0.777 0.197 phone_freq_diff_friend ""
## 54 -0.0591 0.824 phone_freq_diff_friend ""
## 55 -0.253 0.638 phone_freq_diff_friend ""
## 56 0.148 0.866 elec_freq_diff_friend "*"
## 57 -2.20 -1.48 elec_freq_diff_friend "*"
## 58 -1.10 -0.180 elec_freq_diff_friend "*"
## 59 -1.43 -0.416 elec_freq_diff_friend "*"
## 60 -0.574 0.145 elec_freq_diff_friend ""
## 61 -0.344 0.553 elec_freq_diff_friend ""
## 62 -0.530 0.515 elec_freq_diff_friend ""
## 63 -0.570 0.310 elec_freq_diff_friend ""
## 64 -0.767 0.199 elec_freq_diff_friend ""
## 65 -0.442 0.541 elec_freq_diff_friend ""
## 66 -0.472 0.667 elec_freq_diff_friend ""
Conduct FDR correction for multiple comparisons
h7_phq_fdr_results_df <- h7_phq_results_df %>%
group_by(predictor) %>%
mutate(
FDRp = p.adjust(p, method = "fdr"),
FDRp.sig = ifelse(FDRp < 0.05, "*", "")
) %>%
ungroup()
print(h7_phq_fdr_results_df)
## # A tibble: 66 × 11
## Variable B SE t p
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z 0.0440 0.184 0.240 8.11e- 1
## 2 t1_age_z -1.59 0.194 -8.20 6.42e-16
## 3 time2 -0.541 0.241 -2.24 2.53e- 2
## 4 time3 -0.655 0.262 -2.50 1.26e- 2
## 5 t1_f2f_freq_diff_family_z:t1_age_z -0.103 0.192 -0.535 5.93e- 1
## 6 t1_f2f_freq_diff_family_z:time2 0.948 0.262 3.62 3.21e- 4
## 7 t1_f2f_freq_diff_family_z:time3 0.630 0.291 2.16 3.10e- 2
## 8 t1_age_z:time2 -0.158 0.237 -0.666 5.06e- 1
## 9 t1_age_z:time3 -0.323 0.269 -1.20 2.31e- 1
## 10 t1_f2f_freq_diff_family_z:t1_age_z:time2 -0.123 0.238 -0.516 6.06e- 1
## 11 t1_f2f_freq_diff_family_z:t1_age_z:time3 -0.202 0.304 -0.664 5.07e- 1
## 12 t1_phone_freq_diff_family_z 0.601 0.185 3.25 1.18e- 3
## 13 t1_age_z -1.59 0.194 -8.20 6.27e-16
## 14 time2 -0.543 0.241 -2.25 2.46e- 2
## 15 time3 -0.625 0.262 -2.39 1.71e- 2
## 16 t1_phone_freq_diff_family_z:t1_age_z -0.136 0.190 -0.714 4.76e- 1
## 17 t1_phone_freq_diff_family_z:time2 -0.160 0.242 -0.664 5.07e- 1
## 18 t1_phone_freq_diff_family_z:time3 0.722 0.280 2.58 1.02e- 2
## 19 t1_age_z:time2 -0.242 0.237 -1.02 3.07e- 1
## 20 t1_age_z:time3 -0.322 0.262 -1.23 2.21e- 1
## 21 t1_phone_freq_diff_family_z:t1_age_z:time2 0.0201 0.230 0.0873 9.30e- 1
## 22 t1_phone_freq_diff_family_z:t1_age_z:time3 -0.801 0.311 -2.58 1.02e- 2
## 23 t1_elec_freq_diff_family_z 0.388 0.183 2.12 3.46e- 2
## 24 t1_age_z -1.58 0.194 -8.15 9.10e-16
## 25 time2 -0.558 0.243 -2.30 2.18e- 2
## 26 time3 -0.650 0.263 -2.47 1.37e- 2
## 27 t1_elec_freq_diff_family_z:t1_age_z -0.286 0.207 -1.38 1.67e- 1
## 28 t1_elec_freq_diff_family_z:time2 0.290 0.244 1.19 2.33e- 1
## 29 t1_elec_freq_diff_family_z:time3 -0.0580 0.323 -0.180 8.58e- 1
## 30 t1_age_z:time2 -0.206 0.238 -0.864 3.88e- 1
## 31 t1_age_z:time3 -0.263 0.265 -0.992 3.22e- 1
## 32 t1_elec_freq_diff_family_z:t1_age_z:time2 0.175 0.228 0.769 4.42e- 1
## 33 t1_elec_freq_diff_family_z:t1_age_z:time3 -0.425 0.319 -1.33 1.83e- 1
## 34 t1_f2f_freq_diff_friend_z 0.374 0.190 1.97 4.87e- 2
## 35 t1_age_z -1.82 0.188 -9.71 1.61e-21
## 36 time2 -0.656 0.241 -2.73 6.56e- 3
## 37 time3 -1.01 0.266 -3.81 1.52e- 4
## 38 t1_f2f_freq_diff_friend_z:t1_age_z -0.416 0.170 -2.44 1.48e- 2
## 39 t1_f2f_freq_diff_friend_z:time2 0.0473 0.244 0.194 8.47e- 1
## 40 t1_f2f_freq_diff_friend_z:time3 0.129 0.270 0.479 6.32e- 1
## 41 t1_age_z:time2 -0.0964 0.233 -0.413 6.80e- 1
## 42 t1_age_z:time3 -0.224 0.256 -0.872 3.83e- 1
## 43 t1_f2f_freq_diff_friend_z:t1_age_z:time2 -0.00479 0.213 -0.0225 9.82e- 1
## 44 t1_f2f_freq_diff_friend_z:t1_age_z:time3 -0.157 0.229 -0.684 4.94e- 1
## 45 t1_phone_freq_diff_friend_z 0.443 0.183 2.41 1.59e- 2
## 46 t1_age_z -1.84 0.184 -10.0 8.04e-23
## 47 time2 -0.576 0.236 -2.44 1.51e- 2
## 48 time3 -0.914 0.257 -3.55 4.14e- 4
## 49 t1_phone_freq_diff_friend_z:t1_age_z -0.368 0.185 -1.99 4.72e- 2
## 50 t1_phone_freq_diff_friend_z:time2 -0.0795 0.241 -0.330 7.42e- 1
## 51 t1_phone_freq_diff_friend_z:time3 -0.263 0.249 -1.06 2.91e- 1
## 52 t1_age_z:time2 -0.0708 0.230 -0.308 7.58e- 1
## 53 t1_age_z:time3 -0.290 0.248 -1.17 2.42e- 1
## 54 t1_phone_freq_diff_friend_z:t1_age_z:time2 0.382 0.225 1.70 8.95e- 2
## 55 t1_phone_freq_diff_friend_z:t1_age_z:time3 0.192 0.227 0.849 3.96e- 1
## 56 t1_elec_freq_diff_friend_z 0.507 0.183 2.77 5.72e- 3
## 57 t1_age_z -1.84 0.182 -10.1 4.85e-23
## 58 time2 -0.638 0.233 -2.74 6.37e- 3
## 59 time3 -0.921 0.257 -3.58 3.64e- 4
## 60 t1_elec_freq_diff_friend_z:t1_age_z -0.214 0.183 -1.17 2.42e- 1
## 61 t1_elec_freq_diff_friend_z:time2 0.105 0.228 0.458 6.47e- 1
## 62 t1_elec_freq_diff_friend_z:time3 -0.00760 0.266 -0.0285 9.77e- 1
## 63 t1_age_z:time2 -0.130 0.224 -0.578 5.63e- 1
## 64 t1_age_z:time3 -0.284 0.246 -1.15 2.49e- 1
## 65 t1_elec_freq_diff_friend_z:t1_age_z:time2 0.0499 0.250 0.199 8.42e- 1
## 66 t1_elec_freq_diff_friend_z:t1_age_z:time3 0.0977 0.290 0.337 7.36e- 1
## CI_Lower CI_Upper predictor p.sig FDRp FDRp.sig
## <dbl> <dbl> <chr> <chr> <dbl> <chr>
## 1 -0.316 0.404 f2f_freq_diff_family "" 8.11e- 1 ""
## 2 -1.97 -1.21 f2f_freq_diff_family "*" 7.06e-15 "*"
## 3 -1.01 -0.0671 f2f_freq_diff_family "*" 6.82e- 2 ""
## 4 -1.17 -0.141 f2f_freq_diff_family "*" 4.64e- 2 "*"
## 5 -0.479 0.274 f2f_freq_diff_family "" 6.67e- 1 ""
## 6 0.433 1.46 f2f_freq_diff_family "*" 1.76e- 3 "*"
## 7 0.0577 1.20 f2f_freq_diff_family "*" 6.82e- 2 ""
## 8 -0.623 0.308 f2f_freq_diff_family "" 6.67e- 1 ""
## 9 -0.852 0.206 f2f_freq_diff_family "" 4.23e- 1 ""
## 10 -0.590 0.345 f2f_freq_diff_family "" 6.67e- 1 ""
## 11 -0.799 0.395 f2f_freq_diff_family "" 6.67e- 1 ""
## 12 0.239 0.964 phone_freq_diff_family "*" 6.50e- 3 "*"
## 13 -1.97 -1.21 phone_freq_diff_family "*" 6.90e-15 "*"
## 14 -1.02 -0.0699 phone_freq_diff_family "*" 4.50e- 2 "*"
## 15 -1.14 -0.112 phone_freq_diff_family "*" 3.76e- 2 "*"
## 16 -0.509 0.237 phone_freq_diff_family "" 5.58e- 1 ""
## 17 -0.635 0.314 phone_freq_diff_family "" 5.58e- 1 ""
## 18 0.171 1.27 phone_freq_diff_family "*" 2.82e- 2 "*"
## 19 -0.707 0.223 phone_freq_diff_family "" 4.23e- 1 ""
## 20 -0.837 0.194 phone_freq_diff_family "" 3.47e- 1 ""
## 21 -0.432 0.473 phone_freq_diff_family "" 9.30e- 1 ""
## 22 -1.41 -0.191 phone_freq_diff_family "*" 2.82e- 2 "*"
## 23 0.0281 0.747 elec_freq_diff_family "*" 9.52e- 2 ""
## 24 -1.96 -1.20 elec_freq_diff_family "*" 1.00e-14 "*"
## 25 -1.03 -0.0814 elec_freq_diff_family "*" 8.01e- 2 ""
## 26 -1.17 -0.134 elec_freq_diff_family "*" 7.52e- 2 ""
## 27 -0.692 0.120 elec_freq_diff_family "" 3.35e- 1 ""
## 28 -0.188 0.769 elec_freq_diff_family "" 3.67e- 1 ""
## 29 -0.692 0.576 elec_freq_diff_family "" 8.58e- 1 ""
## 30 -0.673 0.262 elec_freq_diff_family "" 4.74e- 1 ""
## 31 -0.785 0.258 elec_freq_diff_family "" 4.42e- 1 ""
## 32 -0.273 0.624 elec_freq_diff_family "" 4.87e- 1 ""
## 33 -1.05 0.201 elec_freq_diff_family "" 3.35e- 1 ""
## 34 0.00222 0.746 f2f_freq_diff_friend "*" 1.07e- 1 ""
## 35 -2.19 -1.45 f2f_freq_diff_friend "*" 1.77e-20 "*"
## 36 -1.13 -0.184 f2f_freq_diff_friend "*" 2.41e- 2 "*"
## 37 -1.53 -0.491 f2f_freq_diff_friend "*" 8.38e- 4 "*"
## 38 -0.750 -0.0816 f2f_freq_diff_friend "*" 4.07e- 2 "*"
## 39 -0.433 0.527 f2f_freq_diff_friend "" 9.31e- 1 ""
## 40 -0.400 0.659 f2f_freq_diff_friend "" 8.31e- 1 ""
## 41 -0.555 0.362 f2f_freq_diff_friend "" 8.31e- 1 ""
## 42 -0.727 0.280 f2f_freq_diff_friend "" 7.03e- 1 ""
## 43 -0.422 0.413 f2f_freq_diff_friend "" 9.82e- 1 ""
## 44 -0.607 0.293 f2f_freq_diff_friend "" 7.76e- 1 ""
## 45 0.0829 0.802 phone_freq_diff_friend "*" 4.38e- 2 "*"
## 46 -2.21 -1.48 phone_freq_diff_friend "*" 8.85e-22 "*"
## 47 -1.04 -0.112 phone_freq_diff_friend "*" 4.38e- 2 "*"
## 48 -1.42 -0.408 phone_freq_diff_friend "*" 2.28e- 3 "*"
## 49 -0.731 -0.00462 phone_freq_diff_friend "*" 1.04e- 1 ""
## 50 -0.553 0.394 phone_freq_diff_friend "" 7.58e- 1 ""
## 51 -0.751 0.225 phone_freq_diff_friend "" 3.99e- 1 ""
## 52 -0.522 0.380 phone_freq_diff_friend "" 7.58e- 1 ""
## 53 -0.777 0.197 phone_freq_diff_friend "" 3.81e- 1 ""
## 54 -0.0591 0.824 phone_freq_diff_friend "" 1.64e- 1 ""
## 55 -0.253 0.638 phone_freq_diff_friend "" 4.85e- 1 ""
## 56 0.148 0.866 elec_freq_diff_friend "*" 1.75e- 2 "*"
## 57 -2.20 -1.48 elec_freq_diff_friend "*" 5.34e-22 "*"
## 58 -1.10 -0.180 elec_freq_diff_friend "*" 1.75e- 2 "*"
## 59 -1.43 -0.416 elec_freq_diff_friend "*" 2.00e- 3 "*"
## 60 -0.574 0.145 elec_freq_diff_friend "" 4.57e- 1 ""
## 61 -0.344 0.553 elec_freq_diff_friend "" 8.90e- 1 ""
## 62 -0.530 0.515 elec_freq_diff_friend "" 9.77e- 1 ""
## 63 -0.570 0.310 elec_freq_diff_friend "" 8.85e- 1 ""
## 64 -0.767 0.199 elec_freq_diff_friend "" 4.57e- 1 ""
## 65 -0.442 0.541 elec_freq_diff_friend "" 9.26e- 1 ""
## 66 -0.472 0.667 elec_freq_diff_friend "" 9.00e- 1 ""
Examine 3 way interactions
# 3 way interaction
h7_phq_3int_results_df <- h7_phq_fdr_results_df %>%
filter(str_detect(Variable, ":t1_age_z:"))
print(h7_phq_3int_results_df)
## # A tibble: 12 × 11
## Variable B SE t p
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z:t1_age_z:time2 -0.123 0.238 -0.516 0.606
## 2 t1_f2f_freq_diff_family_z:t1_age_z:time3 -0.202 0.304 -0.664 0.507
## 3 t1_phone_freq_diff_family_z:t1_age_z:time2 0.0201 0.230 0.0873 0.930
## 4 t1_phone_freq_diff_family_z:t1_age_z:time3 -0.801 0.311 -2.58 0.0102
## 5 t1_elec_freq_diff_family_z:t1_age_z:time2 0.175 0.228 0.769 0.442
## 6 t1_elec_freq_diff_family_z:t1_age_z:time3 -0.425 0.319 -1.33 0.183
## 7 t1_f2f_freq_diff_friend_z:t1_age_z:time2 -0.00479 0.213 -0.0225 0.982
## 8 t1_f2f_freq_diff_friend_z:t1_age_z:time3 -0.157 0.229 -0.684 0.494
## 9 t1_phone_freq_diff_friend_z:t1_age_z:time2 0.382 0.225 1.70 0.0895
## 10 t1_phone_freq_diff_friend_z:t1_age_z:time3 0.192 0.227 0.849 0.396
## 11 t1_elec_freq_diff_friend_z:t1_age_z:time2 0.0499 0.250 0.199 0.842
## 12 t1_elec_freq_diff_friend_z:t1_age_z:time3 0.0977 0.290 0.337 0.736
## CI_Lower CI_Upper predictor p.sig FDRp FDRp.sig
## <dbl> <dbl> <chr> <chr> <dbl> <chr>
## 1 -0.590 0.345 f2f_freq_diff_family "" 0.667 ""
## 2 -0.799 0.395 f2f_freq_diff_family "" 0.667 ""
## 3 -0.432 0.473 phone_freq_diff_family "" 0.930 ""
## 4 -1.41 -0.191 phone_freq_diff_family "*" 0.0282 "*"
## 5 -0.273 0.624 elec_freq_diff_family "" 0.487 ""
## 6 -1.05 0.201 elec_freq_diff_family "" 0.335 ""
## 7 -0.422 0.413 f2f_freq_diff_friend "" 0.982 ""
## 8 -0.607 0.293 f2f_freq_diff_friend "" 0.776 ""
## 9 -0.0591 0.824 phone_freq_diff_friend "" 0.164 ""
## 10 -0.253 0.638 phone_freq_diff_friend "" 0.485 ""
## 11 -0.442 0.541 elec_freq_diff_friend "" 0.926 ""
## 12 -0.472 0.667 elec_freq_diff_friend "" 0.900 ""
# significant 3 way interaction: t1_phone_freq_diff_family_z:t1_age_z:time3
# plot interaction
plot_h7_phq_phone_freq_diff_fam <- interact_plot(h7_phq_models[["phone_freq_diff_family"]],
pred ="t1_phone_freq_diff_family_z",
modx = "t1_age_z",
mod2 = "time",
x.label = "Frequency of Contact (Phone) Difference (Standardised)",
y.label = "PHQ-9 Total Score",
legend.main = "Age (Standardised)",
modx.labels = c("-1 SD", "Mean", "+1 SD"),
interval = TRUE)
plot_h7_phq_phone_freq_diff_fam

2 way interaction
h7_phq_2int_results_df <- h7_phq_fdr_results_df %>%
filter(str_detect(Variable, ":") & !str_detect(Variable, ":t1_age_z:"))
print(h7_phq_2int_results_df)
## # A tibble: 30 × 11
## Variable B SE t p CI_Lower
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z:t1_age_z -0.103 0.192 -0.535 0.593 -0.479
## 2 t1_f2f_freq_diff_family_z:time2 0.948 0.262 3.62 0.000321 0.433
## 3 t1_f2f_freq_diff_family_z:time3 0.630 0.291 2.16 0.0310 0.0577
## 4 t1_age_z:time2 -0.158 0.237 -0.666 0.506 -0.623
## 5 t1_age_z:time3 -0.323 0.269 -1.20 0.231 -0.852
## 6 t1_phone_freq_diff_family_z:t1_age_z -0.136 0.190 -0.714 0.476 -0.509
## 7 t1_phone_freq_diff_family_z:time2 -0.160 0.242 -0.664 0.507 -0.635
## 8 t1_phone_freq_diff_family_z:time3 0.722 0.280 2.58 0.0102 0.171
## 9 t1_age_z:time2 -0.242 0.237 -1.02 0.307 -0.707
## 10 t1_age_z:time3 -0.322 0.262 -1.23 0.221 -0.837
## 11 t1_elec_freq_diff_family_z:t1_age_z -0.286 0.207 -1.38 0.167 -0.692
## 12 t1_elec_freq_diff_family_z:time2 0.290 0.244 1.19 0.233 -0.188
## 13 t1_elec_freq_diff_family_z:time3 -0.0580 0.323 -0.180 0.858 -0.692
## 14 t1_age_z:time2 -0.206 0.238 -0.864 0.388 -0.673
## 15 t1_age_z:time3 -0.263 0.265 -0.992 0.322 -0.785
## 16 t1_f2f_freq_diff_friend_z:t1_age_z -0.416 0.170 -2.44 0.0148 -0.750
## 17 t1_f2f_freq_diff_friend_z:time2 0.0473 0.244 0.194 0.847 -0.433
## 18 t1_f2f_freq_diff_friend_z:time3 0.129 0.270 0.479 0.632 -0.400
## 19 t1_age_z:time2 -0.0964 0.233 -0.413 0.680 -0.555
## 20 t1_age_z:time3 -0.224 0.256 -0.872 0.383 -0.727
## 21 t1_phone_freq_diff_friend_z:t1_age_z -0.368 0.185 -1.99 0.0472 -0.731
## 22 t1_phone_freq_diff_friend_z:time2 -0.0795 0.241 -0.330 0.742 -0.553
## 23 t1_phone_freq_diff_friend_z:time3 -0.263 0.249 -1.06 0.291 -0.751
## 24 t1_age_z:time2 -0.0708 0.230 -0.308 0.758 -0.522
## 25 t1_age_z:time3 -0.290 0.248 -1.17 0.242 -0.777
## 26 t1_elec_freq_diff_friend_z:t1_age_z -0.214 0.183 -1.17 0.242 -0.574
## 27 t1_elec_freq_diff_friend_z:time2 0.105 0.228 0.458 0.647 -0.344
## 28 t1_elec_freq_diff_friend_z:time3 -0.00760 0.266 -0.0285 0.977 -0.530
## 29 t1_age_z:time2 -0.130 0.224 -0.578 0.563 -0.570
## 30 t1_age_z:time3 -0.284 0.246 -1.15 0.249 -0.767
## CI_Upper predictor p.sig FDRp FDRp.sig
## <dbl> <chr> <chr> <dbl> <chr>
## 1 0.274 f2f_freq_diff_family "" 0.667 ""
## 2 1.46 f2f_freq_diff_family "*" 0.00176 "*"
## 3 1.20 f2f_freq_diff_family "*" 0.0682 ""
## 4 0.308 f2f_freq_diff_family "" 0.667 ""
## 5 0.206 f2f_freq_diff_family "" 0.423 ""
## 6 0.237 phone_freq_diff_family "" 0.558 ""
## 7 0.314 phone_freq_diff_family "" 0.558 ""
## 8 1.27 phone_freq_diff_family "*" 0.0282 "*"
## 9 0.223 phone_freq_diff_family "" 0.423 ""
## 10 0.194 phone_freq_diff_family "" 0.347 ""
## 11 0.120 elec_freq_diff_family "" 0.335 ""
## 12 0.769 elec_freq_diff_family "" 0.367 ""
## 13 0.576 elec_freq_diff_family "" 0.858 ""
## 14 0.262 elec_freq_diff_family "" 0.474 ""
## 15 0.258 elec_freq_diff_family "" 0.442 ""
## 16 -0.0816 f2f_freq_diff_friend "*" 0.0407 "*"
## 17 0.527 f2f_freq_diff_friend "" 0.931 ""
## 18 0.659 f2f_freq_diff_friend "" 0.831 ""
## 19 0.362 f2f_freq_diff_friend "" 0.831 ""
## 20 0.280 f2f_freq_diff_friend "" 0.703 ""
## 21 -0.00462 phone_freq_diff_friend "*" 0.104 ""
## 22 0.394 phone_freq_diff_friend "" 0.758 ""
## 23 0.225 phone_freq_diff_friend "" 0.399 ""
## 24 0.380 phone_freq_diff_friend "" 0.758 ""
## 25 0.197 phone_freq_diff_friend "" 0.381 ""
## 26 0.145 elec_freq_diff_friend "" 0.457 ""
## 27 0.553 elec_freq_diff_friend "" 0.890 ""
## 28 0.515 elec_freq_diff_friend "" 0.977 ""
## 29 0.310 elec_freq_diff_friend "" 0.885 ""
## 30 0.199 elec_freq_diff_friend "" 0.457 ""
# significant 2 way interactions:
# t1_f2f_freq_diff_family_z:time2
# t1_phone_freq_diff_family_z:time3,
# t1_f2f_freq_diff_friend_z:t1_age_z
# plot interaction
plot_h7_phq_f2f_freq_diff_fam_time <- interact_plot(h7_phq_models[["f2f_freq_diff_family"]],
pred = "t1_f2f_freq_diff_family_z",
modx = "time",
x.label = "Frequency of Contact (F2F) Difference (Standardised)",
y.label = "PHQ-9 Total Score",
legend.main = "Wave",
modx.labels = c("T1", "T2", "T3"),
interval = TRUE)
plot_h7_phq_f2f_freq_diff_fam_time

plot_h7_phq_phone_freq_diff_fam_time <- interact_plot(h7_phq_models[["phone_freq_diff_family"]],
pred = "t1_phone_freq_diff_family_z",
modx = "time",
x.label = "Frequency of Contact (Phone) Difference (Standardised)",
y.label = "PHQ-9 Total Score",
legend.main = "Wave",
modx.labels = c("T1", "T2", "T3"),
interval = TRUE)
plot_h7_phq_phone_freq_diff_fam_time

plot_h7_phq_f2f_freq_diff_friend <- interact_plot(h7_phq_models[["f2f_freq_diff_friend"]],
pred = "t1_f2f_freq_diff_friend_z",
modx = "t1_age_z",
x.label = "Frequency of Contact (F2F) Difference (Standardised)",
y.label = "PHQ-9 Total Score",
legend.main = "Age (Standardised)",
modx.labels = c("-1 SD", "Mean", "+1 SD"),
interval = TRUE)
plot_h7_phq_f2f_freq_diff_friend

Loneliness
h7_ucla_results <- list()
h7_ucla_models <- list()
for (h7_diff_var in h7_diff_vars) {
formula <- as.formula(paste0("ucla_loneliness_total ~ t1_", h7_diff_var, "_z * t1_age_z * time + (1 | ID)"))
model <- lmer(formula, data = full_data_long)
tidy_model <- tidy(model, conf.int = TRUE) %>%
filter(!str_detect(term, "Intercept|sd")) %>% # removes intercept terms
rename(
Variable = term,
B = estimate,
SE = std.error,
t = statistic,
p = p.value,
CI_Lower = conf.low,
CI_Upper = conf.high
) %>%
dplyr::select(Variable, B, SE, t, p, CI_Lower, CI_Upper) %>%
mutate(predictor = h7_diff_var)
h7_ucla_results[[h7_diff_var]] <- tidy_model
h7_ucla_models[[h7_diff_var]] <- model
h7_ucla_results_df <- bind_rows(h7_ucla_results) %>%
mutate(p.sig = ifelse(p < 0.05, "*", ""))
}
h7_ucla_fdr_results_df <- h7_ucla_results_df %>%
group_by(predictor) %>%
mutate(
FDRp = p.adjust(p, method = "fdr"),
FDRp.sig = ifelse(FDRp < 0.05, "*", "")
) %>%
ungroup()
print(h7_ucla_fdr_results_df)
## # A tibble: 66 × 11
## Variable B SE t p
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z 0.192 0.472 0.406 6.85e- 1
## 2 t1_age_z -3.16 0.499 -6.32 3.68e-10
## 3 time2 0.320 0.654 0.490 6.24e- 1
## 4 time3 2.29 0.709 3.23 1.31e- 3
## 5 t1_f2f_freq_diff_family_z:t1_age_z -0.0838 0.494 -0.170 8.65e- 1
## 6 t1_f2f_freq_diff_family_z:time2 0.513 0.709 0.723 4.70e- 1
## 7 t1_f2f_freq_diff_family_z:time3 0.100 0.789 0.127 8.99e- 1
## 8 t1_age_z:time2 -0.288 0.642 -0.449 6.53e- 1
## 9 t1_age_z:time3 -0.380 0.730 -0.521 6.02e- 1
## 10 t1_f2f_freq_diff_family_z:t1_age_z:time2 0.364 0.645 0.565 5.72e- 1
## 11 t1_f2f_freq_diff_family_z:t1_age_z:time3 0.452 0.823 0.549 5.83e- 1
## 12 t1_phone_freq_diff_family_z 0.972 0.477 2.04 4.16e- 2
## 13 t1_age_z -3.15 0.499 -6.31 3.99e-10
## 14 time2 0.260 0.648 0.401 6.89e- 1
## 15 time3 2.24 0.703 3.19 1.52e- 3
## 16 t1_phone_freq_diff_family_z:t1_age_z -0.177 0.490 -0.360 7.19e- 1
## 17 t1_phone_freq_diff_family_z:time2 0.621 0.650 0.956 3.39e- 1
## 18 t1_phone_freq_diff_family_z:time3 2.27 0.754 3.01 2.71e- 3
## 19 t1_age_z:time2 -0.338 0.637 -0.530 5.96e- 1
## 20 t1_age_z:time3 -0.526 0.706 -0.746 4.56e- 1
## 21 t1_phone_freq_diff_family_z:t1_age_z:time2 -0.554 0.620 -0.894 3.72e- 1
## 22 t1_phone_freq_diff_family_z:t1_age_z:time3 -1.16 0.835 -1.39 1.64e- 1
## 23 t1_elec_freq_diff_family_z 0.340 0.472 0.721 4.71e- 1
## 24 t1_age_z -3.15 0.498 -6.32 3.77e-10
## 25 time2 0.183 0.647 0.282 7.78e- 1
## 26 time3 2.28 0.702 3.25 1.21e- 3
## 27 t1_elec_freq_diff_family_z:t1_age_z -0.470 0.532 -0.883 3.77e- 1
## 28 t1_elec_freq_diff_family_z:time2 1.70 0.650 2.62 9.06e- 3
## 29 t1_elec_freq_diff_family_z:time3 0.976 0.862 1.13 2.58e- 1
## 30 t1_age_z:time2 -0.285 0.635 -0.448 6.54e- 1
## 31 t1_age_z:time3 -0.351 0.708 -0.495 6.21e- 1
## 32 t1_elec_freq_diff_family_z:t1_age_z:time2 -0.883 0.609 -1.45 1.48e- 1
## 33 t1_elec_freq_diff_family_z:t1_age_z:time3 -1.13 0.851 -1.33 1.85e- 1
## 34 t1_f2f_freq_diff_friend_z 1.09 0.488 2.23 2.59e- 2
## 35 t1_age_z -3.48 0.482 -7.21 9.65e-13
## 36 time2 0.184 0.615 0.299 7.65e- 1
## 37 time3 1.89 0.679 2.78 5.65e- 3
## 38 t1_f2f_freq_diff_friend_z:t1_age_z -1.15 0.438 -2.63 8.72e- 3
## 39 t1_f2f_freq_diff_friend_z:time2 0.195 0.625 0.313 7.55e- 1
## 40 t1_f2f_freq_diff_friend_z:time3 0.263 0.689 0.382 7.02e- 1
## 41 t1_age_z:time2 -0.839 0.596 -1.41 1.60e- 1
## 42 t1_age_z:time3 -0.647 0.655 -0.988 3.24e- 1
## 43 t1_f2f_freq_diff_friend_z:t1_age_z:time2 -0.167 0.543 -0.307 7.59e- 1
## 44 t1_f2f_freq_diff_friend_z:t1_age_z:time3 0.00658 0.585 0.0112 9.91e- 1
## 45 t1_phone_freq_diff_friend_z 1.17 0.472 2.47 1.35e- 2
## 46 t1_age_z -3.51 0.474 -7.40 2.50e-13
## 47 time2 0.152 0.605 0.251 8.02e- 1
## 48 time3 1.97 0.659 3.00 2.84e- 3
## 49 t1_phone_freq_diff_friend_z:t1_age_z -0.0230 0.477 -0.0482 9.62e- 1
## 50 t1_phone_freq_diff_friend_z:time2 0.399 0.617 0.647 5.18e- 1
## 51 t1_phone_freq_diff_friend_z:time3 -0.121 0.636 -0.191 8.49e- 1
## 52 t1_age_z:time2 -0.916 0.588 -1.56 1.20e- 1
## 53 t1_age_z:time3 -0.768 0.634 -1.21 2.26e- 1
## 54 t1_phone_freq_diff_friend_z:t1_age_z:time2 -0.658 0.575 -1.14 2.53e- 1
## 55 t1_phone_freq_diff_friend_z:t1_age_z:time3 0.0474 0.580 0.0816 9.35e- 1
## 56 t1_elec_freq_diff_friend_z 1.30 0.472 2.76 5.80e- 3
## 57 t1_age_z -3.56 0.469 -7.59 6.34e-14
## 58 time2 0.305 0.594 0.513 6.08e- 1
## 59 time3 2.02 0.655 3.08 2.16e- 3
## 60 t1_elec_freq_diff_friend_z:t1_age_z -0.232 0.472 -0.492 6.23e- 1
## 61 t1_elec_freq_diff_friend_z:time2 0.100 0.581 0.172 8.63e- 1
## 62 t1_elec_freq_diff_friend_z:time3 0.318 0.678 0.468 6.40e- 1
## 63 t1_age_z:time2 -0.941 0.571 -1.65 9.98e- 2
## 64 t1_age_z:time3 -0.773 0.627 -1.23 2.18e- 1
## 65 t1_elec_freq_diff_friend_z:t1_age_z:time2 -0.803 0.638 -1.26 2.08e- 1
## 66 t1_elec_freq_diff_friend_z:t1_age_z:time3 0.348 0.739 0.471 6.38e- 1
## CI_Lower CI_Upper predictor p.sig FDRp FDRp.sig
## <dbl> <dbl> <chr> <chr> <dbl> <chr>
## 1 -0.735 1.12 f2f_freq_diff_family "" 8.37e- 1 ""
## 2 -4.14 -2.18 f2f_freq_diff_family "*" 4.05e- 9 "*"
## 3 -0.963 1.60 f2f_freq_diff_family "" 8.37e- 1 ""
## 4 0.897 3.68 f2f_freq_diff_family "*" 7.18e- 3 "*"
## 5 -1.05 0.885 f2f_freq_diff_family "" 8.99e- 1 ""
## 6 -0.880 1.91 f2f_freq_diff_family "" 8.37e- 1 ""
## 7 -1.45 1.65 f2f_freq_diff_family "" 8.99e- 1 ""
## 8 -1.55 0.972 f2f_freq_diff_family "" 8.37e- 1 ""
## 9 -1.81 1.05 f2f_freq_diff_family "" 8.37e- 1 ""
## 10 -0.902 1.63 f2f_freq_diff_family "" 8.37e- 1 ""
## 11 -1.16 2.07 f2f_freq_diff_family "" 8.37e- 1 ""
## 12 0.0371 1.91 phone_freq_diff_family "*" 1.14e- 1 ""
## 13 -4.13 -2.17 phone_freq_diff_family "*" 4.39e- 9 "*"
## 14 -1.01 1.53 phone_freq_diff_family "" 7.19e- 1 ""
## 15 0.859 3.62 phone_freq_diff_family "*" 8.35e- 3 "*"
## 16 -1.14 0.785 phone_freq_diff_family "" 7.19e- 1 ""
## 17 -0.655 1.90 phone_freq_diff_family "" 5.84e- 1 ""
## 18 0.789 3.75 phone_freq_diff_family "*" 9.95e- 3 "*"
## 19 -1.59 0.913 phone_freq_diff_family "" 7.19e- 1 ""
## 20 -1.91 0.860 phone_freq_diff_family "" 6.27e- 1 ""
## 21 -1.77 0.663 phone_freq_diff_family "" 5.84e- 1 ""
## 22 -2.81 0.476 phone_freq_diff_family "" 3.61e- 1 ""
## 23 -0.586 1.27 elec_freq_diff_family "" 6.48e- 1 ""
## 24 -4.13 -2.17 elec_freq_diff_family "*" 4.14e- 9 "*"
## 25 -1.09 1.45 elec_freq_diff_family "" 7.78e- 1 ""
## 26 0.904 3.66 elec_freq_diff_family "*" 6.64e- 3 "*"
## 27 -1.51 0.574 elec_freq_diff_family "" 5.93e- 1 ""
## 28 0.425 2.98 elec_freq_diff_family "*" 3.32e- 2 "*"
## 29 -0.716 2.67 elec_freq_diff_family "" 4.72e- 1 ""
## 30 -1.53 0.963 elec_freq_diff_family "" 7.20e- 1 ""
## 31 -1.74 1.04 elec_freq_diff_family "" 7.20e- 1 ""
## 32 -2.08 0.313 elec_freq_diff_family "" 4.06e- 1 ""
## 33 -2.80 0.542 elec_freq_diff_family "" 4.07e- 1 ""
## 34 0.131 2.04 f2f_freq_diff_friend "*" 7.12e- 2 ""
## 35 -4.43 -2.53 f2f_freq_diff_friend "*" 1.06e-11 "*"
## 36 -1.02 1.39 f2f_freq_diff_friend "" 8.42e- 1 ""
## 37 0.552 3.22 f2f_freq_diff_friend "*" 3.11e- 2 "*"
## 38 -2.01 -0.292 f2f_freq_diff_friend "*" 3.20e- 2 "*"
## 39 -1.03 1.42 f2f_freq_diff_friend "" 8.42e- 1 ""
## 40 -1.09 1.62 f2f_freq_diff_friend "" 8.42e- 1 ""
## 41 -2.01 0.332 f2f_freq_diff_friend "" 3.52e- 1 ""
## 42 -1.93 0.639 f2f_freq_diff_friend "" 5.93e- 1 ""
## 43 -1.23 0.900 f2f_freq_diff_friend "" 8.42e- 1 ""
## 44 -1.14 1.16 f2f_freq_diff_friend "" 9.91e- 1 ""
## 45 0.241 2.10 phone_freq_diff_friend "*" 4.96e- 2 "*"
## 46 -4.44 -2.58 phone_freq_diff_friend "*" 2.75e-12 "*"
## 47 -1.04 1.34 phone_freq_diff_friend "" 9.62e- 1 ""
## 48 0.680 3.27 phone_freq_diff_friend "*" 1.56e- 2 "*"
## 49 -0.959 0.913 phone_freq_diff_friend "" 9.62e- 1 ""
## 50 -0.812 1.61 phone_freq_diff_friend "" 8.14e- 1 ""
## 51 -1.37 1.13 phone_freq_diff_friend "" 9.62e- 1 ""
## 52 -2.07 0.238 phone_freq_diff_friend "" 3.29e- 1 ""
## 53 -2.01 0.477 phone_freq_diff_friend "" 4.64e- 1 ""
## 54 -1.79 0.472 phone_freq_diff_friend "" 4.64e- 1 ""
## 55 -1.09 1.19 phone_freq_diff_friend "" 9.62e- 1 ""
## 56 0.378 2.23 elec_freq_diff_friend "*" 2.13e- 2 "*"
## 57 -4.48 -2.64 elec_freq_diff_friend "*" 6.97e-13 "*"
## 58 -0.862 1.47 elec_freq_diff_friend "" 7.04e- 1 ""
## 59 0.731 3.30 elec_freq_diff_friend "*" 1.19e- 2 "*"
## 60 -1.16 0.694 elec_freq_diff_friend "" 7.04e- 1 ""
## 61 -1.04 1.24 elec_freq_diff_friend "" 8.63e- 1 ""
## 62 -1.01 1.65 elec_freq_diff_friend "" 7.04e- 1 ""
## 63 -2.06 0.180 elec_freq_diff_friend "" 2.74e- 1 ""
## 64 -2.00 0.459 elec_freq_diff_friend "" 4.00e- 1 ""
## 65 -2.06 0.449 elec_freq_diff_friend "" 4.00e- 1 ""
## 66 -1.10 1.80 elec_freq_diff_friend "" 7.04e- 1 ""
3-way interaction
h7_ucla_3int_results_df <- h7_ucla_fdr_results_df %>%
filter(str_detect(Variable, ":t1_age_z:"))
print(h7_ucla_3int_results_df) # no significant 3-way interactions
## # A tibble: 12 × 11
## Variable B SE t p
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z:t1_age_z:time2 0.364 0.645 0.565 0.572
## 2 t1_f2f_freq_diff_family_z:t1_age_z:time3 0.452 0.823 0.549 0.583
## 3 t1_phone_freq_diff_family_z:t1_age_z:time2 -0.554 0.620 -0.894 0.372
## 4 t1_phone_freq_diff_family_z:t1_age_z:time3 -1.16 0.835 -1.39 0.164
## 5 t1_elec_freq_diff_family_z:t1_age_z:time2 -0.883 0.609 -1.45 0.148
## 6 t1_elec_freq_diff_family_z:t1_age_z:time3 -1.13 0.851 -1.33 0.185
## 7 t1_f2f_freq_diff_friend_z:t1_age_z:time2 -0.167 0.543 -0.307 0.759
## 8 t1_f2f_freq_diff_friend_z:t1_age_z:time3 0.00658 0.585 0.0112 0.991
## 9 t1_phone_freq_diff_friend_z:t1_age_z:time2 -0.658 0.575 -1.14 0.253
## 10 t1_phone_freq_diff_friend_z:t1_age_z:time3 0.0474 0.580 0.0816 0.935
## 11 t1_elec_freq_diff_friend_z:t1_age_z:time2 -0.803 0.638 -1.26 0.208
## 12 t1_elec_freq_diff_friend_z:t1_age_z:time3 0.348 0.739 0.471 0.638
## CI_Lower CI_Upper predictor p.sig FDRp FDRp.sig
## <dbl> <dbl> <chr> <chr> <dbl> <chr>
## 1 -0.902 1.63 f2f_freq_diff_family "" 0.837 ""
## 2 -1.16 2.07 f2f_freq_diff_family "" 0.837 ""
## 3 -1.77 0.663 phone_freq_diff_family "" 0.584 ""
## 4 -2.81 0.476 phone_freq_diff_family "" 0.361 ""
## 5 -2.08 0.313 elec_freq_diff_family "" 0.406 ""
## 6 -2.80 0.542 elec_freq_diff_family "" 0.407 ""
## 7 -1.23 0.900 f2f_freq_diff_friend "" 0.842 ""
## 8 -1.14 1.16 f2f_freq_diff_friend "" 0.991 ""
## 9 -1.79 0.472 phone_freq_diff_friend "" 0.464 ""
## 10 -1.09 1.19 phone_freq_diff_friend "" 0.962 ""
## 11 -2.06 0.449 elec_freq_diff_friend "" 0.400 ""
## 12 -1.10 1.80 elec_freq_diff_friend "" 0.704 ""
2-way interaction
h7_ucla_2int_results_df <- h7_ucla_fdr_results_df %>%
filter(str_detect(Variable, ":") & !str_detect(Variable, ":t1_age_z:"))
# sig 2-way interaction:
# t1_phone_freq_diff_family_z:time3,
# t1_elec_freq_diff_family_z:time2
# t1_f2f_freq_diff_friend_z:t1_age_z
# plot interaction
plot_h7_ucla_phone_freq_diff_fam <- interact_plot(h7_ucla_models[["phone_freq_diff_family"]],
pred = "t1_phone_freq_diff_family_z",
modx = "time",
x.label = "Frequency of Contact (Phone) Difference (Standardised)",
y.label = "UCLA Loneliness Total Score",
legend.main = "Wave",
modx.labels = c("T1", "T2", "T3"),
interval = TRUE)
plot_h7_ucla_phone_freq_diff_fam

plot_h7_ucla_elec_freq_diff_fam <- interact_plot(h7_ucla_models[["elec_freq_diff_family"]],
pred = "t1_elec_freq_diff_family_z",
modx = "time",
x.label = "Frequency of Contact (Electronic) Difference (Standardised)",
y.label = "UCLA Loneliness Total Score",
legend.main = "Wave",
modx.labels = c("T1", "T2", "T3"),
interval = TRUE)
plot_h7_ucla_elec_freq_diff_fam

plot_h7_ucla_f2f_freq_diff_friend <- interact_plot(h7_ucla_models[["f2f_freq_diff_friend"]],
pred = "t1_f2f_freq_diff_friend_z",
modx = "t1_age_z",
x.label = "Frequency of Contact (F2F) Difference (Standardised)",
y.label = "UCLA Loneliness Total Score",
legend.main = "Age (Standardised)",
modx.labels = c("-1 SD", "Mean", "+1 SD"),
interval = TRUE)
plot_h7_ucla_f2f_freq_diff_friend

Sensitivity analyses
Depression
h7_phq_senresults <- list()
h7_phq_senmodels <- list()
for (h7_diff_var in h7_diff_vars) {
formula <- as.formula(paste0("phq_total ~ t1_", h7_diff_var, "_z * t1_age_z * time + country + (1 | ID)"))
model <- lmer(formula, data = full_data_long)
tidy_model <- tidy(model, conf.int = TRUE) %>%
filter(!str_detect(term, "Intercept|sd")) %>% # removes intercept terms
rename(
Variable = term,
B = estimate,
SE = std.error,
t = statistic,
p = p.value,
CI_Lower = conf.low,
CI_Upper = conf.high
) %>%
dplyr::select(Variable, B, SE, t, p, CI_Lower, CI_Upper) %>%
mutate(predictor = h7_diff_var)
h7_phq_senresults[[h7_diff_var]] <- tidy_model
h7_phq_senmodels[[h7_diff_var]] <- model
h7_phq_senresults_df <- bind_rows(h7_phq_senresults) %>%
mutate(p.sig = ifelse(p < 0.05, "*", ""))
}
h7_phq_fdr_senresults_df <- h7_phq_senresults_df %>%
group_by(predictor) %>%
mutate(
FDRp = p.adjust(p, method = "fdr"),
FDRp.sig = ifelse(FDRp < 0.05, "*", "")
) %>%
ungroup() #%>%
#filter(str_detect(Variable, ":"))
print(h7_phq_fdr_senresults_df)
## # A tibble: 78 × 11
## Variable B SE t p
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z -0.0292 0.185 -0.158 8.74e- 1
## 2 t1_age_z -1.58 0.193 -8.20 6.61e-16
## 3 time2 -0.525 0.242 -2.17 3.01e- 2
## 4 time3 -0.647 0.262 -2.47 1.38e- 2
## 5 countryuk 1.44 0.436 3.31 9.70e- 4
## 6 countryusa 0.789 0.465 1.70 9.00e- 2
## 7 t1_f2f_freq_diff_family_z:t1_age_z -0.0932 0.191 -0.487 6.26e- 1
## 8 t1_f2f_freq_diff_family_z:time2 0.947 0.262 3.61 3.26e- 4
## 9 t1_f2f_freq_diff_family_z:time3 0.633 0.291 2.17 3.03e- 2
## 10 t1_age_z:time2 -0.163 0.237 -0.688 4.92e- 1
## 11 t1_age_z:time3 -0.315 0.269 -1.17 2.43e- 1
## 12 t1_f2f_freq_diff_family_z:t1_age_z:time2 -0.115 0.238 -0.482 6.30e- 1
## 13 t1_f2f_freq_diff_family_z:t1_age_z:time3 -0.180 0.304 -0.593 5.53e- 1
## 14 t1_phone_freq_diff_family_z 0.566 0.185 3.07 2.21e- 3
## 15 t1_age_z -1.58 0.193 -8.21 5.90e-16
## 16 time2 -0.526 0.241 -2.18 2.95e- 2
## 17 time3 -0.617 0.262 -2.36 1.87e- 2
## 18 countryuk 1.39 0.433 3.21 1.37e- 3
## 19 countryusa 0.693 0.464 1.49 1.36e- 1
## 20 t1_phone_freq_diff_family_z:t1_age_z -0.146 0.189 -0.770 4.42e- 1
## 21 t1_phone_freq_diff_family_z:time2 -0.168 0.242 -0.696 4.87e- 1
## 22 t1_phone_freq_diff_family_z:time3 0.715 0.281 2.55 1.10e- 2
## 23 t1_age_z:time2 -0.245 0.237 -1.04 3.01e- 1
## 24 t1_age_z:time3 -0.315 0.263 -1.20 2.31e- 1
## 25 t1_phone_freq_diff_family_z:t1_age_z:time2 0.0307 0.231 0.133 8.94e- 1
## 26 t1_phone_freq_diff_family_z:t1_age_z:time3 -0.794 0.311 -2.55 1.09e- 2
## 27 t1_elec_freq_diff_family_z 0.372 0.183 2.04 4.17e- 2
## 28 t1_age_z -1.57 0.193 -8.16 8.53e-16
## 29 time2 -0.540 0.243 -2.22 2.65e- 2
## 30 time3 -0.642 0.263 -2.44 1.49e- 2
## 31 countryuk 1.44 0.432 3.33 8.88e- 4
## 32 countryusa 0.710 0.464 1.53 1.26e- 1
## 33 t1_elec_freq_diff_family_z:t1_age_z -0.280 0.206 -1.36 1.74e- 1
## 34 t1_elec_freq_diff_family_z:time2 0.281 0.244 1.15 2.50e- 1
## 35 t1_elec_freq_diff_family_z:time3 -0.0672 0.323 -0.208 8.35e- 1
## 36 t1_age_z:time2 -0.210 0.238 -0.881 3.79e- 1
## 37 t1_age_z:time3 -0.259 0.266 -0.977 3.29e- 1
## 38 t1_elec_freq_diff_family_z:t1_age_z:time2 0.187 0.228 0.820 4.13e- 1
## 39 t1_elec_freq_diff_family_z:t1_age_z:time3 -0.403 0.319 -1.26 2.07e- 1
## 40 t1_f2f_freq_diff_friend_z 0.297 0.190 1.56 1.19e- 1
## 41 t1_age_z -1.82 0.187 -9.73 1.37e-21
## 42 time2 -0.634 0.241 -2.63 8.66e- 3
## 43 time3 -1.00 0.266 -3.76 1.82e- 4
## 44 countryuk 1.49 0.428 3.47 5.41e- 4
## 45 countryusa 0.679 0.469 1.45 1.48e- 1
## 46 t1_f2f_freq_diff_friend_z:t1_age_z -0.429 0.170 -2.53 1.16e- 2
## 47 t1_f2f_freq_diff_friend_z:time2 0.0518 0.245 0.212 8.32e- 1
## 48 t1_f2f_freq_diff_friend_z:time3 0.133 0.270 0.492 6.23e- 1
## 49 t1_age_z:time2 -0.0980 0.234 -0.420 6.75e- 1
## 50 t1_age_z:time3 -0.220 0.256 -0.859 3.91e- 1
## 51 t1_f2f_freq_diff_friend_z:t1_age_z:time2 -0.00817 0.213 -0.0384 9.69e- 1
## 52 t1_f2f_freq_diff_friend_z:t1_age_z:time3 -0.157 0.229 -0.684 4.94e- 1
## 53 t1_phone_freq_diff_friend_z 0.450 0.183 2.46 1.39e- 2
## 54 t1_age_z -1.82 0.183 -9.91 2.53e-22
## 55 time2 -0.555 0.237 -2.35 1.93e- 2
## 56 time3 -0.906 0.258 -3.52 4.67e- 4
## 57 countryuk 1.57 0.426 3.68 2.40e- 4
## 58 countryusa 0.782 0.468 1.67 9.51e- 2
## 59 t1_phone_freq_diff_friend_z:t1_age_z -0.335 0.184 -1.82 6.94e- 2
## 60 t1_phone_freq_diff_friend_z:time2 -0.0730 0.241 -0.303 7.62e- 1
## 61 t1_phone_freq_diff_friend_z:time3 -0.257 0.249 -1.03 3.03e- 1
## 62 t1_age_z:time2 -0.0722 0.230 -0.314 7.53e- 1
## 63 t1_age_z:time3 -0.287 0.248 -1.16 2.48e- 1
## 64 t1_phone_freq_diff_friend_z:t1_age_z:time2 0.377 0.225 1.67 9.44e- 2
## 65 t1_phone_freq_diff_friend_z:t1_age_z:time3 0.187 0.227 0.823 4.11e- 1
## 66 t1_elec_freq_diff_friend_z 0.513 0.182 2.81 4.97e- 3
## 67 t1_age_z -1.81 0.182 -9.99 1.27e-22
## 68 time2 -0.616 0.234 -2.64 8.54e- 3
## 69 time3 -0.912 0.257 -3.54 4.24e- 4
## 70 countryuk 1.59 0.425 3.75 1.90e- 4
## 71 countryusa 0.761 0.468 1.63 1.04e- 1
## 72 t1_elec_freq_diff_friend_z:t1_age_z -0.199 0.182 -1.09 2.75e- 1
## 73 t1_elec_freq_diff_friend_z:time2 0.106 0.228 0.466 6.41e- 1
## 74 t1_elec_freq_diff_friend_z:time3 0.00273 0.266 0.0102 9.92e- 1
## 75 t1_age_z:time2 -0.132 0.224 -0.591 5.55e- 1
## 76 t1_age_z:time3 -0.281 0.246 -1.14 2.54e- 1
## 77 t1_elec_freq_diff_friend_z:t1_age_z:time2 0.0518 0.250 0.207 8.36e- 1
## 78 t1_elec_freq_diff_friend_z:t1_age_z:time3 0.113 0.290 0.388 6.98e- 1
## CI_Lower CI_Upper predictor p.sig FDRp FDRp.sig
## <dbl> <dbl> <chr> <chr> <dbl> <chr>
## 1 -0.392 0.334 f2f_freq_diff_family "" 8.74e- 1 ""
## 2 -1.96 -1.20 f2f_freq_diff_family "*" 8.59e-15 "*"
## 3 -1.00 -0.0507 f2f_freq_diff_family "*" 6.56e- 2 ""
## 4 -1.16 -0.132 f2f_freq_diff_family "*" 4.50e- 2 "*"
## 5 0.587 2.30 f2f_freq_diff_family "*" 4.20e- 3 "*"
## 6 -0.123 1.70 f2f_freq_diff_family "" 1.67e- 1 ""
## 7 -0.469 0.282 f2f_freq_diff_family "" 6.83e- 1 ""
## 8 0.433 1.46 f2f_freq_diff_family "*" 2.12e- 3 "*"
## 9 0.0605 1.21 f2f_freq_diff_family "*" 6.56e- 2 ""
## 10 -0.629 0.302 f2f_freq_diff_family "" 6.83e- 1 ""
## 11 -0.844 0.214 f2f_freq_diff_family "" 3.95e- 1 ""
## 12 -0.582 0.353 f2f_freq_diff_family "" 6.83e- 1 ""
## 13 -0.778 0.417 f2f_freq_diff_family "" 6.83e- 1 ""
## 14 0.204 0.928 phone_freq_diff_family "*" 9.58e- 3 "*"
## 15 -1.96 -1.21 phone_freq_diff_family "*" 7.68e-15 "*"
## 16 -1.00 -0.0525 phone_freq_diff_family "*" 5.48e- 2 ""
## 17 -1.13 -0.103 phone_freq_diff_family "*" 4.05e- 2 "*"
## 18 0.540 2.24 phone_freq_diff_family "*" 8.92e- 3 "*"
## 19 -0.218 1.60 phone_freq_diff_family "" 2.21e- 1 ""
## 20 -0.517 0.226 phone_freq_diff_family "" 5.22e- 1 ""
## 21 -0.643 0.307 phone_freq_diff_family "" 5.27e- 1 ""
## 22 0.164 1.27 phone_freq_diff_family "*" 2.87e- 2 "*"
## 23 -0.710 0.220 phone_freq_diff_family "" 3.91e- 1 ""
## 24 -0.831 0.201 phone_freq_diff_family "" 3.33e- 1 ""
## 25 -0.422 0.484 phone_freq_diff_family "" 8.94e- 1 ""
## 26 -1.40 -0.183 phone_freq_diff_family "*" 2.87e- 2 "*"
## 27 0.0141 0.730 elec_freq_diff_family "*" 1.08e- 1 ""
## 28 -1.95 -1.20 elec_freq_diff_family "*" 1.11e-14 "*"
## 29 -1.02 -0.0634 elec_freq_diff_family "*" 8.60e- 2 ""
## 30 -1.16 -0.125 elec_freq_diff_family "*" 6.48e- 2 ""
## 31 0.593 2.29 elec_freq_diff_family "*" 5.77e- 3 "*"
## 32 -0.200 1.62 elec_freq_diff_family "" 2.74e- 1 ""
## 33 -0.685 0.124 elec_freq_diff_family "" 3.22e- 1 ""
## 34 -0.198 0.759 elec_freq_diff_family "" 3.61e- 1 ""
## 35 -0.702 0.568 elec_freq_diff_family "" 8.35e- 1 ""
## 36 -0.678 0.258 elec_freq_diff_family "" 4.47e- 1 ""
## 37 -0.781 0.262 elec_freq_diff_family "" 4.28e- 1 ""
## 38 -0.261 0.636 elec_freq_diff_family "" 4.47e- 1 ""
## 39 -1.03 0.224 elec_freq_diff_family "" 3.37e- 1 ""
## 40 -0.0762 0.670 f2f_freq_diff_friend "" 2.57e- 1 ""
## 41 -2.19 -1.45 f2f_freq_diff_friend "*" 1.78e-20 "*"
## 42 -1.11 -0.161 f2f_freq_diff_friend "*" 2.81e- 2 "*"
## 43 -1.52 -0.479 f2f_freq_diff_friend "*" 1.18e- 3 "*"
## 44 0.646 2.33 f2f_freq_diff_friend "*" 2.35e- 3 "*"
## 45 -0.242 1.60 f2f_freq_diff_friend "" 2.75e- 1 ""
## 46 -0.763 -0.0960 f2f_freq_diff_friend "*" 3.03e- 2 "*"
## 47 -0.429 0.532 f2f_freq_diff_friend "" 9.02e- 1 ""
## 48 -0.397 0.663 f2f_freq_diff_friend "" 7.98e- 1 ""
## 49 -0.557 0.361 f2f_freq_diff_friend "" 7.98e- 1 ""
## 50 -0.724 0.283 f2f_freq_diff_friend "" 6.35e- 1 ""
## 51 -0.426 0.409 f2f_freq_diff_friend "" 9.69e- 1 ""
## 52 -0.607 0.293 f2f_freq_diff_friend "" 7.14e- 1 ""
## 53 0.0915 0.808 phone_freq_diff_friend "*" 4.52e- 2 "*"
## 54 -2.18 -1.46 phone_freq_diff_friend "*" 3.28e-21 "*"
## 55 -1.02 -0.0903 phone_freq_diff_friend "*" 5.02e- 2 ""
## 56 -1.41 -0.400 phone_freq_diff_friend "*" 2.02e- 3 "*"
## 57 0.735 2.41 phone_freq_diff_friend "*" 1.56e- 3 "*"
## 58 -0.136 1.70 phone_freq_diff_friend "" 1.55e- 1 ""
## 59 -0.697 0.0266 phone_freq_diff_friend "" 1.50e- 1 ""
## 60 -0.547 0.401 phone_freq_diff_friend "" 7.62e- 1 ""
## 61 -0.745 0.232 phone_freq_diff_friend "" 3.94e- 1 ""
## 62 -0.523 0.379 phone_freq_diff_friend "" 7.62e- 1 ""
## 63 -0.774 0.200 phone_freq_diff_friend "" 3.59e- 1 ""
## 64 -0.0650 0.819 phone_freq_diff_friend "" 1.55e- 1 ""
## 65 -0.259 0.632 phone_freq_diff_friend "" 4.86e- 1 ""
## 66 0.155 0.871 elec_freq_diff_friend "*" 1.61e- 2 "*"
## 67 -2.17 -1.46 elec_freq_diff_friend "*" 1.65e-21 "*"
## 68 -1.07 -0.157 elec_freq_diff_friend "*" 2.22e- 2 "*"
## 69 -1.42 -0.406 elec_freq_diff_friend "*" 1.84e- 3 "*"
## 70 0.759 2.43 elec_freq_diff_friend "*" 1.23e- 3 "*"
## 71 -0.157 1.68 elec_freq_diff_friend "" 2.26e- 1 ""
## 72 -0.557 0.158 elec_freq_diff_friend "" 4.46e- 1 ""
## 73 -0.342 0.555 elec_freq_diff_friend "" 8.25e- 1 ""
## 74 -0.520 0.526 elec_freq_diff_friend "" 9.92e- 1 ""
## 75 -0.573 0.308 elec_freq_diff_friend "" 8.02e- 1 ""
## 76 -0.765 0.202 elec_freq_diff_friend "" 4.46e- 1 ""
## 77 -0.440 0.544 elec_freq_diff_friend "" 9.06e- 1 ""
## 78 -0.457 0.683 elec_freq_diff_friend "" 8.25e- 1 ""
Loneliness
h7_ucla_senresults <- list()
h7_ucla_senmodels <- list()
for (h7_diff_var in h7_diff_vars) {
formula <- as.formula(paste0("ucla_loneliness_total ~ t1_", h7_diff_var, "_z * t1_age_z * time + country + (1 | ID)"))
model <- lmer(formula, data = full_data_long)
tidy_model <- tidy(model, conf.int = TRUE) %>%
filter(!str_detect(term, "Intercept|sd")) %>% # removes intercept terms
rename(
Variable = term,
B = estimate,
SE = std.error,
t = statistic,
p = p.value,
CI_Lower = conf.low,
CI_Upper = conf.high
) %>%
dplyr::select(Variable, B, SE, t, p, CI_Lower, CI_Upper) %>%
mutate(predictor = h7_diff_var)
h7_ucla_senresults[[h7_diff_var]] <- tidy_model
h7_ucla_senmodels[[h7_diff_var]] <- model
h7_ucla_senresults_df <- bind_rows(h7_ucla_senresults) %>%
mutate(p.sig = ifelse(p < 0.05, "*", ""))
}
h7_ucla_fdr_senresults_df <- h7_ucla_senresults_df %>%
group_by(predictor) %>%
mutate(
FDRp = p.adjust(p, method = "fdr"),
FDRp.sig = ifelse(FDRp < 0.05, "*", "")
) %>%
ungroup() #%>%
#filter(str_detect(Variable, ":"))
print(h7_ucla_fdr_senresults_df)
## # A tibble: 78 × 11
## Variable B SE t p
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 t1_f2f_freq_diff_family_z 0.107 0.477 0.225 8.22e- 1
## 2 t1_age_z -3.13 0.498 -6.29 4.38e-10
## 3 time2 0.345 0.654 0.528 5.98e- 1
## 4 time3 2.30 0.710 3.24 1.26e- 3
## 5 countryuk 2.91 1.12 2.60 9.35e- 3
## 6 countryusa 2.79 1.19 2.34 1.95e- 2
## 7 t1_f2f_freq_diff_family_z:t1_age_z -0.109 0.494 -0.221 8.25e- 1
## 8 t1_f2f_freq_diff_family_z:time2 0.506 0.710 0.712 4.76e- 1
## 9 t1_f2f_freq_diff_family_z:time3 0.111 0.789 0.140 8.88e- 1
## 10 t1_age_z:time2 -0.295 0.642 -0.459 6.47e- 1
## 11 t1_age_z:time3 -0.359 0.730 -0.492 6.23e- 1
## 12 t1_f2f_freq_diff_family_z:t1_age_z:time2 0.392 0.645 0.608 5.44e- 1
## 13 t1_f2f_freq_diff_family_z:t1_age_z:time3 0.500 0.824 0.607 5.44e- 1
## 14 t1_phone_freq_diff_family_z 0.887 0.477 1.86 6.31e- 2
## 15 t1_age_z -3.13 0.498 -6.28 4.64e-10
## 16 time2 0.285 0.648 0.439 6.61e- 1
## 17 time3 2.25 0.704 3.19 1.47e- 3
## 18 countryuk 2.73 1.11 2.45 1.43e- 2
## 19 countryusa 2.60 1.19 2.18 2.98e- 2
## 20 t1_phone_freq_diff_family_z:t1_age_z -0.202 0.489 -0.414 6.79e- 1
## 21 t1_phone_freq_diff_family_z:time2 0.607 0.650 0.933 3.51e- 1
## 22 t1_phone_freq_diff_family_z:time3 2.26 0.754 2.99 2.87e- 3
## 23 t1_age_z:time2 -0.342 0.637 -0.537 5.92e- 1
## 24 t1_age_z:time3 -0.510 0.706 -0.722 4.71e- 1
## 25 t1_phone_freq_diff_family_z:t1_age_z:time2 -0.541 0.620 -0.872 3.83e- 1
## 26 t1_phone_freq_diff_family_z:t1_age_z:time3 -1.16 0.836 -1.38 1.67e- 1
## 27 t1_elec_freq_diff_family_z 0.303 0.471 0.644 5.20e- 1
## 28 t1_age_z -3.13 0.497 -6.29 4.47e-10
## 29 time2 0.209 0.648 0.323 7.47e- 1
## 30 time3 2.29 0.702 3.26 1.17e- 3
## 31 countryuk 2.80 1.11 2.52 1.20e- 2
## 32 countryusa 2.59 1.19 2.17 3.02e- 2
## 33 t1_elec_freq_diff_family_z:t1_age_z -0.434 0.531 -0.817 4.14e- 1
## 34 t1_elec_freq_diff_family_z:time2 1.67 0.650 2.57 1.03e- 2
## 35 t1_elec_freq_diff_family_z:time3 0.958 0.862 1.11 2.67e- 1
## 36 t1_age_z:time2 -0.288 0.636 -0.454 6.50e- 1
## 37 t1_age_z:time3 -0.341 0.709 -0.480 6.31e- 1
## 38 t1_elec_freq_diff_family_z:t1_age_z:time2 -0.854 0.610 -1.40 1.62e- 1
## 39 t1_elec_freq_diff_family_z:t1_age_z:time3 -1.08 0.852 -1.27 2.06e- 1
## 40 t1_f2f_freq_diff_friend_z 0.933 0.490 1.90 5.70e- 2
## 41 t1_age_z -3.44 0.482 -7.14 1.56e-12
## 42 time2 0.211 0.615 0.343 7.32e- 1
## 43 time3 1.89 0.679 2.78 5.59e- 3
## 44 countryuk 2.88 1.11 2.60 9.36e- 3
## 45 countryusa 2.46 1.21 2.03 4.22e- 2
## 46 t1_f2f_freq_diff_friend_z:t1_age_z -1.15 0.438 -2.62 8.97e- 3
## 47 t1_f2f_freq_diff_friend_z:time2 0.196 0.625 0.313 7.54e- 1
## 48 t1_f2f_freq_diff_friend_z:time3 0.270 0.689 0.392 6.95e- 1
## 49 t1_age_z:time2 -0.836 0.597 -1.40 1.62e- 1
## 50 t1_age_z:time3 -0.628 0.655 -0.959 3.38e- 1
## 51 t1_f2f_freq_diff_friend_z:t1_age_z:time2 -0.168 0.543 -0.310 7.57e- 1
## 52 t1_f2f_freq_diff_friend_z:t1_age_z:time3 0.00342 0.586 0.00583 9.95e- 1
## 53 t1_phone_freq_diff_friend_z 1.15 0.471 2.45 1.45e- 2
## 54 t1_age_z -3.42 0.473 -7.23 8.37e-13
## 55 time2 0.180 0.606 0.298 7.66e- 1
## 56 time3 1.97 0.659 2.99 2.86e- 3
## 57 countryuk 3.22 1.10 2.93 3.50e- 3
## 58 countryusa 2.77 1.21 2.29 2.25e- 2
## 59 t1_phone_freq_diff_friend_z:t1_age_z 0.0317 0.476 0.0667 9.47e- 1
## 60 t1_phone_freq_diff_friend_z:time2 0.420 0.617 0.680 4.97e- 1
## 61 t1_phone_freq_diff_friend_z:time3 -0.105 0.636 -0.165 8.69e- 1
## 62 t1_age_z:time2 -0.907 0.588 -1.54 1.23e- 1
## 63 t1_age_z:time3 -0.747 0.635 -1.18 2.39e- 1
## 64 t1_phone_freq_diff_friend_z:t1_age_z:time2 -0.657 0.576 -1.14 2.54e- 1
## 65 t1_phone_freq_diff_friend_z:t1_age_z:time3 0.0466 0.581 0.0803 9.36e- 1
## 66 t1_elec_freq_diff_friend_z 1.28 0.471 2.73 6.47e- 3
## 67 t1_age_z -3.48 0.469 -7.43 2.05e-13
## 68 time2 0.331 0.595 0.557 5.78e- 1
## 69 time3 2.02 0.655 3.08 2.16e- 3
## 70 countryuk 3.20 1.10 2.91 3.67e- 3
## 71 countryusa 2.71 1.21 2.24 2.51e- 2
## 72 t1_elec_freq_diff_friend_z:t1_age_z -0.210 0.471 -0.445 6.56e- 1
## 73 t1_elec_freq_diff_friend_z:time2 0.106 0.581 0.182 8.55e- 1
## 74 t1_elec_freq_diff_friend_z:time3 0.338 0.678 0.498 6.19e- 1
## 75 t1_age_z:time2 -0.936 0.571 -1.64 1.01e- 1
## 76 t1_age_z:time3 -0.753 0.627 -1.20 2.30e- 1
## 77 t1_elec_freq_diff_friend_z:t1_age_z:time2 -0.789 0.638 -1.24 2.16e- 1
## 78 t1_elec_freq_diff_friend_z:t1_age_z:time3 0.382 0.739 0.517 6.06e- 1
## CI_Lower CI_Upper predictor p.sig FDRp FDRp.sig
## <dbl> <dbl> <chr> <chr> <dbl> <chr>
## 1 -0.829 1.04 f2f_freq_diff_family "" 8.88e- 1 ""
## 2 -4.11 -2.16 f2f_freq_diff_family "*" 5.69e- 9 "*"
## 3 -0.940 1.63 f2f_freq_diff_family "" 8.41e- 1 ""
## 4 0.905 3.69 f2f_freq_diff_family "*" 8.22e- 3 "*"
## 5 0.718 5.11 f2f_freq_diff_family "*" 4.05e- 2 "*"
## 6 0.450 5.13 f2f_freq_diff_family "*" 6.34e- 2 ""
## 7 -1.08 0.859 f2f_freq_diff_family "" 8.88e- 1 ""
## 8 -0.888 1.90 f2f_freq_diff_family "" 8.41e- 1 ""
## 9 -1.44 1.66 f2f_freq_diff_family "" 8.88e- 1 ""
## 10 -1.56 0.967 f2f_freq_diff_family "" 8.41e- 1 ""
## 11 -1.79 1.07 f2f_freq_diff_family "" 8.41e- 1 ""
## 12 -0.875 1.66 f2f_freq_diff_family "" 8.41e- 1 ""
## 13 -1.12 2.12 f2f_freq_diff_family "" 8.41e- 1 ""
## 14 -0.0485 1.82 phone_freq_diff_family "" 1.37e- 1 ""
## 15 -4.11 -2.15 phone_freq_diff_family "*" 6.04e- 9 "*"
## 16 -0.989 1.56 phone_freq_diff_family "" 6.79e- 1 ""
## 17 0.867 3.63 phone_freq_diff_family "*" 9.55e- 3 "*"
## 18 0.547 4.92 phone_freq_diff_family "*" 4.66e- 2 "*"
## 19 0.254 4.94 phone_freq_diff_family "*" 7.76e- 2 ""
## 20 -1.16 0.757 phone_freq_diff_family "" 6.79e- 1 ""
## 21 -0.670 1.88 phone_freq_diff_family "" 5.54e- 1 ""
## 22 0.776 3.74 phone_freq_diff_family "*" 1.25e- 2 "*"
## 23 -1.59 0.909 phone_freq_diff_family "" 6.79e- 1 ""
## 24 -1.90 0.877 phone_freq_diff_family "" 6.12e- 1 ""
## 25 -1.76 0.677 phone_freq_diff_family "" 5.54e- 1 ""
## 26 -2.80 0.485 phone_freq_diff_family "" 3.10e- 1 ""
## 27 -0.621 1.23 elec_freq_diff_family "" 6.76e- 1 ""
## 28 -4.10 -2.15 elec_freq_diff_family "*" 5.81e- 9 "*"
## 29 -1.06 1.48 elec_freq_diff_family "" 7.47e- 1 ""
## 30 0.910 3.67 elec_freq_diff_family "*" 7.63e- 3 "*"
## 31 0.617 4.98 elec_freq_diff_family "*" 3.89e- 2 "*"
## 32 0.248 4.93 elec_freq_diff_family "*" 7.86e- 2 ""
## 33 -1.48 0.608 elec_freq_diff_family "" 5.98e- 1 ""
## 34 0.397 2.95 elec_freq_diff_family "*" 3.89e- 2 "*"
## 35 -0.735 2.65 elec_freq_diff_family "" 4.34e- 1 ""
## 36 -1.54 0.960 elec_freq_diff_family "" 7.04e- 1 ""
## 37 -1.73 1.05 elec_freq_diff_family "" 7.04e- 1 ""
## 38 -2.05 0.344 elec_freq_diff_family "" 3.51e- 1 ""
## 39 -2.75 0.594 elec_freq_diff_family "" 3.82e- 1 ""
## 40 -0.0280 1.89 f2f_freq_diff_friend "" 1.24e- 1 ""
## 41 -4.39 -2.50 f2f_freq_diff_friend "*" 2.03e-11 "*"
## 42 -0.997 1.42 f2f_freq_diff_friend "" 8.20e- 1 ""
## 43 0.555 3.22 f2f_freq_diff_friend "*" 3.04e- 2 "*"
## 44 0.709 5.05 f2f_freq_diff_friend "*" 3.04e- 2 "*"
## 45 0.0872 4.84 f2f_freq_diff_friend "*" 1.10e- 1 ""
## 46 -2.01 -0.287 f2f_freq_diff_friend "*" 3.04e- 2 "*"
## 47 -1.03 1.42 f2f_freq_diff_friend "" 8.20e- 1 ""
## 48 -1.08 1.62 f2f_freq_diff_friend "" 8.20e- 1 ""
## 49 -2.01 0.336 f2f_freq_diff_friend "" 3.01e- 1 ""
## 50 -1.91 0.658 f2f_freq_diff_friend "" 5.49e- 1 ""
## 51 -1.24 0.899 f2f_freq_diff_friend "" 8.20e- 1 ""
## 52 -1.15 1.15 f2f_freq_diff_friend "" 9.95e- 1 ""
## 53 0.230 2.08 phone_freq_diff_friend "*" 4.70e- 2 "*"
## 54 -4.35 -2.49 phone_freq_diff_friend "*" 1.09e-11 "*"
## 55 -1.01 1.37 phone_freq_diff_friend "" 9.47e- 1 ""
## 56 0.679 3.27 phone_freq_diff_friend "*" 1.52e- 2 "*"
## 57 1.06 5.39 phone_freq_diff_friend "*" 1.52e- 2 "*"
## 58 0.392 5.14 phone_freq_diff_friend "*" 5.84e- 2 ""
## 59 -0.903 0.966 phone_freq_diff_friend "" 9.47e- 1 ""
## 60 -0.792 1.63 phone_freq_diff_friend "" 7.17e- 1 ""
## 61 -1.35 1.14 phone_freq_diff_friend "" 9.47e- 1 ""
## 62 -2.06 0.247 phone_freq_diff_friend "" 2.67e- 1 ""
## 63 -1.99 0.499 phone_freq_diff_friend "" 4.13e- 1 ""
## 64 -1.79 0.474 phone_freq_diff_friend "" 4.13e- 1 ""
## 65 -1.09 1.19 phone_freq_diff_friend "" 9.47e- 1 ""
## 66 0.361 2.21 elec_freq_diff_friend "*" 2.10e- 2 "*"
## 67 -4.40 -2.56 elec_freq_diff_friend "*" 2.66e-12 "*"
## 68 -0.836 1.50 elec_freq_diff_friend "" 7.11e- 1 ""
## 69 0.731 3.30 elec_freq_diff_friend "*" 1.40e- 2 "*"
## 70 1.04 5.36 elec_freq_diff_friend "*" 1.59e- 2 "*"
## 71 0.339 5.09 elec_freq_diff_friend "*" 6.53e- 2 ""
## 72 -1.13 0.714 elec_freq_diff_friend "" 7.11e- 1 ""
## 73 -1.04 1.25 elec_freq_diff_friend "" 8.55e- 1 ""
## 74 -0.994 1.67 elec_freq_diff_friend "" 7.11e- 1 ""
## 75 -2.06 0.185 elec_freq_diff_friend "" 2.20e- 1 ""
## 76 -1.99 0.478 elec_freq_diff_friend "" 3.74e- 1 ""
## 77 -2.04 0.463 elec_freq_diff_friend "" 3.74e- 1 ""
## 78 -1.07 1.83 elec_freq_diff_friend "" 7.11e- 1 ""