Crithidia vs. logical chemicals
Models
## ---------------- AZAGUARD ----------------
mod_azaguard <- glm(
cbind(crith_pos, crith_neg) ~ azaguard_L,
data = gh6,
family = binomial("logit")
)
summary(mod_azaguard)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ azaguard_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.2955 0.1167 -2.531 0.0114 *
## azaguard_LTRUE 0.8890 0.1449 6.137 8.4e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 204.05 on 4 degrees of freedom
## AIC: 237.95
##
## Number of Fisher Scoring iterations: 4
Anova(mod_azaguard)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## azaguard_L 38.371 1 5.849e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ma.s <- emmeans(mod_azaguard, pairwise ~ azaguard_L, type = "response")
ma.s.df <- as.data.frame(ma.s$emmeans)
ma.s
## $emmeans
## azaguard_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- BOTANIGARD 22WP ----------------
mod_botanigard_22wp <- glm(
cbind(crith_pos, crith_neg) ~ botanigard_22wp_L,
data = gh6,
family = binomial("logit")
)
summary(mod_botanigard_22wp)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ botanigard_22wp_L,
## family = binomial("logit"), data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.2955 0.1167 -2.531 0.0114 *
## botanigard_22wp_LTRUE 0.8890 0.1449 6.137 8.4e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 204.05 on 4 degrees of freedom
## AIC: 237.95
##
## Number of Fisher Scoring iterations: 4
Anova(mod_botanigard_22wp)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## botanigard_22wp_L 38.371 1 5.849e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mb22.s <- emmeans(mod_botanigard_22wp, pairwise ~ botanigard_22wp_L, type = "response")
mb22.s.df <- as.data.frame(mb22.s$emmeans)
mb22.s
## $emmeans
## botanigard_22wp_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- BOTANIGARD ES ----------------
mod_botanigard_es <- glm(
cbind(crith_pos, crith_neg) ~ botanigard_es_L,
data = gh6,
family = binomial("logit")
)
summary(mod_botanigard_es)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ botanigard_es_L,
## family = binomial("logit"), data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.50013 0.09244 5.410 6.29e-08 ***
## botanigard_es_LTRUE -0.47481 0.13665 -3.475 0.000511 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 230.29 on 4 degrees of freedom
## AIC: 264.19
##
## Number of Fisher Scoring iterations: 3
Anova(mod_botanigard_es)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## botanigard_es_L 12.134 1 0.000495 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mbe.s <- emmeans(mod_botanigard_es, pairwise ~ botanigard_es_L, type = "response")
mbe.s.df <- as.data.frame(mbe.s$emmeans)
mbe.s
## $emmeans
## botanigard_es_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.622 0.0217 Inf 0.579 0.664
## TRUE 0.506 0.0252 Inf 0.457 0.555
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.61 0.22 Inf 1 3.475 0.0005
##
## Tests are performed on the log odds ratio scale
## ---------------- CAPTIVA PRIME ----------------
mod_captiva_prime <- glm(
cbind(crith_pos, crith_neg) ~ captiva_prime_L,
data = gh6,
family = binomial("logit")
)
summary(mod_captiva_prime)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ captiva_prime_L,
## family = binomial("logit"), data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.2955 0.1167 -2.531 0.0114 *
## captiva_prime_LTRUE 0.8890 0.1449 6.137 8.4e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 204.05 on 4 degrees of freedom
## AIC: 237.95
##
## Number of Fisher Scoring iterations: 4
Anova(mod_captiva_prime)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## captiva_prime_L 38.371 1 5.849e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcp.s <- emmeans(mod_captiva_prime, pairwise ~ captiva_prime_L, type = "response")
mcp.s.df <- as.data.frame(mcp.s$emmeans)
mcp.s
## $emmeans
## captiva_prime_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- NOFLY ----------------
mod_nofly <- glm(
cbind(crith_pos, crith_neg) ~ nofly_L,
data = gh6,
family = binomial("logit")
)
summary(mod_nofly)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ nofly_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.2955 0.1167 -2.531 0.0114 *
## nofly_LTRUE 0.8890 0.1449 6.137 8.4e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 204.05 on 4 degrees of freedom
## AIC: 237.95
##
## Number of Fisher Scoring iterations: 4
Anova(mod_nofly)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## nofly_L 38.371 1 5.849e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mnf.s <- emmeans(mod_nofly, pairwise ~ nofly_L, type = "response")
mnf.s.df <- as.data.frame(mnf.s$emmeans)
mnf.s
## $emmeans
## nofly_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- VENERATE CG ----------------
mod_venerate_cg <- glm(
cbind(crith_pos, crith_neg) ~ venerate_cg_L,
data = gh6,
family = binomial("logit")
)
summary(mod_venerate_cg)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ venerate_cg_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.14425 0.08968 -1.609 0.108
## venerate_cg_LTRUE 1.02694 0.14260 7.202 5.95e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 188.26 on 4 degrees of freedom
## AIC: 222.16
##
## Number of Fisher Scoring iterations: 4
Anova(mod_venerate_cg)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## venerate_cg_L 54.162 1 1.846e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mvc.s <- emmeans(mod_venerate_cg, pairwise ~ venerate_cg_L, type = "response")
mvc.s.df <- as.data.frame(mvc.s$emmeans)
mvc.s
## $emmeans
## venerate_cg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.464 0.0223 Inf 0.421 0.508
## TRUE 0.707 0.0229 Inf 0.660 0.750
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.358 0.0511 Inf 1 -7.202 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- M-PEDE ----------------
mod_m_pede <- glm(
cbind(crith_pos, crith_neg) ~ m_pede_L,
data = gh6,
family = binomial("logit")
)
summary(mod_m_pede)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ m_pede_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.2955 0.1167 -2.531 0.0114 *
## m_pede_LTRUE 0.8890 0.1449 6.137 8.4e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 204.05 on 4 degrees of freedom
## AIC: 237.95
##
## Number of Fisher Scoring iterations: 4
Anova(mod_m_pede)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## m_pede_L 38.371 1 5.849e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mmp.s <- emmeans(mod_m_pede, pairwise ~ m_pede_L, type = "response")
mmp.s.df <- as.data.frame(mmp.s$emmeans)
mmp.s
## $emmeans
## m_pede_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- ROOTSHIELD PLUS ----------------
mod_rootshield_plus <- glm(
cbind(crith_pos, crith_neg) ~ rootshield_plus_L,
data = gh6,
family = binomial("logit")
)
summary(mod_rootshield_plus)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ rootshield_plus_L,
## family = binomial("logit"), data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.25497 0.07654 3.331 0.000864 ***
## rootshield_plus_LTRUE 0.14212 0.16359 0.869 0.384969
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 241.66 on 4 degrees of freedom
## AIC: 275.56
##
## Number of Fisher Scoring iterations: 4
Anova(mod_rootshield_plus)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## rootshield_plus_L 0.75854 1 0.3838
mrs.s <- emmeans(mod_rootshield_plus, pairwise ~ rootshield_plus_L, type = "response")
mrs.s.df <- as.data.frame(mrs.s$emmeans)
mrs.s
## $emmeans
## rootshield_plus_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.563 0.0188 Inf 0.526 0.600
## TRUE 0.598 0.0348 Inf 0.528 0.664
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.868 0.142 Inf 1 -0.869 0.3850
##
## Tests are performed on the log odds ratio scale
## ---------------- LALSTOP K61 ----------------
mod_lalstop_k61 <- glm(
cbind(crith_pos, crith_neg) ~ lalstop_k61_L,
data = gh6,
family = binomial("logit")
)
summary(mod_lalstop_k61)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ lalstop_k61_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.63304 0.07987 7.926 2.26e-15 ***
## lalstop_k61_LTRUE -1.53544 0.17497 -8.776 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 157.44 on 4 degrees of freedom
## AIC: 191.34
##
## Number of Fisher Scoring iterations: 4
Anova(mod_lalstop_k61)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## lalstop_k61_L 84.983 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mlk.s <- emmeans(mod_lalstop_k61, pairwise ~ lalstop_k61_L, type = "response")
mlk.s.df <- as.data.frame(mlk.s$emmeans)
mlk.s
## $emmeans
## lalstop_k61_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.653 0.0181 Inf 0.617 0.688
## TRUE 0.289 0.0320 Inf 0.230 0.355
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 4.64 0.812 Inf 1 8.776 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- BELEAF 50SG ----------------
mod_beleaf_50sg <- glm(
cbind(crith_pos, crith_neg) ~ beleaf_50sg_L,
data = gh6,
family = binomial("logit")
)
summary(mod_beleaf_50sg)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ beleaf_50sg_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.11235 0.07598 -1.479 0.139
## beleaf_50sg_LTRUE 2.54377 0.27160 9.366 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.420 on 5 degrees of freedom
## Residual deviance: 95.039 on 4 degrees of freedom
## AIC: 128.94
##
## Number of Fisher Scoring iterations: 4
Anova(mod_beleaf_50sg)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## beleaf_50sg_L 147.38 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mbf.s <- emmeans(mod_beleaf_50sg, pairwise ~ beleaf_50sg_L, type = "response")
mbf.s.df <- as.data.frame(mbf.s$emmeans)
mbf.s
## $emmeans
## beleaf_50sg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.472 0.0189 Inf 0.435 0.509
## TRUE 0.919 0.0194 Inf 0.872 0.950
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.0786 0.0213 Inf 1 -9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- CORAGEN ----------------
mod_coragen <- glm(
cbind(crith_pos, crith_neg) ~ coragen_L,
data = gh6,
family = binomial("logit")
)
summary(mod_coragen)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ coragen_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.08004 0.14153 0.566 0.5717
## coragen_LTRUE 0.26683 0.16118 1.655 0.0978 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 239.69 on 4 degrees of freedom
## AIC: 273.59
##
## Number of Fisher Scoring iterations: 4
Anova(mod_coragen)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## coragen_L 2.7331 1 0.09829 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcg.s <- emmeans(mod_coragen, pairwise ~ coragen_L, type = "response")
mcg.s.df <- as.data.frame(mcg.s$emmeans)
mcg.s
## $emmeans
## coragen_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.520 0.0353 Inf 0.451 0.588
## TRUE 0.586 0.0187 Inf 0.549 0.622
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.766 0.123 Inf 1 -1.655 0.0978
##
## Tests are performed on the log odds ratio scale
## ---------------- ENTRUST SC ----------------
mod_entrust_sc <- glm(
cbind(crith_pos, crith_neg) ~ entrust_sc_L,
data = gh6,
family = binomial("logit")
)
summary(mod_entrust_sc)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ entrust_sc_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.11235 0.07598 -1.479 0.139
## entrust_sc_LTRUE 2.54377 0.27160 9.366 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.420 on 5 degrees of freedom
## Residual deviance: 95.039 on 4 degrees of freedom
## AIC: 128.94
##
## Number of Fisher Scoring iterations: 4
Anova(mod_entrust_sc)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## entrust_sc_L 147.38 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mes.s <- emmeans(mod_entrust_sc, pairwise ~ entrust_sc_L, type = "response")
mes.s.df <- as.data.frame(mes.s$emmeans)
mes.s
## $emmeans
## entrust_sc_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.472 0.0189 Inf 0.435 0.509
## TRUE 0.919 0.0194 Inf 0.872 0.950
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.0786 0.0213 Inf 1 -9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- PYLON ----------------
mod_pylon <- glm(
cbind(crith_pos, crith_neg) ~ pylon_L,
data = gh6,
family = binomial("logit")
)
summary(mod_pylon)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ pylon_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.59356 0.08577 6.920 4.51e-12 ***
## pylon_LTRUE -0.88903 0.14486 -6.137 8.40e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 204.05 on 4 degrees of freedom
## AIC: 237.95
##
## Number of Fisher Scoring iterations: 4
Anova(mod_pylon)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## pylon_L 38.371 1 5.849e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mpy.s <- emmeans(mod_pylon, pairwise ~ pylon_L, type = "response")
mpy.s.df <- as.data.frame(mpy.s$emmeans)
mpy.s
## $emmeans
## pylon_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.644 0.0197 Inf 0.605 0.682
## TRUE 0.427 0.0286 Inf 0.372 0.483
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 2.43 0.352 Inf 1 6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- GROTTO ----------------
mod_grotto <- glm(
cbind(crith_pos, crith_neg) ~ grotto_L,
data = gh6,
family = binomial("logit")
)
summary(mod_grotto)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ grotto_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.11235 0.07598 -1.479 0.139
## grotto_LTRUE 2.54377 0.27160 9.366 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.420 on 5 degrees of freedom
## Residual deviance: 95.039 on 4 degrees of freedom
## AIC: 128.94
##
## Number of Fisher Scoring iterations: 4
Anova(mod_grotto)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## grotto_L 147.38 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mgr.s <- emmeans(mod_grotto, pairwise ~ grotto_L, type = "response")
mgr.s.df <- as.data.frame(mgr.s$emmeans)
mgr.s
## $emmeans
## grotto_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.472 0.0189 Inf 0.435 0.509
## TRUE 0.919 0.0194 Inf 0.872 0.950
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.0786 0.0213 Inf 1 -9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- LUNA TRANQUILITY ----------------
mod_luna_tranquility <- glm(
cbind(crith_pos, crith_neg) ~ luna_tranquility_L,
data = gh6,
family = binomial("logit")
)
summary(mod_luna_tranquility)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ luna_tranquility_L,
## family = binomial("logit"), data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.26415 0.08278 3.191 0.00142 **
## luna_tranquility_LTRUE 0.06659 0.14352 0.464 0.64267
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 242.20 on 4 degrees of freedom
## AIC: 276.11
##
## Number of Fisher Scoring iterations: 4
Anova(mod_luna_tranquility)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## luna_tranquility_L 0.21553 1 0.6425
mlt.s <- emmeans(mod_luna_tranquility, pairwise ~ luna_tranquility_L, type = "response")
mlt.s.df <- as.data.frame(mlt.s$emmeans)
mlt.s
## $emmeans
## luna_tranquility_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.566 0.0203 Inf 0.525 0.605
## TRUE 0.582 0.0285 Inf 0.525 0.637
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.936 0.134 Inf 1 -0.464 0.6427
##
## Tests are performed on the log odds ratio scale
## ---------------- PREVICUR FLEX ----------------
mod_previcur_flex <- glm(
cbind(crith_pos, crith_neg) ~ previcur_flex_L,
data = gh6,
family = binomial("logit")
)
summary(mod_previcur_flex)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ previcur_flex_L,
## family = binomial("logit"), data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 2.4314 0.2608 9.324 <2e-16 ***
## previcur_flex_LTRUE -2.5438 0.2716 -9.366 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.420 on 5 degrees of freedom
## Residual deviance: 95.039 on 4 degrees of freedom
## AIC: 128.94
##
## Number of Fisher Scoring iterations: 4
Anova(mod_previcur_flex)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## previcur_flex_L 147.38 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mpf.s <- emmeans(mod_previcur_flex, pairwise ~ previcur_flex_L, type = "response")
mpf.s.df <- as.data.frame(mpf.s$emmeans)
mpf.s
## $emmeans
## previcur_flex_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.919 0.0194 Inf 0.872 0.950
## TRUE 0.472 0.0189 Inf 0.435 0.509
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 12.7 3.46 Inf 1 9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
## ---------------- FONTELIS ----------------
mod_fontelis <- glm(
cbind(crith_pos, crith_neg) ~ fontelis_L,
data = gh6,
family = binomial("logit")
)
summary(mod_fontelis)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ fontelis_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.34687 0.07712 4.498 6.86e-06 ***
## fontelis_LTRUE -0.26683 0.16118 -1.655 0.0978 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 239.69 on 4 degrees of freedom
## AIC: 273.59
##
## Number of Fisher Scoring iterations: 4
Anova(mod_fontelis)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## fontelis_L 2.7331 1 0.09829 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mfo.s <- emmeans(mod_fontelis, pairwise ~ fontelis_L, type = "response")
mfo.s.df <- as.data.frame(mfo.s$emmeans)
mfo.s
## $emmeans
## fontelis_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.586 0.0187 Inf 0.549 0.622
## TRUE 0.520 0.0353 Inf 0.451 0.588
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.31 0.21 Inf 1 1.655 0.0978
##
## Tests are performed on the log odds ratio scale
## ---------------- QUADRISTOP ----------------
mod_quadristop <- glm(
cbind(crith_pos, crith_neg) ~ quadristop_L,
data = gh6,
family = binomial("logit")
)
summary(mod_quadristop)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ quadristop_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.2175 0.0714 3.046 0.00232 **
## quadristop_LTRUE 0.6637 0.2321 2.860 0.00424 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 233.71 on 4 degrees of freedom
## AIC: 267.61
##
## Number of Fisher Scoring iterations: 4
Anova(mod_quadristop)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## quadristop_L 8.7121 1 0.003161 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mqs.s <- emmeans(mod_quadristop, pairwise ~ quadristop_L, type = "response")
mqs.s.df <- as.data.frame(mqs.s$emmeans)
mqs.s
## $emmeans
## quadristop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.554 0.0176 Inf 0.519 0.588
## TRUE 0.707 0.0457 Inf 0.610 0.788
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.515 0.12 Inf 1 -2.860 0.0042
##
## Tests are performed on the log odds ratio scale
## ---------------- MILSTOP ----------------
mod_milstop <- glm(
cbind(crith_pos, crith_neg) ~ milstop_L,
data = gh6,
family = binomial("logit")
)
summary(mod_milstop)
##
## Call:
## glm(formula = cbind(crith_pos, crith_neg) ~ milstop_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.2175 0.0714 3.046 0.00232 **
## milstop_LTRUE 0.6637 0.2321 2.860 0.00424 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 242.42 on 5 degrees of freedom
## Residual deviance: 233.71 on 4 degrees of freedom
## AIC: 267.61
##
## Number of Fisher Scoring iterations: 4
Anova(mod_milstop)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(crith_pos, crith_neg)
## LR Chisq Df Pr(>Chisq)
## milstop_L 8.7121 1 0.003161 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mms.s <- emmeans(mod_milstop, pairwise ~ milstop_L, type = "response")
mms.s.df <- as.data.frame(mms.s$emmeans)
mms.s
## $emmeans
## milstop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.554 0.0176 Inf 0.519 0.588
## TRUE 0.707 0.0457 Inf 0.610 0.788
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.515 0.12 Inf 1 -2.860 0.0042
##
## Tests are performed on the log odds ratio scale
ma.s
## $emmeans
## azaguard_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
mb22.s
## $emmeans
## botanigard_22wp_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
mbe.s
## $emmeans
## botanigard_es_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.622 0.0217 Inf 0.579 0.664
## TRUE 0.506 0.0252 Inf 0.457 0.555
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.61 0.22 Inf 1 3.475 0.0005
##
## Tests are performed on the log odds ratio scale
mcp.s
## $emmeans
## captiva_prime_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
mnf.s
## $emmeans
## nofly_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
mvc.s
## $emmeans
## venerate_cg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.464 0.0223 Inf 0.421 0.508
## TRUE 0.707 0.0229 Inf 0.660 0.750
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.358 0.0511 Inf 1 -7.202 <0.0001
##
## Tests are performed on the log odds ratio scale
mmp.s
## $emmeans
## m_pede_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.427 0.0286 Inf 0.372 0.483
## TRUE 0.644 0.0197 Inf 0.605 0.682
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.411 0.0595 Inf 1 -6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
mrs.s
## $emmeans
## rootshield_plus_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.563 0.0188 Inf 0.526 0.600
## TRUE 0.598 0.0348 Inf 0.528 0.664
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.868 0.142 Inf 1 -0.869 0.3850
##
## Tests are performed on the log odds ratio scale
mlk.s
## $emmeans
## lalstop_k61_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.653 0.0181 Inf 0.617 0.688
## TRUE 0.289 0.0320 Inf 0.230 0.355
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 4.64 0.812 Inf 1 8.776 <0.0001
##
## Tests are performed on the log odds ratio scale
mbf.s
## $emmeans
## beleaf_50sg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.472 0.0189 Inf 0.435 0.509
## TRUE 0.919 0.0194 Inf 0.872 0.950
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.0786 0.0213 Inf 1 -9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
mcg.s
## $emmeans
## coragen_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.520 0.0353 Inf 0.451 0.588
## TRUE 0.586 0.0187 Inf 0.549 0.622
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.766 0.123 Inf 1 -1.655 0.0978
##
## Tests are performed on the log odds ratio scale
mes.s
## $emmeans
## entrust_sc_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.472 0.0189 Inf 0.435 0.509
## TRUE 0.919 0.0194 Inf 0.872 0.950
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.0786 0.0213 Inf 1 -9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
mpy.s
## $emmeans
## pylon_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.644 0.0197 Inf 0.605 0.682
## TRUE 0.427 0.0286 Inf 0.372 0.483
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 2.43 0.352 Inf 1 6.137 <0.0001
##
## Tests are performed on the log odds ratio scale
mgr.s
## $emmeans
## grotto_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.472 0.0189 Inf 0.435 0.509
## TRUE 0.919 0.0194 Inf 0.872 0.950
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.0786 0.0213 Inf 1 -9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
mlt.s
## $emmeans
## luna_tranquility_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.566 0.0203 Inf 0.525 0.605
## TRUE 0.582 0.0285 Inf 0.525 0.637
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.936 0.134 Inf 1 -0.464 0.6427
##
## Tests are performed on the log odds ratio scale
mpf.s
## $emmeans
## previcur_flex_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.919 0.0194 Inf 0.872 0.950
## TRUE 0.472 0.0189 Inf 0.435 0.509
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 12.7 3.46 Inf 1 9.366 <0.0001
##
## Tests are performed on the log odds ratio scale
mfo.s
## $emmeans
## fontelis_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.586 0.0187 Inf 0.549 0.622
## TRUE 0.520 0.0353 Inf 0.451 0.588
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.31 0.21 Inf 1 1.655 0.0978
##
## Tests are performed on the log odds ratio scale
mqs.s
## $emmeans
## quadristop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.554 0.0176 Inf 0.519 0.588
## TRUE 0.707 0.0457 Inf 0.610 0.788
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.515 0.12 Inf 1 -2.860 0.0042
##
## Tests are performed on the log odds ratio scale
mms.s
## $emmeans
## milstop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.554 0.0176 Inf 0.519 0.588
## TRUE 0.707 0.0457 Inf 0.610 0.788
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.515 0.12 Inf 1 -2.860 0.0042
##
## Tests are performed on the log odds ratio scale
Summary Table
mod.crith <- list(
azaguard = mod_azaguard,
botanigard_22wp = mod_botanigard_22wp,
botanigard_es = mod_botanigard_es,
captiva_prime = mod_captiva_prime,
nofly = mod_nofly,
venerate_cg = mod_venerate_cg,
m_pede = mod_m_pede,
rootshield_plus = mod_rootshield_plus,
lalstop_k61 = mod_lalstop_k61,
beleaf_50sg = mod_beleaf_50sg,
coragen = mod_coragen,
entrust_sc = mod_entrust_sc,
pylon = mod_pylon,
grotto = mod_grotto,
luna_tranquility = mod_luna_tranquility,
previcur_flex = mod_previcur_flex,
fontelis = mod_fontelis,
quadristop = mod_quadristop,
milstop = mod_milstop
)
summary_table.crith <- do.call(rbind, lapply(names(mod.crith), function(name){
m <- mod.crith[[name]]
s <- summary(m)
coef_row <- s$coefficients[2, ]
lr <- Anova(m)$`LR Chisq`[1]
LRT <- Anova(m)$'Pr(>Chisq)'[1]
data.frame(
Pesticide = name,
Estimate = coef_row["Estimate"],
SE = coef_row["Std. Error"],
z = coef_row["z value"],
p = coef_row["Pr(>|z|)"],
LR_ChiSq = lr,
LRT = LRT,
AIC = AIC(m),
Residual_Deviance = m$deviance
)
}))
summary_table.crith$Estimate <- signif(summary_table.crith$Estimate,3)
summary_table.crith$SE <- signif(summary_table.crith$SE,3)
summary_table.crith$z <- round(summary_table.crith$z,2)
summary_table.crith$p <- signif(summary_table.crith$p,3)
summary_table.crith$LR_ChiSq <- round(summary_table.crith$LR_ChiSq,2)
summary_table.crith$AIC <- round(summary_table.crith$AIC,2)
summary_table.crith$Residual_Deviance <- round(summary_table.crith$Residual_Deviance,2)
summary_table.crith$LRT <- signif(summary_table.crith$LRT,2)
summary_table.crith
## Pesticide Estimate SE z p LR_ChiSq LRT
## Estimate azaguard 0.8890 0.145 6.14 8.40e-10 38.37 5.8e-10
## Estimate1 botanigard_22wp 0.8890 0.145 6.14 8.40e-10 38.37 5.8e-10
## Estimate2 botanigard_es -0.4750 0.137 -3.47 5.11e-04 12.13 5.0e-04
## Estimate3 captiva_prime 0.8890 0.145 6.14 8.40e-10 38.37 5.8e-10
## Estimate4 nofly 0.8890 0.145 6.14 8.40e-10 38.37 5.8e-10
## Estimate5 venerate_cg 1.0300 0.143 7.20 5.95e-13 54.16 1.8e-13
## Estimate6 m_pede 0.8890 0.145 6.14 8.40e-10 38.37 5.8e-10
## Estimate7 rootshield_plus 0.1420 0.164 0.87 3.85e-01 0.76 3.8e-01
## Estimate8 lalstop_k61 -1.5400 0.175 -8.78 1.70e-18 84.98 3.0e-20
## Estimate9 beleaf_50sg 2.5400 0.272 9.37 7.55e-21 147.38 6.5e-34
## Estimate10 coragen 0.2670 0.161 1.66 9.78e-02 2.73 9.8e-02
## Estimate11 entrust_sc 2.5400 0.272 9.37 7.55e-21 147.38 6.5e-34
## Estimate12 pylon -0.8890 0.145 -6.14 8.40e-10 38.37 5.8e-10
## Estimate13 grotto 2.5400 0.272 9.37 7.55e-21 147.38 6.5e-34
## Estimate14 luna_tranquility 0.0666 0.144 0.46 6.43e-01 0.22 6.4e-01
## Estimate15 previcur_flex -2.5400 0.272 -9.37 7.55e-21 147.38 6.5e-34
## Estimate16 fontelis -0.2670 0.161 -1.66 9.78e-02 2.73 9.8e-02
## Estimate17 quadristop 0.6640 0.232 2.86 4.24e-03 8.71 3.2e-03
## Estimate18 milstop 0.6640 0.232 2.86 4.24e-03 8.71 3.2e-03
## AIC Residual_Deviance
## Estimate 237.95 204.05
## Estimate1 237.95 204.05
## Estimate2 264.19 230.29
## Estimate3 237.95 204.05
## Estimate4 237.95 204.05
## Estimate5 222.16 188.26
## Estimate6 237.95 204.05
## Estimate7 275.56 241.66
## Estimate8 191.34 157.44
## Estimate9 128.94 95.04
## Estimate10 273.59 239.69
## Estimate11 128.94 95.04
## Estimate12 237.95 204.05
## Estimate13 128.94 95.04
## Estimate14 276.11 242.20
## Estimate15 128.94 95.04
## Estimate16 273.59 239.69
## Estimate17 267.61 233.71
## Estimate18 267.61 233.71
csdf <- as.data.frame(summary_table.crith)
Plots
basic plot
predictors <- c(
"azaguard_L",
"botanigard_22wp_L",
"botanigard_es_L",
"captiva_prime_L",
"nofly_L",
"venerate_cg_L",
"m_pede_L",
"rootshield_plus_L",
"lalstop_k61_L",
"beleaf_50sg_L",
"coragen_L",
"entrust_sc_L",
"pylon_L",
"grotto_L",
"luna_tranquility_L",
"previcur_flex_L",
"fontelis_L",
"quadristop_L",
"milstop_L"
)
gh6$crith_prop <- gh6$crith_pos /
(gh6$crith_pos + gh6$crith_neg)
gh6$crith_se <- sqrt(
(gh6$crith_prop * (1 - gh6$crith_prop)) /
(gh6$crith_pos + gh6$crith_neg)
)
plots <- lapply(predictors, function(var){
ggplot(gh6, aes_string(x = var, y = "crith_prop")) +
geom_boxplot(width = 0.5, alpha = 0.6, outlier.shape = NA) +
geom_jitter(width = 0.1, size = 3, alpha = 0.8) +
labs(
x = var,
y = "Crithidia infection proportion"
) +
theme_classic(base_size = 14) +
scale_x_discrete(labels = c("FALSE" = "Not applied", "TRUE" = "Applied"))
})
## 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 per session.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
wrap_plots(plots, ncol = 4)

updated crithidia vs. logical chems
plots <- lapply(1:nrow(chem_info), function(i){
chem_type <- chem_info$type[i]
chem_mode <- chem_info$mode[i]
var <- chem_info$pesticide[i]
chem_class <- chem_info$class[i]
var2 <- gsub("_L$", "", chem_info$pesticide[i]) # remove _L at end
stats_row <- summary_table.crith[summary_table.crith$Pesticide == var2, ]
annot_text <- paste0(
"LRT χ² = ", round(stats_row$LR_ChiSq, 2),
", p = ", signif(stats_row$LRT, 3))
ggplot(gh6, aes_string(x = var, y = "crith_prop")) +
geom_boxplot(fill = chem_colors[chem_class],
alpha = 0.6,
width = 0.5,
outlier.shape = NA) +
geom_jitter(width = 0.1, size = 3, alpha = 0.8) +
labs(
title = gsub("_L","",var),
subtitle = paste(chem_type, "-", chem_mode),
x = "",
y = "Probability of *Crithidia* Detection"
) +
scale_x_discrete(labels = c("FALSE"="Not applied","TRUE"="Applied")) +
theme_classic(base_size = 16) +
theme(axis.title.y = ggtext::element_markdown()) +
annotate(
"text",
x = 2, y = 1.05, # top-right
label = annot_text,
hjust = 1, vjust = 1, size = 4
) +
coord_cartesian(ylim = c(0, 1))
})
wrap_plots(plots, ncol = 4)

