Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- qf6 = readRDS("~/Downloads/qf6.rds") %>%
- na.omit %>% arrange(name, month) %>% group_by(name) %>% mutate(lmonth=lead(month), lperc=lead(perc), lntot=lead(ntot))
- avg = function(x, leadx, frac) ifelse(frac == 0, x, frac * leadx + (1-frac) * x)
- steps = tibble(day=0:31)
- qf6b = crossing(qf6, steps) %>% mutate(x=month + day) %>% filter(day == 0 | (!is.na(lmonth) & x < lmonth)) %>%
- group_by(name, month) %>% mutate(frac=day/(max(day)+1)) %>% mutate(perc2=avg(perc, lperc, frac), ntot2=avg(ntot, lntot, frac))
- ggplot(data=qf6b, aes(x=x, y=perc2, group=name, colour=factor(name), alpha=ifelse(ntot2>100, 100, ntot2))) +
- geom_line(size=1.2)+
- scale_x_date(date_breaks = "2 month") +
- xlab("Datum") +
- ylab("% Bron in het nieuws")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement