generate_report.Rmd

---
title: "MPGQuestionnaireReport"
author: "Kleemeyer"
date: "`r format(Sys.time(), '%d %B, %Y')`"
output:
#  pdf_document: default
  html_document: default
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, message = FALSE)
library(knitr)
library(kableExtra)
library(ggplot2)
library(dplyr)
library(corrr)
library(umx)
library(matrixStats)

source("get_data_exported.R")
```
# Scale consistency

## SCALE A: Group atmosphere

```{r echo = FALSE, results='asis'}
interitem_A <- responses_coded %>%
  select(A03:A05) %>%
  correlate() 

kable(interitem_A, caption = "Interitem correlation scale A")

scale_A <- responses_coded %>%
  select(A03:A05) 

items <- colnames(scale_A)

rely_A <- psych::alpha(scale_A, use = "complete.obs", check.keys=TRUE)

output_A <- as.data.frame(bind_cols(rely_A$item.stats$n, rely_A$item.stats$mean, rely_A$item.stats$sd, rely_A$item.stats$r.drop, rely_A$alpha.drop$raw_alpha))
colnames(output_A) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_A) <- c(items)

kable(output_A, caption = "Item statistics scale A")

```

Cronbach's alpha for Scale A is `r rely_A$total$std.alpha`. 

## SCALE B1: Leadership leadstyle

```{r echo = FALSE, results='asis'}
interitem_B1 <- responses_coded %>%
  select(B02:B07) %>%
  correlate() 

kable(interitem_B1, caption = "Interitem correlation scale B1")

scale_B1 <- responses_coded %>%
  select(B02:B07) 

items <- colnames(scale_B1)

rely_B1 <- psych::alpha(scale_B1, use = "complete.obs", check.keys=TRUE)

output_B1 <- as.data.frame(bind_cols(rely_B1$item.stats$n, rely_B1$item.stats$mean, rely_B1$item.stats$sd, rely_B1$item.stats$r.drop, rely_B1$alpha.drop$raw_alpha))
colnames(output_B1) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_B1) <- c(items)

kable(output_B1, caption = "Item statistics scale B1")

```

Cronbach's alpha for Scale B1 is `r rely_B1$total$std.alpha`. 

## SCALE B2: Leadership development

```{r echo = FALSE, results='asis'}
interitem_B2 <- responses_coded %>%
  select(B09:B13) %>%
  correlate() 

kable(interitem_B2, caption = "Interitem correlation scale B2")

scale_B2 <- responses_coded %>%
  select(B09:B13) 

items <- colnames(scale_B2)

rely_B2 <- psych::alpha(scale_B2, use = "complete.obs", check.keys=TRUE)

output_B2 <- as.data.frame(bind_cols(rely_B2$item.stats$n, rely_B2$item.stats$mean, rely_B2$item.stats$sd, rely_B2$item.stats$r.drop, rely_B2$alpha.drop$raw_alpha))
colnames(output_B2) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_B2) <- c(items)

kable(output_B2, caption = "Item statistics scale B2")

```

Cronbach's alpha for Scale B2 is `r rely_B2$total$std.alpha`. 

## SCALE C: Commitment

```{r echo = FALSE, results='asis'}
interitem_C <- responses_coded %>%
  select(C02:C04) %>%
  correlate() 

kable(interitem_C, caption = "Interitem correlation scale C")

scale_C <- responses_coded %>%
  select(C02:C04) 

items <- colnames(scale_C)

rely_C <- psych::alpha(scale_C, use = "complete.obs", check.keys=TRUE)

output_C <- as.data.frame(bind_cols(rely_C$item.stats$n, rely_C$item.stats$mean, rely_C$item.stats$sd, rely_C$item.stats$r.drop, rely_C$alpha.drop$raw_alpha))
colnames(output_C) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_C) <- c(items)

kable(output_C, caption = "Item statistics scale C")

```

Cronbach's alpha for Scale C is `r rely_C$total$std.alpha`. 

## SCALE D: Bullying

```{r echo = FALSE, results='asis'}
interitem_D <- responses_coded %>%
  select(D02:D06) %>%
  correlate() 

kable(interitem_D, caption = "Interitem correlation scale D")

scale_D <- responses_coded %>%
  select(D02:D06) 

items <- colnames(scale_D)

rely_D <- psych::alpha(scale_D, use = "complete.obs", check.keys=TRUE)

output_D <- as.data.frame(bind_cols(rely_D$item.stats$n, rely_D$item.stats$mean, rely_D$item.stats$sd, rely_D$item.stats$r.drop, rely_D$alpha.drop$raw_alpha))
colnames(output_D) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_D) <- c(items)

kable(output_D, caption = "Item statistics scale D")

```

Cronbach's alpha for Scale D is `r rely_D$total$std.alpha`. 

## SCALE E: Gender discrimination

```{r echo = FALSE, results='asis'}
interitem_E <- responses_coded %>%
  select(E03:E05) %>%
  correlate() 

kable(interitem_E, caption = "Interitem correlation scale E")

scale_E <- responses_coded %>%
  select(E03:E05) 

items <- colnames(scale_E)

rely_E <- psych::alpha(scale_E, use = "complete.obs", check.keys=TRUE)

output_E <- as.data.frame(bind_cols(rely_E$item.stats$n, rely_E$item.stats$mean, rely_E$item.stats$sd, rely_E$item.stats$r.drop, rely_E$alpha.drop$raw_alpha))
colnames(output_E) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_E) <- c(items)

kable(output_E, caption = "Item statistics scale E")

```

Cronbach's alpha for Scale E is `r rely_E$total$std.alpha`. 

## SCALE F: Sexual harassment

```{r echo = FALSE, results='asis'}
interitem_F <- responses_coded %>%
  select(F02:F08) %>%
  correlate() 

kable(interitem_F, caption = "Interitem correlation scale F")

scale_F <- responses_coded %>%
  select(F02:F08) 

items <- colnames(scale_F)

rely_F <- psych::alpha(scale_F, use = "complete.obs", check.keys=TRUE)

output_F <- as.data.frame(bind_cols(rely_F$item.stats$n, rely_F$item.stats$mean, rely_F$item.stats$sd, rely_F$item.stats$r.drop, rely_F$alpha.drop$raw_alpha))
colnames(output_F) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_F) <- c(items)

kable(output_F, caption = "Item statistics scale F")

```

Cronbach's alpha for Scale F is `r rely_F$total$std.alpha`. 

## SCALE H1: Work-life Balance

```{r echo = FALSE, results='asis'}
interitem_H1 <- responses_coded %>%
  select(H02:H06) %>%
  correlate() 

kable(interitem_H1, caption = "Interitem correlation scale H1")

scale_H1 <- responses_coded %>%
  select(H02:H06) 

items <- colnames(scale_H1)

rely_H1 <- psych::alpha(scale_H1, use = "complete.obs", check.keys=TRUE)

output_H1 <- as.data.frame(bind_cols(rely_H1$item.stats$n, rely_H1$item.stats$mean, rely_H1$item.stats$sd, rely_H1$item.stats$r.drop, rely_H1$alpha.drop$raw_alpha))
colnames(output_H1) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_H1) <- c(items)

kable(output_H1, caption = "Item statistics scale H1")

```

Cronbach's alpha for Scale H1 is `r rely_H1$total$std.alpha`. 

## SCALE H2: Work-family Balance

```{r echo = FALSE, results='asis'}
interitem_H2 <- responses_coded %>%
  select(H09:H13) %>%
  correlate() 

kable(interitem_H2, caption = "Interitem correlation scale H2")

scale_H2 <- responses_coded %>%
  select(H09:H13) 

items <- colnames(scale_H2)

rely_H2 <- psych::alpha(scale_H2, use = "complete.obs", check.keys=TRUE)

output_H2 <- as.data.frame(bind_cols(rely_H2$item.stats$n, rely_H2$item.stats$mean, rely_H2$item.stats$sd, rely_H2$item.stats$r.drop, rely_H2$alpha.drop$raw_alpha))
colnames(output_H2) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_H2) <- c(items)

kable(output_H2, caption = "Item statistics scale H2")

```

Cronbach's alpha for Scale H2 is `r rely_H2$total$std.alpha`. 

## SCALE H3: Parental leave - yes

```{r echo = FALSE, results='asis'}
interitem_H3 <- responses_coded %>%
  select(H17:H21) %>%
  correlate() 

kable(interitem_H3, caption = "Interitem correlation scale H3")

scale_H3 <- responses_coded %>%
  select(H17:H21) 

items <- colnames(scale_H3)

rely_H3 <- psych::alpha(scale_H3, use = "complete.obs", check.keys=TRUE)

output_H3 <- as.data.frame(bind_cols(rely_H3$item.stats$n, rely_H3$item.stats$mean, rely_H3$item.stats$sd, rely_H3$item.stats$r.drop, rely_H3$alpha.drop$raw_alpha))
colnames(output_H3) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_H3) <- c(items)

kable(output_H3, caption = "Item statistics scale H3")

```

Cronbach's alpha for Scale H3 is `r rely_H3$total$std.alpha`. 

## SCALE H4: Parental leave - no

```{r echo = FALSE, results='asis'}
interitem_H4 <- responses_coded %>%
  select(H23:H25) %>%
  correlate() 

kable(interitem_H4, caption = "Interitem correlation scale H4")

scale_H4 <- responses_coded %>%
  select(H23:H25) 

items <- colnames(scale_H4)

rely_H4 <- psych::alpha(scale_H4, use = "complete.obs", check.keys=TRUE)

output_H4 <- as.data.frame(bind_cols(rely_H4$item.stats$n, rely_H4$item.stats$mean, rely_H4$item.stats$sd, rely_H4$item.stats$r.drop, rely_H4$alpha.drop$raw_alpha))
colnames(output_H4) <- c("total", "mean", "SD", "Corrected Item total r", "Alpha if item deleted")
rownames(output_H4) <- c(items)

kable(output_H4, caption = "Item statistics scale H4")

```

Cronbach's alpha for Scale H4 is `r rely_H4$total$std.alpha`. 

# Factor Analyses

```{r echo = FALSE, message = FALSE, include = FALSE}
likert <- responses_coded %>%
  select(A03:B13,C02:D06,E03:E05,F02:F08,H02:H06)

efa <- umxEFA(likert, factors = 8, minManifests = 5)
factors <- as.data.frame(loadings(efa))

color.me <- rowMaxs(abs(loadings(efa)))
```

```{r echo = FALSE, message = FALSE}

factors %>% 
  mutate(items = row.names(.)) %>%
  relocate(items) %>%
  mutate_all(~cell_spec(
    .x, 
    color = ifelse(abs(.x) %in% color.me, "red", "black"))) %>%
  kbl(format = "html", escape = F) %>%
  kable_paper("striped") %>%
  pack_rows("Scale A", 1, 3) %>%
  pack_rows("Scale B", 4, 14) %>%
  pack_rows("Scale C", 15, 17) %>%
  pack_rows("Scale D", 18, 22) %>%
  pack_rows("Scale E", 23, 25) %>%
  pack_rows("Scale F", 26, 32) %>%
  pack_rows("Scale H", 33, 37) %>%

  kable_styling()

```

# General Analyses

In total, `r nrow(responses)` participants (partly) completed the questionnaire.

# Binary Questions

* Do you work in a group?
```{r echo = FALSE, results='asis'}
kable(table(responses$A01, useNA = "ifany"))
```

* Did you have a personal conversation with your supervisor?
```{r echo = FALSE, results='asis'}
kable(table(responses$B14, useNA = "ifany"))
```

* Have you been subjected to bullying during the past 12 months?
```{r echo = FALSE, results='asis'}
kable(table(responses$D12, useNA = "ifany"))
```

* Have you been subjected to gender discrimination during the past 12 months?
```{r echo = FALSE, results='asis'}
kable(table(responses$E05, useNA = "ifany"))
```

* Have you been subjected to sexual harassment during the past 12 months?
```{r echo = FALSE, results='asis'}
kable(table(responses$F11, useNA = "ifany"))
```

* Do children under the age of 18 live in your household?
```{r echo = FALSE, results='asis'}
kable(table(responses$H07, useNA = "ifany"))
```

* Do people with care needs live in your household?
```{r echo = FALSE, results='asis'}
kable(table(responses$H30, useNA = "ifany"))
```

* Have you taken parental leave?
```{r echo = FALSE, results='asis'}
kable(table(responses$H14, useNA = "ifany"))
```

* What is your gender?
```{r echo = FALSE, results='asis'}
kable(table(responses$I02, useNA = "ifany"))
```

* What is your position in the MPG?
```{r echo = FALSE, results='asis'}
kable(table(responses$I03, useNA = "ifany"))
```