IAT: 0.15, 0.35, and 0.65 are considered small, medium, and large level of bias for individual scores.
Positive means bias towards arts / against Math.
iat = read_csv(here::here(params$arquivo_dados), col_types = "cccdc")
iat = iat %>%
mutate(sex = factor(sex, levels = c("m", "f"), ordered = TRUE))
glimpse(iat)
## Observations: 91
## Variables: 5
## $ session_id <chr> "2399787", "2399902", "2400037", "2400210", "2400235…
## $ referrer <chr> "mcdaniel", "mcdaniel", "mcdaniel", "mcdaniel", "mcd…
## $ sex <ord> f, m, f, f, f, f, f, m, f, m, f, f, f, f, f, f, f, f…
## $ d_art <dbl> 0.994230746, 0.845950103, 1.130822093, 0.408574905, …
## $ iat_exclude <chr> "Include", "Include", "Include", "Include", "Include…
iat %>%
ggplot(aes(x = d_art, fill = sex, color = sex)) +
geom_histogram(binwidth = .2, alpha = .4) +
geom_rug() +
facet_grid(sex ~ ., scales = "free_y") +
theme(legend.position = "None")
iat %>%
ggplot(aes(x = sex, y = d_art)) +
geom_quasirandom(width = .1)
iat %>%
ggplot(aes(x = sex, y = d_art)) +
geom_quasirandom(width = .1) +
stat_summary(geom = "point", fun.y = "mean", color = "red", size = 5)
#cont_m = iat %>%
# filter(sex == "m")
#View(cont_m)
#cont_f = iat %>%
# filter(sex == "f")
#View(cont_f)
iat %>%
group_by(sex) %>%
summarise(media = mean(d_art), desvio = sd(d_art))
## # A tibble: 2 x 3
## sex media desvio
## <ord> <dbl> <dbl>
## 1 m 0.251 0.538
## 2 f 0.553 0.458
agrupado = iat %>%
group_by(sex) %>%
summarise(media = mean(d_art))
m = agrupado %>% filter(sex == "m") %>% pull(media)
f = agrupado %>% filter(sex == "f") %>% pull(media)
m - f
## [1] -0.3028058
library(boot)
theta <- function(d, i) {
agrupado = d %>%
slice(i) %>%
group_by(sex) %>%
summarise(media = mean(d_art))
m = agrupado %>% filter(sex == "m") %>% pull(media)
f = agrupado %>% filter(sex == "f") %>% pull(media)
m - f
}
booted <- boot(data = iat,
statistic = theta,
R = 2000)
ci = tidy(booted,
conf.level = .95,
conf.method = "bca",
conf.int = TRUE)
glimpse(ci)
## Observations: 1
## Variables: 5
## $ statistic <dbl> -0.3028058
## $ bias <dbl> -0.0001595605
## $ std.error <dbl> 0.1308302
## $ conf.low <dbl> -0.5495447
## $ conf.high <dbl> -0.02722118
ci %>%
ggplot(aes(
x = "",
y = statistic,
ymin = conf.low,
ymax = conf.high
)) +
geom_pointrange() +
geom_point(size = 3) +
labs(x = "Diferença",
y = "IAT homens - mulheres")
p1 = iat %>%
ggplot(aes(x = sex, y = d_art)) +
geom_quasirandom(width = .1) +
stat_summary(geom = "point", fun.y = "mean", color = "red", size = 5)
p2 = ci %>%
ggplot(aes(
x = "",
y = statistic,
ymin = conf.low,
ymax = conf.high
)) +
geom_pointrange() +
geom_point(size = 3) +
ylim(-1, 1) +
labs(x = "Diferença",
y = "IAT homens - mulheres")
grid.arrange(p1, p2, ncol = 2)
Em média, as mulheres que participaram do experimento tiveram uma associação implícita (medida pelo IAT) com a matemática negativa e média (média 0.553345, desv. padrão 0.4577992, N = 71). Homens tiveram uma associação negativa com a matemática e menor que a das mulheres (0.250539 , desv. padrão 0.5379578, N = 20). Houve portanto uma considerável diferença entre homens e mulheres (media homens - média mulheres = -0.3028058, 95% CI [-0.5535792, -0.02694593]). A partir desta amostra, estimamos que mulheres têm uma associação negativa, com uma diferença que provavelmente está entre 0.55 e 0.03 pontos na escala IAT, em que a partir dos dados não é possível dizer se a distancia é grande ou pequena. É necessário coletar mais dados para determinar se a diferença é relevante ou negligenciável.
Exemplos de possíveis conclusões para completar