Average of DVs

# Average pairs for IN1,IN2 
climate_labels$Inj <- rowMeans(climate_labels[, c("IN1", "IN2")], na.rm = TRUE)

# Reverse coding for Harm Attribution 
climate_labels$HA2_r <- NA
climate_labels$HA2_r <- 7 - climate_labels$HA2  
climate_labels$HA4_r <- NA
climate_labels$HA4_r <- 7 - climate_labels$HA4  

# Average pairs for HA 
climate_labels$HA_12 <- rowMeans(climate_labels[, c("HA1", "HA2_r")], na.rm = TRUE)
climate_labels$HA_34 <- rowMeans(climate_labels[, c("HA3", "HA4_r")], na.rm = TRUE)


# Reverse coding for Efficacy Beliefs 
climate_labels$EB2_r <- NA
climate_labels$EB2_r <- 6 - climate_labels$EB2  
climate_labels$EB4_r <- NA
climate_labels$EB4_r <- 6 - climate_labels$EB4  
climate_labels$EB6_r <- NA
climate_labels$EB6_r <- 6 - climate_labels$EB6

# Average pairs for EB
climate_labels$EB_12 <- rowMeans(climate_labels[, c("EB1", "EB2_r")], na.rm = TRUE)
climate_labels$EB_34 <- rowMeans(climate_labels[, c("EB3", "EB4_r")], na.rm = TRUE)
climate_labels$EB_56 <- rowMeans(climate_labels[, c("EB5", "EB6_r")], na.rm = TRUE)


# Reverse coding for Systems Responsibility 
climate_labels$SR2_r <- NA 
climate_labels$SR2_r <- 7 - climate_labels$SR2
climate_labels$SR4_r <- NA 
climate_labels$SR4_r <- 7 - climate_labels$SR4

# Average pairs for SR 
climate_labels$SR_12 <- rowMeans(climate_labels[, c("SR1", "SR2_r")], na.rm = TRUE)
climate_labels$SR_34 <- rowMeans(climate_labels[, c("SR3", "SR4_r")], na.rm = TRUE)


# Reverse coding for Integral Solutions Mindset 
climate_labels$ISM2_r <- NA 
climate_labels$ISM2_r <- 6 - climate_labels$ISM2
climate_labels$ISM4_r <- NA 
climate_labels$ISM4_r <- 6 - climate_labels$ISM4
climate_labels$ISM6_r <- NA
climate_labels$ISM6_r <- 6 - climate_labels$ISM6

# Average pairs for ISM
climate_labels$ISM_12 <- rowMeans(climate_labels[, c("ISM1", "ISM2_r")], na.rm = TRUE)
climate_labels$ISM_34 <- rowMeans(climate_labels[, c("ISM3", "ISM4_r")], na.rm = TRUE)
climate_labels$ISM_56 <- rowMeans(climate_labels[, c("ISM5", "ISM6_r")], na.rm = TRUE)

Descriptive

variables <- c("Inj", "DN1", "DN2", "DescN", 
               "HA_12", "HA_34", 
               "EB_12", "EB_34", "EB_56",
               "SR_12", "SR_34",
               "ISM_12", "ISM_34", "ISM_56")

summary_table <- climate_labels %>%
  group_by(Condition) %>%
  summarise(across(all_of(variables), list(mean = ~ mean(.x, na.rm = TRUE),
                                           sd = ~ sd(.x, na.rm = TRUE)))) %>%
  pivot_longer(-Condition, names_to = c("Variable", ".value"), 
               names_pattern = "(.*)_(mean|sd)") %>%
  arrange(Variable, Condition)

summary_table <- summary_table %>%
  mutate(Mean_SD = paste0(round(mean, 2), " (", round(sd, 2), ")")) %>%
  select(Variable, Condition, Mean_SD)

final_table <- summary_table %>%
  pivot_wider(names_from = Condition, values_from = Mean_SD)

final_table <- final_table[, colSums(is.na(final_table)) < nrow(final_table)]
final_table <- final_table %>% select(-`NA`)

final_table %>%
  kbl(caption = "Means and Standard Deviation by Condition", 
      format = "html", 
      align = "c", 
      digits = 2) %>%
  kable_styling(full_width = FALSE, 
                position = "center", 
                font_size = 12) %>%
  row_spec(0, bold = TRUE)  
Means and Standard Deviation by Condition
Variable Direct With Societal Change Direct With Support Direct With Support With Societal Change Direct Without Support Indirect With Societal Change Indirect With Support Indirect With Support With Societal Change Indirect Without Support No Label
DN1 4.51 (1.26) 4.49 (1.28) 4.55 (1.2) 4.49 (1.23) 4.55 (1.23) 4.58 (1.24) 4.54 (1.28) 4.5 (1.15) 4.47 (1.23)
DN2 3.4 (1.41) 3.44 (1.42) 3.44 (1.42) 3.42 (1.32) 3.47 (1.46) 3.37 (1.42) 3.43 (1.42) 3.37 (1.46) 3.3 (1.42)
DescN 43.43 (23.82) 44.16 (24.18) 44.8 (23.16) 41.3 (23.68) 43.74 (24.64) 46.83 (24.43) 43.67 (24.45) 39.91 (25.22) 33.98 (24.84)
EB_12 2.95 (1.05) 2.76 (1.05) 2.82 (1.04) 2.87 (1.01) 2.86 (1.08) 2.86 (1.11) 2.82 (1.04) 2.83 (1.06) 2.87 (1.06)
EB_34 3.18 (0.93) 3.06 (0.96) 3.07 (0.94) 3.08 (0.83) 3.17 (0.95) 3.15 (0.97) 3.1 (0.96) 3.1 (0.93) 3.16 (0.96)
EB_56 2.61 (1.08) 2.58 (1.12) 2.5 (1.07) 2.6 (1.08) 2.6 (1.09) 2.62 (1.07) 2.59 (1.09) 2.64 (1.08) 2.63 (1.04)
HA_12 4.26 (1.22) 4.16 (1.32) 4.08 (1.29) 4.26 (1.22) 4.19 (1.27) 4.22 (1.19) 4.18 (1.26) 4.18 (1.26) 4.03 (1.23)
HA_34 4.31 (1.3) 4.22 (1.36) 4.19 (1.33) 4.3 (1.29) 4.27 (1.32) 4.27 (1.26) 4.2 (1.31) 4.22 (1.3) 4.13 (1.34)
ISM_12 3.36 (1.04) 3.21 (1.09) 3.13 (1.1) 3.29 (1.09) 3.17 (1.1) 3.27 (1.08) 3.18 (1.09) 3.29 (1.08) 3.22 (1.13)
ISM_34 3.65 (0.89) 3.53 (0.91) 3.51 (0.91) 3.51 (0.92) 3.47 (0.94) 3.63 (0.89) 3.51 (0.94) 3.55 (0.89) 3.53 (0.91)
ISM_56 3.58 (1.1) 3.49 (1.1) 3.37 (1.13) 3.48 (1.07) 3.4 (1.13) 3.52 (1.05) 3.46 (1.14) 3.51 (1.1) 3.5 (1.09)
Inj 2.71 (0.99) 2.89 (1.04) 2.72 (0.95) 2.73 (0.99) 2.67 (0.98) 2.87 (1) 2.69 (0.99) 2.63 (0.99) 2.56 (0.96)
SR_12 3.78 (1.32) 3.71 (1.35) 3.66 (1.37) 3.8 (1.33) 3.7 (1.34) 3.8 (1.35) 3.75 (1.32) 3.8 (1.41) 3.69 (1.25)
SR_34 3.8 (1.3) 3.73 (1.32) 3.75 (1.3) 3.81 (1.38) 3.67 (1.3) 3.88 (1.32) 3.78 (1.26) 3.9 (1.37) 3.76 (1.27)

Correlation

dvs <- climate_labels[, c("Inj", "DN1", "DN2", "DescN", 
                          "HA_12", "HA_34", 
                          "EB_12", "EB_34", "EB_56",
                          "SR_12", "SR_34",
                          "ISM_12", "ISM_34", "ISM_56")]

cor_matrix <- cor(dvs, use = "complete.obs")
print(cor_matrix)
##              Inj       DN1       DN2     DescN     HA_12     HA_34     EB_12
## Inj    1.0000000 0.3034393 0.5356997 0.5816107 0.3680490 0.3528971 0.4383451
## DN1    0.3034393 1.0000000 0.2575024 0.2305618 0.2127825 0.2201628 0.1683195
## DN2    0.5356997 0.2575024 1.0000000 0.4013786 0.3388615 0.3574229 0.4381992
## DescN  0.5816107 0.2305618 0.4013786 1.0000000 0.3446094 0.3223164 0.3272923
## HA_12  0.3680490 0.2127825 0.3388615 0.3446094 1.0000000 0.8512745 0.4426437
## HA_34  0.3528971 0.2201628 0.3574229 0.3223164 0.8512745 1.0000000 0.4609926
## EB_12  0.4383451 0.1683195 0.4381992 0.3272923 0.4426437 0.4609926 1.0000000
## EB_34  0.3878705 0.2050196 0.4131477 0.2944261 0.4690087 0.4902185 0.6148327
## EB_56  0.4728845 0.1696008 0.4230804 0.3442003 0.3832823 0.3808996 0.6487695
## SR_12  0.2760276 0.1266947 0.2842176 0.2179726 0.4517847 0.4561996 0.5102184
## SR_34  0.2673525 0.1547755 0.2739165 0.2291372 0.4472835 0.4582385 0.5145910
## ISM_12 0.4216142 0.2122344 0.3904008 0.3291507 0.6110485 0.6217566 0.5688091
## ISM_34 0.2215406 0.1582244 0.1720420 0.1890209 0.3831671 0.3929884 0.3377211
## ISM_56 0.4012560 0.2295021 0.4021979 0.3278693 0.6431643 0.6624370 0.5922310
##            EB_34     EB_56     SR_12     SR_34    ISM_12    ISM_34    ISM_56
## Inj    0.3878705 0.4728845 0.2760276 0.2673525 0.4216142 0.2215406 0.4012560
## DN1    0.2050196 0.1696008 0.1266947 0.1547755 0.2122344 0.1582244 0.2295021
## DN2    0.4131477 0.4230804 0.2842176 0.2739165 0.3904008 0.1720420 0.4021979
## DescN  0.2944261 0.3442003 0.2179726 0.2291372 0.3291507 0.1890209 0.3278693
## HA_12  0.4690087 0.3832823 0.4517847 0.4472835 0.6110485 0.3831671 0.6431643
## HA_34  0.4902185 0.3808996 0.4561996 0.4582385 0.6217566 0.3929884 0.6624370
## EB_12  0.6148327 0.6487695 0.5102184 0.5145910 0.5688091 0.3377211 0.5922310
## EB_34  1.0000000 0.4561410 0.4422934 0.4533827 0.5337360 0.3514100 0.5849420
## EB_56  0.4561410 1.0000000 0.3946176 0.4018261 0.5103852 0.2463809 0.4769093
## SR_12  0.4422934 0.3946176 1.0000000 0.8089257 0.5094293 0.3727764 0.5302526
## SR_34  0.4533827 0.4018261 0.8089257 1.0000000 0.5112970 0.3790889 0.5401448
## ISM_12 0.5337360 0.5103852 0.5094293 0.5112970 1.0000000 0.4607437 0.7831017
## ISM_34 0.3514100 0.2463809 0.3727764 0.3790889 0.4607437 1.0000000 0.5138200
## ISM_56 0.5849420 0.4769093 0.5302526 0.5401448 0.7831017 0.5138200 1.0000000
corrplot(cor_matrix, method = "circle", type = "upper", 
         tl.col = "black", tl.srt = 45, 
         title = "Correlation Plot", mar = c(0, 0, 1, 0))

Point plots for DVs

climate_labels$Condition <- factor(climate_labels$Condition, 
                                   levels = c("No Label", 
                                              "Direct With Support", 
                                              "Indirect With Support", 
                                              "Direct Without Support", 
                                              "Indirect Without Support", 
                                              "Direct With Support With Societal Change",
                                              "Indirect With Support With Societal Change",
                                              "Direct With Societal Change",
                                              "Indirect With Societal Change"),
                                   labels = c("NoLab", 
                                              "DirNorm", 
                                              "IndNorm", 
                                              "Direct", 
                                              "Indirect", 
                                              "DirNormBE",
                                              "IndNormBE",
                                              "DirBE",
                                              "IndBE"))

variables <- c("Inj", "DN1", "DN2", "DescN", 
               "HA_12", "HA_34", 
               "EB_12", "EB_34", "EB_56",
               "SR_12", "SR_34",
               "ISM_12", "ISM_34", "ISM_56")

for (var in variables) {
  plot <- ggplot(climate_labels, aes_string(x = "Condition", y = var)) +
    stat_summary(fun.data = "mean_cl_normal", geom = "pointrange", 
                 color = "blue", size = 1) +
    theme_minimal() +
    labs(title = paste("Point Plot with 95% CIs for", var))
  
  print(plot)
}
## Warning: `aes_string()` was deprecated in ggplot2 3.0.0.
## ℹ Please use tidy evaluation idioms with `aes()`.
## ℹ See also `vignette("ggplot2-in-packages")` for more information.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 3 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 7 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 3 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 3 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 3 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 3 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

## Warning: Removed 2 rows containing non-finite outside the scale range
## (`stat_summary()`).

Pairwise

climate_labels$Condition <- relevel(climate_labels$Condition, ref = "NoLab")

variables <- c("Inj", "DN1", "DN2", "DescN", 
               "HA_12", "HA_34", 
               "EB_12", "EB_34", "EB_56",
               "SR_12", "SR_34",
               "ISM_12", "ISM_34", "ISM_56")

results <- lapply(variables, function(var) {
  formula <- as.formula(paste(var, "~ Condition"))
  model <- lm(formula, data = climate_labels)
  
  pairwise <- emmeans(model, pairwise ~ Condition)
  
  list(
    variable = var,
    model_summary = summary(model),
    pairwise_summary = summary(pairwise)
  )
})

for (res in results) {
  cat("Results for ", res$variable, "\n\n")
  cat("Model Summary\n")
  print(res$model_summary)
  cat("\n\n Pairwise Comparisons\n")
  print(res$pairwise_summary)
  cat("\n\n--------------------------------------------------------------------------------\n\n")
}
## Results for  Inj 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8948 -0.7138 -0.1275  0.7684  2.4440 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.55597    0.04934  51.802  < 2e-16 ***
## ConditionDirNorm    0.33883    0.06969   4.862 1.21e-06 ***
## ConditionIndNorm    0.31658    0.06952   4.554 5.44e-06 ***
## ConditionDirect     0.17565    0.06952   2.527   0.0116 *  
## ConditionIndirect   0.07151    0.06969   1.026   0.3050    
## ConditionDirNormBE  0.16584    0.06952   2.385   0.0171 *  
## ConditionIndNormBE  0.12985    0.06948   1.869   0.0617 .  
## ConditionDirBE      0.15779    0.06956   2.268   0.0234 *  
## ConditionIndBE      0.11195    0.06991   1.601   0.1094    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9893 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.01048,    Adjusted R-squared:  0.008309 
## F-statistic: 4.821 on 8 and 3640 DF,  p-value: 6.361e-06
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       2.56 0.0493 3640     2.46     2.65
##  DirNorm     2.89 0.0492 3640     2.80     2.99
##  IndNorm     2.87 0.0490 3640     2.78     2.97
##  Direct      2.73 0.0490 3640     2.64     2.83
##  Indirect    2.63 0.0492 3640     2.53     2.72
##  DirNormBE   2.72 0.0490 3640     2.63     2.82
##  IndNormBE   2.69 0.0489 3640     2.59     2.78
##  DirBE       2.71 0.0490 3640     2.62     2.81
##  IndBE       2.67 0.0495 3640     2.57     2.77
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm       -0.33883 0.0697 3640  -4.862  <.0001
##  NoLab - IndNorm       -0.31658 0.0695 3640  -4.554  0.0002
##  NoLab - Direct        -0.17565 0.0695 3640  -2.527  0.2194
##  NoLab - Indirect      -0.07151 0.0697 3640  -1.026  0.9835
##  NoLab - DirNormBE     -0.16584 0.0695 3640  -2.385  0.2921
##  NoLab - IndNormBE     -0.12985 0.0695 3640  -1.869  0.6354
##  NoLab - DirBE         -0.15779 0.0696 3640  -2.268  0.3618
##  NoLab - IndBE         -0.11195 0.0699 3640  -1.601  0.8044
##  DirNorm - IndNorm      0.02225 0.0694 3640   0.320  1.0000
##  DirNorm - Direct       0.16318 0.0694 3640   2.350  0.3123
##  DirNorm - Indirect     0.26733 0.0696 3640   3.841  0.0040
##  DirNorm - DirNormBE    0.17299 0.0694 3640   2.491  0.2363
##  DirNorm - IndNormBE    0.20898 0.0694 3640   3.012  0.0653
##  DirNorm - DirBE        0.18104 0.0695 3640   2.606  0.1843
##  DirNorm - IndBE        0.22688 0.0698 3640   3.249  0.0320
##  IndNorm - Direct       0.14093 0.0693 3640   2.035  0.5190
##  IndNorm - Indirect     0.24507 0.0694 3640   3.530  0.0125
##  IndNorm - DirNormBE    0.15074 0.0693 3640   2.176  0.4215
##  IndNorm - IndNormBE    0.18673 0.0692 3640   2.698  0.1489
##  IndNorm - DirBE        0.15879 0.0693 3640   2.291  0.3476
##  IndNorm - IndBE        0.20463 0.0697 3640   2.938  0.0803
##  Direct - Indirect      0.10414 0.0694 3640   1.500  0.8560
##  Direct - DirNormBE     0.00980 0.0693 3640   0.142  1.0000
##  Direct - IndNormBE     0.04580 0.0692 3640   0.662  0.9992
##  Direct - DirBE         0.01786 0.0693 3640   0.258  1.0000
##  Direct - IndBE         0.06370 0.0697 3640   0.915  0.9922
##  Indirect - DirNormBE  -0.09434 0.0694 3640  -1.359  0.9131
##  Indirect - IndNormBE  -0.05834 0.0694 3640  -0.841  0.9956
##  Indirect - DirBE      -0.08628 0.0695 3640  -1.242  0.9472
##  Indirect - IndBE      -0.04044 0.0698 3640  -0.579  0.9997
##  DirNormBE - IndNormBE  0.03599 0.0692 3640   0.520  0.9999
##  DirNormBE - DirBE      0.00805 0.0693 3640   0.116  1.0000
##  DirNormBE - IndBE      0.05389 0.0697 3640   0.774  0.9975
##  IndNormBE - DirBE     -0.02794 0.0693 3640  -0.403  1.0000
##  IndNormBE - IndBE      0.01790 0.0696 3640   0.257  1.0000
##  DirBE - IndBE          0.04584 0.0697 3640   0.658  0.9992
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  DN1 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5833 -0.5490  0.4167  0.5124  2.5348 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.46517    0.06154  72.554   <2e-16 ***
## ConditionDirNorm    0.02245    0.08693   0.258    0.796    
## ConditionIndNorm    0.11816    0.08671   1.363    0.173    
## ConditionDirect     0.02747    0.08671   0.317    0.751    
## ConditionIndirect   0.02988    0.08693   0.344    0.731    
## ConditionDirNormBE  0.08385    0.08671   0.967    0.334    
## ConditionIndNormBE  0.07517    0.08666   0.867    0.386    
## ConditionDirBE      0.04097    0.08677   0.472    0.637    
## ConditionIndBE      0.08620    0.08720   0.989    0.323    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.234 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.0008664,  Adjusted R-squared:  -0.00133 
## F-statistic: 0.3945 on 8 and 3640 DF,  p-value: 0.9241
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       4.47 0.0615 3640     4.34     4.59
##  DirNorm     4.49 0.0614 3640     4.37     4.61
##  IndNorm     4.58 0.0611 3640     4.46     4.70
##  Direct      4.49 0.0611 3640     4.37     4.61
##  Indirect    4.50 0.0614 3640     4.37     4.62
##  DirNormBE   4.55 0.0611 3640     4.43     4.67
##  IndNormBE   4.54 0.0610 3640     4.42     4.66
##  DirBE       4.51 0.0612 3640     4.39     4.63
##  IndBE       4.55 0.0618 3640     4.43     4.67
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm       -0.02245 0.0869 3640  -0.258  1.0000
##  NoLab - IndNorm       -0.11816 0.0867 3640  -1.363  0.9117
##  NoLab - Direct        -0.02747 0.0867 3640  -0.317  1.0000
##  NoLab - Indirect      -0.02988 0.0869 3640  -0.344  1.0000
##  NoLab - DirNormBE     -0.08385 0.0867 3640  -0.967  0.9887
##  NoLab - IndNormBE     -0.07517 0.0867 3640  -0.867  0.9946
##  NoLab - DirBE         -0.04097 0.0868 3640  -0.472  0.9999
##  NoLab - IndBE         -0.08620 0.0872 3640  -0.989  0.9870
##  DirNorm - IndNorm     -0.09571 0.0866 3640  -1.105  0.9736
##  DirNorm - Direct      -0.00502 0.0866 3640  -0.058  1.0000
##  DirNorm - Indirect    -0.00743 0.0868 3640  -0.086  1.0000
##  DirNorm - DirNormBE   -0.06140 0.0866 3640  -0.709  0.9987
##  DirNorm - IndNormBE   -0.05272 0.0866 3640  -0.609  0.9996
##  DirNorm - DirBE       -0.01852 0.0867 3640  -0.214  1.0000
##  DirNorm - IndBE       -0.06375 0.0871 3640  -0.732  0.9983
##  IndNorm - Direct       0.09069 0.0864 3640   1.050  0.9809
##  IndNorm - Indirect     0.08828 0.0866 3640   1.019  0.9841
##  IndNorm - DirNormBE    0.03431 0.0864 3640   0.397  1.0000
##  IndNorm - IndNormBE    0.04299 0.0863 3640   0.498  0.9999
##  IndNorm - DirBE        0.07719 0.0864 3640   0.893  0.9934
##  IndNorm - IndBE        0.03195 0.0869 3640   0.368  1.0000
##  Direct - Indirect     -0.00240 0.0866 3640  -0.028  1.0000
##  Direct - DirNormBE    -0.05637 0.0864 3640  -0.653  0.9993
##  Direct - IndNormBE    -0.04770 0.0863 3640  -0.552  0.9998
##  Direct - DirBE        -0.01350 0.0864 3640  -0.156  1.0000
##  Direct - IndBE        -0.05873 0.0869 3640  -0.676  0.9991
##  Indirect - DirNormBE  -0.05397 0.0866 3640  -0.623  0.9995
##  Indirect - IndNormBE  -0.04529 0.0866 3640  -0.523  0.9999
##  Indirect - DirBE      -0.01109 0.0867 3640  -0.128  1.0000
##  Indirect - IndBE      -0.05633 0.0871 3640  -0.647  0.9993
##  DirNormBE - IndNormBE  0.00868 0.0863 3640   0.101  1.0000
##  DirNormBE - DirBE      0.04288 0.0864 3640   0.496  0.9999
##  DirNormBE - IndBE     -0.00236 0.0869 3640  -0.027  1.0000
##  IndNormBE - DirBE      0.03420 0.0864 3640   0.396  1.0000
##  IndNormBE - IndBE     -0.01104 0.0868 3640  -0.127  1.0000
##  DirBE - IndBE         -0.04524 0.0869 3640  -0.520  0.9999
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  DN2 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4662 -1.3652  0.5338  0.6990  2.6990 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.30100    0.07071  46.683   <2e-16 ***
## ConditionDirNorm    0.14207    0.09988   1.423   0.1550    
## ConditionIndNorm    0.06420    0.09963   0.644   0.5194    
## ConditionDirect     0.11812    0.09963   1.186   0.2359    
## ConditionIndirect   0.07277    0.09988   0.729   0.4663    
## ConditionDirNormBE  0.14018    0.09963   1.407   0.1595    
## ConditionIndNormBE  0.13038    0.09963   1.309   0.1908    
## ConditionDirBE      0.09950    0.09969   0.998   0.3183    
## ConditionIndBE      0.16517    0.10019   1.649   0.0993 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.418 on 3639 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.001144,   Adjusted R-squared:  -0.001052 
## F-statistic: 0.5208 on 8 and 3639 DF,  p-value: 0.8417
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       3.30 0.0707 3639     3.16     3.44
##  DirNorm     3.44 0.0705 3639     3.30     3.58
##  IndNorm     3.37 0.0702 3639     3.23     3.50
##  Direct      3.42 0.0702 3639     3.28     3.56
##  Indirect    3.37 0.0705 3639     3.24     3.51
##  DirNormBE   3.44 0.0702 3639     3.30     3.58
##  IndNormBE   3.43 0.0702 3639     3.29     3.57
##  DirBE       3.40 0.0703 3639     3.26     3.54
##  IndBE       3.47 0.0710 3639     3.33     3.61
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm       -0.14207 0.0999 3639  -1.423  0.8894
##  NoLab - IndNorm       -0.06420 0.0996 3639  -0.644  0.9993
##  NoLab - Direct        -0.11812 0.0996 3639  -1.186  0.9597
##  NoLab - Indirect      -0.07277 0.0999 3639  -0.729  0.9984
##  NoLab - DirNormBE     -0.14018 0.0996 3639  -1.407  0.8955
##  NoLab - IndNormBE     -0.13038 0.0996 3639  -1.309  0.9291
##  NoLab - DirBE         -0.09950 0.0997 3639  -0.998  0.9862
##  NoLab - IndBE         -0.16517 0.1002 3639  -1.649  0.7776
##  DirNorm - IndNorm      0.07787 0.0995 3639   0.783  0.9973
##  DirNorm - Direct       0.02395 0.0995 3639   0.241  1.0000
##  DirNorm - Indirect     0.06931 0.0998 3639   0.695  0.9989
##  DirNorm - DirNormBE    0.00189 0.0995 3639   0.019  1.0000
##  DirNorm - IndNormBE    0.01170 0.0995 3639   0.118  1.0000
##  DirNorm - DirBE        0.04258 0.0996 3639   0.428  1.0000
##  DirNorm - IndBE       -0.02310 0.1001 3639  -0.231  1.0000
##  IndNorm - Direct      -0.05392 0.0993 3639  -0.543  0.9998
##  IndNorm - Indirect    -0.00857 0.0995 3639  -0.086  1.0000
##  IndNorm - DirNormBE   -0.07598 0.0993 3639  -0.765  0.9977
##  IndNorm - IndNormBE   -0.06618 0.0993 3639  -0.667  0.9992
##  IndNorm - DirBE       -0.03530 0.0993 3639  -0.355  1.0000
##  IndNorm - IndBE       -0.10097 0.0998 3639  -1.012  0.9849
##  Direct - Indirect      0.04536 0.0995 3639   0.456  1.0000
##  Direct - DirNormBE    -0.02206 0.0993 3639  -0.222  1.0000
##  Direct - IndNormBE    -0.01225 0.0993 3639  -0.123  1.0000
##  Direct - DirBE         0.01863 0.0993 3639   0.188  1.0000
##  Direct - IndBE        -0.04705 0.0998 3639  -0.471  0.9999
##  Indirect - DirNormBE  -0.06741 0.0995 3639  -0.677  0.9991
##  Indirect - IndNormBE  -0.05761 0.0995 3639  -0.579  0.9997
##  Indirect - DirBE      -0.02673 0.0996 3639  -0.268  1.0000
##  Indirect - IndBE      -0.09240 0.1001 3639  -0.923  0.9917
##  DirNormBE - IndNormBE  0.00980 0.0993 3639   0.099  1.0000
##  DirNormBE - DirBE      0.04069 0.0993 3639   0.410  1.0000
##  DirNormBE - IndBE     -0.02499 0.0998 3639  -0.250  1.0000
##  IndNormBE - DirBE      0.03088 0.0993 3639   0.311  1.0000
##  IndNormBE - IndBE     -0.03479 0.0998 3639  -0.349  1.0000
##  DirBE - IndBE         -0.06567 0.0999 3639  -0.658  0.9992
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  DescN 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -46.833 -19.908   0.092  18.167  60.092 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          33.980      1.211  28.067  < 2e-16 ***
## ConditionDirNorm     10.176      1.710   5.951 2.92e-09 ***
## ConditionIndNorm     12.853      1.707   7.530 6.36e-14 ***
## ConditionDirect       7.321      1.706   4.292 1.82e-05 ***
## ConditionIndirect     5.928      1.711   3.465 0.000537 ***
## ConditionDirNormBE   10.821      1.706   6.344 2.52e-10 ***
## ConditionIndNormBE    9.693      1.707   5.679 1.46e-08 ***
## ConditionDirBE        9.453      1.708   5.535 3.33e-08 ***
## ConditionIndBE        9.757      1.715   5.688 1.39e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.27 on 3635 degrees of freedom
##   (7 observations deleted due to missingness)
## Multiple R-squared:  0.02053,    Adjusted R-squared:  0.01837 
## F-statistic: 9.524 on 8 and 3635 DF,  p-value: 3.98e-13
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean   SE   df lower.CL upper.CL
##  NoLab       34.0 1.21 3635     31.6     36.4
##  DirNorm     44.2 1.21 3635     41.8     46.5
##  IndNorm     46.8 1.20 3635     44.5     49.2
##  Direct      41.3 1.20 3635     38.9     43.7
##  Indirect    39.9 1.21 3635     37.5     42.3
##  DirNormBE   44.8 1.20 3635     42.4     47.2
##  IndNormBE   43.7 1.20 3635     41.3     46.0
##  DirBE       43.4 1.20 3635     41.1     45.8
##  IndBE       43.7 1.22 3635     41.4     46.1
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate   SE   df t.ratio p.value
##  NoLab - DirNorm       -10.1758 1.71 3635  -5.951  <.0001
##  NoLab - IndNorm       -12.8528 1.71 3635  -7.530  <.0001
##  NoLab - Direct         -7.3214 1.71 3635  -4.292  0.0006
##  NoLab - Indirect       -5.9281 1.71 3635  -3.465  0.0157
##  NoLab - DirNormBE     -10.8214 1.71 3635  -6.344  <.0001
##  NoLab - IndNormBE      -9.6931 1.71 3635  -5.679  <.0001
##  NoLab - DirBE          -9.4532 1.71 3635  -5.535  <.0001
##  NoLab - IndBE          -9.7567 1.72 3635  -5.688  <.0001
##  DirNorm - IndNorm      -2.6770 1.70 3635  -1.570  0.8211
##  DirNorm - Direct        2.8545 1.70 3635   1.675  0.7617
##  DirNorm - Indirect      4.2478 1.71 3635   2.486  0.2392
##  DirNorm - DirNormBE    -0.6455 1.70 3635  -0.379  1.0000
##  DirNorm - IndNormBE     0.4827 1.70 3635   0.283  1.0000
##  DirNorm - DirBE         0.7227 1.71 3635   0.424  1.0000
##  DirNorm - IndBE         0.4191 1.71 3635   0.245  1.0000
##  IndNorm - Direct        5.5315 1.70 3635   3.253  0.0316
##  IndNorm - Indirect      6.9247 1.71 3635   4.060  0.0017
##  IndNorm - DirNormBE     2.0315 1.70 3635   1.195  0.9578
##  IndNorm - IndNormBE     3.1597 1.70 3635   1.857  0.6437
##  IndNorm - DirBE         3.3997 1.70 3635   1.997  0.5457
##  IndNorm - IndBE         3.0961 1.71 3635   1.810  0.6753
##  Direct - Indirect       1.3933 1.70 3635   0.817  0.9964
##  Direct - DirNormBE     -3.5000 1.70 3635  -2.059  0.5016
##  Direct - IndNormBE     -2.3717 1.70 3635  -1.395  0.9002
##  Direct - DirBE         -2.1318 1.70 3635  -1.253  0.9444
##  Direct - IndBE         -2.4354 1.71 3635  -1.425  0.8884
##  Indirect - DirNormBE   -4.8933 1.70 3635  -2.870  0.0963
##  Indirect - IndNormBE   -3.7650 1.71 3635  -2.207  0.4010
##  Indirect - DirBE       -3.5251 1.71 3635  -2.065  0.4975
##  Indirect - IndBE       -3.8287 1.71 3635  -2.233  0.3840
##  DirNormBE - IndNormBE   1.1283 1.70 3635   0.663  0.9992
##  DirNormBE - DirBE       1.3682 1.70 3635   0.804  0.9968
##  DirNormBE - IndBE       1.0646 1.71 3635   0.623  0.9995
##  IndNormBE - DirBE       0.2400 1.70 3635   0.141  1.0000
##  IndNormBE - IndBE      -0.0636 1.71 3635  -0.037  1.0000
##  DirBE - IndBE          -0.3036 1.71 3635  -0.177  1.0000
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  HA_12 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2617 -0.6942 -0.0348  0.8203  1.9652 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.03483    0.06234  64.721  < 2e-16 ***
## ConditionDirNorm    0.12646    0.08811   1.435  0.15129    
## ConditionIndNorm    0.18454    0.08784   2.101  0.03572 *  
## ConditionDirect     0.22130    0.08784   2.519  0.01180 *  
## ConditionIndirect   0.14834    0.08806   1.685  0.09214 .  
## ConditionDirNormBE  0.04238    0.08784   0.482  0.62950    
## ConditionIndNormBE  0.14488    0.08779   1.650  0.09895 .  
## ConditionDirBE      0.22684    0.08789   2.581  0.00989 ** 
## ConditionIndBE      0.15941    0.08833   1.805  0.07120 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.25 on 3639 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.003262,   Adjusted R-squared:  0.00107 
## F-statistic: 1.488 on 8 and 3639 DF,  p-value: 0.1558
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       4.03 0.0623 3639     3.91     4.16
##  DirNorm     4.16 0.0623 3639     4.04     4.28
##  IndNorm     4.22 0.0619 3639     4.10     4.34
##  Direct      4.26 0.0619 3639     4.13     4.38
##  Indirect    4.18 0.0622 3639     4.06     4.31
##  DirNormBE   4.08 0.0619 3639     3.96     4.20
##  IndNormBE   4.18 0.0618 3639     4.06     4.30
##  DirBE       4.26 0.0620 3639     4.14     4.38
##  IndBE       4.19 0.0626 3639     4.07     4.32
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm       -0.12646 0.0881 3639  -1.435  0.8843
##  NoLab - IndNorm       -0.18454 0.0878 3639  -2.101  0.4728
##  NoLab - Direct        -0.22130 0.0878 3639  -2.519  0.2228
##  NoLab - Indirect      -0.14834 0.0881 3639  -1.685  0.7562
##  NoLab - DirNormBE     -0.04238 0.0878 3639  -0.482  0.9999
##  NoLab - IndNormBE     -0.14488 0.0878 3639  -1.650  0.7766
##  NoLab - DirBE         -0.22684 0.0879 3639  -2.581  0.1949
##  NoLab - IndBE         -0.15941 0.0883 3639  -1.805  0.6792
##  DirNorm - IndNorm     -0.05807 0.0878 3639  -0.662  0.9992
##  DirNorm - Direct      -0.09484 0.0878 3639  -1.080  0.9771
##  DirNorm - Indirect    -0.02188 0.0880 3639  -0.249  1.0000
##  DirNorm - DirNormBE    0.08408 0.0878 3639   0.958  0.9894
##  DirNorm - IndNormBE   -0.01842 0.0877 3639  -0.210  1.0000
##  DirNorm - DirBE       -0.10038 0.0878 3639  -1.143  0.9676
##  DirNorm - IndBE       -0.03295 0.0883 3639  -0.373  1.0000
##  IndNorm - Direct      -0.03676 0.0875 3639  -0.420  1.0000
##  IndNorm - Indirect     0.03619 0.0877 3639   0.413  1.0000
##  IndNorm - DirNormBE    0.14216 0.0875 3639   1.624  0.7916
##  IndNorm - IndNormBE    0.03966 0.0875 3639   0.453  1.0000
##  IndNorm - DirBE       -0.04231 0.0876 3639  -0.483  0.9999
##  IndNorm - IndBE        0.02513 0.0880 3639   0.286  1.0000
##  Direct - Indirect      0.07296 0.0877 3639   0.832  0.9959
##  Direct - DirNormBE     0.17892 0.0875 3639   2.044  0.5121
##  Direct - IndNormBE     0.07642 0.0875 3639   0.874  0.9943
##  Direct - DirBE        -0.00554 0.0876 3639  -0.063  1.0000
##  Direct - IndBE         0.06189 0.0880 3639   0.703  0.9988
##  Indirect - DirNormBE   0.10596 0.0877 3639   1.208  0.9550
##  Indirect - IndNormBE   0.00346 0.0877 3639   0.039  1.0000
##  Indirect - DirBE      -0.07850 0.0878 3639  -0.894  0.9933
##  Indirect - IndBE      -0.01107 0.0882 3639  -0.125  1.0000
##  DirNormBE - IndNormBE -0.10250 0.0875 3639  -1.172  0.9623
##  DirNormBE - DirBE     -0.18446 0.0876 3639  -2.107  0.4689
##  DirNormBE - IndBE     -0.11703 0.0880 3639  -1.330  0.9226
##  IndNormBE - DirBE     -0.08196 0.0875 3639  -0.937  0.9909
##  IndNormBE - IndBE     -0.01453 0.0880 3639  -0.165  1.0000
##  DirBE - IndBE          0.06744 0.0881 3639   0.766  0.9977
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  HA_34 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3096 -0.7221  0.2281  0.8064  1.8731 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.12687    0.06536  63.142   <2e-16 ***
## ConditionDirNorm    0.09522    0.09237   1.031   0.3027    
## ConditionIndNorm    0.14152    0.09209   1.537   0.1245    
## ConditionDirect     0.17583    0.09209   1.909   0.0563 .  
## ConditionIndirect   0.08972    0.09232   0.972   0.3312    
## ConditionDirNormBE  0.06676    0.09209   0.725   0.4685    
## ConditionIndNormBE  0.07607    0.09203   0.827   0.4086    
## ConditionDirBE      0.18272    0.09215   1.983   0.0475 *  
## ConditionIndBE      0.14506    0.09260   1.566   0.1173    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.31 on 3639 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.001782,   Adjusted R-squared:  -0.0004124 
## F-statistic: 0.8121 on 8 and 3639 DF,  p-value: 0.5918
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       4.13 0.0654 3639     4.00     4.26
##  DirNorm     4.22 0.0653 3639     4.09     4.35
##  IndNorm     4.27 0.0649 3639     4.14     4.40
##  Direct      4.30 0.0649 3639     4.18     4.43
##  Indirect    4.22 0.0652 3639     4.09     4.34
##  DirNormBE   4.19 0.0649 3639     4.07     4.32
##  IndNormBE   4.20 0.0648 3639     4.08     4.33
##  DirBE       4.31 0.0650 3639     4.18     4.44
##  IndBE       4.27 0.0656 3639     4.14     4.40
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm       -0.09522 0.0924 3639  -1.031  0.9830
##  NoLab - IndNorm       -0.14152 0.0921 3639  -1.537  0.8382
##  NoLab - Direct        -0.17583 0.0921 3639  -1.909  0.6073
##  NoLab - Indirect      -0.08972 0.0923 3639  -0.972  0.9884
##  NoLab - DirNormBE     -0.06676 0.0921 3639  -0.725  0.9985
##  NoLab - IndNormBE     -0.07607 0.0920 3639  -0.827  0.9961
##  NoLab - DirBE         -0.18272 0.0921 3639  -1.983  0.5555
##  NoLab - IndBE         -0.14506 0.0926 3639  -1.566  0.8231
##  DirNorm - IndNorm     -0.04630 0.0920 3639  -0.503  0.9999
##  DirNorm - Direct      -0.08061 0.0920 3639  -0.876  0.9942
##  DirNorm - Indirect     0.00550 0.0923 3639   0.060  1.0000
##  DirNorm - DirNormBE    0.02846 0.0920 3639   0.309  1.0000
##  DirNorm - IndNormBE    0.01915 0.0920 3639   0.208  1.0000
##  DirNorm - DirBE       -0.08750 0.0921 3639  -0.950  0.9900
##  DirNorm - IndBE       -0.04985 0.0925 3639  -0.539  0.9998
##  IndNorm - Direct      -0.03431 0.0917 3639  -0.374  1.0000
##  IndNorm - Indirect     0.05180 0.0920 3639   0.563  0.9998
##  IndNorm - DirNormBE    0.07475 0.0917 3639   0.815  0.9965
##  IndNorm - IndNormBE    0.06545 0.0917 3639   0.714  0.9986
##  IndNorm - DirBE       -0.04120 0.0918 3639  -0.449  1.0000
##  IndNorm - IndBE       -0.00355 0.0923 3639  -0.038  1.0000
##  Direct - Indirect      0.08611 0.0920 3639   0.936  0.9909
##  Direct - DirNormBE     0.10907 0.0917 3639   1.189  0.9590
##  Direct - IndNormBE     0.09976 0.0917 3639   1.088  0.9760
##  Direct - DirBE        -0.00689 0.0918 3639  -0.075  1.0000
##  Direct - IndBE         0.03077 0.0923 3639   0.333  1.0000
##  Indirect - DirNormBE   0.02296 0.0920 3639   0.250  1.0000
##  Indirect - IndNormBE   0.01365 0.0919 3639   0.149  1.0000
##  Indirect - DirBE      -0.09300 0.0920 3639  -1.011  0.9850
##  Indirect - IndBE      -0.05535 0.0925 3639  -0.598  0.9996
##  DirNormBE - IndNormBE -0.00931 0.0917 3639  -0.101  1.0000
##  DirNormBE - DirBE     -0.11595 0.0918 3639  -1.263  0.9418
##  DirNormBE - IndBE     -0.07830 0.0923 3639  -0.849  0.9953
##  IndNormBE - DirBE     -0.10665 0.0917 3639  -1.162  0.9641
##  IndNormBE - IndBE     -0.06900 0.0922 3639  -0.748  0.9981
##  DirBE - IndBE          0.03765 0.0923 3639   0.408  1.0000
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  EB_12 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.9484 -0.8247  0.1287  0.6683  2.2364 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.874378   0.052639  54.605   <2e-16 ***
## ConditionDirNorm   -0.110764   0.074351  -1.490    0.136    
## ConditionIndNorm   -0.015309   0.074169  -0.206    0.836    
## ConditionDirect    -0.003055   0.074169  -0.041    0.967    
## ConditionIndirect  -0.042695   0.074351  -0.574    0.566    
## ConditionDirNormBE -0.053300   0.074169  -0.719    0.472    
## ConditionIndNormBE -0.049623   0.074169  -0.669    0.504    
## ConditionDirBE      0.074025   0.074214   0.997    0.319    
## ConditionIndBE     -0.014729   0.074583  -0.197    0.843    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.055 on 3639 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.002016,   Adjusted R-squared:  -0.0001783 
## F-statistic: 0.9187 on 8 and 3639 DF,  p-value: 0.4995
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       2.87 0.0526 3639     2.77     2.98
##  DirNorm     2.76 0.0525 3639     2.66     2.87
##  IndNorm     2.86 0.0523 3639     2.76     2.96
##  Direct      2.87 0.0523 3639     2.77     2.97
##  Indirect    2.83 0.0525 3639     2.73     2.93
##  DirNormBE   2.82 0.0523 3639     2.72     2.92
##  IndNormBE   2.82 0.0523 3639     2.72     2.93
##  DirBE       2.95 0.0523 3639     2.85     3.05
##  IndBE       2.86 0.0528 3639     2.76     2.96
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm        0.11076 0.0744 3639   1.490  0.8607
##  NoLab - IndNorm        0.01531 0.0742 3639   0.206  1.0000
##  NoLab - Direct         0.00305 0.0742 3639   0.041  1.0000
##  NoLab - Indirect       0.04269 0.0744 3639   0.574  0.9997
##  NoLab - DirNormBE      0.05330 0.0742 3639   0.719  0.9985
##  NoLab - IndNormBE      0.04962 0.0742 3639   0.669  0.9991
##  NoLab - DirBE         -0.07402 0.0742 3639  -0.997  0.9862
##  NoLab - IndBE          0.01473 0.0746 3639   0.197  1.0000
##  DirNorm - IndNorm     -0.09545 0.0741 3639  -1.289  0.9349
##  DirNorm - Direct      -0.10771 0.0741 3639  -1.454  0.8764
##  DirNorm - Indirect    -0.06807 0.0743 3639  -0.917  0.9921
##  DirNorm - DirNormBE   -0.05746 0.0741 3639  -0.776  0.9975
##  DirNorm - IndNormBE   -0.06114 0.0741 3639  -0.825  0.9961
##  DirNorm - DirBE       -0.18479 0.0741 3639  -2.493  0.2355
##  DirNorm - IndBE       -0.09604 0.0745 3639  -1.289  0.9347
##  IndNorm - Direct      -0.01225 0.0739 3639  -0.166  1.0000
##  IndNorm - Indirect     0.02739 0.0741 3639   0.370  1.0000
##  IndNorm - DirNormBE    0.03799 0.0739 3639   0.514  0.9999
##  IndNorm - IndNormBE    0.03431 0.0739 3639   0.464  0.9999
##  IndNorm - DirBE       -0.08933 0.0739 3639  -1.208  0.9549
##  IndNorm - IndBE       -0.00058 0.0743 3639  -0.008  1.0000
##  Direct - Indirect      0.03964 0.0741 3639   0.535  0.9998
##  Direct - DirNormBE     0.05025 0.0739 3639   0.680  0.9990
##  Direct - IndNormBE     0.04657 0.0739 3639   0.630  0.9994
##  Direct - DirBE        -0.07708 0.0739 3639  -1.042  0.9817
##  Direct - IndBE         0.01167 0.0743 3639   0.157  1.0000
##  Indirect - DirNormBE   0.01060 0.0741 3639   0.143  1.0000
##  Indirect - IndNormBE   0.00693 0.0741 3639   0.094  1.0000
##  Indirect - DirBE      -0.11672 0.0741 3639  -1.575  0.8188
##  Indirect - IndBE      -0.02797 0.0745 3639  -0.375  1.0000
##  DirNormBE - IndNormBE -0.00368 0.0739 3639  -0.050  1.0000
##  DirNormBE - DirBE     -0.12732 0.0739 3639  -1.722  0.7330
##  DirNormBE - IndBE     -0.03857 0.0743 3639  -0.519  0.9999
##  IndNormBE - DirBE     -0.12365 0.0739 3639  -1.672  0.7636
##  IndNormBE - IndBE     -0.03489 0.0743 3639  -0.470  0.9999
##  DirBE - IndBE          0.08875 0.0744 3639   1.194  0.9580
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  EB_34 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.17813 -0.58333 -0.08333  0.43521  1.93936 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.159204   0.046739  67.593   <2e-16 ***
## ConditionDirNorm   -0.098560   0.066017  -1.493    0.136    
## ConditionIndNorm   -0.009694   0.065855  -0.147    0.883    
## ConditionDirect    -0.075871   0.065855  -1.152    0.249    
## ConditionIndirect  -0.058956   0.066017  -0.893    0.372    
## ConditionDirNormBE -0.093028   0.065855  -1.413    0.158    
## ConditionIndNormBE -0.059939   0.065855  -0.910    0.363    
## ConditionDirBE      0.018929   0.065895   0.287    0.774    
## ConditionIndBE      0.006210   0.066223   0.094    0.925    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9371 on 3639 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.002074,   Adjusted R-squared:  -0.0001198 
## F-statistic: 0.9454 on 8 and 3639 DF,  p-value: 0.4774
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       3.16 0.0467 3639     3.07     3.25
##  DirNorm     3.06 0.0466 3639     2.97     3.15
##  IndNorm     3.15 0.0464 3639     3.06     3.24
##  Direct      3.08 0.0464 3639     2.99     3.17
##  Indirect    3.10 0.0466 3639     3.01     3.19
##  DirNormBE   3.07 0.0464 3639     2.98     3.16
##  IndNormBE   3.10 0.0464 3639     3.01     3.19
##  DirBE       3.18 0.0465 3639     3.09     3.27
##  IndBE       3.17 0.0469 3639     3.07     3.26
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast               estimate     SE   df t.ratio p.value
##  NoLab - DirNorm        0.098560 0.0660 3639   1.493  0.8592
##  NoLab - IndNorm        0.009694 0.0659 3639   0.147  1.0000
##  NoLab - Direct         0.075871 0.0659 3639   1.152  0.9660
##  NoLab - Indirect       0.058956 0.0660 3639   0.893  0.9934
##  NoLab - DirNormBE      0.093028 0.0659 3639   1.413  0.8933
##  NoLab - IndNormBE      0.059939 0.0659 3639   0.910  0.9925
##  NoLab - DirBE         -0.018929 0.0659 3639  -0.287  1.0000
##  NoLab - IndBE         -0.006210 0.0662 3639  -0.094  1.0000
##  DirNorm - IndNorm     -0.088866 0.0658 3639  -1.351  0.9156
##  DirNorm - Direct      -0.022690 0.0658 3639  -0.345  1.0000
##  DirNorm - Indirect    -0.039604 0.0659 3639  -0.601  0.9996
##  DirNorm - DirNormBE   -0.005533 0.0658 3639  -0.084  1.0000
##  DirNorm - IndNormBE   -0.038621 0.0658 3639  -0.587  0.9997
##  DirNorm - DirBE       -0.117489 0.0658 3639  -1.785  0.6922
##  DirNorm - IndBE       -0.104770 0.0661 3639  -1.584  0.8138
##  IndNorm - Direct       0.066176 0.0656 3639   1.009  0.9852
##  IndNorm - Indirect     0.049262 0.0658 3639   0.749  0.9980
##  IndNorm - DirNormBE    0.083333 0.0656 3639   1.270  0.9400
##  IndNorm - IndNormBE    0.050245 0.0656 3639   0.766  0.9977
##  IndNorm - DirBE       -0.028623 0.0657 3639  -0.436  1.0000
##  IndNorm - IndBE       -0.015904 0.0660 3639  -0.241  1.0000
##  Direct - Indirect     -0.016914 0.0658 3639  -0.257  1.0000
##  Direct - DirNormBE     0.017157 0.0656 3639   0.261  1.0000
##  Direct - IndNormBE    -0.015931 0.0656 3639  -0.243  1.0000
##  Direct - DirBE        -0.094799 0.0657 3639  -1.444  0.8807
##  Direct - IndBE        -0.082080 0.0660 3639  -1.244  0.9466
##  Indirect - DirNormBE   0.034071 0.0658 3639   0.518  0.9999
##  Indirect - IndNormBE   0.000983 0.0658 3639   0.015  1.0000
##  Indirect - DirBE      -0.077885 0.0658 3639  -1.183  0.9601
##  Indirect - IndBE      -0.065166 0.0661 3639  -0.985  0.9873
##  DirNormBE - IndNormBE -0.033088 0.0656 3639  -0.504  0.9999
##  DirNormBE - DirBE     -0.111956 0.0657 3639  -1.705  0.7435
##  DirNormBE - IndBE     -0.099237 0.0660 3639  -1.504  0.8540
##  IndNormBE - DirBE     -0.078868 0.0657 3639  -1.201  0.9564
##  IndNormBE - IndBE     -0.066149 0.0660 3639  -1.003  0.9857
##  DirBE - IndBE          0.012719 0.0660 3639   0.193  1.0000
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  EB_56 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.64109 -0.99632 -0.09559  0.87438  2.50368 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.625622   0.053819  48.786   <2e-16 ***
## ConditionDirNorm   -0.042701   0.076018  -0.562   0.5743    
## ConditionIndNorm   -0.007975   0.075832  -0.105   0.9163    
## ConditionDirect    -0.030034   0.075832  -0.396   0.6921    
## ConditionIndirect   0.015467   0.076018   0.203   0.8388    
## ConditionDirNormBE -0.129298   0.075832  -1.705   0.0883 .  
## ConditionIndNormBE -0.035157   0.075786  -0.464   0.6427    
## ConditionDirBE     -0.011371   0.075878  -0.150   0.8809    
## ConditionIndBE     -0.024118   0.076255  -0.316   0.7518    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.079 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.001331,   Adjusted R-squared:  -0.000864 
## F-statistic: 0.6063 on 8 and 3640 DF,  p-value: 0.7733
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       2.63 0.0538 3640     2.52     2.73
##  DirNorm     2.58 0.0537 3640     2.48     2.69
##  IndNorm     2.62 0.0534 3640     2.51     2.72
##  Direct      2.60 0.0534 3640     2.49     2.70
##  Indirect    2.64 0.0537 3640     2.54     2.75
##  DirNormBE   2.50 0.0534 3640     2.39     2.60
##  IndNormBE   2.59 0.0534 3640     2.49     2.70
##  DirBE       2.61 0.0535 3640     2.51     2.72
##  IndBE       2.60 0.0540 3640     2.50     2.71
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm        0.04270 0.0760 3640   0.562  0.9998
##  NoLab - IndNorm        0.00797 0.0758 3640   0.105  1.0000
##  NoLab - Direct         0.03003 0.0758 3640   0.396  1.0000
##  NoLab - Indirect      -0.01547 0.0760 3640  -0.203  1.0000
##  NoLab - DirNormBE      0.12930 0.0758 3640   1.705  0.7436
##  NoLab - IndNormBE      0.03516 0.0758 3640   0.464  0.9999
##  NoLab - DirBE          0.01137 0.0759 3640   0.150  1.0000
##  NoLab - IndBE          0.02412 0.0763 3640   0.316  1.0000
##  DirNorm - IndNorm     -0.03473 0.0757 3640  -0.459  0.9999
##  DirNorm - Direct      -0.01267 0.0757 3640  -0.167  1.0000
##  DirNorm - Indirect    -0.05817 0.0759 3640  -0.766  0.9977
##  DirNorm - DirNormBE    0.08660 0.0757 3640   1.143  0.9675
##  DirNorm - IndNormBE   -0.00754 0.0757 3640  -0.100  1.0000
##  DirNorm - DirBE       -0.03133 0.0758 3640  -0.413  1.0000
##  DirNorm - IndBE       -0.01858 0.0762 3640  -0.244  1.0000
##  IndNorm - Direct       0.02206 0.0756 3640   0.292  1.0000
##  IndNorm - Indirect    -0.02344 0.0757 3640  -0.310  1.0000
##  IndNorm - DirNormBE    0.12132 0.0756 3640   1.606  0.8019
##  IndNorm - IndNormBE    0.02718 0.0755 3640   0.360  1.0000
##  IndNorm - DirBE        0.00340 0.0756 3640   0.045  1.0000
##  IndNorm - IndBE        0.01614 0.0760 3640   0.212  1.0000
##  Direct - Indirect     -0.04550 0.0757 3640  -0.601  0.9996
##  Direct - DirNormBE     0.09926 0.0756 3640   1.314  0.9275
##  Direct - IndNormBE     0.00512 0.0755 3640   0.068  1.0000
##  Direct - DirBE        -0.01866 0.0756 3640  -0.247  1.0000
##  Direct - IndBE        -0.00592 0.0760 3640  -0.078  1.0000
##  Indirect - DirNormBE   0.14477 0.0757 3640   1.911  0.6058
##  Indirect - IndNormBE   0.05062 0.0757 3640   0.669  0.9991
##  Indirect - DirBE       0.02684 0.0758 3640   0.354  1.0000
##  Indirect - IndBE       0.03959 0.0762 3640   0.520  0.9999
##  DirNormBE - IndNormBE -0.09414 0.0755 3640  -1.247  0.9459
##  DirNormBE - DirBE     -0.11793 0.0756 3640  -1.560  0.8265
##  DirNormBE - IndBE     -0.10518 0.0760 3640  -1.384  0.9040
##  IndNormBE - DirBE     -0.02379 0.0756 3640  -0.315  1.0000
##  IndNormBE - IndBE     -0.01104 0.0759 3640  -0.145  1.0000
##  DirBE - IndBE          0.01275 0.0760 3640   0.168  1.0000
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  SR_12 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2995 -0.7494 -0.1940  0.7506  2.8431 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.69403    0.06664  55.429   <2e-16 ***
## ConditionDirNorm    0.01513    0.09413   0.161    0.872    
## ConditionIndNorm    0.10254    0.09390   1.092    0.275    
## ConditionDirect     0.10254    0.09390   1.092    0.275    
## ConditionIndirect   0.10548    0.09413   1.120    0.263    
## ConditionDirNormBE -0.03717    0.09390  -0.396    0.692    
## ConditionIndNormBE  0.05536    0.09385   0.590    0.555    
## ConditionDirBE      0.08607    0.09396   0.916    0.360    
## ConditionIndBE      0.01023    0.09443   0.108    0.914    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.336 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.001423,   Adjusted R-squared:  -0.0007717 
## F-statistic: 0.6484 on 8 and 3640 DF,  p-value: 0.7374
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       3.69 0.0666 3640     3.56     3.82
##  DirNorm     3.71 0.0665 3640     3.58     3.84
##  IndNorm     3.80 0.0662 3640     3.67     3.93
##  Direct      3.80 0.0662 3640     3.67     3.93
##  Indirect    3.80 0.0665 3640     3.67     3.93
##  DirNormBE   3.66 0.0662 3640     3.53     3.79
##  IndNormBE   3.75 0.0661 3640     3.62     3.88
##  DirBE       3.78 0.0662 3640     3.65     3.91
##  IndBE       3.70 0.0669 3640     3.57     3.84
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm       -0.01513 0.0941 3640  -0.161  1.0000
##  NoLab - IndNorm       -0.10254 0.0939 3640  -1.092  0.9755
##  NoLab - Direct        -0.10254 0.0939 3640  -1.092  0.9755
##  NoLab - Indirect      -0.10548 0.0941 3640  -1.120  0.9713
##  NoLab - DirNormBE      0.03717 0.0939 3640   0.396  1.0000
##  NoLab - IndNormBE     -0.05536 0.0938 3640  -0.590  0.9997
##  NoLab - DirBE         -0.08607 0.0940 3640  -0.916  0.9921
##  NoLab - IndBE         -0.01023 0.0944 3640  -0.108  1.0000
##  DirNorm - IndNorm     -0.08741 0.0938 3640  -0.932  0.9912
##  DirNorm - Direct      -0.08741 0.0938 3640  -0.932  0.9912
##  DirNorm - Indirect    -0.09035 0.0940 3640  -0.961  0.9892
##  DirNorm - DirNormBE    0.05230 0.0938 3640   0.558  0.9998
##  DirNorm - IndNormBE   -0.04023 0.0937 3640  -0.429  1.0000
##  DirNorm - DirBE       -0.07094 0.0938 3640  -0.756  0.9979
##  DirNorm - IndBE        0.00490 0.0943 3640   0.052  1.0000
##  IndNorm - Direct       0.00000 0.0936 3640   0.000  1.0000
##  IndNorm - Indirect    -0.00294 0.0938 3640  -0.031  1.0000
##  IndNorm - DirNormBE    0.13971 0.0936 3640   1.493  0.8590
##  IndNorm - IndNormBE    0.04718 0.0935 3640   0.505  0.9999
##  IndNorm - DirBE        0.01647 0.0936 3640   0.176  1.0000
##  IndNorm - IndBE        0.09231 0.0941 3640   0.981  0.9876
##  Direct - Indirect     -0.00294 0.0938 3640  -0.031  1.0000
##  Direct - DirNormBE     0.13971 0.0936 3640   1.493  0.8590
##  Direct - IndNormBE     0.04718 0.0935 3640   0.505  0.9999
##  Direct - DirBE         0.01647 0.0936 3640   0.176  1.0000
##  Direct - IndBE         0.09231 0.0941 3640   0.981  0.9876
##  Indirect - DirNormBE   0.14264 0.0938 3640   1.521  0.8460
##  Indirect - IndNormBE   0.05012 0.0937 3640   0.535  0.9998
##  Indirect - DirBE       0.01941 0.0938 3640   0.207  1.0000
##  Indirect - IndBE       0.09524 0.0943 3640   1.010  0.9851
##  DirNormBE - IndNormBE -0.09253 0.0935 3640  -0.990  0.9869
##  DirNormBE - DirBE     -0.12324 0.0936 3640  -1.316  0.9267
##  DirNormBE - IndBE     -0.04740 0.0941 3640  -0.504  0.9999
##  IndNormBE - DirBE     -0.03071 0.0936 3640  -0.328  1.0000
##  IndNormBE - IndBE      0.04513 0.0940 3640   0.480  0.9999
##  DirBE - IndBE          0.07584 0.0941 3640   0.806  0.9967
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  SR_34 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3973 -0.7575 -0.2302  0.7225  2.8346 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.757463   0.065543  57.328   <2e-16 ***
## ConditionDirNorm   -0.027265   0.092577  -0.295    0.768    
## ConditionIndNorm    0.119988   0.092351   1.299    0.194    
## ConditionDirect     0.057488   0.092351   0.622    0.534    
## ConditionIndirect   0.139815   0.092577   1.510    0.131    
## ConditionDirNormBE -0.006237   0.092351  -0.068    0.946    
## ConditionIndNormBE  0.020043   0.092295   0.217    0.828    
## ConditionDirBE      0.041063   0.092407   0.444    0.657    
## ConditionIndBE     -0.092049   0.092866  -0.991    0.322    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.314 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.002654,   Adjusted R-squared:  0.0004625 
## F-statistic: 1.211 on 8 and 3640 DF,  p-value: 0.288
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       3.76 0.0655 3640     3.63     3.89
##  DirNorm     3.73 0.0654 3640     3.60     3.86
##  IndNorm     3.88 0.0651 3640     3.75     4.01
##  Direct      3.81 0.0651 3640     3.69     3.94
##  Indirect    3.90 0.0654 3640     3.77     4.03
##  DirNormBE   3.75 0.0651 3640     3.62     3.88
##  IndNormBE   3.78 0.0650 3640     3.65     3.90
##  DirBE       3.80 0.0651 3640     3.67     3.93
##  IndBE       3.67 0.0658 3640     3.54     3.79
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm        0.02726 0.0926 3640   0.295  1.0000
##  NoLab - IndNorm       -0.11999 0.0924 3640  -1.299  0.9318
##  NoLab - Direct        -0.05749 0.0924 3640  -0.622  0.9995
##  NoLab - Indirect      -0.13981 0.0926 3640  -1.510  0.8511
##  NoLab - DirNormBE      0.00624 0.0924 3640   0.068  1.0000
##  NoLab - IndNormBE     -0.02004 0.0923 3640  -0.217  1.0000
##  NoLab - DirBE         -0.04106 0.0924 3640  -0.444  1.0000
##  NoLab - IndBE          0.09205 0.0929 3640   0.991  0.9868
##  DirNorm - IndNorm     -0.14725 0.0922 3640  -1.596  0.8071
##  DirNorm - Direct      -0.08475 0.0922 3640  -0.919  0.9920
##  DirNorm - Indirect    -0.16708 0.0925 3640  -1.807  0.6776
##  DirNorm - DirNormBE   -0.02103 0.0922 3640  -0.228  1.0000
##  DirNorm - IndNormBE   -0.04731 0.0922 3640  -0.513  0.9999
##  DirNorm - DirBE       -0.06833 0.0923 3640  -0.740  0.9982
##  DirNorm - IndBE        0.06478 0.0928 3640   0.698  0.9988
##  IndNorm - Direct       0.06250 0.0920 3640   0.679  0.9990
##  IndNorm - Indirect    -0.01983 0.0922 3640  -0.215  1.0000
##  IndNorm - DirNormBE    0.12623 0.0920 3640   1.372  0.9085
##  IndNorm - IndNormBE    0.09994 0.0920 3640   1.087  0.9762
##  IndNorm - DirBE        0.07893 0.0921 3640   0.857  0.9950
##  IndNorm - IndBE        0.21204 0.0925 3640   2.292  0.3473
##  Direct - Indirect     -0.08233 0.0922 3640  -0.893  0.9934
##  Direct - DirNormBE     0.06373 0.0920 3640   0.693  0.9989
##  Direct - IndNormBE     0.03744 0.0920 3640   0.407  1.0000
##  Direct - DirBE         0.01643 0.0921 3640   0.178  1.0000
##  Direct - IndBE         0.14954 0.0925 3640   1.616  0.7962
##  Indirect - DirNormBE   0.14605 0.0922 3640   1.583  0.8141
##  Indirect - IndNormBE   0.11977 0.0922 3640   1.299  0.9318
##  Indirect - DirBE       0.09875 0.0923 3640   1.070  0.9784
##  Indirect - IndBE       0.23186 0.0928 3640   2.500  0.2322
##  DirNormBE - IndNormBE -0.02628 0.0920 3640  -0.286  1.0000
##  DirNormBE - DirBE     -0.04730 0.0921 3640  -0.514  0.9999
##  DirNormBE - IndBE      0.08581 0.0925 3640   0.927  0.9915
##  IndNormBE - DirBE     -0.02102 0.0920 3640  -0.228  1.0000
##  IndNormBE - IndBE      0.11209 0.0925 3640   1.212  0.9541
##  DirBE - IndBE          0.13311 0.0926 3640   1.438  0.8832
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  ISM_12 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.3563 -0.6785  0.1437  0.7908  1.8738 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.218905   0.054314  59.265   <2e-16 ***
## ConditionDirNorm   -0.009747   0.076716  -0.127   0.8989    
## ConditionIndNorm    0.047026   0.076528   0.614   0.5389    
## ConditionDirect     0.066634   0.076528   0.871   0.3840    
## ConditionIndirect   0.073174   0.076716   0.954   0.3402    
## ConditionDirNormBE -0.092680   0.076528  -1.211   0.2260    
## ConditionIndNormBE -0.040421   0.076482  -0.529   0.5972    
## ConditionDirBE      0.137360   0.076575   1.794   0.0729 .  
## ConditionIndBE     -0.049733   0.076955  -0.646   0.5182    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.089 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.003914,   Adjusted R-squared:  0.001725 
## F-statistic: 1.788 on 8 and 3640 DF,  p-value: 0.07455
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       3.22 0.0543 3640     3.11     3.33
##  DirNorm     3.21 0.0542 3640     3.10     3.32
##  IndNorm     3.27 0.0539 3640     3.16     3.37
##  Direct      3.29 0.0539 3640     3.18     3.39
##  Indirect    3.29 0.0542 3640     3.19     3.40
##  DirNormBE   3.13 0.0539 3640     3.02     3.23
##  IndNormBE   3.18 0.0538 3640     3.07     3.28
##  DirBE       3.36 0.0540 3640     3.25     3.46
##  IndBE       3.17 0.0545 3640     3.06     3.28
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm        0.00975 0.0767 3640   0.127  1.0000
##  NoLab - IndNorm       -0.04703 0.0765 3640  -0.614  0.9995
##  NoLab - Direct        -0.06663 0.0765 3640  -0.871  0.9944
##  NoLab - Indirect      -0.07317 0.0767 3640  -0.954  0.9897
##  NoLab - DirNormBE      0.09268 0.0765 3640   1.211  0.9543
##  NoLab - IndNormBE      0.04042 0.0765 3640   0.529  0.9998
##  NoLab - DirBE         -0.13736 0.0766 3640  -1.794  0.6865
##  NoLab - IndBE          0.04973 0.0770 3640   0.646  0.9993
##  DirNorm - IndNorm     -0.05677 0.0764 3640  -0.743  0.9982
##  DirNorm - Direct      -0.07638 0.0764 3640  -0.999  0.9860
##  DirNorm - Indirect    -0.08292 0.0766 3640  -1.082  0.9768
##  DirNorm - DirNormBE    0.08293 0.0764 3640   1.085  0.9765
##  DirNorm - IndNormBE    0.03067 0.0764 3640   0.402  1.0000
##  DirNorm - DirBE       -0.14711 0.0765 3640  -1.923  0.5974
##  DirNorm - IndBE        0.03999 0.0769 3640   0.520  0.9999
##  IndNorm - Direct      -0.01961 0.0762 3640  -0.257  1.0000
##  IndNorm - Indirect    -0.02615 0.0764 3640  -0.342  1.0000
##  IndNorm - DirNormBE    0.13971 0.0762 3640   1.832  0.6605
##  IndNorm - IndNormBE    0.08745 0.0762 3640   1.148  0.9668
##  IndNorm - DirBE       -0.09033 0.0763 3640  -1.184  0.9600
##  IndNorm - IndBE        0.09676 0.0767 3640   1.262  0.9421
##  Direct - Indirect     -0.00654 0.0764 3640  -0.086  1.0000
##  Direct - DirNormBE     0.15931 0.0762 3640   2.090  0.4806
##  Direct - IndNormBE     0.10706 0.0762 3640   1.405  0.8963
##  Direct - DirBE        -0.07073 0.0763 3640  -0.927  0.9915
##  Direct - IndBE         0.11637 0.0767 3640   1.518  0.8475
##  Indirect - DirNormBE   0.16585 0.0764 3640   2.170  0.4257
##  Indirect - IndNormBE   0.11360 0.0764 3640   1.487  0.8619
##  Indirect - DirBE      -0.06419 0.0765 3640  -0.839  0.9957
##  Indirect - IndBE       0.12291 0.0769 3640   1.599  0.8056
##  DirNormBE - IndNormBE -0.05226 0.0762 3640  -0.686  0.9990
##  DirNormBE - DirBE     -0.23004 0.0763 3640  -3.015  0.0646
##  DirNormBE - IndBE     -0.04295 0.0767 3640  -0.560  0.9998
##  IndNormBE - DirBE     -0.17778 0.0762 3640  -2.332  0.3231
##  IndNormBE - IndBE      0.00931 0.0766 3640   0.122  1.0000
##  DirBE - IndBE          0.18709 0.0767 3640   2.439  0.2632
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  ISM_34 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.65233 -0.54579 -0.03465  0.84767  1.52632 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.532338   0.045384  77.832   <2e-16 ***
## ConditionDirNorm    0.002315   0.064103   0.036   0.9712    
## ConditionIndNorm    0.093887   0.063946   1.468   0.1421    
## ConditionDirect    -0.022534   0.063946  -0.352   0.7246    
## ConditionIndirect   0.013454   0.064103   0.210   0.8338    
## ConditionDirNormBE -0.026211   0.063946  -0.410   0.6819    
## ConditionIndNormBE -0.022558   0.063907  -0.353   0.7241    
## ConditionDirBE      0.119996   0.063985   1.875   0.0608 .  
## ConditionIndBE     -0.058654   0.064303  -0.912   0.3618    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9099 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.003682,   Adjusted R-squared:  0.001493 
## F-statistic: 1.682 on 8 and 3640 DF,  p-value: 0.09758
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       3.53 0.0454 3640     3.44     3.62
##  DirNorm     3.53 0.0453 3640     3.45     3.62
##  IndNorm     3.63 0.0450 3640     3.54     3.71
##  Direct      3.51 0.0450 3640     3.42     3.60
##  Indirect    3.55 0.0453 3640     3.46     3.63
##  DirNormBE   3.51 0.0450 3640     3.42     3.59
##  IndNormBE   3.51 0.0450 3640     3.42     3.60
##  DirBE       3.65 0.0451 3640     3.56     3.74
##  IndBE       3.47 0.0456 3640     3.38     3.56
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast               estimate     SE   df t.ratio p.value
##  NoLab - DirNorm       -0.002315 0.0641 3640  -0.036  1.0000
##  NoLab - IndNorm       -0.093887 0.0639 3640  -1.468  0.8703
##  NoLab - Direct         0.022534 0.0639 3640   0.352  1.0000
##  NoLab - Indirect      -0.013454 0.0641 3640  -0.210  1.0000
##  NoLab - DirNormBE      0.026211 0.0639 3640   0.410  1.0000
##  NoLab - IndNormBE      0.022558 0.0639 3640   0.353  1.0000
##  NoLab - DirBE         -0.119996 0.0640 3640  -1.875  0.6309
##  NoLab - IndBE          0.058654 0.0643 3640   0.912  0.9924
##  DirNorm - IndNorm     -0.091572 0.0639 3640  -1.434  0.8849
##  DirNorm - Direct       0.024850 0.0639 3640   0.389  1.0000
##  DirNorm - Indirect    -0.011139 0.0640 3640  -0.174  1.0000
##  DirNorm - DirNormBE    0.028526 0.0639 3640   0.447  1.0000
##  DirNorm - IndNormBE    0.024874 0.0638 3640   0.390  1.0000
##  DirNorm - DirBE       -0.117681 0.0639 3640  -1.841  0.6543
##  DirNorm - IndBE        0.060969 0.0642 3640   0.949  0.9900
##  IndNorm - Direct       0.116422 0.0637 3640   1.827  0.6639
##  IndNorm - Indirect     0.080433 0.0639 3640   1.259  0.9428
##  IndNorm - DirNormBE    0.120098 0.0637 3640   1.885  0.6242
##  IndNorm - IndNormBE    0.116446 0.0637 3640   1.829  0.6628
##  IndNorm - DirBE       -0.026109 0.0637 3640  -0.410  1.0000
##  IndNorm - IndBE        0.152541 0.0641 3640   2.381  0.2947
##  Direct - Indirect     -0.035988 0.0639 3640  -0.563  0.9998
##  Direct - DirNormBE     0.003676 0.0637 3640   0.058  1.0000
##  Direct - IndNormBE     0.000024 0.0637 3640   0.000  1.0000
##  Direct - DirBE        -0.142530 0.0637 3640  -2.236  0.3824
##  Direct - IndBE         0.036120 0.0641 3640   0.564  0.9998
##  Indirect - DirNormBE   0.039665 0.0639 3640   0.621  0.9995
##  Indirect - IndNormBE   0.036012 0.0638 3640   0.564  0.9998
##  Indirect - DirBE      -0.106542 0.0639 3640  -1.667  0.7667
##  Indirect - IndBE       0.072108 0.0642 3640   1.123  0.9709
##  DirNormBE - IndNormBE -0.003652 0.0637 3640  -0.057  1.0000
##  DirNormBE - DirBE     -0.146207 0.0637 3640  -2.294  0.3462
##  DirNormBE - IndBE      0.032443 0.0641 3640   0.506  0.9999
##  IndNormBE - DirBE     -0.142554 0.0637 3640  -2.238  0.3813
##  IndNormBE - IndBE      0.036096 0.0640 3640   0.564  0.9998
##  DirBE - IndBE          0.178650 0.0641 3640   2.787  0.1195
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------
## 
## Results for  ISM_56 
## 
## Model Summary
## 
## Call:
## lm(formula = formula, data = climate_labels)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.57985 -0.52206  0.01716  0.99381  1.63235 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.502488   0.054935  63.757   <2e-16 ***
## ConditionDirNorm   -0.011151   0.077594  -0.144   0.8857    
## ConditionIndNorm    0.019571   0.077404   0.253   0.8004    
## ConditionDirect    -0.019644   0.077404  -0.254   0.7997    
## ConditionIndirect   0.003701   0.077594   0.048   0.9620    
## ConditionDirNormBE -0.134841   0.077404  -1.742   0.0816 .  
## ConditionIndNormBE -0.045275   0.077357  -0.585   0.5584    
## ConditionDirBE      0.077365   0.077451   0.999   0.3179    
## ConditionIndBE     -0.097726   0.077836  -1.256   0.2094    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.101 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.002921,   Adjusted R-squared:  0.0007293 
## F-statistic: 1.333 on 8 and 3640 DF,  p-value: 0.222
## 
## 
## 
##  Pairwise Comparisons
## $emmeans
##  Condition emmean     SE   df lower.CL upper.CL
##  NoLab       3.50 0.0549 3640     3.39     3.61
##  DirNorm     3.49 0.0548 3640     3.38     3.60
##  IndNorm     3.52 0.0545 3640     3.42     3.63
##  Direct      3.48 0.0545 3640     3.38     3.59
##  Indirect    3.51 0.0548 3640     3.40     3.61
##  DirNormBE   3.37 0.0545 3640     3.26     3.47
##  IndNormBE   3.46 0.0545 3640     3.35     3.56
##  DirBE       3.58 0.0546 3640     3.47     3.69
##  IndBE       3.40 0.0551 3640     3.30     3.51
## 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast              estimate     SE   df t.ratio p.value
##  NoLab - DirNorm        0.01115 0.0776 3640   0.144  1.0000
##  NoLab - IndNorm       -0.01957 0.0774 3640  -0.253  1.0000
##  NoLab - Direct         0.01964 0.0774 3640   0.254  1.0000
##  NoLab - Indirect      -0.00370 0.0776 3640  -0.048  1.0000
##  NoLab - DirNormBE      0.13484 0.0774 3640   1.742  0.7203
##  NoLab - IndNormBE      0.04527 0.0774 3640   0.585  0.9997
##  NoLab - DirBE         -0.07737 0.0775 3640  -0.999  0.9861
##  NoLab - IndBE          0.09773 0.0778 3640   1.256  0.9438
##  DirNorm - IndNorm     -0.03072 0.0773 3640  -0.397  1.0000
##  DirNorm - Direct       0.00849 0.0773 3640   0.110  1.0000
##  DirNorm - Indirect    -0.01485 0.0775 3640  -0.192  1.0000
##  DirNorm - DirNormBE    0.12369 0.0773 3640   1.600  0.8052
##  DirNorm - IndNormBE    0.03412 0.0773 3640   0.442  1.0000
##  DirNorm - DirBE       -0.08852 0.0774 3640  -1.144  0.9674
##  DirNorm - IndBE        0.08657 0.0777 3640   1.114  0.9723
##  IndNorm - Direct       0.03922 0.0771 3640   0.509  0.9999
##  IndNorm - Indirect     0.01587 0.0773 3640   0.205  1.0000
##  IndNorm - DirNormBE    0.15441 0.0771 3640   2.002  0.5418
##  IndNorm - IndNormBE    0.06485 0.0771 3640   0.841  0.9956
##  IndNorm - DirBE       -0.05779 0.0772 3640  -0.749  0.9980
##  IndNorm - IndBE        0.11730 0.0776 3640   1.513  0.8500
##  Direct - Indirect     -0.02334 0.0773 3640  -0.302  1.0000
##  Direct - DirNormBE     0.11520 0.0771 3640   1.494  0.8588
##  Direct - IndNormBE     0.02563 0.0771 3640   0.333  1.0000
##  Direct - DirBE        -0.09701 0.0772 3640  -1.257  0.9433
##  Direct - IndBE         0.07808 0.0776 3640   1.007  0.9853
##  Indirect - DirNormBE   0.13854 0.0773 3640   1.792  0.6876
##  Indirect - IndNormBE   0.04898 0.0773 3640   0.634  0.9994
##  Indirect - DirBE      -0.07366 0.0774 3640  -0.952  0.9898
##  Indirect - IndBE       0.10143 0.0777 3640   1.305  0.9302
##  DirNormBE - IndNormBE -0.08957 0.0771 3640  -1.162  0.9642
##  DirNormBE - DirBE     -0.21221 0.0772 3640  -2.750  0.1310
##  DirNormBE - IndBE     -0.03711 0.0776 3640  -0.479  0.9999
##  IndNormBE - DirBE     -0.12264 0.0771 3640  -1.590  0.8104
##  IndNormBE - IndBE      0.05245 0.0775 3640   0.677  0.9991
##  DirBE - IndBE          0.17509 0.0776 3640   2.256  0.3693
## 
## P value adjustment: tukey method for comparing a family of 9 estimates 
## 
## 
## 
## --------------------------------------------------------------------------------

3-way

InteractAnalyses <- subset(climate_labels, Condition != "NoLab")
InteractAnalyses$Condition <- as.factor(InteractAnalyses$Condition)


# Harm
InteractAnalyses$DN <- ifelse (InteractAnalyses$Condition=="DirNorm", 1,0)
InteractAnalyses$D <- ifelse(InteractAnalyses$Condition=="Direct", 1,0)
InteractAnalyses$DBE <- ifelse(InteractAnalyses$Condition=="DirBE", 1,0)

# Norm
InteractAnalyses$IN <- ifelse (InteractAnalyses$Condition== "IndNorm", 1,0)
InteractAnalyses$DNBE <- ifelse(InteractAnalyses$Condition=="DirNormBE", 1,0)
InteractAnalyses$INBE <- ifelse (InteractAnalyses$Condition== "IndNormBE", 1,0)

# Broader effort
InteractAnalyses$IBE <- ifelse(InteractAnalyses$Condition=="IndBE", 1,0)

InteractAnalyses$Harm<-InteractAnalyses$D+InteractAnalyses$DN+InteractAnalyses$DBE+InteractAnalyses$DNBE # Direct Harm

InteractAnalyses$Norm<-InteractAnalyses$DN+InteractAnalyses$IN+InteractAnalyses$DNBE+InteractAnalyses$INBE # Norm

InteractAnalyses$BroaderEffort<-InteractAnalyses$DNBE+InteractAnalyses$INBE+InteractAnalyses$DBE+InteractAnalyses$IBE # Broader effort


outcome_vars <- c("Inj", "DN1", "DN2", "DescN", 
                  "HA_12", "HA_34", 
                  "EB_12", "EB_34", "EB_56",
                  "SR_12", "SR_34",
                  "ISM_12", "ISM_34", "ISM_56")

for (outcome in outcome_vars) {
  model <- lm(as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), data = InteractAnalyses)
  print(summary(model))
}
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8948 -0.7218 -0.1275  0.7684  2.3725 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              2.62748    0.04939  53.200  < 2e-16 ***
## Norm                     0.24507    0.06967   3.517 0.000442 ***
## Harm                     0.10414    0.06967   1.495 0.135093    
## BroaderEffort            0.04044    0.07006   0.577 0.563814    
## Norm:Harm               -0.08189    0.09854  -0.831 0.405997    
## Norm:BroaderEffort      -0.22717    0.09866  -2.303 0.021365 *  
## Harm:BroaderEffort      -0.05830    0.09872  -0.591 0.554838    
## Norm:Harm:BroaderEffort  0.07204    0.13937   0.517 0.605250    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9927 on 3239 degrees of freedom
## Multiple R-squared:  0.007982,   Adjusted R-squared:  0.005838 
## F-statistic: 3.723 on 7 and 3239 DF,  p-value: 0.0005085
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5833 -0.5490  0.4167  0.5074  2.5124 
## 
## Coefficients:
##                          Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              4.495050   0.061431  73.173   <2e-16 ***
## Norm                     0.088284   0.086663   1.019    0.308    
## Harm                    -0.002402   0.086663  -0.028    0.978    
## BroaderEffort            0.056329   0.087148   0.646    0.518    
## Norm:Harm               -0.093307   0.122560  -0.761    0.447    
## Norm:BroaderEffort      -0.099320   0.122715  -0.809    0.418    
## Harm:BroaderEffort      -0.042833   0.122790  -0.349    0.727    
## Norm:Harm:BroaderEffort  0.147220   0.173355   0.849    0.396    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.235 on 3239 degrees of freedom
## Multiple R-squared:  0.0007072,  Adjusted R-squared:  -0.001452 
## F-statistic: 0.3274 on 7 and 3239 DF,  p-value: 0.9419
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4662 -1.3652  0.5338  0.6348  2.6348 
## 
## Coefficients:
##                          Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              3.373762   0.070522  47.840   <2e-16 ***
## Norm                    -0.008566   0.099488  -0.086    0.931    
## Harm                     0.045355   0.099488   0.456    0.649    
## BroaderEffort            0.092403   0.100045   0.924    0.356    
## Norm:Harm                0.032518   0.140697   0.231    0.817    
## Norm:BroaderEffort      -0.026227   0.140919  -0.186    0.852    
## Harm:BroaderEffort      -0.111029   0.140962  -0.788    0.431    
## Norm:Harm:BroaderEffort  0.042960   0.199041   0.216    0.829    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.417 on 3238 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0005414,  Adjusted R-squared:  -0.001619 
## F-statistic: 0.2506 on 7 and 3238 DF,  p-value: 0.9722
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -46.833 -19.801   0.263  18.167  60.092 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)               39.908      1.206  33.102  < 2e-16 ***
## Norm                       6.925      1.701   4.071 4.78e-05 ***
## Harm                       1.393      1.700   0.820  0.41246    
## BroaderEffort              3.829      1.709   2.240  0.02516 *  
## Norm:Harm                 -4.070      2.404  -1.693  0.09051 .  
## Norm:BroaderEffort        -6.988      2.408  -2.902  0.00374 ** 
## Harm:BroaderEffort        -1.697      2.408  -0.705  0.48112    
## Norm:Harm:BroaderEffort    5.502      3.401   1.618  0.10577    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.2 on 3234 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.006603,   Adjusted R-squared:  0.004453 
## F-statistic: 3.071 on 7 and 3234 DF,  p-value: 0.003167
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2617 -0.6942 -0.0772  0.8203  1.9228 
## 
## Coefficients:
##                          Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              4.183168   0.062335  67.108   <2e-16 ***
## Norm                     0.036194   0.087938   0.412    0.681    
## Harm                     0.072959   0.087938   0.830    0.407    
## BroaderEffort            0.011067   0.088431   0.125    0.900    
## Norm:Harm               -0.131032   0.124403  -1.053    0.292    
## Norm:BroaderEffort      -0.050723   0.124522  -0.407    0.684    
## Harm:BroaderEffort      -0.005524   0.124597  -0.044    0.965    
## Norm:Harm:BroaderEffort -0.038904   0.175935  -0.221    0.825    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.253 on 3238 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.001926,   Adjusted R-squared:  -0.0002314 
## F-statistic: 0.8928 on 7 and 3238 DF,  p-value: 0.5111
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3096 -0.7221  0.2281  0.7971  1.8064 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              4.21658    0.06503  64.841   <2e-16 ***
## Norm                     0.05180    0.09174   0.565    0.572    
## Harm                     0.08611    0.09174   0.939    0.348    
## BroaderEffort            0.05535    0.09225   0.600    0.549    
## Norm:Harm               -0.13241    0.12978  -1.020    0.308    
## Norm:BroaderEffort      -0.12079    0.12990  -0.930    0.353    
## Harm:BroaderEffort      -0.04846    0.12998  -0.373    0.709    
## Norm:Harm:BroaderEffort  0.08545    0.18354   0.466    0.642    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.307 on 3238 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.00106,    Adjusted R-squared:  -0.0011 
## F-statistic: 0.4906 on 7 and 3238 DF,  p-value: 0.842
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.9484 -0.8247  0.1287  0.6683  2.2364 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              2.83168    0.05249  53.948   <2e-16 ***
## Norm                     0.02739    0.07405   0.370    0.712    
## Harm                     0.03964    0.07405   0.535    0.592    
## BroaderEffort            0.02797    0.07446   0.376    0.707    
## Norm:Harm               -0.13510    0.10472  -1.290    0.197    
## Norm:BroaderEffort      -0.06228    0.10489  -0.594    0.553    
## Harm:BroaderEffort       0.04911    0.10492   0.468    0.640    
## Norm:Harm:BroaderEffort  0.04266    0.14815   0.288    0.773    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.055 on 3238 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.002195,   Adjusted R-squared:  3.803e-05 
## F-statistic: 1.018 on 7 and 3238 DF,  p-value: 0.4164
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.17813 -0.58333 -0.08333  0.43797  1.93936 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              3.10025    0.04646  66.727   <2e-16 ***
## Norm                     0.04926    0.06555   0.752    0.452    
## Harm                    -0.01691    0.06555  -0.258    0.796    
## BroaderEffort            0.06517    0.06591   0.989    0.323    
## Norm:Harm               -0.07195    0.09269  -0.776    0.438    
## Norm:BroaderEffort      -0.11541    0.09284  -1.243    0.214    
## Harm:BroaderEffort       0.02963    0.09287   0.319    0.750    
## Norm:Harm:BroaderEffort  0.02614    0.13113   0.199    0.842    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9339 on 3238 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.002074,   Adjusted R-squared:  -8.291e-05 
## F-statistic: 0.9616 on 7 and 3238 DF,  p-value: 0.4576
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.64109 -0.99632 -0.09559  0.88235  2.50368 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              2.64109    0.05394  48.965   <2e-16 ***
## Norm                    -0.02344    0.07609  -0.308    0.758    
## Harm                    -0.04550    0.07609  -0.598    0.550    
## BroaderEffort           -0.03959    0.07652  -0.517    0.605    
## Norm:Harm                0.01077    0.10761   0.100    0.920    
## Norm:BroaderEffort       0.01240    0.10775   0.115    0.908    
## Harm:BroaderEffort       0.05825    0.10781   0.540    0.589    
## Norm:Harm:BroaderEffort -0.11766    0.15221  -0.773    0.440    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3239 degrees of freedom
## Multiple R-squared:  0.001378,   Adjusted R-squared:  -0.0007802 
## F-statistic: 0.6385 on 7 and 3239 DF,  p-value: 0.7243
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2995 -0.7494 -0.2043  0.7506  2.8431 
## 
## Coefficients:
##                          Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              3.799505   0.067009  56.701   <2e-16 ***
## Norm                    -0.002936   0.094533  -0.031    0.975    
## Harm                    -0.002936   0.094533  -0.031    0.975    
## BroaderEffort           -0.095244   0.095062  -1.002    0.316    
## Norm:Harm               -0.084474   0.133690  -0.632    0.528    
## Norm:BroaderEffort       0.048064   0.133860   0.359    0.720    
## Harm:BroaderEffort       0.078774   0.133941   0.588    0.556    
## Norm:Harm:BroaderEffort -0.083890   0.189099  -0.444    0.657    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.347 on 3239 degrees of freedom
## Multiple R-squared:  0.001389,   Adjusted R-squared:  -0.0007689 
## F-statistic: 0.6437 on 7 and 3239 DF,  p-value: 0.72
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3973 -0.7512 -0.2302  0.7225  2.8346 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              3.89728    0.06566  59.353   <2e-16 ***
## Norm                    -0.01983    0.09263  -0.214   0.8305    
## Harm                    -0.08233    0.09263  -0.889   0.3742    
## BroaderEffort           -0.23186    0.09315  -2.489   0.0129 *  
## Norm:Harm               -0.06493    0.13100  -0.496   0.6202    
## Norm:BroaderEffort       0.13192    0.13117   1.006   0.3146    
## Harm:BroaderEffort       0.21544    0.13125   1.641   0.1008    
## Norm:Harm:BroaderEffort -0.09447    0.18530  -0.510   0.6102    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.32 on 3239 degrees of freedom
## Multiple R-squared:  0.002893,   Adjusted R-squared:  0.0007376 
## F-statistic: 1.342 on 7 and 3239 DF,  p-value: 0.2259
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.3563 -0.6785  0.1437  0.7908  1.8738 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              3.29208    0.05393  61.039   <2e-16 ***
## Norm                    -0.02615    0.07609  -0.344   0.7311    
## Harm                    -0.00654    0.07609  -0.086   0.9315    
## BroaderEffort           -0.12291    0.07651  -1.606   0.1083    
## Norm:Harm               -0.05023    0.10760  -0.467   0.6406    
## Norm:BroaderEffort       0.03546    0.10774   0.329   0.7421    
## Harm:BroaderEffort       0.19363    0.10780   1.796   0.0726 .  
## Norm:Harm:BroaderEffort -0.18912    0.15220  -1.243   0.2141    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3239 degrees of freedom
## Multiple R-squared:  0.004411,   Adjusted R-squared:  0.00226 
## F-statistic:  2.05 on 7 and 3239 DF,  p-value: 0.0456
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.65233 -0.54579 -0.03465  0.68699  1.52632 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              3.54579    0.04529  78.287   <2e-16 ***
## Norm                     0.08043    0.06390   1.259   0.2082    
## Harm                    -0.03599    0.06390  -0.563   0.5733    
## BroaderEffort           -0.07211    0.06425  -1.122   0.2618    
## Norm:Harm               -0.05558    0.09036  -0.615   0.5385    
## Norm:BroaderEffort      -0.04434    0.09048  -0.490   0.6241    
## Harm:BroaderEffort       0.21464    0.09053   2.371   0.0178 *  
## Norm:Harm:BroaderEffort -0.12672    0.12781  -0.991   0.3215    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9104 on 3239 degrees of freedom
## Multiple R-squared:  0.004112,   Adjusted R-squared:  0.001959 
## F-statistic:  1.91 on 7 and 3239 DF,  p-value: 0.06391
## 
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm*Harm*BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.57985 -0.52206  0.01716  0.99381  1.63235 
## 
## Coefficients:
##                          Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              3.506188   0.054844  63.930   <2e-16 ***
## Norm                     0.015871   0.077371   0.205   0.8375    
## Harm                    -0.023345   0.077371  -0.302   0.7629    
## BroaderEffort           -0.101426   0.077804  -1.304   0.1925    
## Norm:Harm               -0.007377   0.109420  -0.067   0.9463    
## Norm:BroaderEffort       0.036580   0.109559   0.334   0.7385    
## Harm:BroaderEffort       0.198436   0.109625   1.810   0.0704 .  
## Norm:Harm:BroaderEffort -0.257279   0.154770  -1.662   0.0965 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.102 on 3239 degrees of freedom
## Multiple R-squared:  0.003215,   Adjusted R-squared:  0.001061 
## F-statistic: 1.492 on 7 and 3239 DF,  p-value: 0.165

