Created
November 9, 2025 17:42
-
-
Save bbolker/d6a61d4e8a05616aa3b754583c1603cf to your computer and use it in GitHub Desktop.
scraped data from Edwards et al 2021 PNAS on aerosol/COVID-19
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| ## scraped from Edwards et al 2021 PNAS "Exhaled aerosol increases with COVID-19 infection, age, and obesity" doi: 10.1073/pnas.2021830118 | |
| library(ggplot2); theme_set(theme_bw()) | |
| dd <- read.table(header = TRUE, text = " | |
| BMI_years exhaled_particles super_spreader | |
| 200.035325376 50.0311243148 n | |
| 353.718561745 26.9914433584 n | |
| 362.042435728 75.3702213313 n | |
| 374.750834553 7.02162110504 n | |
| 386.246395482 30.3918564146 n | |
| 393.241199689 0.387986576661 n | |
| 399.646108108 40.4246277275 n | |
| 411.173575829 110.489387877 n | |
| 416.220926266 30.4518411814 n | |
| 432.148872224 7.13648555207 n | |
| 446.179503655 7.16456352802 n | |
| 469.137579195 5.54284770539 n | |
| 486.363827666 15.5832766481 n | |
| 488.97524656 103.974446611 n | |
| 497.849132842 23.9445575447 n | |
| 505.512460286 38.9688272173 n | |
| 518.883684706 7.31005849426 n | |
| 526.548151678 24.0019897683 n | |
| 535.539409369 115.74124565 n | |
| 536.740851981 7.34579410001 n | |
| 548.236412911 30.7160294096 n | |
| 560.335543967 4.05769294758 n | |
| 569.917836975 27.4240994422 n | |
| 575.686128948 69.1271259229 n | |
| 589.709923209 59.1492342902 n | |
| 594.140029181 9.12832014849 n | |
| 581.388328283 14.1057795201 n | |
| 619.035302737 42.531326771 n | |
| 647.709251951 5.90020376284 n | |
| 659.196836183 17.5968078623 n | |
| 659.223045333 55.9530246954 n | |
| 693.635658786 17.6657265305 n | |
| 711.492826062 17.7014621363 n | |
| 726.206794597 84.4372804414 n | |
| 731.288330881 54.4295817886 n | |
| 746.576241809 27.7776266848 n | |
| 756.766663056 7.78610781365 n | |
| 770.808689769 24.490801804 n | |
| 781.033296863 54.529130976 n | |
| 786.102298337 6.17715470737 n | |
| 816.112154497 57.9346491187 n | |
| 832.661903033 11.2733073409 n | |
| 846.690255408 7.96606211401 n | |
| 870.323691354 61.3784554105 n | |
| 872.853064214 29.6979900697 n | |
| 893.307976044 98.1129564209 n | |
| 910.500038668 58.1235373205 n | |
| 930.313775939 121.534242596 n | |
| 924.524972458 49.8133072893 n | |
| 934.705137951 14.8128340052 n | |
| 944.943419384 64.8631023945 n | |
| 951.949618873 51.535848571 n | |
| 976.798171771 16.5647295345 n | |
| 979.402753495 94.9499298886 n | |
| 982.584696196 84.9503416382 n | |
| 995.307908888 38.2813422306 n | |
| 1012.55125028 73.3366951949 n | |
| 1016.39943716 105.029923252 n | |
| 1033.56187205 21.6813025142 n | |
| 1032.29661586 36.6877043839 n | |
| 1056.55071486 65.0864499304 n | |
| 1064.81191479 21.7438398243 n | |
| 1085.23833841 48.4672661396 n | |
| 1090.95193303 10.1225357512 n | |
| 1115.20831108 41.8566045006 n | |
| 1171.977709 55.3114854875 n | |
| 1196.18166876 10.3331205707 n | |
| 1198.74522747 28.6825032731 n | |
| 1245.95854153 57.1271945985 n | |
| 1253.58882266 23.7892778293 n | |
| 1281.09209652 140.58046961 n | |
| 1318.62853673 7.24284152155 n | |
| 1327.57763188 37.2786181503 n | |
| 1332.65802864 5.60325789606 n | |
| 1337.14055291 32.2947774205 n | |
| 1354.38845242 74.0207767906 n | |
| 1377.31348164 24.0368745262 n | |
| 1389.47300769 85.7646029405 n | |
| 1401.56416205 47.4326352685 n | |
| 1421.96209747 32.4645215478 n | |
| 1426.42752881 34.1411170507 n | |
| 1469.79037694 27.557257116 n | |
| 1499.13170987 34.2866120169 n | |
| 1514.46634145 76.0087825721 n | |
| 1576.3389267 91.1415353335 n | |
| 1624.77533441 42.876345536 n | |
| 1708.32364608 46.3788603229 n | |
| 1758.11419319 113.184873568 n | |
| 1777.22522138 81.5375912895 n | |
| 1861.39761468 65.0294431308 n | |
| 1918.83097773 116.841817222 n | |
| 2064.90102589 152.154977057 n | |
| 2083.33555415 63.8059240579 n | |
| 2653.47572068 44.9349716812 n | |
| 2306.18929448 469.493112008 y | |
| 2032.74012025 685.741599664 y | |
| 1954.56354473 143.595874057 y | |
| 1839.5570365 768.737969431 y | |
| 1698.59131518 736.770342972 y | |
| 1552.14864128 156.131839468 y | |
| 1492.86582389 197.704686242 y | |
| 1192.45427187 155.412022266 y | |
| 1173.51759151 442.211529564 y | |
| 1152.914541 157.000555026 y | |
| 1132.53255898 195.315931167 y | |
| 1122.31023094 168.612925198 y | |
| 947.797557793 508.466193464 y | |
| 892.195416635 336.586012883 y | |
| 885.749485199 236.513554081 y | |
| 761.964431221 147.879892508 y | |
| 726.886713116 146.142035966 y | |
| 718.004850136 214.498293822 y | |
| 630.582142438 140.946334145 y | |
| 477.51842959 137.30470575 y | |
| ") | |
| ggplot(dd, aes(BMI_years,exhaled_particles, colour = super_spreader, | |
| shape = super_spreader)) + | |
| geom_point() + | |
| scale_y_log10() + | |
| scale_colour_manual(values = c("black", "red")) + | |
| geom_smooth(method = "lm") + | |
| geom_hline(yintercept = 156, lty = 2) | |
| ggsave("pnas_aerosol.png") |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment