Get data - global.R

setwd("/Documents/GRADUATE_SCHOOL/Projects/metalab/metalab/write-ups/synthesis_paper/paper_analyses")
ld.df = data.frame(dataset = datasets$name,
                   domain = c("Prosody", "Words", "Communication", "Sounds",
                              "Sounds", "Sounds", "Sounds", "Sounds", "Words",
                              "Words", "Communication", "Words"))

ld.df$domain = factor(ld.df$domain, levels = c("Prosody","Sounds", "Words", "Communication"))

# Fit models using raw d_calc
raw_data = all_data %>%
          filter(dataset != "Gaze following" | d_calc > 0) %>% # remove negative ES for gaze following
          filter(dataset != "Pointing and vocabulary") %>%
          left_join(ld.df)

Unweighted logarithmic

un.l = ggplot(raw_data,
            aes(x = mean_age/365, y = d_calc, col = dataset)) +
  facet_grid(~ domain) +
  geom_point(aes(size = n, y = d_calc), alpha = .1) + 
  geom_smooth(method="lm", se = FALSE, size = 1, formula = y ~ log(x)) + 
  ylim(c(-.3, 2))+
  xlim(0,3.1) +
  geom_hline(yintercept = 0, linetype = "dashed") +
  xlab("Age (years)") + 
  ylab("Effect size (d)") + 
  theme_bw() + 
  theme(legend.position = "none",
        legend.key = element_blank(),
        legend.background = element_rect(fill = "transparent"),
        axis.line= element_line(size = 3),
        axis.text = element_text(colour = "black", size = 18),
        strip.text.x = element_text(size = 14),
        axis.title = element_text(colour = "black", size = 18))

direct.label(un.l, list("top.bumptwice", dl.trans(y = y + 0.1), cex=.7))

Weighted logarithmic

weight = 1/n

wl = ggplot(raw_data,
            aes(x = mean_age/365, y = d_calc, col = dataset, weight = 1/n)) +
  facet_grid(~ domain) +
  geom_point(aes(size = n, y = d_calc), alpha = .1) + 
  geom_smooth(method="lm", se = FALSE, size = 1, formula = y ~ log(x)) + 
  ylim(c(-.3, 2))+
  xlim(0,3.1) +
  geom_hline(yintercept = 0, linetype = "dashed") +
  xlab("Age (years)") + 
  ylab("Effect size (d)") + 
  theme_bw() + 
  theme(legend.position = "none",
        legend.key = element_blank(),
        legend.background = element_rect(fill = "transparent"),
        axis.line= element_line(size = 3),
        axis.text = element_text(colour = "black", size = 18),
        strip.text.x = element_text(size = 14),
        axis.title = element_text(colour = "black", size = 18))

direct.label(wl, list("top.bumptwice", dl.trans(y = y + 0.1), cex=.7))

Unweighted loess

uwloess = ggplot(raw_data,
            aes(x = mean_age/365, y = d_calc, col = dataset)) +
  facet_grid(~ domain) +
  geom_point(aes(size = n, y = d_calc), alpha = .1) + 
  geom_smooth(method="loess", se = FALSE, size = 1) + 
  ylim(c(-.3, 2))+
  xlim(0,3.1) +
  geom_hline(yintercept = 0, linetype = "dashed") +
  xlab("Age (years)") + 
  ylab("Effect size (d)") + 
  theme_bw() + 
  theme(legend.position = "none",
        legend.key = element_blank(),
        legend.background = element_rect(fill = "transparent"),
        axis.line= element_line(size = 3),
        axis.text = element_text(colour = "black", size = 18),
        strip.text.x = element_text(size = 14),
        axis.title = element_text(colour = "black", size = 18))

direct.label(uwloess, list("top.bumptwice", dl.trans(y = y + 0.1), cex=.7))

Weighted loess

weight = 1/n

wloess = ggplot(raw_data,
            aes(x = mean_age/365, y = d_calc, col = dataset, weight = 1/n)) +
  facet_grid(~ domain) +
  geom_point(aes(size = n, y = d_calc), alpha = .1) + 
  geom_smooth(method="loess", se = FALSE, size = 1) + 
  ylim(c(-.3, 2))+
  xlim(0,3.1) +
  geom_hline(yintercept = 0, linetype = "dashed") +
  xlab("Age (years)") + 
  ylab("Effect size (d)") + 
  theme_bw() + 
  theme(legend.position = "none",
        legend.key = element_blank(),
        legend.background = element_rect(fill = "transparent"),
        axis.line= element_line(size = 3),
        axis.text = element_text(colour = "black", size = 18),
        strip.text.x = element_text(size = 14),
        axis.title = element_text(colour = "black", size = 18))

direct.label(wloess, list("top.bumptwice", dl.trans(y = y + 0.1), cex=.7))

Unweighted loess with medium span

span = 1

loess.med = ggplot(raw_data,
            aes(x = mean_age/365, y = d_calc, col = dataset)) +
  facet_grid(~ domain) +
  geom_point(aes(size = n, y = d_calc), alpha = .1) + 
  geom_smooth(method="loess",span = 1, se = FALSE, size = 1) + 
  ylim(c(-.3, 2))+
  xlim(0,3.1) +
  geom_hline(yintercept = 0, linetype = "dashed") +
  xlab("Age (years)") + 
  ylab("Effect size (d)") + 
  theme_bw() + 
  theme(legend.position = "none",
        legend.key = element_blank(),
        legend.background = element_rect(fill = "transparent"),
        axis.line= element_line(size = 3),
        axis.text = element_text(colour = "black", size = 18),
        strip.text.x = element_text(size = 14),
        axis.title = element_text(colour = "black", size = 18))

direct.label(loess.med, list("top.bumptwice", dl.trans(y = y + 0.1), cex=.7))

Unweighted loess with giant span

span = 100

loess.giant = ggplot(raw_data,
            aes(x = mean_age/365, y = d_calc, col = dataset)) +
  facet_grid(~ domain) +
  geom_point(aes(size = n, y = d_calc), alpha = .1) + 
  geom_smooth(method="loess",span = 100, se = FALSE, size = 1) + 
  ylim(c(-.3, 2))+
  xlim(0,3.1) +
  geom_hline(yintercept = 0, linetype = "dashed") +
  xlab("Age (years)") + 
  ylab("Effect size (d)") + 
  theme_bw() + 
  theme(legend.position = "none",
        legend.key = element_blank(),
        legend.background = element_rect(fill = "transparent"),
        axis.line= element_line(size = 3),
        axis.text = element_text(colour = "black", size = 18),
        strip.text.x = element_text(size = 14),
        axis.title = element_text(colour = "black", size = 18))

direct.label(loess.giant, list("top.bumptwice", dl.trans(y = y + 0.1), cex=.7))

Weighted loess with giant span

span = 100 weight = 1/n

loess.giant = ggplot(raw_data,
            aes(x = mean_age/365, y = d_calc, col = dataset, weight = 1/n)) +
  facet_grid(~ domain) +
  geom_point(aes(size = n, y = d_calc), alpha = .1) + 
  geom_smooth(method="loess",span = 100, se = FALSE, size = 1) + 
  ylim(c(-.3, 2))+
  xlim(0,3.1) +
  geom_hline(yintercept = 0, linetype = "dashed") +
  xlab("Age (years)") + 
  ylab("Effect size (d)") + 
  theme_bw() + 
  theme(legend.position = "none",
        legend.key = element_blank(),
        legend.background = element_rect(fill = "transparent"),
        axis.line= element_line(size = 3),
        axis.text = element_text(colour = "black", size = 18),
        strip.text.x = element_text(size = 14),
        axis.title = element_text(colour = "black", size = 18))

direct.label(loess.giant, list("top.bumptwice", dl.trans(y = y + 0.1), cex=.7))