Apicystis vs. logical chemicals
mod.a_azaguard <- glm(
cbind(api_pos, api_neg) ~ azaguard_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_azaguard)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ azaguard_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.7346 0.1617 -10.728 <2e-16 ***
## azaguard_LTRUE -0.1236 0.2015 -0.614 0.539
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 213.41 on 4 degrees of freedom
## AIC: 235.79
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_azaguard)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## azaguard_L 0.37344 1 0.5411
ma.a <- emmeans(mod.a_azaguard, pairwise ~ azaguard_L, type = "response")
ma.a.df <- as.data.frame(ma.a$emmeans)
ma.a
## $emmeans
## azaguard_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mod.a_botanigard_22wp <- glm(
cbind(api_pos, api_neg) ~ botanigard_22wp_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_botanigard_22wp)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ botanigard_22wp_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.7346 0.1617 -10.728 <2e-16 ***
## botanigard_22wp_LTRUE -0.1236 0.2015 -0.614 0.539
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 213.41 on 4 degrees of freedom
## AIC: 235.79
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_botanigard_22wp)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## botanigard_22wp_L 0.37344 1 0.5411
mb22.a <- emmeans(mod.a_botanigard_22wp, pairwise ~ botanigard_22wp_L, type = "response")
mb22.a.df <- as.data.frame(mb22.a$emmeans)
mb22.a
## $emmeans
## botanigard_22wp_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mod.a_botanigard_es <- glm(
cbind(api_pos, api_neg) ~ botanigard_es_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_botanigard_es)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ botanigard_es_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -2.3092 0.1563 -14.774 < 2e-16 ***
## botanigard_es_LTRUE 0.9387 0.2003 4.687 2.77e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 190.83 on 4 degrees of freedom
## AIC: 213.21
##
## Number of Fisher Scoring iterations: 5
Anova(mod.a_botanigard_es)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## botanigard_es_L 22.958 1 1.656e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mbe.a <- emmeans(mod.a_botanigard_es, pairwise ~ botanigard_es_L, type = "response")
mbe.a.df <- as.data.frame(mbe.a$emmeans)
mbe.a
## $emmeans
## botanigard_es_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.0904 0.0128 Inf 0.0681 0.119
## TRUE 0.2025 0.0202 Inf 0.1658 0.245
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.391 0.0783 Inf 1 -4.687 <0.0001
##
## Tests are performed on the log odds ratio scale
mod.a_captiva_prime <- glm(
cbind(api_pos, api_neg) ~ captiva_prime_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_captiva_prime)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ captiva_prime_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.7346 0.1617 -10.728 <2e-16 ***
## captiva_prime_LTRUE -0.1236 0.2015 -0.614 0.539
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 213.41 on 4 degrees of freedom
## AIC: 235.79
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_captiva_prime)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## captiva_prime_L 0.37344 1 0.5411
mcp.a <- emmeans(mod.a_captiva_prime, pairwise ~ captiva_prime_L, type = "response")
mcp.a.df <- as.data.frame(mcp.a$emmeans)
mcp.a
## $emmeans
## captiva_prime_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mod.a_nofly <- glm(
cbind(api_pos, api_neg) ~ nofly_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_nofly)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ nofly_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.7346 0.1617 -10.728 <2e-16 ***
## nofly_LTRUE -0.1236 0.2015 -0.614 0.539
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 213.41 on 4 degrees of freedom
## AIC: 235.79
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_nofly)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## nofly_L 0.37344 1 0.5411
mnf.a <- emmeans(mod.a_nofly, pairwise ~ nofly_L, type = "response")
mnf.a.df <- as.data.frame(mnf.a$emmeans)
mnf.a
## $emmeans
## nofly_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mod.a_venerate_cg <- glm(
cbind(api_pos, api_neg) ~ venerate_cg_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_venerate_cg)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ venerate_cg_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.0986 0.1033 -10.637 <2e-16 ***
## venerate_cg_LTRUE -20.8817 1813.9497 -0.012 0.991
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.783 on 5 degrees of freedom
## Residual deviance: 52.926 on 4 degrees of freedom
## AIC: 75.307
##
## Number of Fisher Scoring iterations: 16
Anova(mod.a_venerate_cg)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## venerate_cg_L 160.86 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mvc.a <- emmeans(mod.a_venerate_cg, pairwise ~ venerate_cg_L, type = "response")
mvc.a.df <- as.data.frame(mvc.a$emmeans)
mvc.a
## $emmeans
## venerate_cg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.25 1.94e-02 Inf 0.214 0.29
## TRUE 0.00 5.16e-07 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.17e+09 2.13e+12 Inf 1 0.012 0.9908
##
## Tests are performed on the log odds ratio scale
mod.a_m_pede <- glm(
cbind(api_pos, api_neg) ~ m_pede_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_m_pede)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ m_pede_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.7346 0.1617 -10.728 <2e-16 ***
## m_pede_LTRUE -0.1236 0.2015 -0.614 0.539
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 213.41 on 4 degrees of freedom
## AIC: 235.79
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_m_pede)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## m_pede_L 0.37344 1 0.5411
mmp.a <- emmeans(mod.a_m_pede, pairwise ~ m_pede_L, type = "response")
mmp.a.df <- as.data.frame(mmp.a$emmeans)
mmp.a
## $emmeans
## m_pede_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mod.a_rootshield_plus <- glm(
cbind(api_pos, api_neg) ~ rootshield_plus_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_rootshield_plus)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ rootshield_plus_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.7817 0.1081 -16.484 <2e-16 ***
## rootshield_plus_LTRUE -0.1585 0.2396 -0.661 0.508
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 213.34 on 4 degrees of freedom
## AIC: 235.72
##
## Number of Fisher Scoring iterations: 5
Anova(mod.a_rootshield_plus)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## rootshield_plus_L 0.44747 1 0.5035
mrs.a <- emmeans(mod.a_rootshield_plus, pairwise ~ rootshield_plus_L, type = "response")
mrs.a.df <- as.data.frame(mrs.a$emmeans)
mrs.a
## $emmeans
## rootshield_plus_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.144 0.0133 Inf 0.1199 0.172
## TRUE 0.126 0.0235 Inf 0.0863 0.179
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.17 0.281 Inf 1 0.661 0.5084
##
## Tests are performed on the log odds ratio scale
mod.a_lalstop_k61 <- glm(
cbind(api_pos, api_neg) ~ lalstop_k61_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_lalstop_k61)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ lalstop_k61_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.9659 0.1158 -16.975 < 2e-16 ***
## lalstop_k61_LTRUE 0.5734 0.2112 2.714 0.00664 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 206.78 on 4 degrees of freedom
## AIC: 229.16
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_lalstop_k61)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## lalstop_k61_L 7.0039 1 0.008133 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mlk.a <- emmeans(mod.a_lalstop_k61, pairwise ~ lalstop_k61_L, type = "response")
mlk.a.df <- as.data.frame(mlk.a$emmeans)
mlk.a
## $emmeans
## lalstop_k61_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.123 0.0125 Inf 0.100 0.149
## TRUE 0.199 0.0282 Inf 0.149 0.260
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.564 0.119 Inf 1 -2.714 0.0066
##
## Tests are performed on the log odds ratio scale
mod.a_beleaf_50sg <- glm(
cbind(api_pos, api_neg) ~ beleaf_50sg_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_beleaf_50sg)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ beleaf_50sg_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.51732 0.09876 -15.36 <2e-16 ***
## beleaf_50sg_LTRUE -18.47059 943.71566 -0.02 0.984
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 145.52 on 4 degrees of freedom
## AIC: 167.91
##
## Number of Fisher Scoring iterations: 14
Anova(mod.a_beleaf_50sg)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## beleaf_50sg_L 68.259 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mbf.a <- emmeans(mod.a_beleaf_50sg, pairwise ~ beleaf_50sg_L, type = "response")
mbf.a.df <- as.data.frame(mbf.a$emmeans)
mbf.a
## $emmeans
## beleaf_50sg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.18 1.46e-02 Inf 0.153 0.21
## TRUE 0.00 1.97e-06 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.05e+08 9.92e+10 Inf 1 0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mod.a_coragen <- glm(
cbind(api_pos, api_neg) ~ coragen_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_coragen)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ coragen_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.4055 0.1443 -2.809 0.00497 **
## coragen_LTRUE -2.2618 0.2112 -10.710 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.783 on 5 degrees of freedom
## Residual deviance: 92.901 on 4 degrees of freedom
## AIC: 115.28
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_coragen)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## coragen_L 120.88 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcg.a <- emmeans(mod.a_coragen, pairwise ~ coragen_L, type = "response")
mcg.a.df <- as.data.frame(mcg.a$emmeans)
mcg.a
## $emmeans
## coragen_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.4000 0.03460 Inf 0.3344 0.4694
## TRUE 0.0649 0.00936 Inf 0.0488 0.0859
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 9.6 2.03 Inf 1 10.710 <0.0001
##
## Tests are performed on the log odds ratio scale
mod.a_entrust_sc <- glm(
cbind(api_pos, api_neg) ~ entrust_sc_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_entrust_sc)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ entrust_sc_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.51732 0.09876 -15.36 <2e-16 ***
## entrust_sc_LTRUE -18.47059 943.71566 -0.02 0.984
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 145.52 on 4 degrees of freedom
## AIC: 167.91
##
## Number of Fisher Scoring iterations: 14
Anova(mod.a_entrust_sc)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## entrust_sc_L 68.259 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mes.a <- emmeans(mod.a_entrust_sc, pairwise ~ entrust_sc_L, type = "response")
mes.a.df <- as.data.frame(mes.a$emmeans)
mes.a
## $emmeans
## entrust_sc_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.18 1.46e-02 Inf 0.153 0.21
## TRUE 0.00 1.97e-06 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.05e+08 9.92e+10 Inf 1 0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mod.a_pylon <- glm(
cbind(api_pos, api_neg) ~ pylon_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_pylon)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ pylon_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.8582 0.1202 -15.459 <2e-16 ***
## pylon_LTRUE 0.1236 0.2015 0.614 0.539
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 213.41 on 4 degrees of freedom
## AIC: 235.79
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_pylon)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## pylon_L 0.37344 1 0.5411
mpy.a <- emmeans(mod.a_pylon, pairwise ~ pylon_L, type = "response")
mpy.a.df <- as.data.frame(mpy.a$emmeans)
mpy.a
## $emmeans
## pylon_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.135 0.0140 Inf 0.110 0.165
## TRUE 0.150 0.0206 Inf 0.114 0.195
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.884 0.178 Inf 1 -0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mod.a_grotto <- glm(
cbind(api_pos, api_neg) ~ grotto_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_grotto)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ grotto_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.51732 0.09876 -15.36 <2e-16 ***
## grotto_LTRUE -18.47059 943.71566 -0.02 0.984
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 145.52 on 4 degrees of freedom
## AIC: 167.91
##
## Number of Fisher Scoring iterations: 14
Anova(mod.a_grotto)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## grotto_L 68.259 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mgr.a <- emmeans(mod.a_grotto, pairwise ~ grotto_L, type = "response")
mgr.a.df <- as.data.frame(mgr.a$emmeans)
mgr.a
## $emmeans
## grotto_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.18 1.46e-02 Inf 0.153 0.21
## TRUE 0.00 1.97e-06 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.05e+08 9.92e+10 Inf 1 0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mod.a_luna_tranquility <- glm(
cbind(api_pos, api_neg) ~ luna_tranquility_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_luna_tranquility)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ luna_tranquility_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -2.6283 0.1637 -16.053 < 2e-16 ***
## luna_tranquility_LTRUE 1.7050 0.2079 8.199 2.42e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 140.65 on 4 degrees of freedom
## AIC: 163.03
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_luna_tranquility)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## luna_tranquility_L 73.132 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mlt.a <- emmeans(mod.a_luna_tranquility, pairwise ~ luna_tranquility_L, type = "response")
mlt.a.df <- as.data.frame(mlt.a$emmeans)
mlt.a
## $emmeans
## luna_tranquility_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.0673 0.0103 Inf 0.0498 0.0905
## TRUE 0.2843 0.0261 Inf 0.2360 0.3380
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.182 0.0378 Inf 1 -8.199 <0.0001
##
## Tests are performed on the log odds ratio scale
mod.a_previcur_flex <- glm(
cbind(api_pos, api_neg) ~ previcur_flex_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_previcur_flex)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ previcur_flex_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -19.99 943.72 -0.021 0.983
## previcur_flex_LTRUE 18.47 943.72 0.020 0.984
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 145.52 on 4 degrees of freedom
## AIC: 167.91
##
## Number of Fisher Scoring iterations: 14
Anova(mod.a_previcur_flex)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## previcur_flex_L 68.259 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mpf.a <- emmeans(mod.a_previcur_flex, pairwise ~ previcur_flex_L, type = "response")
mpf.a.df <- as.data.frame(mpf.a$emmeans)
mpf.a
## $emmeans
## previcur_flex_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.00 1.97e-06 Inf 0.000 1.00
## TRUE 0.18 1.46e-02 Inf 0.153 0.21
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 9.51e-09 8.98e-06 Inf 1 -0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mod.a_fontelis <- glm(
cbind(api_pos, api_neg) ~ fontelis_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_fontelis)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ fontelis_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -2.6672 0.1542 -17.30 <2e-16 ***
## fontelis_LTRUE 2.2618 0.2112 10.71 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.783 on 5 degrees of freedom
## Residual deviance: 92.901 on 4 degrees of freedom
## AIC: 115.28
##
## Number of Fisher Scoring iterations: 6
Anova(mod.a_fontelis)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## fontelis_L 120.88 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mfo.a <- emmeans(mod.a_fontelis, pairwise ~ fontelis_L, type = "response")
mfo.a.df <- as.data.frame(mfo.a$emmeans)
mfo.a
## $emmeans
## fontelis_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.0649 0.00936 Inf 0.0488 0.0859
## TRUE 0.4000 0.03460 Inf 0.3344 0.4694
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.104 0.022 Inf 1 -10.710 <0.0001
##
## Tests are performed on the log odds ratio scale
mod.a_quadristop <- glm(
cbind(api_pos, api_neg) ~ quadristop_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_quadristop)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ quadristop_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.72574 0.09908 -17.418 <2e-16 ***
## quadristop_LTRUE -1.20812 0.46953 -2.573 0.0101 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 204.57 on 4 degrees of freedom
## AIC: 226.95
##
## Number of Fisher Scoring iterations: 5
Anova(mod.a_quadristop)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## quadristop_L 9.2159 1 0.002399 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mqs.a <- emmeans(mod.a_quadristop, pairwise ~ quadristop_L, type = "response")
mqs.a.df <- as.data.frame(mqs.a$emmeans)
mqs.a
## $emmeans
## quadristop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.1511 0.0127 Inf 0.1279 0.178
## TRUE 0.0505 0.0220 Inf 0.0212 0.116
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 3.35 1.57 Inf 1 2.573 0.0101
##
## Tests are performed on the log odds ratio scale
mod.a_milstop <- glm(
cbind(api_pos, api_neg) ~ milstop_L,
data = gh6,
family = binomial("logit"))
summary(mod.a_milstop)
##
## Call:
## glm(formula = cbind(api_pos, api_neg) ~ milstop_L, family = binomial("logit"),
## data = gh6)
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.72574 0.09908 -17.418 <2e-16 ***
## milstop_LTRUE -1.20812 0.46953 -2.573 0.0101 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 213.78 on 5 degrees of freedom
## Residual deviance: 204.57 on 4 degrees of freedom
## AIC: 226.95
##
## Number of Fisher Scoring iterations: 5
Anova(mod.a_milstop)
## Analysis of Deviance Table (Type II tests)
##
## Response: cbind(api_pos, api_neg)
## LR Chisq Df Pr(>Chisq)
## milstop_L 9.2159 1 0.002399 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mms.a <- emmeans(mod.a_milstop, pairwise ~ milstop_L, type = "response")
mms.a.df <- as.data.frame(mms.a$emmeans)
mms.a
## $emmeans
## milstop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.1511 0.0127 Inf 0.1279 0.178
## TRUE 0.0505 0.0220 Inf 0.0212 0.116
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 3.35 1.57 Inf 1 2.573 0.0101
##
## Tests are performed on the log odds ratio scale
ma.a
## $emmeans
## azaguard_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mb22.a
## $emmeans
## botanigard_22wp_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mbe.a
## $emmeans
## botanigard_es_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.0904 0.0128 Inf 0.0681 0.119
## TRUE 0.2025 0.0202 Inf 0.1658 0.245
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.391 0.0783 Inf 1 -4.687 <0.0001
##
## Tests are performed on the log odds ratio scale
mcp.a
## $emmeans
## captiva_prime_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mnf.a
## $emmeans
## nofly_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mvc.a
## $emmeans
## venerate_cg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.25 1.94e-02 Inf 0.214 0.29
## TRUE 0.00 5.16e-07 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.17e+09 2.13e+12 Inf 1 0.012 0.9908
##
## Tests are performed on the log odds ratio scale
mmp.a
## $emmeans
## m_pede_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.150 0.0206 Inf 0.114 0.195
## TRUE 0.135 0.0140 Inf 0.110 0.165
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.13 0.228 Inf 1 0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mrs.a
## $emmeans
## rootshield_plus_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.144 0.0133 Inf 0.1199 0.172
## TRUE 0.126 0.0235 Inf 0.0863 0.179
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.17 0.281 Inf 1 0.661 0.5084
##
## Tests are performed on the log odds ratio scale
mlk.a
## $emmeans
## lalstop_k61_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.123 0.0125 Inf 0.100 0.149
## TRUE 0.199 0.0282 Inf 0.149 0.260
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.564 0.119 Inf 1 -2.714 0.0066
##
## Tests are performed on the log odds ratio scale
mbf.a
## $emmeans
## beleaf_50sg_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.18 1.46e-02 Inf 0.153 0.21
## TRUE 0.00 1.97e-06 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.05e+08 9.92e+10 Inf 1 0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mcg.a
## $emmeans
## coragen_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.4000 0.03460 Inf 0.3344 0.4694
## TRUE 0.0649 0.00936 Inf 0.0488 0.0859
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 9.6 2.03 Inf 1 10.710 <0.0001
##
## Tests are performed on the log odds ratio scale
mes.a
## $emmeans
## entrust_sc_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.18 1.46e-02 Inf 0.153 0.21
## TRUE 0.00 1.97e-06 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.05e+08 9.92e+10 Inf 1 0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mpy.a
## $emmeans
## pylon_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.135 0.0140 Inf 0.110 0.165
## TRUE 0.150 0.0206 Inf 0.114 0.195
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.884 0.178 Inf 1 -0.614 0.5394
##
## Tests are performed on the log odds ratio scale
mgr.a
## $emmeans
## grotto_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.18 1.46e-02 Inf 0.153 0.21
## TRUE 0.00 1.97e-06 Inf 0.000 1.00
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 1.05e+08 9.92e+10 Inf 1 0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mlt.a
## $emmeans
## luna_tranquility_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.0673 0.0103 Inf 0.0498 0.0905
## TRUE 0.2843 0.0261 Inf 0.2360 0.3380
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.182 0.0378 Inf 1 -8.199 <0.0001
##
## Tests are performed on the log odds ratio scale
mpf.a
## $emmeans
## previcur_flex_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.00 1.97e-06 Inf 0.000 1.00
## TRUE 0.18 1.46e-02 Inf 0.153 0.21
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 9.51e-09 8.98e-06 Inf 1 -0.020 0.9844
##
## Tests are performed on the log odds ratio scale
mfo.a
## $emmeans
## fontelis_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.0649 0.00936 Inf 0.0488 0.0859
## TRUE 0.4000 0.03460 Inf 0.3344 0.4694
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 0.104 0.022 Inf 1 -10.710 <0.0001
##
## Tests are performed on the log odds ratio scale
mqs.a
## $emmeans
## quadristop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.1511 0.0127 Inf 0.1279 0.178
## TRUE 0.0505 0.0220 Inf 0.0212 0.116
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 3.35 1.57 Inf 1 2.573 0.0101
##
## Tests are performed on the log odds ratio scale
mms.a
## $emmeans
## milstop_L prob SE df asymp.LCL asymp.UCL
## FALSE 0.1511 0.0127 Inf 0.1279 0.178
## TRUE 0.0505 0.0220 Inf 0.0212 0.116
##
## Confidence level used: 0.95
## Intervals are back-transformed from the logit scale
##
## $contrasts
## contrast odds.ratio SE df null z.ratio p.value
## FALSE / TRUE 3.35 1.57 Inf 1 2.573 0.0101
##
## Tests are performed on the log odds ratio scale
summary table
mod.a.api <- list(
azaguard = mod.a_azaguard,
botanigard_22wp = mod.a_botanigard_22wp,
botanigard_es = mod.a_botanigard_es,
captiva_prime = mod.a_captiva_prime,
nofly = mod.a_nofly,
venerate_cg = mod.a_venerate_cg,
m_pede = mod.a_m_pede,
rootshield_plus = mod.a_rootshield_plus,
lalstop_k61 = mod.a_lalstop_k61,
beleaf_50sg = mod.a_beleaf_50sg,
coragen = mod.a_coragen,
entrust_sc = mod.a_entrust_sc,
pylon = mod.a_pylon,
grotto = mod.a_grotto,
luna_tranquility = mod.a_luna_tranquility,
previcur_flex = mod.a_previcur_flex,
fontelis = mod.a_fontelis,
quadristop = mod.a_quadristop,
milstop = mod.a_milstop
)
summary_table.api <- do.call(rbind, lapply(names(mod.a.api), function(name){
m <- mod.a.api[[name]]
s <- summary(m)
coef_row <- s$coefficients[2, ]
lr <- Anova(m)$`LR Chisq`[1]
LRT <- Anova(m)$'Pr(>Chisq)'[1]
data.frame(
Pesticide = name,
Estimate = coef_row["Estimate"],
SE = coef_row["Std. Error"],
z = coef_row["z value"],
p = coef_row["Pr(>|z|)"],
LR_ChiSq = lr,
LRT = LRT,
AIC = AIC(m),
Residual_Deviance = m$deviance
)
}))
summary_table.api$Estimate <- signif(summary_table.api$Estimate,3)
summary_table.api$SE <- signif(summary_table.api$SE,3)
summary_table.api$z <- round(summary_table.api$z,2)
summary_table.api$p <- signif(summary_table.api$p,3)
summary_table.api$LR_ChiSq <- round(summary_table.api$LR_ChiSq,2)
summary_table.api$AIC <- round(summary_table.api$AIC,2)
summary_table.api$Residual_Deviance <- round(summary_table.api$Residual_Deviance,2)
summary_table.api$LRT <- signif(summary_table.api$LRT,2)
summary_table.api
## Pesticide Estimate SE z p LR_ChiSq LRT
## Estimate azaguard -0.124 0.201 -0.61 5.39e-01 0.37 5.4e-01
## Estimate1 botanigard_22wp -0.124 0.201 -0.61 5.39e-01 0.37 5.4e-01
## Estimate2 botanigard_es 0.939 0.200 4.69 2.77e-06 22.96 1.7e-06
## Estimate3 captiva_prime -0.124 0.201 -0.61 5.39e-01 0.37 5.4e-01
## Estimate4 nofly -0.124 0.201 -0.61 5.39e-01 0.37 5.4e-01
## Estimate5 venerate_cg -20.900 1810.000 -0.01 9.91e-01 160.86 7.4e-37
## Estimate6 m_pede -0.124 0.201 -0.61 5.39e-01 0.37 5.4e-01
## Estimate7 rootshield_plus -0.158 0.240 -0.66 5.08e-01 0.45 5.0e-01
## Estimate8 lalstop_k61 0.573 0.211 2.71 6.64e-03 7.00 8.1e-03
## Estimate9 beleaf_50sg -18.500 944.000 -0.02 9.84e-01 68.26 1.4e-16
## Estimate10 coragen -2.260 0.211 -10.71 9.15e-27 120.88 4.1e-28
## Estimate11 entrust_sc -18.500 944.000 -0.02 9.84e-01 68.26 1.4e-16
## Estimate12 pylon 0.124 0.201 0.61 5.39e-01 0.37 5.4e-01
## Estimate13 grotto -18.500 944.000 -0.02 9.84e-01 68.26 1.4e-16
## Estimate14 luna_tranquility 1.700 0.208 8.20 2.42e-16 73.13 1.2e-17
## Estimate15 previcur_flex 18.500 944.000 0.02 9.84e-01 68.26 1.4e-16
## Estimate16 fontelis 2.260 0.211 10.71 9.15e-27 120.88 4.1e-28
## Estimate17 quadristop -1.210 0.470 -2.57 1.01e-02 9.22 2.4e-03
## Estimate18 milstop -1.210 0.470 -2.57 1.01e-02 9.22 2.4e-03
## AIC Residual_Deviance
## Estimate 235.79 213.41
## Estimate1 235.79 213.41
## Estimate2 213.21 190.83
## Estimate3 235.79 213.41
## Estimate4 235.79 213.41
## Estimate5 75.31 52.93
## Estimate6 235.79 213.41
## Estimate7 235.72 213.34
## Estimate8 229.16 206.78
## Estimate9 167.91 145.52
## Estimate10 115.28 92.90
## Estimate11 167.91 145.52
## Estimate12 235.79 213.41
## Estimate13 167.91 145.52
## Estimate14 163.03 140.65
## Estimate15 167.91 145.52
## Estimate16 115.28 92.90
## Estimate17 226.95 204.57
## Estimate18 226.95 204.57
asdf <- as.data.frame(summary_table.api)
plots
plots <- lapply(1:nrow(chem_info), function(i){
chem_type <- chem_info$type[i]
chem_mode <- chem_info$mode[i]
chem_class <- chem_info$class[i]
var <- chem_info$pesticide[i]
var2 <- gsub("_L$", "", chem_info$pesticide[i]) # remove _L at end
stats_row <- summary_table.api[summary_table.api$Pesticide == var2, ]
annot_text <- paste0(
"LRT χ² = ", round(stats_row$LR_ChiSq, 2),
", p = ", signif(stats_row$LRT, 3))
ggplot(gh6, aes_string(x = var, y = "api_prop")) +
geom_boxplot(fill = chem_colors[chem_class],
alpha = 0.6,
width = 0.5,
outlier.shape = NA) +
geom_jitter(width = 0.1, size = 3, alpha = 0.8) +
labs(
title = gsub("_L","",var),
subtitle = paste(chem_type, "-", chem_mode),
x = "",
y = "Probability of *Apicystis* detection"
) +
scale_x_discrete(labels = c("FALSE"="Not applied","TRUE"="Applied")) +
theme_classic(base_size = 16) +
theme(axis.title.y = ggtext::element_markdown()) +
annotate(
"text",
x = 2, y = 1.05, # top-right
label = annot_text,
hjust = 1, vjust = 1, size = 4
) +
coord_cartesian(ylim = c(0, 1))
})
wrap_plots(plots, ncol = 4)

anther bruising vs. specific chemicals presence absence
anth_chem_m1 <- glm(logbruise ~ azaguard_L, data = gh6)
Anova(anth_chem_m1)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## azaguard_L 0.65877 1 0.417
summary(anth_chem_m1)
##
## Call:
## glm(formula = logbruise ~ azaguard_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.82271 0.03698 22.247 2.42e-05 ***
## azaguard_LTRUE -0.04245 0.05230 -0.812 0.463
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.00410272)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.016411 on 4 degrees of freedom
## AIC: -12.382
##
## Number of Fisher Scoring iterations: 2
anth_chem_m2 <- glm(logbruise ~ botanigard_22wp_L, data = gh6)
Anova(anth_chem_m2)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## botanigard_22wp_L 0.65877 1 0.417
summary(anth_chem_m2)
##
## Call:
## glm(formula = logbruise ~ botanigard_22wp_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.82271 0.03698 22.247 2.42e-05 ***
## botanigard_22wp_LTRUE -0.04245 0.05230 -0.812 0.463
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.00410272)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.016411 on 4 degrees of freedom
## AIC: -12.382
##
## Number of Fisher Scoring iterations: 2
anth_chem_m3 <- glm(logbruise ~ botanigard_es_L, data = gh6)
Anova(anth_chem_m3)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## botanigard_es_L 0.12131 1 0.7276
summary(anth_chem_m3)
##
## Call:
## glm(formula = logbruise ~ botanigard_es_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.79464 0.03405 23.337 2e-05 ***
## botanigard_es_LTRUE 0.02054 0.05898 0.348 0.745
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.004637755)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.018551 on 4 degrees of freedom
## AIC: -11.647
##
## Number of Fisher Scoring iterations: 2
anth_chem_m4 <- glm(logbruise ~ captiva_prime_L, data = gh6)
Anova(anth_chem_m4)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## captiva_prime_L 0.65877 1 0.417
summary(anth_chem_m4)
##
## Call:
## glm(formula = logbruise ~ captiva_prime_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.82271 0.03698 22.247 2.42e-05 ***
## captiva_prime_LTRUE -0.04245 0.05230 -0.812 0.463
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.00410272)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.016411 on 4 degrees of freedom
## AIC: -12.382
##
## Number of Fisher Scoring iterations: 2
anth_chem_m5 <- glm(logbruise ~ nofly_L, data = gh6)
Anova(anth_chem_m5)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## nofly_L 0.65877 1 0.417
summary(anth_chem_m5)
##
## Call:
## glm(formula = logbruise ~ nofly_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.82271 0.03698 22.247 2.42e-05 ***
## nofly_LTRUE -0.04245 0.05230 -0.812 0.463
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.00410272)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.016411 on 4 degrees of freedom
## AIC: -12.382
##
## Number of Fisher Scoring iterations: 2
anth_chem_m6 <- glm(logbruise ~ venerate_cg_L, data = gh6)
Anova(anth_chem_m6)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## venerate_cg_L 6.6774 1 0.009765 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(anth_chem_m6)
##
## Call:
## glm(formula = logbruise ~ venerate_cg_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.83305 0.02115 39.379 2.48e-06 ***
## venerate_cg_LTRUE -0.09468 0.03664 -2.584 0.0611 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.001790108)
##
## Null deviance: 0.0191136 on 5 degrees of freedom
## Residual deviance: 0.0071604 on 4 degrees of freedom
## AIC: -17.358
##
## Number of Fisher Scoring iterations: 2
anth_chem_m7 <- glm(logbruise ~ m_pede_L, data = gh6)
Anova(anth_chem_m7)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## m_pede_L 0.65877 1 0.417
summary(anth_chem_m7)
##
## Call:
## glm(formula = logbruise ~ m_pede_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.82271 0.03698 22.247 2.42e-05 ***
## m_pede_LTRUE -0.04245 0.05230 -0.812 0.463
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.00410272)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.016411 on 4 degrees of freedom
## AIC: -12.382
##
## Number of Fisher Scoring iterations: 2
anth_chem_m8 <- glm(logbruise ~ rootshield_plus_L, data = gh6)
Anova(anth_chem_m8)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## rootshield_plus_L 0.0039631 1 0.9498
summary(anth_chem_m8)
##
## Call:
## glm(formula = logbruise ~ rootshield_plus_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.800232 0.034546 23.164 2.06e-05 ***
## rootshield_plus_LTRUE 0.003767 0.059835 0.063 0.953
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.004773677)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.019095 on 4 degrees of freedom
## AIC: -11.473
##
## Number of Fisher Scoring iterations: 2
anth_chem_m9 <- glm(logbruise ~ lalstop_k61_L, data = gh6)
Anova(anth_chem_m9)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## lalstop_k61_L 2.1541 1 0.1422
summary(anth_chem_m9)
##
## Call:
## glm(formula = logbruise ~ lalstop_k61_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.77788 0.02787 27.916 9.8e-06 ***
## lalstop_k61_LTRUE 0.07084 0.04826 1.468 0.216
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.003105838)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.012423 on 4 degrees of freedom
## AIC: -14.052
##
## Number of Fisher Scoring iterations: 2
anth_chem_m10 <- glm(logbruise ~ beleaf_50sg_L, data = gh6)
Anova(anth_chem_m10)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## beleaf_50sg_L 4.3437 1 0.03715 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(anth_chem_m10)
##
## Call:
## glm(formula = logbruise ~ beleaf_50sg_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.81970 0.02140 38.295 2.78e-06 ***
## beleaf_50sg_LTRUE -0.10927 0.05243 -2.084 0.106
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.002290794)
##
## Null deviance: 0.0191136 on 5 degrees of freedom
## Residual deviance: 0.0091632 on 4 degrees of freedom
## AIC: -15.879
##
## Number of Fisher Scoring iterations: 2
anth_chem_m11 <- glm(logbruise ~ coragen_L, data = gh6)
Anova(anth_chem_m11)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## coragen_L 1.3037 1 0.2535
summary(anth_chem_m11)
##
## Call:
## glm(formula = logbruise ~ coragen_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.86406 0.06003 14.393 0.000135 ***
## coragen_LTRUE -0.07509 0.06576 -1.142 0.317246
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.003603822)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.014415 on 4 degrees of freedom
## AIC: -13.16
##
## Number of Fisher Scoring iterations: 2
anth_chem_m12 <- glm(logbruise ~ entrust_sc_L, data = gh6)
Anova(anth_chem_m12)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## entrust_sc_L 4.3437 1 0.03715 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(anth_chem_m12)
##
## Call:
## glm(formula = logbruise ~ entrust_sc_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.81970 0.02140 38.295 2.78e-06 ***
## entrust_sc_LTRUE -0.10927 0.05243 -2.084 0.106
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.002290794)
##
## Null deviance: 0.0191136 on 5 degrees of freedom
## Residual deviance: 0.0091632 on 4 degrees of freedom
## AIC: -15.879
##
## Number of Fisher Scoring iterations: 2
anth_chem_m13 <- glm(logbruise ~ pylon_L, data = gh6)
Anova(anth_chem_m13)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## pylon_L 0.65877 1 0.417
summary(anth_chem_m13)
##
## Call:
## glm(formula = logbruise ~ pylon_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.78026 0.03698 21.099 2.98e-05 ***
## pylon_LTRUE 0.04245 0.05230 0.812 0.463
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.00410272)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.016411 on 4 degrees of freedom
## AIC: -12.382
##
## Number of Fisher Scoring iterations: 2
anth_chem_m14 <- glm(logbruise ~ grotto_L, data = gh6)
Anova(anth_chem_m14)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## grotto_L 4.3437 1 0.03715 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(anth_chem_m14)
##
## Call:
## glm(formula = logbruise ~ grotto_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.81970 0.02140 38.295 2.78e-06 ***
## grotto_LTRUE -0.10927 0.05243 -2.084 0.106
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.002290794)
##
## Null deviance: 0.0191136 on 5 degrees of freedom
## Residual deviance: 0.0091632 on 4 degrees of freedom
## AIC: -15.879
##
## Number of Fisher Scoring iterations: 2
anth_chem_m15 <- glm(logbruise ~ luna_tranquility_L, data = gh6)
Anova(anth_chem_m15)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## luna_tranquility_L 0.16524 1 0.6844
summary(anth_chem_m15)
##
## Call:
## glm(formula = logbruise ~ luna_tranquility_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.79354 0.03387 23.429 1.97e-05 ***
## luna_tranquility_LTRUE 0.02385 0.05867 0.406 0.705
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.004588841)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.018355 on 4 degrees of freedom
## AIC: -11.71
##
## Number of Fisher Scoring iterations: 2
anth_chem_m16 <- glm(logbruise ~ previcur_flex_L, data = gh6)
Anova(anth_chem_m16)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## previcur_flex_L 4.3437 1 0.03715 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(anth_chem_m16)
##
## Call:
## glm(formula = logbruise ~ previcur_flex_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.71043 0.04786 14.843 0.00012 ***
## previcur_flex_LTRUE 0.10927 0.05243 2.084 0.10553
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.002290794)
##
## Null deviance: 0.0191136 on 5 degrees of freedom
## Residual deviance: 0.0091632 on 4 degrees of freedom
## AIC: -15.879
##
## Number of Fisher Scoring iterations: 2
anth_chem_m17 <- glm(logbruise ~ fontelis_L, data = gh6)
Anova(anth_chem_m17)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## fontelis_L 1.3037 1 0.2535
summary(anth_chem_m17)
##
## Call:
## glm(formula = logbruise ~ fontelis_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.78897 0.02685 29.388 7.98e-06 ***
## fontelis_LTRUE 0.07509 0.06576 1.142 0.317
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.003603822)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.014415 on 4 degrees of freedom
## AIC: -13.16
##
## Number of Fisher Scoring iterations: 2
anth_chem_m18 <- glm(logbruise ~ quadristop_L, data = gh6)
Anova(anth_chem_m18)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## quadristop_L 0.25289 1 0.615
summary(anth_chem_m18)
##
## Call:
## glm(formula = logbruise ~ quadristop_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.80764 0.02998 26.939 1.13e-05 ***
## quadristop_LTRUE -0.03693 0.07344 -0.503 0.641
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.004494262)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.017977 on 4 degrees of freedom
## AIC: -11.835
##
## Number of Fisher Scoring iterations: 2
anth_chem_m19 <- glm(logbruise ~ milstop_L, data = gh6)
Anova(anth_chem_m19)
## Analysis of Deviance Table (Type II tests)
##
## Response: logbruise
## LR Chisq Df Pr(>Chisq)
## milstop_L 0.25289 1 0.615
summary(anth_chem_m19)
##
## Call:
## glm(formula = logbruise ~ milstop_L, data = gh6)
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.80764 0.02998 26.939 1.13e-05 ***
## milstop_LTRUE -0.03693 0.07344 -0.503 0.641
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.004494262)
##
## Null deviance: 0.019114 on 5 degrees of freedom
## Residual deviance: 0.017977 on 4 degrees of freedom
## AIC: -11.835
##
## Number of Fisher Scoring iterations: 2
summary table
anth_models <- list(
azaguard_L = anth_chem_m1,
botanigard_22wp_L = anth_chem_m2,
botanigard_es_L = anth_chem_m3,
captiva_prime_L = anth_chem_m4,
nofly_L = anth_chem_m5,
venerate_cg_L = anth_chem_m6,
m_pede_L = anth_chem_m7,
rootshield_plus_L = anth_chem_m8,
lalstop_k61_L = anth_chem_m9,
beleaf_50sg_L = anth_chem_m10,
coragen_L = anth_chem_m11,
entrust_sc_L = anth_chem_m12,
pylon_L = anth_chem_m13,
grotto_L = anth_chem_m14,
luna_tranquility_L = anth_chem_m15,
previcur_flex_L = anth_chem_m16,
fontelis_L = anth_chem_m17,
quadristop_L = anth_chem_m18,
milstop_L = anth_chem_m19
)
summary_table.anth <- do.call(rbind, lapply(names(anth_models), function(name){
m <- anth_models[[name]]
s <- summary(m)
coef_row <- s$coefficients[2, ]
lr <- Anova(m)$`LR Chisq`[1]
LRT <- Anova(m)$`Pr(>Chisq)`[1]
data.frame(
Predictor = name,
Estimate = coef_row["Estimate"],
SE = coef_row["Std. Error"],
z = coef_row["t value"],
p = coef_row["Pr(>|t|)"],
LR_ChiSq = lr,
AIC = AIC(m),
LRT = LRT,
Residual_Deviance = m$deviance
)
}))
summary_table.anth$SE <- signif(summary_table.anth$SE, 3)
summary_table.anth$z <- round(summary_table.anth$z, 2)
summary_table.anth$p <- signif(summary_table.anth$p, 3)
summary_table.anth$LR_ChiSq <- round(summary_table.anth$LR_ChiSq, 2)
summary_table.anth$AIC <- round(summary_table.anth$AIC, 2)
summary_table.anth$Residual_Deviance <- round(summary_table.anth$Residual_Deviance, 2)
summary_table.anth$LRT <- signif(summary_table.anth$LRT, 3)
anth_table <- as.data.frame(summary_table.anth)
anth_table
## Predictor Estimate SE z p LR_ChiSq AIC
## Estimate azaguard_L -0.04244796 0.0523 -0.81 0.4630 0.66 -12.38
## Estimate1 botanigard_22wp_L -0.04244796 0.0523 -0.81 0.4630 0.66 -12.38
## Estimate2 botanigard_es_L 0.02054152 0.0590 0.35 0.7450 0.12 -11.65
## Estimate3 captiva_prime_L -0.04244796 0.0523 -0.81 0.4630 0.66 -12.38
## Estimate4 nofly_L -0.04244796 0.0523 -0.81 0.4630 0.66 -12.38
## Estimate5 venerate_cg_L -0.09468312 0.0366 -2.58 0.0611 6.68 -17.36
## Estimate6 m_pede_L -0.04244796 0.0523 -0.81 0.4630 0.66 -12.38
## Estimate7 rootshield_plus_L 0.00376682 0.0598 0.06 0.9530 0.00 -11.47
## Estimate8 lalstop_k61_L 0.07083577 0.0483 1.47 0.2160 2.15 -14.05
## Estimate9 beleaf_50sg_L -0.10927277 0.0524 -2.08 0.1060 4.34 -15.88
## Estimate10 coragen_L -0.07508666 0.0658 -1.14 0.3170 1.30 -13.16
## Estimate11 entrust_sc_L -0.10927277 0.0524 -2.08 0.1060 4.34 -15.88
## Estimate12 pylon_L 0.04244796 0.0523 0.81 0.4630 0.66 -12.38
## Estimate13 grotto_L -0.10927277 0.0524 -2.08 0.1060 4.34 -15.88
## Estimate14 luna_tranquility_L 0.02384735 0.0587 0.41 0.7050 0.17 -11.71
## Estimate15 previcur_flex_L 0.10927277 0.0524 2.08 0.1060 4.34 -15.88
## Estimate16 fontelis_L 0.07508666 0.0658 1.14 0.3170 1.30 -13.16
## Estimate17 quadristop_L -0.03693090 0.0734 -0.50 0.6410 0.25 -11.84
## Estimate18 milstop_L -0.03693090 0.0734 -0.50 0.6410 0.25 -11.84
## LRT Residual_Deviance
## Estimate 0.41700 0.02
## Estimate1 0.41700 0.02
## Estimate2 0.72800 0.02
## Estimate3 0.41700 0.02
## Estimate4 0.41700 0.02
## Estimate5 0.00976 0.01
## Estimate6 0.41700 0.02
## Estimate7 0.95000 0.02
## Estimate8 0.14200 0.01
## Estimate9 0.03710 0.01
## Estimate10 0.25400 0.01
## Estimate11 0.03710 0.01
## Estimate12 0.41700 0.02
## Estimate13 0.03710 0.01
## Estimate14 0.68400 0.02
## Estimate15 0.03710 0.01
## Estimate16 0.25400 0.01
## Estimate17 0.61500 0.02
## Estimate18 0.61500 0.02
plot of anther bruise vs pesticides applied/not applied
plots <- lapply(1:nrow(chem_info), function(i){
chem_type <- chem_info$type[i]
chem_mode <- chem_info$mode[i]
chem_class <- chem_info$class[i]
var <- chem_info$pesticide[i]
var2 <- gsub("_L$", "", chem_info$pesticide[i]) # remove _L at end
ggplot(gh6, aes_string(x = var, y = "average_bruise")) +
geom_boxplot(fill = chem_colors[chem_class],
alpha = 0.6,
width = 0.5,
outlier.shape = NA) +
geom_jitter(width = 0.1, size = 3, alpha = 0.8) +
labs(
title = gsub("_L","",var),
subtitle = paste(chem_type, "-", chem_mode),
x = "",
y = "Average Anther Score"
) +
scale_x_discrete(labels = c("FALSE"="Not applied","TRUE"="Applied")) +
theme_classic(base_size = 16)
})
wrap_plots(plots, ncol = 4)
