31:T60a,graphic = { const plot = addTooltips( Plot.plot({ height: 600, width: width, marginLeft: 160, marginRight: 80, x: { inset: 30, domain: [0, 1], axis: "top", ticks: 5, tickPadding: 10, tickSize: 5, tickFormat: "%", label: null }, y: { axis: null, label: null }, facet: { data: data, y: "region", marginLeft: 190, grid: true }, marks: [ Plot.dot(data, { x: "proportion", y: "x", fill: "#cb1a1a", fillOpacity: 0.2, stroke: "#cb1a1a", title : (d) => d.msoa + "\n" + d3.format(".0%")(d["proportion"]) }) ], style: { fontFamily: "'Source Sans Pro', sans-serif;", fontSize: "14", fontWeight: 600, color: "#707070", overflow: "visible" } }) ) return wrap( html`

Population with at least one COVID-19 vaccine dose

MSOAs in England, 8 December 2020 to 20 June 2021

`, plot, html`

Source: Public Health England

` ); }32:T635,```r # NIMS vaccination counts and populations by MSOA # Source: Public Health England # URL: https://www.england.nhs.uk/statistics/statistical-work-areas/covid-19-vaccinations # load libraries ----------------------- library(tidyverse) ; library(httr) ; library(readxl) ; library(ggforce) # retrieve data ----------------------- tmp <- tempfile(fileext = ".xlsx") GET(url = url <- "https://www.england.nhs.uk/statistics/wp-content/uploads/sites/2/2021/06/COVID-19-weekly-announced-vaccinations-24-June-2021.xlsx", write_disk(tmp)) vaccinations <- read_xlsx(tmp, sheet = "MSOA", skip = 12) %>% select(msoa_code = 5, msoa = 6, region = 2, 7:19) %>% filter(!is.na(region)) %>% rowwise() %>% mutate(first_dose = sum(c_across(4:16))) %>% select(msoa_code, msoa, region, first_dose) population <- read_xlsx(tmp, sheet = "Population estimates (NIMS)", range = "V15:AK6806") %>% mutate(population = rowSums(.[3:16])) %>% select(msoa_code = 1, population) # join data ----------------------- df <- left_join(vaccinations, population, by = "msoa_code") %>% mutate(proportion = round(first_dose/population,3)) # plot ----------------------- p <- ggplot(df, aes(x = region, y = proportion)) + geom_sina() + coord_flip() + facet_wrap(~region, scales = "free_y", ncol = 1) p # extract coordinates from geom_sina() ----------------------- left_join(rowid_to_column(ggplot_build(p)$plot$data, "id"), rowid_to_column(ggplot_build(p)$data[[1]], "id"), by = "id") %>% select(msoa, region, proportion, x) %>% write_csv("sina_data.csv") ```33:[{"id":"466da02ba7b8f9b98214dacbb8d88741a0fb28ad25b290bfe0d26b1fd3b396740480f0ce5aa3d077d8f936e9e1490d353c4d9d2859e4c6e181a84ba366e4bbeb","href":"https://static.observableusercontent.com/files/466da02ba7b8f9b98214dacbb8d88741a0fb28ad25b290bfe0d26b1fd3b396740480f0ce5aa3d077d8f936e9e1490d353c4d9d2859e4c6e181a84ba366e4bbeb","name":"sina_data.csv","type":"text/csv","lastModified":1625056523666.06,"size":379140}] 2f:{"title":"COVID-19 vaccine take-up","cells":[{"id":0,"value":"# COVID-19 vaccine take-up\n\nRecreation of a graphic published in [the Guardian](https://www.theguardian.com/world/2021/jun/14/covid-vaccines-are-working-but-charts-show-right-delay-reopening-england) newspaper as a [SinaPlot](https://www.tandfonline.com/doi/abs/10.1080/10618600.2017.1366914). The width of the jitter was calculated using the `geom_sina()` function in the [ggforce](https://ggforce.data-imaginist.com/index.html) R package.\n\n
","pinned":false,"hidden":false,"mode":"md"},{"id":5,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":79,"value":"---\n## Appendix","pinned":false,"hidden":false,"mode":"md"},{"id":2,"value":"import {Plot, d3} from \"@observablehq/plot\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":41,"value":"import { addTooltips } from \"@mkfreeman/plot-tooltip\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":132,"value":"wrap = (...elems) => {\n const div = document.createElement(\"div\");\n elems.forEach(el => div.appendChild(el));\n return div;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":231,"value":"html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":358,"value":"$32","pinned":false,"hidden":false,"mode":"md"},{"id":7,"value":"data = d3.csvParse(await FileAttachment(\"sina_data.csv\").text(), d3.autoType)","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}