2-way

for (outcome in outcome_vars) {
  
  # Norm x Harm
  model1 <- lm(as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
  cat("\n--- Norm x Harm Interaction for", outcome, "---\n")
  print(summary(model1))
  
  # Norm x Broader Effort
  model2 <- lm(as.formula(paste(outcome, "~ Norm * BroaderEffort")), data = InteractAnalyses)
  cat("\n--- Norm x BroaderEffort Interaction for", outcome, "---\n")
  print(summary(model2))
  
  # Harm x Broader Effort
  model3 <- lm(as.formula(paste(outcome, "~ Harm * BroaderEffort")), data = InteractAnalyses)
  cat("\n--- Harm x BroaderEffort Interaction for", outcome, "---\n")
  print(summary(model3))
}
## 
## --- Norm x Harm Interaction for Inj ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8079 -0.7791 -0.1476  0.7209  2.3524 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.64757    0.03508  75.463  < 2e-16 ***
## Norm         0.13150    0.04940   2.662  0.00781 ** 
## Harm         0.07513    0.04943   1.520  0.12867    
## Norm:Harm   -0.04632    0.06979  -0.664  0.50697    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9942 on 3243 degrees of freedom
## Multiple R-squared:  0.003758,   Adjusted R-squared:  0.002837 
## F-statistic: 4.078 on 3 and 3243 DF,  p-value: 0.006689
## 
## 
## --- Norm x BroaderEffort Interaction for Inj ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8836 -0.7038 -0.1798  0.7962  2.3202 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.67980    0.03483  76.935  < 2e-16 ***
## Norm                0.20382    0.04926   4.138  3.6e-05 ***
## BroaderEffort       0.01126    0.04935   0.228  0.81947    
## Norm:BroaderEffort -0.19109    0.06968  -2.743  0.00613 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9926 on 3243 degrees of freedom
## Multiple R-squared:  0.007052,   Adjusted R-squared:  0.006134 
## F-statistic: 7.678 on 3 and 3243 DF,  p-value: 4.125e-05
## 
## 
## --- Harm x BroaderEffort Interaction for Inj ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8128 -0.7506 -0.1770  0.7494  2.3230 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.75062    0.03491  78.787   <2e-16 ***
## Harm                0.06219    0.04937   1.260    0.208    
## BroaderEffort      -0.07364    0.04943  -1.490    0.136    
## Harm:BroaderEffort -0.02138    0.06983  -0.306    0.759    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9948 on 3243 degrees of freedom
## Multiple R-squared:  0.002488,   Adjusted R-squared:  0.001566 
## F-statistic: 2.697 on 3 and 3243 DF,  p-value: 0.04437
## 
## 
## --- Norm x Harm Interaction for DN1 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5618 -0.5230  0.4382  0.5006  2.5006 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.52304    0.04355 103.849   <2e-16 ***
## Norm         0.03877    0.06133   0.632    0.527    
## Harm        -0.02365    0.06137  -0.385    0.700    
## Norm:Harm   -0.01969    0.08664  -0.227    0.820    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.234 on 3243 degrees of freedom
## Multiple R-squared:  0.0003395,  Adjusted R-squared:  -0.0005852 
## F-statistic: 0.3672 on 3 and 3243 DF,  p-value: 0.7767
## 
## 
## --- Norm x BroaderEffort Interaction for DN1 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5447 -0.5357  0.4553  0.5062  2.5062 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.49384    0.04331 103.750   <2e-16 ***
## Norm                0.04187    0.06126   0.684    0.494    
## BroaderEffort       0.03469    0.06137   0.565    0.572    
## Norm:BroaderEffort -0.02573    0.08664  -0.297    0.766    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.234 on 3243 degrees of freedom
## Multiple R-squared:  0.0002441,  Adjusted R-squared:  -0.0006808 
## F-statistic: 0.2639 on 3 and 3243 DF,  p-value: 0.8514
## 
## 
## --- Harm x BroaderEffort Interaction for DN1 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5458 -0.5394  0.4542  0.5099  2.5099 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.539409   0.043313 104.805   <2e-16 ***
## Harm               -0.049261   0.061254  -0.804    0.421    
## BroaderEffort       0.006383   0.061329   0.104    0.917    
## Harm:BroaderEffort  0.031076   0.086639   0.359    0.720    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.234 on 3243 degrees of freedom
## Multiple R-squared:  0.0003052,  Adjusted R-squared:  -0.0006196 
## F-statistic:  0.33 on 3 and 3243 DF,  p-value: 0.8037
## 
## 
## --- Norm x Harm Interaction for DN2 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4421 -1.3983  0.5579  0.6017  2.6017 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.41968    0.05000  68.392   <2e-16 ***
## Norm        -0.02139    0.07043  -0.304    0.761    
## Harm        -0.00986    0.07045  -0.140    0.889    
## Norm:Harm    0.05369    0.09948   0.540    0.589    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.417 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0001299,  Adjusted R-squared:  -0.0007954 
## F-statistic: 0.1404 on 3 and 3242 DF,  p-value: 0.9358
## 
## 
## --- Norm x BroaderEffort Interaction for DN2 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4363 -1.3965  0.5637  0.6035  2.6035 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.396552   0.049722  68.310   <2e-16 ***
## Norm                0.007389   0.070318   0.105    0.916    
## BroaderEffort       0.036451   0.070449   0.517    0.605    
## Norm:BroaderEffort -0.004117   0.099476  -0.041    0.967    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.417 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0001516,  Adjusted R-squared:  -0.0007736 
## F-statistic: 0.1639 on 3 and 3242 DF,  p-value: 0.9207
## 
## 
## --- Harm x BroaderEffort Interaction for DN2 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4486 -1.3695  0.5514  0.6305  2.6305 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.36946    0.04972  67.775   <2e-16 ***
## Harm                0.06158    0.07031   0.876    0.381    
## BroaderEffort       0.07912    0.07042   1.124    0.261    
## Harm:BroaderEffort -0.08929    0.09946  -0.898    0.369    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.417 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0004319,  Adjusted R-squared:  -0.0004931 
## F-statistic: 0.4669 on 3 and 3242 DF,  p-value: 0.7054
## 
## 
## --- Norm x Harm Interaction for DescN ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -45.25 -20.25   0.52  18.19  58.19 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  41.8130     0.8554  48.878  < 2e-16 ***
## Norm          3.4401     1.2053   2.854  0.00434 ** 
## Harm          0.5518     1.2053   0.458  0.64714    
## Norm:Harm    -1.3246     1.7019  -0.778  0.43647    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.23 on 3238 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.003468,   Adjusted R-squared:  0.002544 
## F-statistic: 3.756 on 3 and 3238 DF,  p-value: 0.01045
## 
## 
## --- Norm x BroaderEffort Interaction for DescN ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -45.499 -20.499   0.076  19.391  59.391 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         40.6091     0.8498  47.785  < 2e-16 ***
## Norm                 4.8903     1.2018   4.069 4.83e-05 ***
## BroaderEffort        2.9746     1.2041   2.470   0.0135 *  
## Norm:BroaderEffort  -4.2359     1.7002  -2.491   0.0128 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.2 on 3238 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.005489,   Adjusted R-squared:  0.004568 
## F-statistic: 5.958 on 3 and 3238 DF,  p-value: 0.0004776
## 
## 
## --- Harm x BroaderEffort Interaction for DescN ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -44.119 -19.119   0.881  17.278  57.278 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         43.3877     0.8525  50.894   <2e-16 ***
## Harm                -0.6660     1.2049  -0.553    0.580    
## BroaderEffort        0.3171     1.2071   0.263    0.793    
## Harm:BroaderEffort   1.0803     1.7045   0.634    0.526    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.26 on 3238 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.0004441,  Adjusted R-squared:  -0.000482 
## F-statistic: 0.4796 on 3 and 3238 DF,  p-value: 0.6965
## 
## 
## --- Norm x Harm Interaction for HA_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2589 -0.6995 -0.1190  0.8113  1.8810 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.18867    0.04419  94.777   <2e-16 ***
## Norm         0.01084    0.06223   0.174   0.8617    
## Harm         0.07023    0.06227   1.128   0.2595    
## Norm:Harm   -0.15075    0.08793  -1.714   0.0865 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.252 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.001576,   Adjusted R-squared:  0.0006518 
## F-statistic: 1.706 on 3 and 3242 DF,  p-value: 0.1637
## 
## 
## --- Norm x BroaderEffort Interaction for HA_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2283 -0.7198 -0.1285  0.8095  1.8715 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.21983    0.04396  95.987   <2e-16 ***
## Norm               -0.02932    0.06219  -0.471    0.637    
## BroaderEffort       0.00846    0.06229   0.136    0.892    
## Norm:BroaderEffort -0.07045    0.08795  -0.801    0.423    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.253 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0009786,  Adjusted R-squared:  5.416e-05 
## F-statistic: 1.059 on 3 and 3242 DF,  p-value: 0.3655
## 
## 
## --- Harm x BroaderEffort Interaction for HA_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2090 -0.7014 -0.1693  0.8131  1.8307 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.201355   0.043981  95.527   <2e-16 ***
## Harm                0.007647   0.062217   0.123    0.902    
## BroaderEffort      -0.014473   0.062275  -0.232    0.816    
## Harm:BroaderEffort -0.025203   0.087989  -0.286    0.775    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.253 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0001463,  Adjusted R-squared:  -0.0007789 
## F-statistic: 0.1582 on 3 and 3242 DF,  p-value: 0.9245
## 
## 
## --- Norm x Harm Interaction for HA_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3061 -0.7356  0.2559  0.7922  1.7922 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.244085   0.046104  92.054   <2e-16 ***
## Norm        -0.008467   0.064922  -0.130    0.896    
## Harm         0.062050   0.064961   0.955    0.340    
## Norm:Harm   -0.089900   0.091727  -0.980    0.327    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.306 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0007591,  Adjusted R-squared:  -0.0001655 
## F-statistic: 0.821 on 3 and 3242 DF,  p-value: 0.4821
## 
## 
## --- Norm x BroaderEffort Interaction for HA_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2909 -0.7454  0.2401  0.8017  1.8017 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.25985    0.04585  92.907   <2e-16 ***
## Norm               -0.01448    0.06486  -0.223    0.823    
## BroaderEffort       0.03109    0.06496   0.479    0.632    
## Norm:BroaderEffort -0.07818    0.09173  -0.852    0.394    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.307 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0006545,  Adjusted R-squared:  -0.0002703 
## F-statistic: 0.7077 on 3 and 3242 DF,  p-value: 0.5473
## 
## 
## --- Harm x BroaderEffort Interaction for HA_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2626 -0.7426  0.2485  0.7630  1.7630 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.242611   0.045864  92.504   <2e-16 ***
## Harm                0.020028   0.064882   0.309    0.758    
## BroaderEffort      -0.005606   0.064942  -0.086    0.931    
## Harm:BroaderEffort -0.005499   0.091757  -0.060    0.952    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.307 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  5.499e-05,  Adjusted R-squared:  -0.0008703 
## F-statistic: 0.05943 on 3 and 3242 DF,  p-value: 0.981
## 
## 
## --- Norm x Harm Interaction for EB_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.90982 -0.84191  0.09018  0.65809  2.20751 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.845579   0.037219  76.454   <2e-16 ***
## Norm        -0.003667   0.052426  -0.070    0.944    
## Harm         0.064237   0.052442   1.225    0.221    
## Norm:Harm   -0.113661   0.074050  -1.535    0.125    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.055 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.001564,   Adjusted R-squared:  0.0006405 
## F-statistic: 1.693 on 3 and 3242 DF,  p-value: 0.1663
## 
## 
## --- Norm x BroaderEffort Interaction for EB_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.9045 -0.8229  0.1484  0.6771  2.1884 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.85160    0.03702  77.028   <2e-16 ***
## Norm               -0.04002    0.05235  -0.764    0.445    
## BroaderEffort       0.05287    0.05245   1.008    0.314    
## Norm:BroaderEffort -0.04153    0.07406  -0.561    0.575    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.055 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.001155,   Adjusted R-squared:  0.0002303 
## F-statistic: 1.249 on 3 and 3242 DF,  p-value: 0.2902
## 
## 
## --- Harm x BroaderEffort Interaction for EB_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8847 -0.8420  0.1546  0.6580  2.1823 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.845443   0.037032  76.837   <2e-16 ***
## Harm               -0.027709   0.052371  -0.529    0.597    
## BroaderEffort      -0.003436   0.052452  -0.066    0.948    
## Harm:BroaderEffort  0.070365   0.074087   0.950    0.342    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.055 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0005185,  Adjusted R-squared:  -0.0004064 
## F-statistic: 0.5606 on 3 and 3242 DF,  p-value: 0.641
## 
## 
## --- Norm x Harm Interaction for EB_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.1326 -0.6244 -0.1244  0.4366  1.9366 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.132628   0.032954  95.062   <2e-16 ***
## Norm        -0.008240   0.046418  -0.178    0.859    
## Harm        -0.001953   0.046432  -0.042    0.966    
## Norm:Harm   -0.059011   0.065563  -0.900    0.368    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9338 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0009419,  Adjusted R-squared:  1.741e-05 
## F-statistic: 1.019 on 3 and 3242 DF,  p-value: 0.3831
## 
## 
## --- Norm x BroaderEffort Interaction for EB_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.17184 -0.59175 -0.09175  0.41728  1.91728 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.09175    0.03276  94.367   <2e-16 ***
## Norm                0.01355    0.04633   0.292   0.7700    
## BroaderEffort       0.08009    0.04642   1.725   0.0846 .  
## Norm:BroaderEffort -0.10266    0.06555  -1.566   0.1174    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9336 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.001397,   Adjusted R-squared:  0.0004729 
## F-statistic: 1.512 on 3 and 3242 DF,  p-value: 0.2094
## 
## 
## --- Harm x BroaderEffort Interaction for EB_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.1320 -0.6221 -0.1221  0.4280  1.9280 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.12500    0.03277  95.346   <2e-16 ***
## Harm               -0.05296    0.04635  -1.142    0.253    
## BroaderEffort       0.00697    0.04642   0.150    0.881    
## Harm:BroaderEffort  0.04307    0.06557   0.657    0.511    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.934 on 3242 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.000649,   Adjusted R-squared:  -0.0002758 
## F-statistic: 0.7018 on 3 and 3242 DF,  p-value: 0.5509
## 
## 
## --- Norm x Harm Interaction for EB_56 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.6214 -1.0394 -0.1040  0.8786  2.4606 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.62142    0.03825  68.542   <2e-16 ***
## Norm        -0.01738    0.05385  -0.323    0.747    
## Harm        -0.01651    0.05389  -0.306    0.759    
## Norm:Harm   -0.04812    0.07608  -0.632    0.527    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3243 degrees of freedom
## Multiple R-squared:  0.000838,   Adjusted R-squared:  -8.633e-05 
## F-statistic: 0.9066 on 3 and 3243 DF,  p-value: 0.437
## 
## 
## --- Norm x BroaderEffort Interaction for EB_56 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.6182 -1.0434 -0.1004  0.8818  2.4566 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.61823    0.03804  68.837   <2e-16 ***
## Norm               -0.01786    0.05379  -0.332    0.740    
## BroaderEffort      -0.01029    0.05389  -0.191    0.849    
## Norm:BroaderEffort -0.04663    0.07608  -0.613    0.540    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3243 degrees of freedom
## Multiple R-squared:  0.0007201,  Adjusted R-squared:  -0.0002043 
## F-statistic: 0.779 on 3 and 3243 DF,  p-value: 0.5056
## 
## 
## --- Harm x BroaderEffort Interaction for EB_56 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.62931 -1.05521 -0.08929  0.87069  2.44479 
## 
## Coefficients:
##                      Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.6293103  0.0380376  69.124   <2e-16 ***
## Harm               -0.0400246  0.0537932  -0.744    0.457    
## BroaderEffort      -0.0333945  0.0538598  -0.620    0.535    
## Harm:BroaderEffort -0.0006765  0.0760872  -0.009    0.993    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3243 degrees of freedom
## Multiple R-squared:  0.0005905,  Adjusted R-squared:  -0.0003341 
## F-statistic: 0.6387 on 3 and 3243 DF,  p-value: 0.5901
## 
## 
## --- Norm x Harm Interaction for SR_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2883 -0.7522 -0.1829  0.7478  2.8171 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.75218    0.04751  78.973   <2e-16 ***
## Norm         0.02077    0.06690   0.310    0.756    
## Harm         0.03616    0.06694   0.540    0.589    
## Norm:Harm   -0.12623    0.09451  -1.336    0.182    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.346 on 3243 degrees of freedom
## Multiple R-squared:  0.0008987,  Adjusted R-squared:  -2.552e-05 
## F-statistic: 0.9724 on 3 and 3243 DF,  p-value: 0.4047
## 
## 
## --- Norm x BroaderEffort Interaction for SR_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2980 -0.7426 -0.2032  0.7574  2.7968 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.798030   0.047255  80.374   <2e-16 ***
## Norm               -0.044951   0.066828  -0.673    0.501    
## BroaderEffort      -0.055474   0.066953  -0.829    0.407    
## Norm:BroaderEffort  0.005577   0.094525   0.059    0.953    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.347 on 3243 degrees of freedom
## Multiple R-squared:  0.0006312,  Adjusted R-squared:  -0.0002933 
## F-statistic: 0.6828 on 3 and 3243 DF,  p-value: 0.5625
## 
## 
## --- Harm x BroaderEffort Interaction for SR_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2980 -0.7271 -0.2184  0.7729  2.7816 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.79803    0.04726  80.369   <2e-16 ***
## Harm               -0.04495    0.06683  -0.673    0.501    
## BroaderEffort      -0.07093    0.06691  -1.060    0.289    
## Harm:BroaderEffort  0.03625    0.09453   0.383    0.701    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.347 on 3243 degrees of freedom
## Multiple R-squared:  0.0005295,  Adjusted R-squared:  -0.0003951 
## F-statistic: 0.5727 on 3 and 3243 DF,  p-value: 0.6329
## 
## 
## --- Norm x Harm Interaction for SR_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3274 -0.7821 -0.2408  0.7179  2.7592 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.78207    0.04660  81.160   <2e-16 ***
## Norm         0.04535    0.06562   0.691    0.490    
## Harm         0.02468    0.06566   0.376    0.707    
## Norm:Harm   -0.11134    0.09270  -1.201    0.230    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.321 on 3243 degrees of freedom
## Multiple R-squared:  0.0005991,  Adjusted R-squared:  -0.0003254 
## F-statistic: 0.648 on 3 and 3243 DF,  p-value: 0.5841
## 
## 
## --- Norm x BroaderEffort Interaction for SR_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3559 -0.7644 -0.2326  0.7356  2.7674 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.85591    0.04633  83.233   <2e-16 ***
## Norm               -0.05172    0.06552  -0.789   0.4299    
## BroaderEffort      -0.12328    0.06564  -1.878   0.0604 .  
## Norm:BroaderEffort  0.08348    0.09267   0.901   0.3678    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.32 on 3243 degrees of freedom
## Multiple R-squared:  0.001216,   Adjusted R-squared:  0.0002917 
## F-statistic: 1.316 on 3 and 3243 DF,  p-value: 0.2674
## 
## 
## --- Harm x BroaderEffort Interaction for SR_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3873 -0.7728 -0.2222  0.7272  2.7778 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.88732    0.04631  83.948   <2e-16 ***
## Harm               -0.11453    0.06549  -1.749   0.0804 .  
## BroaderEffort      -0.16516    0.06557  -2.519   0.0118 *  
## Harm:BroaderEffort  0.16723    0.09263   1.805   0.0711 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.32 on 3243 degrees of freedom
## Multiple R-squared:  0.002092,   Adjusted R-squared:  0.001169 
## F-statistic: 2.266 on 3 and 3243 DF,  p-value: 0.07885
## 
## 
## --- Norm x Harm Interaction for ISM_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.3209 -0.7221  0.1791  0.7779  1.8325 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.231009   0.038267  84.433   <2e-16 ***
## Norm        -0.008854   0.053885  -0.164   0.8695    
## Harm         0.089850   0.053918   1.666   0.0957 .  
## Norm:Harm   -0.144517   0.076122  -1.898   0.0577 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3243 degrees of freedom
## Multiple R-squared:  0.002578,   Adjusted R-squared:  0.001655 
## F-statistic: 2.794 on 3 and 3243 DF,  p-value: 0.03891
## 
## 
## --- Norm x BroaderEffort Interaction for ISM_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.2888 -0.7377  0.2112  0.7623  1.8476 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.28879    0.03806  86.409   <2e-16 ***
## Norm               -0.05111    0.05383  -0.950    0.342    
## BroaderEffort      -0.02515    0.05393  -0.466    0.641    
## Norm:BroaderEffort -0.06015    0.07613  -0.790    0.430    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.085 on 3243 degrees of freedom
## Multiple R-squared:  0.002247,   Adjusted R-squared:  0.001324 
## F-statistic: 2.435 on 3 and 3243 DF,  p-value: 0.06299
## 
## 
## --- Harm x BroaderEffort Interaction for ISM_12 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.2789 -0.7411  0.2211  0.7589  1.8261 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.27894    0.03808  86.107   <2e-16 ***
## Harm               -0.03140    0.05385  -0.583   0.5598    
## BroaderEffort      -0.10505    0.05392  -1.948   0.0515 .  
## Harm:BroaderEffort  0.09862    0.07617   1.295   0.1955    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.085 on 3243 degrees of freedom
## Multiple R-squared:  0.001241,   Adjusted R-squared:  0.0003166 
## F-statistic: 1.343 on 3 and 3243 DF,  p-value: 0.2586
## 
## 
## --- Norm x Harm Interaction for ISM_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.58098 -0.56793 -0.06793  0.70453  1.49004 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.50996    0.03215 109.159   <2e-16 ***
## Norm         0.05797    0.04528   1.280   0.2005    
## Harm         0.07102    0.04531   1.568   0.1171    
## Norm:Harm   -0.11863    0.06396  -1.855   0.0637 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9112 on 3243 degrees of freedom
## Multiple R-squared:  0.0011, Adjusted R-squared:  0.000176 
## F-statistic:  1.19 on 3 and 3243 DF,  p-value: 0.3118
## 
## 
## --- Norm x BroaderEffort Interaction for ISM_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.5807 -0.5639 -0.0639  0.7057  1.4920 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.52771    0.03198 110.318   <2e-16 ***
## Norm                0.05296    0.04522   1.171   0.2417    
## BroaderEffort       0.03619    0.04531   0.799   0.4245    
## Norm:BroaderEffort -0.10890    0.06397  -1.702   0.0888 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9112 on 3243 degrees of freedom
## Multiple R-squared:  0.000996,   Adjusted R-squared:  7.181e-05 
## F-statistic: 1.078 on 3 and 3243 DF,  p-value: 0.3572
## 
## 
## --- Harm x BroaderEffort Interaction for ISM_34 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.58621 -0.57914 -0.07914  0.71092  1.50804 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.58621    0.03196 112.196   <2e-16 ***
## Harm               -0.06404    0.04520  -1.417   0.1567    
## BroaderEffort      -0.09425    0.04526  -2.082   0.0374 *  
## Harm:BroaderEffort  0.15123    0.06394   2.365   0.0181 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9108 on 3243 degrees of freedom
## Multiple R-squared:  0.001864,   Adjusted R-squared:  0.0009411 
## F-statistic: 2.019 on 3 and 3243 DF,  p-value: 0.1091
## 
## 
## --- Norm x Harm Interaction for ISM_56 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * Harm")), data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.53129 -0.53129  0.04421  1.01040  1.57081 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.45579    0.03892  88.800   <2e-16 ***
## Norm         0.03381    0.05480   0.617   0.5374    
## Harm         0.07550    0.05483   1.377   0.1687    
## Norm:Harm   -0.13591    0.07741  -1.756   0.0793 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.103 on 3243 degrees of freedom
## Multiple R-squared:  0.001203,   Adjusted R-squared:  0.0002786 
## F-statistic: 1.302 on 3 and 3243 DF,  p-value: 0.2721
## 
## 
## --- Norm x BroaderEffort Interaction for ISM_56 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Norm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.50677 -0.50677  0.00682  1.00554  1.58752 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.494458   0.038701  90.293   <2e-16 ***
## Norm                0.012315   0.054732   0.225    0.822    
## BroaderEffort      -0.001282   0.054834  -0.023    0.981    
## Norm:BroaderEffort -0.093007   0.077415  -1.201    0.230    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.103 on 3243 degrees of freedom
## Multiple R-squared:  0.001159,   Adjusted R-squared:  0.0002354 
## F-statistic: 1.255 on 3 and 3243 DF,  p-value: 0.2883
## 
## 
## --- Harm x BroaderEffort Interaction for ISM_56 ---
## 
## Call:
## lm(formula = as.formula(paste(outcome, "~ Harm * BroaderEffort")), 
##     data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.51416 -0.51416  0.02638  1.01293  1.56869 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         3.51416    0.03871  90.783   <2e-16 ***
## Harm               -0.02709    0.05474  -0.495    0.621    
## BroaderEffort      -0.08285    0.05481  -1.512    0.131    
## Harm:BroaderEffort  0.06940    0.07743   0.896    0.370    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.103 on 3243 degrees of freedom
## Multiple R-squared:  0.0007342,  Adjusted R-squared:  -0.0001902 
## F-statistic: 0.7943 on 3 and 3243 DF,  p-value: 0.4969

Causalmodel with collapsed condition variables

#HA with harm and EB
summary(lm(HA_12~Harm, data = InteractAnalyses))
## 
## Call:
## lm(formula = HA_12 ~ Harm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.1941 -0.6941 -0.1891  0.8109  1.8109 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.194136   0.031130 134.729   <2e-16 ***
## Harm        -0.005021   0.043984  -0.114    0.909    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.253 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  4.018e-06,  Adjusted R-squared:  -0.0003042 
## F-statistic: 0.01303 on 1 and 3244 DF,  p-value: 0.9091
summary(lm(HA_34~Harm, data = InteractAnalyses))
## 
## Call:
## lm(formula = HA_34 ~ Harm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2571 -0.7398  0.2429  0.7602  1.7602 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.23981    0.03246 130.612   <2e-16 ***
## Harm         0.01726    0.04586   0.376    0.707    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.307 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  4.364e-05,  Adjusted R-squared:  -0.0002646 
## F-statistic: 0.1416 on 1 and 3244 DF,  p-value: 0.7067
summary(lm(EB_56~Harm, data = InteractAnalyses)) 
## 
## Call:
## lm(formula = EB_56 ~ Harm, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.61265 -1.07222 -0.07222  0.88735  2.42778 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.61265    0.02692  97.035   <2e-16 ***
## Harm        -0.04044    0.03804  -1.063    0.288    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3245 degrees of freedom
## Multiple R-squared:  0.0003481,  Adjusted R-squared:  4.009e-05 
## F-statistic:  1.13 on 1 and 3245 DF,  p-value: 0.2878
summary(lm(HA_12~Condition, data = climate_labels))
## 
## Call:
## lm(formula = HA_12 ~ Condition, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2617 -0.6942 -0.0348  0.8203  1.9652 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.03483    0.06234  64.721  < 2e-16 ***
## ConditionDirNorm    0.12646    0.08811   1.435  0.15129    
## ConditionIndNorm    0.18454    0.08784   2.101  0.03572 *  
## ConditionDirect     0.22130    0.08784   2.519  0.01180 *  
## ConditionIndirect   0.14834    0.08806   1.685  0.09214 .  
## ConditionDirNormBE  0.04238    0.08784   0.482  0.62950    
## ConditionIndNormBE  0.14488    0.08779   1.650  0.09895 .  
## ConditionDirBE      0.22684    0.08789   2.581  0.00989 ** 
## ConditionIndBE      0.15941    0.08833   1.805  0.07120 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.25 on 3639 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.003262,   Adjusted R-squared:  0.00107 
## F-statistic: 1.488 on 8 and 3639 DF,  p-value: 0.1558
summary(lm(HA_34~Condition, data = climate_labels))
## 
## Call:
## lm(formula = HA_34 ~ Condition, data = climate_labels)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3096 -0.7221  0.2281  0.8064  1.8731 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         4.12687    0.06536  63.142   <2e-16 ***
## ConditionDirNorm    0.09522    0.09237   1.031   0.3027    
## ConditionIndNorm    0.14152    0.09209   1.537   0.1245    
## ConditionDirect     0.17583    0.09209   1.909   0.0563 .  
## ConditionIndirect   0.08972    0.09232   0.972   0.3312    
## ConditionDirNormBE  0.06676    0.09209   0.725   0.4685    
## ConditionIndNormBE  0.07607    0.09203   0.827   0.4086    
## ConditionDirBE      0.18272    0.09215   1.983   0.0475 *  
## ConditionIndBE      0.14506    0.09260   1.566   0.1173    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.31 on 3639 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.001782,   Adjusted R-squared:  -0.0004124 
## F-statistic: 0.8121 on 8 and 3639 DF,  p-value: 0.5918
summary(lm(EB_56~Condition, data = climate_labels))
## 
## Call:
## lm(formula = EB_56 ~ Condition, data = climate_labels)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.64109 -0.99632 -0.09559  0.87438  2.50368 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         2.625622   0.053819  48.786   <2e-16 ***
## ConditionDirNorm   -0.042701   0.076018  -0.562   0.5743    
## ConditionIndNorm   -0.007975   0.075832  -0.105   0.9163    
## ConditionDirect    -0.030034   0.075832  -0.396   0.6921    
## ConditionIndirect   0.015467   0.076018   0.203   0.8388    
## ConditionDirNormBE -0.129298   0.075832  -1.705   0.0883 .  
## ConditionIndNormBE -0.035157   0.075786  -0.464   0.6427    
## ConditionDirBE     -0.011371   0.075878  -0.150   0.8809    
## ConditionIndBE     -0.024118   0.076255  -0.316   0.7518    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.079 on 3640 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.001331,   Adjusted R-squared:  -0.000864 
## F-statistic: 0.6063 on 8 and 3640 DF,  p-value: 0.7733
#norm items with norms
summary(lm(DescN~Norm, data = InteractAnalyses))
## 
## Call:
## lm(formula = DescN ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -44.867 -19.867   0.133  17.909  57.909 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  42.0909     0.6025  69.859  < 2e-16 ***
## Norm          2.7763     0.8508   3.263  0.00111 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.22 on 3240 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.003276,   Adjusted R-squared:  0.002968 
## F-statistic: 10.65 on 1 and 3240 DF,  p-value: 0.001113
summary(lm(Inj~Norm, data = InteractAnalyses))  
## 
## Call:
## lm(formula = Inj ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.7934 -0.7934 -0.1854  0.7066  2.3146 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.68541    0.02472 108.638  < 2e-16 ***
## Norm         0.10802    0.03490   3.095  0.00198 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9943 on 3245 degrees of freedom
## Multiple R-squared:  0.002944,   Adjusted R-squared:  0.002636 
## F-statistic:  9.58 on 1 and 3245 DF,  p-value: 0.001984
summary(lm(DN1~Norm, data = InteractAnalyses))  
## 
## Call:
## lm(formula = DN1 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5402 -0.5402  0.4598  0.4889  2.4889 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.51112    0.03068 147.054   <2e-16 ***
## Norm         0.02908    0.04331   0.672    0.502    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.234 on 3245 degrees of freedom
## Multiple R-squared:  0.0001389,  Adjusted R-squared:  -0.0001692 
## F-statistic: 0.451 on 1 and 3245 DF,  p-value: 0.5019
summary(lm(DN2~Norm, data = InteractAnalyses))  
## 
## Call:
## lm(formula = DN2 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4202 -1.4147  0.5798  0.5853  2.5853 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 3.414710   0.035216  96.965   <2e-16 ***
## Norm        0.005438   0.049726   0.109    0.913    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.417 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  3.686e-06,  Adjusted R-squared:  -0.0003046 
## F-statistic: 0.01196 on 1 and 3244 DF,  p-value: 0.9129
# dotted lines for norms
summary(lm(HA_12~Norm, data = InteractAnalyses))  
## 
## Call:
## lm(formula = HA_12 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2240 -0.7240 -0.1594  0.8406  1.8406 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.22404    0.03114  135.65   <2e-16 ***
## Norm        -0.06464    0.04397   -1.47    0.142    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.253 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0006659,  Adjusted R-squared:  0.0003578 
## F-statistic: 2.161 on 1 and 3244 DF,  p-value: 0.1416
summary(lm(HA_34~Norm, data = InteractAnalyses))  
## 
## Call:
## lm(formula = HA_34 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2753 -0.7217  0.2247  0.7783  1.7783 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.27534    0.03248 131.650   <2e-16 ***
## Norm        -0.05360    0.04586  -1.169    0.243    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.306 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0004209,  Adjusted R-squared:  0.0001128 
## F-statistic: 1.366 on 1 and 3244 DF,  p-value: 0.2426
summary(lm(EB_12~Norm, data = InteractAnalyses))  
## 
## Call:
## lm(formula = EB_12 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8779 -0.8173  0.1221  0.6827  2.1827 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.87794    0.02622 109.753   <2e-16 ***
## Norm        -0.06068    0.03703  -1.639    0.101    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.055 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0008271,  Adjusted R-squared:  0.0005191 
## F-statistic: 2.685 on 1 and 3244 DF,  p-value: 0.1014
summary(lm(EB_34~Norm, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB_34 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.13164 -0.59398 -0.09398  0.40602  1.90602 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.13164    0.02321 134.902   <2e-16 ***
## Norm        -0.03766    0.03278  -1.149    0.251    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9338 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0004068,  Adjusted R-squared:  9.867e-05 
## F-statistic:  1.32 on 1 and 3244 DF,  p-value: 0.2506
summary(lm(EB_56~Norm, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB_56 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.61310 -1.07182 -0.07182  0.88690  2.42818 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.61310    0.02694  96.993   <2e-16 ***
## Norm        -0.04128    0.03804  -1.085    0.278    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3245 degrees of freedom
## Multiple R-squared:  0.0003628,  Adjusted R-squared:  5.477e-05 
## F-statistic: 1.178 on 1 and 3245 DF,  p-value: 0.2779
#EB for broader effort
summary(lm(EB_12~BroaderEffort, data = InteractAnalyses)) #check the individual items for BE
## 
## Call:
## lm(formula = EB_12 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8634 -0.8316  0.1366  0.6684  2.1684 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.83159    0.02618  108.15   <2e-16 ***
## BroaderEffort  0.03185    0.03704    0.86     0.39    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.055 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0002279,  Adjusted R-squared:  -8.026e-05 
## F-statistic: 0.7396 on 1 and 3244 DF,  p-value: 0.3899
summary(lm(EB_34~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB_34 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.12700 -0.59852 -0.09852  0.40148  1.90148 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.09852    0.02317 133.711   <2e-16 ***
## BroaderEffort  0.02848    0.03278   0.869    0.385    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9339 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0002326,  Adjusted R-squared:  -7.556e-05 
## F-statistic: 0.7548 on 1 and 3244 DF,  p-value: 0.385
summary(lm(EB_56~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB_56 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.60930 -1.07548 -0.07548  0.89070  2.42452 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.60930    0.02689  97.025   <2e-16 ***
## BroaderEffort -0.03382    0.03804  -0.889    0.374    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.084 on 3245 degrees of freedom
## Multiple R-squared:  0.0002436,  Adjusted R-squared:  -6.454e-05 
## F-statistic: 0.7905 on 1 and 3245 DF,  p-value: 0.374
#individual EB items for Broader Effort
summary(lm(EB1~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB1 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.6128 -0.6128 -0.6041  0.3959  2.3959 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.60406    0.03028  86.009   <2e-16 ***
## BroaderEffort  0.00876    0.04283   0.205    0.838    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.22 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  1.289e-05,  Adjusted R-squared:  -0.0002954 
## F-statistic: 0.04183 on 1 and 3244 DF,  p-value: 0.838
summary(lm(EB2~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB2 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.94089 -0.94089  0.05911  1.05911  2.11406 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.94089    0.03166  92.894   <2e-16 ***
## BroaderEffort -0.05494    0.04479  -1.227     0.22    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.276 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0004637,  Adjusted R-squared:  0.0001556 
## F-statistic: 1.505 on 1 and 3244 DF,  p-value: 0.22
summary(lm(EB3~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB3 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.8461 -0.8461  0.1539  1.1539  2.1677 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.84606    0.02865  99.343   <2e-16 ***
## BroaderEffort -0.01375    0.04053  -0.339    0.734    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.155 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  3.55e-05,   Adjusted R-squared:  -0.0002728 
## F-statistic: 0.1152 on 1 and 3244 DF,  p-value: 0.7344
summary(lm(EB4~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB4 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -1.649 -0.649  0.351  1.351  2.422 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.64901    0.03055  86.723   <2e-16 ***
## BroaderEffort -0.07072    0.04321  -1.637    0.102    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.231 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0008249,  Adjusted R-squared:  0.0005169 
## F-statistic: 2.678 on 1 and 3244 DF,  p-value: 0.1018
summary(lm(EB5~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB5 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.6804 -0.6804  0.3196  1.3196  2.3660 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.68042    0.03089  86.779   <2e-16 ***
## BroaderEffort -0.04641    0.04369  -1.062    0.288    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.245 on 3245 degrees of freedom
## Multiple R-squared:  0.0003476,  Adjusted R-squared:  3.954e-05 
## F-statistic: 1.128 on 1 and 3245 DF,  p-value: 0.2882
summary(lm(EB6~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = EB6 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4831 -0.4831  0.5169  1.5169  1.5382 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.46182    0.03140 110.259   <2e-16 ***
## BroaderEffort  0.02123    0.04441   0.478    0.633    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.265 on 3245 degrees of freedom
## Multiple R-squared:  7.045e-05,  Adjusted R-squared:  -0.0002377 
## F-statistic: 0.2286 on 1 and 3245 DF,  p-value: 0.6326
#SR for Broader effort
summary(lm(SR_12~BroaderEffort, data = InteractAnalyses)) 
## 
## Call:
## lm(formula = SR_12 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2756 -0.7227 -0.2227  0.7773  2.7773 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.77555    0.03341 113.014   <2e-16 ***
## BroaderEffort -0.05282    0.04725  -1.118    0.264    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.346 on 3245 degrees of freedom
## Multiple R-squared:  0.0003849,  Adjusted R-squared:  7.683e-05 
## F-statistic: 1.249 on 1 and 3245 DF,  p-value: 0.2637
summary(lm(SR_34~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = SR_34 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3300 -0.7486 -0.2486  0.7514  2.7514 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.83005    0.03275 116.940   <2e-16 ***
## BroaderEffort -0.08144    0.04633  -1.758   0.0789 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.32 on 3245 degrees of freedom
## Multiple R-squared:  0.0009514,  Adjusted R-squared:  0.0006435 
## F-statistic:  3.09 on 1 and 3245 DF,  p-value: 0.07886
#individual SR items for Broader Effort
summary(lm(SR1~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = SR1 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -4.0486 -0.9532  0.0468  1.0468  2.0468 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    5.04865    0.04026 125.402   <2e-16 ***
## BroaderEffort -0.09547    0.05694  -1.677   0.0937 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.622 on 3245 degrees of freedom
## Multiple R-squared:  0.0008655,  Adjusted R-squared:  0.0005576 
## F-statistic: 2.811 on 1 and 3245 DF,  p-value: 0.09372
summary(lm(SR2~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = SR2 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5077 -0.5077  0.4923  1.4923  2.5025 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    4.49754    0.04390 102.452   <2e-16 ***
## BroaderEffort  0.01016    0.06209   0.164     0.87    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.769 on 3245 degrees of freedom
## Multiple R-squared:  8.259e-06,  Adjusted R-squared:  -0.0002999 
## F-statistic: 0.0268 on 1 and 3245 DF,  p-value: 0.87
summary(lm(SR3~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = SR3 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -4.0998 -0.0998 -0.0074  0.9926  1.9926 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    5.09975    0.04059 125.645   <2e-16 ***
## BroaderEffort -0.09236    0.05741  -1.609    0.108    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.636 on 3245 degrees of freedom
## Multiple R-squared:  0.000797,   Adjusted R-squared:  0.000489 
## F-statistic: 2.588 on 1 and 3245 DF,  p-value: 0.1078
summary(lm(SR4~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = SR4 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5102 -0.5102  0.4898  1.4898  2.5603 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    4.43966    0.04396 100.985   <2e-16 ***
## BroaderEffort  0.07051    0.06218   1.134    0.257    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.772 on 3245 degrees of freedom
## Multiple R-squared:  0.0003961,  Adjusted R-squared:  8.803e-05 
## F-statistic: 1.286 on 1 and 3245 DF,  p-value: 0.2569
# Norm and broader effort 
summary(lm(Inj~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = Inj ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.7817 -0.7817 -0.1975  0.7183  2.3025 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.78171    0.02469 112.677   <2e-16 ***
## BroaderEffort -0.08424    0.03492  -2.412   0.0159 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9949 on 3245 degrees of freedom
## Multiple R-squared:  0.00179,    Adjusted R-squared:  0.001483 
## F-statistic:  5.82 on 1 and 3245 DF,  p-value: 0.0159
summary(lm(IN2~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = IN2 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.95012 -0.86128  0.04988  1.04988  2.13872 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    2.95012    0.02775 106.313   <2e-16 ***
## BroaderEffort -0.08884    0.03926  -2.263   0.0237 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.118 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.001576,   Adjusted R-squared:  0.001269 
## F-statistic: 5.122 on 1 and 3244 DF,  p-value: 0.02369
summary(lm(DN1~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = DN1 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5367 -0.5367  0.4633  0.4852  2.4852 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    4.51478    0.03062 147.442   <2e-16 ***
## BroaderEffort  0.02188    0.04331   0.505    0.613    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.234 on 3245 degrees of freedom
## Multiple R-squared:  7.866e-05,  Adjusted R-squared:  -0.0002295 
## F-statistic: 0.2553 on 1 and 3245 DF,  p-value: 0.6134
summary(lm(DN2~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = DN2 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4346 -1.4002  0.5654  0.5998  2.5997 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.40025    0.03515  96.740   <2e-16 ***
## BroaderEffort  0.03440    0.04972   0.692    0.489    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.416 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.0001475,  Adjusted R-squared:  -0.0001607 
## F-statistic: 0.4787 on 1 and 3244 DF,  p-value: 0.4891
summary(lm(DescN~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = DescN ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -43.913 -18.913   1.087  16.946  56.946 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    43.0543     0.6023  71.484   <2e-16 ***
## BroaderEffort   0.8586     0.8520   1.008    0.314    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 24.26 on 3240 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.0003134,  Adjusted R-squared:  4.808e-06 
## F-statistic: 1.016 on 1 and 3240 DF,  p-value: 0.3136
summary(lm(ISM_12~Harm, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ Harm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.2443 -0.7265  0.2557  0.7735  1.7735 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.22654    0.02697 119.647   <2e-16 ***
## Harm         0.01777    0.03810   0.466    0.641    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.085 on 3245 degrees of freedom
## Multiple R-squared:  6.706e-05,  Adjusted R-squared:  -0.0002411 
## F-statistic: 0.2176 on 1 and 3245 DF,  p-value: 0.6409
summary(lm(ISM_34~Harm, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ Harm, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.55071 -0.55071 -0.05071  0.70505  1.46080 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.53920    0.02264  156.30   <2e-16 ***
## Harm         0.01151    0.03199    0.36    0.719    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9114 on 3245 degrees of freedom
## Multiple R-squared:  3.989e-05,  Adjusted R-squared:  -0.0002683 
## F-statistic: 0.1295 on 1 and 3245 DF,  p-value: 0.719
summary(lm(ISM_56~Harm, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ Harm, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.48033 -0.48033  0.02716  1.01967  1.52716 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 3.472840   0.027407 126.714   <2e-16 ***
## Harm        0.007492   0.038718   0.194    0.847    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.103 on 3245 degrees of freedom
## Multiple R-squared:  1.154e-05,  Adjusted R-squared:  -0.0002966 
## F-statistic: 0.03745 on 1 and 3245 DF,  p-value: 0.8466
summary(lm(ISM_12~Norm, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.2763 -0.6949  0.2237  0.8051  1.8051 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.27627    0.02697 121.497   <2e-16 ***
## Norm        -0.08136    0.03807  -2.137   0.0327 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.085 on 3245 degrees of freedom
## Multiple R-squared:  0.001405,   Adjusted R-squared:  0.001098 
## F-statistic: 4.567 on 1 and 3245 DF,  p-value: 0.03266
summary(lm(ISM_34~Norm, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.54574 -0.54574 -0.04574  0.70503  1.45580 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.545735   0.022658 156.491   <2e-16 ***
## Norm        -0.001537   0.031989  -0.048    0.962    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9114 on 3245 degrees of freedom
## Multiple R-squared:  7.11e-07,   Adjusted R-squared:  -0.0003075 
## F-statistic: 0.002307 on 1 and 3245 DF,  p-value: 0.9617
summary(lm(ISM_56~Norm, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ Norm, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.49382 -0.49382  0.04052  1.00618  1.54052 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.49382    0.02742 127.415   <2e-16 ***
## Norm        -0.03434    0.03871  -0.887    0.375    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.103 on 3245 degrees of freedom
## Multiple R-squared:  0.0002423,  Adjusted R-squared:  -6.574e-05 
## F-statistic: 0.7866 on 1 and 3245 DF,  p-value: 0.3752
summary(lm(ISM_12~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.2632 -0.7076  0.2368  0.7924  1.7924 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.26324    0.02693  121.19   <2e-16 ***
## BroaderEffort -0.05560    0.03809   -1.46    0.144    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.085 on 3245 degrees of freedom
## Multiple R-squared:  0.0006563,  Adjusted R-squared:  0.0003484 
## F-statistic: 2.131 on 1 and 3245 DF,  p-value: 0.1444
summary(lm(ISM_34~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.55419 -0.55419 -0.05419  0.70504  1.46426 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.55419    0.02261 157.162   <2e-16 ***
## BroaderEffort -0.01845    0.03199  -0.577    0.564    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9113 on 3245 degrees of freedom
## Multiple R-squared:  0.0001025,  Adjusted R-squared:  -0.0002056 
## F-statistic: 0.3327 on 1 and 3245 DF,  p-value: 0.5641
summary(lm(ISM_56~BroaderEffort, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ BroaderEffort, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.50062 -0.50062  0.04744  0.99938  1.54744 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.50062    0.02737 127.914   <2e-16 ***
## BroaderEffort -0.04806    0.03871  -1.242    0.214    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.103 on 3245 degrees of freedom
## Multiple R-squared:  0.0004748,  Adjusted R-squared:  0.0001668 
## F-statistic: 1.541 on 1 and 3245 DF,  p-value: 0.2145

correlational for control group only: beliefs to coalescent beliefs

DiscreteBeliefs <- climate_labels[, c("Inj", "DN1", "DN2", "DescN", 
                                      "HA_12", "HA_34", 
                                      "EB_12", "EB_34", "EB_56",
                                      "SR_12", "SR_34")]

CoalescentBeliefs <- climate_labels[, c("ISM_12", "ISM_34", "ISM_56")]

cor_matrix_control <- cor(DiscreteBeliefs, CoalescentBeliefs, use = "complete.obs")
print(cor_matrix_control) 
##          ISM_12    ISM_34    ISM_56
## Inj   0.4216142 0.2215406 0.4012560
## DN1   0.2122344 0.1582244 0.2295021
## DN2   0.3904008 0.1720420 0.4021979
## DescN 0.3291507 0.1890209 0.3278693
## HA_12 0.6110485 0.3831671 0.6431643
## HA_34 0.6217566 0.3929884 0.6624370
## EB_12 0.5688091 0.3377211 0.5922310
## EB_34 0.5337360 0.3514100 0.5849420
## EB_56 0.5103852 0.2463809 0.4769093
## SR_12 0.5094293 0.3727764 0.5302526
## SR_34 0.5112970 0.3790889 0.5401448
cor_matrix_control <- corr.test(DiscreteBeliefs, CoalescentBeliefs)
cor_matrix_control
## Call:corr.test(x = DiscreteBeliefs, y = CoalescentBeliefs)
## Correlation matrix 
##       ISM_12 ISM_34 ISM_56
## Inj     0.42   0.22   0.40
## DN1     0.21   0.16   0.23
## DN2     0.39   0.17   0.40
## DescN   0.33   0.19   0.33
## HA_12   0.61   0.38   0.64
## HA_34   0.62   0.39   0.66
## EB_12   0.57   0.34   0.59
## EB_34   0.53   0.35   0.58
## EB_56   0.51   0.25   0.48
## SR_12   0.51   0.37   0.53
## SR_34   0.51   0.38   0.54
## Sample Size 
##       ISM_12 ISM_34 ISM_56
## Inj     3649   3649   3649
## DN1     3649   3649   3649
## DN2     3648   3648   3648
## DescN   3644   3644   3644
## HA_12   3648   3648   3648
## HA_34   3648   3648   3648
## EB_12   3648   3648   3648
## EB_34   3648   3648   3648
## EB_56   3649   3649   3649
## SR_12   3649   3649   3649
## SR_34   3649   3649   3649
## These are the unadjusted probability values.
##   The probability values  adjusted for multiple tests are in the p.adj object. 
##       ISM_12 ISM_34 ISM_56
## Inj        0      0      0
## DN1        0      0      0
## DN2        0      0      0
## DescN      0      0      0
## HA_12      0      0      0
## HA_34      0      0      0
## EB_12      0      0      0
## EB_34      0      0      0
## EB_56      0      0      0
## SR_12      0      0      0
## SR_34      0      0      0
## 
##  To see confidence intervals of the correlations, print with the short=FALSE option
CoalescentBeliefs
## # A tibble: 3,651 × 3
##    ISM_12 ISM_34 ISM_56
##     <dbl>  <dbl>  <dbl>
##  1    2.5    4.5    3  
##  2    3      4.5    3  
##  3    2.5    3      3  
##  4    3      2.5    2.5
##  5    3.5    3.5    3.5
##  6    4      4      4  
##  7    4.5    4.5    4  
##  8    4      3      3.5
##  9    4      5      4.5
## 10    2.5    4      4.5
## # ℹ 3,641 more rows
corrplot(as.matrix(cor_matrix_control$r), method = "circle", p.mat = cor_matrix_control$p, insig = "p-value",
         tl.col = "black", tl.srt = 45, 
         title = "Correlation Plot of Control DB on CB", mar = c(0, 0, 1, 0)) 

#discrete beliefs and coalescent beliefs

#harm attribution
summary(lm(ISM_12~HA_12, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ HA_12, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.1898 -0.6065  0.1019  0.6297  3.4491 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.02314    0.05278   19.39   <2e-16 ***
## HA_12        0.52777    0.01206   43.75   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.861 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.371,  Adjusted R-squared:  0.3709 
## F-statistic:  1914 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~HA_12, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ HA_12, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.5055 -0.5055  0.0635  0.5635  3.3393 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.09176    0.05162   21.15   <2e-16 ***
## HA_12        0.56895    0.01180   48.22   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.842 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.4175, Adjusted R-squared:  0.4173 
## F-statistic:  2325 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_12~HA_34, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ HA_34, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.1366 -0.6075 -0.0639  0.6352  3.4361 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.04939    0.05091   20.61   <2e-16 ***
## HA_34        0.51454    0.01145   44.92   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8524 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.3835, Adjusted R-squared:  0.3833 
## F-statistic:  2018 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~HA_34, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ HA_34, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.4551 -0.5585  0.1036  0.6036  3.3381 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.10326    0.04942   22.33   <2e-16 ***
## HA_34        0.55863    0.01112   50.25   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8273 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.4377, Adjusted R-squared:  0.4375 
## F-statistic:  2525 on 1 and 3244 DF,  p-value: < 2.2e-16
#norms
summary(lm(ISM_12~DN2, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ DN2, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -3.00198 -0.70520  0.09158  0.79480  2.48193 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.22128    0.04589   48.40   <2e-16 ***
## DN2          0.29678    0.01241   23.92   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.001 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:   0.15,  Adjusted R-squared:  0.1497 
## F-statistic: 572.3 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_34~DN2, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ DN2, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.8208 -0.6073 -0.1073  0.7129  1.7129 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.18034    0.04121  77.169   <2e-16 ***
## DN2          0.10675    0.01114   9.581   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8989 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.02752,    Adjusted R-squared:  0.02722 
## F-statistic:  91.8 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~DN2, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ DN2, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2735 -0.6560  0.0352  0.7703  2.2703 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.42095    0.04642   52.15   <2e-16 ***
## DN2          0.30876    0.01255   24.60   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.013 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1572, Adjusted R-squared:  0.157 
## F-statistic: 605.3 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_12~Inj, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ Inj, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2620 -0.6725  0.1005  0.6921  2.5546 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.99121    0.05071   39.27   <2e-16 ***
## Inj          0.45417    0.01740   26.11   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9867 on 3245 degrees of freedom
## Multiple R-squared:  0.1736, Adjusted R-squared:  0.1733 
## F-statistic: 681.6 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_34~Inj, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ Inj, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.70201 -0.59874 -0.09874  0.69473  1.81431 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.97916    0.04563   65.30   <2e-16 ***
## Inj          0.20653    0.01565   13.19   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8879 on 3245 degrees of freedom
## Multiple R-squared:  0.05092,    Adjusted R-squared:  0.05062 
## F-statistic: 174.1 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~Inj, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ Inj, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.4733 -0.6505  0.0700  0.7904  2.2904 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.26864    0.05200   43.63   <2e-16 ***
## Inj          0.44092    0.01784   24.71   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.012 on 3245 degrees of freedom
## Multiple R-squared:  0.1584, Adjusted R-squared:  0.1582 
## F-statistic: 610.8 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_12~DN1, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ DN1, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.68613 -0.63969  0.04245  0.76924  2.40674 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.41111    0.07086   34.03   <2e-16 ***
## DN1          0.18215    0.01511   12.06   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.062 on 3245 degrees of freedom
## Multiple R-squared:  0.04288,    Adjusted R-squared:  0.04259 
## F-statistic: 145.4 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_34~DN1, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ DN1, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.82123 -0.59792 -0.09792  0.73704  1.84869 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  3.03965    0.06012  50.562   <2e-16 ***
## DN1          0.11165    0.01282   8.712   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9009 on 3245 degrees of freedom
## Multiple R-squared:  0.02285,    Adjusted R-squared:  0.02255 
## F-statistic:  75.9 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~DN1, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ DN1, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.9693 -0.6728  0.1281  0.9290  2.2254 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.57545    0.07176   35.89   <2e-16 ***
## DN1          0.19912    0.01530   13.02   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.075 on 3245 degrees of freedom
## Multiple R-squared:  0.04962,    Adjusted R-squared:  0.04932 
## F-statistic: 169.4 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_12~DescN, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ DescN, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.91889 -0.68117  0.03561  0.75588  2.40272 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 2.5972762  0.0369793   70.24   <2e-16 ***
## DescN       0.0146846  0.0007427   19.77   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.026 on 3240 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.1077, Adjusted R-squared:  0.1074 
## F-statistic: 390.9 on 1 and 3240 DF,  p-value: < 2.2e-16
summary(lm(ISM_34~DescN, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ DescN, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.59483 -0.59483 -0.09483  0.70894  1.78489 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 3.2151145  0.0321916   99.87   <2e-16 ***
## DescN       0.0075944  0.0006466   11.75   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8928 on 3240 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.04084,    Adjusted R-squared:  0.04055 
## F-statistic:   138 on 1 and 3240 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~DescN, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ DescN, data = InteractAnalyses)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -3.03335 -0.69404  0.04296  0.81405  2.18747 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 2.8125325  0.0374727   75.06   <2e-16 ***
## DescN       0.0152602  0.0007526   20.28   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.039 on 3240 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.1126, Adjusted R-squared:  0.1123 
## F-statistic: 411.1 on 1 and 3240 DF,  p-value: < 2.2e-16
#efficacy beliefs
summary(lm(ISM_12~EB_12, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ EB_12, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.4887 -0.6154  0.0490  0.5935  2.8401 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.57777    0.04522   34.89   <2e-16 ***
## EB_12        0.58218    0.01489   39.09   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8951 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.3202, Adjusted R-squared:   0.32 
## F-statistic:  1528 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~EB_12, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ EB_12, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.8022 -0.5701  0.1219  0.6219  2.6620 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.72192    0.04502   38.25   <2e-16 ***
## EB_12        0.61605    0.01483   41.55   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.891 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.3474, Adjusted R-squared:  0.3472 
## F-statistic:  1727 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_34~EB_34, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ EB_34, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.1940 -0.5064 -0.0064  0.6498  2.1813 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.47484    0.05211   47.49   <2e-16 ***
## EB_34        0.34384    0.01604   21.44   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.853 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.1241, Adjusted R-squared:  0.1239 
## F-statistic: 459.8 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~EB_34, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ EB_34, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.7725 -0.5552  0.0710  0.6013  2.9751 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.33798    0.05481   24.41   <2e-16 ***
## EB_34        0.68690    0.01687   40.73   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8972 on 3244 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.3383, Adjusted R-squared:  0.3381 
## F-statistic:  1659 on 1 and 3244 DF,  p-value: < 2.2e-16
summary(lm(ISM_12~EB_56, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ EB_56, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.4695 -0.6755  0.0556  0.5808  2.5808 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.90663    0.04243   44.93   <2e-16 ***
## EB_56        0.51258    0.01510   33.94   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9325 on 3245 degrees of freedom
## Multiple R-squared:  0.262,  Adjusted R-squared:  0.2618 
## F-statistic:  1152 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~EB_56, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ EB_56, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.6418 -0.6739  0.0681  0.6002  2.2941 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.22196    0.04416   50.32   <2e-16 ***
## EB_56        0.48397    0.01572   30.79   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9704 on 3245 degrees of freedom
## Multiple R-squared:  0.2261, Adjusted R-squared:  0.2259 
## F-statistic: 948.2 on 1 and 3245 DF,  p-value: < 2.2e-16
#system responsibility beliefs
summary(lm(ISM_34~SR_12, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ SR_12, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2393 -0.6035  0.0131  0.6441  2.2751 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.59867    0.04392   59.17   <2e-16 ***
## SR_12        0.25240    0.01102   22.89   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8457 on 3245 degrees of freedom
## Multiple R-squared:  0.1391, Adjusted R-squared:  0.1388 
## F-statistic: 524.2 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~SR_12, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ SR_12, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.6769 -0.5861  0.0685  0.6321  2.9412 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.84066    0.04848   37.97   <2e-16 ***
## SR_12        0.43635    0.01217   35.85   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9336 on 3245 degrees of freedom
## Multiple R-squared:  0.2837, Adjusted R-squared:  0.2835 
## F-statistic:  1285 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_12~SR_34, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_12 ~ SR_34, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.3864 -0.6126  0.0382  0.6136  3.1612 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.62650    0.04958   32.81   <2e-16 ***
## SR_34        0.42460    0.01236   34.36   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9294 on 3245 degrees of freedom
## Multiple R-squared:  0.2668, Adjusted R-squared:  0.2666 
## F-statistic:  1181 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_34~SR_34, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_34 ~ SR_34, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.2550 -0.4931  0.0069  0.6618  2.3167 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.55232    0.04498   56.74   <2e-16 ***
## SR_34        0.26196    0.01121   23.37   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8432 on 3245 degrees of freedom
## Multiple R-squared:  0.1441, Adjusted R-squared:  0.1438 
## F-statistic: 546.1 on 1 and 3245 DF,  p-value: < 2.2e-16
summary(lm(ISM_56~SR_34, data = InteractAnalyses))
## 
## Call:
## lm(formula = ISM_56 ~ SR_34, data = InteractAnalyses)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.7085 -0.5723  0.1094  0.6549  3.0183 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1.75451    0.04938   35.53   <2e-16 ***
## SR_34        0.45445    0.01230   36.93   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9256 on 3245 degrees of freedom
## Multiple R-squared:  0.296,  Adjusted R-squared:  0.2957 
## F-statistic:  1364 on 1 and 3245 DF,  p-value: < 2.2e-16