# R Dataset / Package HistData / Snow.streets

Webform
Category

Webform
Category

Webform
Category

Webform
Category

Webform
Category

Webform
Category

## Visual Summaries

Embed
<iframe src="https://embed.picostat.com/r-dataset-package-histdata-snowstreets.html" frameBorder="0" width="100%" height="307px" />
Attachment Size
31.99 KB
Documentation

## John Snow's Map and Data on the 1854 London Cholera Outbreak

### Description

The Snow data consists of the relevant 1854 London streets, the location of 578 deaths from cholera, and the position of 13 water pumps (wells) that can be used to re-create John Snow's map showing deaths from cholera in the area surrounding Broad Street, London in the 1854 outbreak. Another data frame provides boundaries of a tesselation of the map into Thiessen (Voronoi) regions which include all cholera deaths nearer to a given pump than to any other.

The apocryphal story of the significance of Snow's map is that, by closing the Broad Street pump (by removing its handle), Dr. Snow stopped the epidemic, and demonstrated that cholera is a water borne disease. The method of contagion of cholera was not previously understood. Snow's map is the most famous and classical example in the field of medical cartography, even if it didn't happen exactly this way. (the apocryphal part is that the epidemic ended when the pump handle was removed.) At any rate, the map, together with various statistical annotations, is compelling because it points to the Broad Street pump as the source of the outbreak.

### Usage

	data(Snow.deaths)
data(Snow.pumps)
data(Snow.streets)
data(Snow.polygons)
data(Snow.dates)


### Format

Snow.deaths: A data frame with 578 observations on the following 3 variables, giving the address of a person who died from cholera. When many points are associated with a single street address, they are "stacked" in a line away from the street so that they are more easily visualized. This is how they are displayed on John Snow's original map. The dates of the deaths are not individually recorded in this data set.

case

Sequential case number, in some arbitrary, randomized order

x

x coordinate

y

y coordinate

Snow.pumps: A data frame with 13 observations on the following 4 variables, giving the locations of water pumps within the boundaries of the map.

pump

pump number

label

pump label: Briddle St Broad St ... Warwick

x

x coordinate

y

y coordinate

Snow.streets: A data frame with 1241 observations on the following 4 variables, giving coordinates used to draw the 528 street segment lines within the boundaries of the map. The map is created by drawing lines connecting the n points in each street segment.

street

street segment number: 1:528

n

number of points in this street line segment

x

x coordinate

y

y coordinate

Snow.polygons: A list of 13 data frames, giving the vertices of Thiessen (Voronoi) polygons containing each pump. Their boundaries define the area that is closest to each pump relative to all other pumps. They are mathematically defined by the perpendicular bisectors of the lines between all pumps. Each data frame contains:

x

x coordinate

y

y coordinate

Snow.deaths2: An alternative version of Snow.deaths correcting some possible duplicate and missing cases, as described in vignette("Snow_deaths-duplicates").

Snow.dates: A data frame of 44 observations and 3 variables from Table 1 of Snow (1855), giving the number of fatal attacks and number of deaths by date from Aug. 19 – Sept. 30, 1854. There are a total of 616 deaths represented in both columns attacks and deaths; of these, the date of the attack is unknown for 45 cases.

### Details

The scale of the source map is approx. 1:2000. The (x, y) coordinate units are 100 meters, with an arbitrary origin.

Of the data in the Snow.dates table, Snow says, “The deaths in the above table are compiled from the sources mentioned above in describing the map; but some deaths which were omitted from the map on account of the number of the house not being known, are included in the table.”

One limitation of these data sets is the lack of exact street addresses. Another is the lack of any data that would serve as a population denominator to allow for a comparison of mortality rates in the Broad Street pump area as opposed to others. See Koch (2000), Koch (2004), Koch \& Denike (2009) and Tufte (1999), p. 27-37, for further discussion.

### Source

Tobler, W. (1994). Snow's Cholera Map, http://www.ncgia.ucsb.edu/pubs/snow/snow.html; data files were obtained from http://ncgia.ucsb.edu/Publications/Software/cholera/, but these sites seem to be down.

The data in these files were first digitized in 1992 by Rusty Dodson of the NCGIA, Santa Barbara, from the map included in the book by John Snow: "Snow on Cholera...", London, Oxford University Press, 1936.

### References

Koch, T. (2000). Cartographies of Disease: Maps, Mapping, and Medicine. ESRI Press. ISBN: 9781589481206.

Koch, T. (2004). The Map as Intent: Variations on the Theme of John Snow Cartographica, 39 (4), 1-14.

Koch, T. and Denike, K. (2009). Crediting his critics' concerns: Remaking John Snow's map of Broad Street cholera, 1854. Social Science \& Medicine 69, 1246-1251.

Snow, J. (1885). On the Mode of Communication of Cholera. London: John Churchill. http://www.ph.ucla.edu/epi/snow/snowbook.html.

Tufte, E. (1997). Visual Explanations. Cheshire, CT: Graphics Press.

SnowMap

### Examples

data(Snow.deaths)
data(Snow.pumps)
data(Snow.streets)
data(Snow.polygons)
data(Snow.deaths)## Plot deaths over time
require(lubridate)
clr <- ifelse(Snow.dates$date < mdy("09/08/1854"), "red", "darkgreen") plot(deaths ~ date, data=Snow.dates, type="h", lwd=2, col=clr) points(deaths ~ date, data=Snow.dates, cex=0.5, pch=16, col=clr) text( mdy("09/08/1854"), 40, "Pump handle\nremoved Sept. 8", pos=4) ## draw Snow's map and dataSnowMap()# add polygons SnowMap(polygons=TRUE, main="Snow's Cholera Map with Pump Polygons")# zoom in a bit, and show density estimate SnowMap(xlim=c(7.5,16.5), ylim=c(7,16), polygons=TRUE, density=TRUE, main="Snow's Cholera Map, Annotated") ## re-do this the sp way... [thx: Stephane Dray]library(sp)# streets slist <- split(Snow.streets[,c("x","y")],as.factor(Snow.streets[,"street"])) Ll1 <- lapply(slist,Line) Lsl1 <- Lines(Ll1,"Street") Snow.streets.sp <- SpatialLines(list(Lsl1)) plot(Snow.streets.sp, col="gray") title(main="Snow's Cholera Map of London (sp)")# deaths Snow.deaths.sp = SpatialPoints(Snow.deaths[,c("x","y")]) plot(Snow.deaths.sp, add=TRUE, col ='red', pch=15, cex=0.6)# pumps spp <- SpatialPoints(Snow.pumps[,c("x","y")]) Snow.pumps.sp <- SpatialPointsDataFrame(spp,Snow.pumps[,c("x","y")]) plot(Snow.pumps.sp, add=TRUE, col='blue', pch=17, cex=1.5) text(Snow.pumps[,c("x","y")], labels=Snow.pumps$label, pos=1, cex=0.8)
--

Dataset imported from https://www.r-project.org.

Picostat Manual
###### How To Register With a Username
1. Go to the user registration page.
4. Click Submit.
5. Click the link that was sent to the email address you registered with.
6. Clicking the link will open another page on Picostat where you can select a password.
7. Click Save and enter any profile details you wish to enter.
###### How To Register With Google Single Sign On (SSO)
1. Go to the user login page.
5. Google will redirect you back to Picostat with your new account created and you will be logged in.
6. Enter any profile details you wish to share.
1. Go to the user login page.
3. Click "Login". You will be redirected to your user homepage authenticated.
1. Go to the user login page.
3. If you already registered with Picostat via Google SSO, you will be redirected to your user homepage authenticated.
###### How To Import a Dataset
1. Create a Picostat account or login with your existing picostat account (see above).
2. Go to the dataset import page.
3. Select a license for the dataset. The default is "No License" but allows Picostat to host a copy of the dataset as per the privacy policy. You may wish to uncheck the "Public" option if you do not wish to share your dataset with others. R Datasets that come by downloading R have a GNU General Public License v3.0 which may also be selected from the Picostat dropdown.
4. Enter a title for the dataset
5. Choose a dataset input methods. Available options include:
• Random data - this populates your dataset with random numbers between 0 and 100. You can specify the number of rows and columns for the random dataset.
• CSV, TSV or TXT file - you will have the option upload a file within the current file size limit and also specify the header and whether or not the dataset is a contingency table. With contingency tables, the first column becomes a label for the rows. Currently with Picostat, there is limited support for contingency tables. Choose "Yes" to the Header option if the first line of the data contains titles for the rows. Also choose the Separator for the dataset. A separator is what breaks the data up. In some cases, a comma would separate data values in a row. You will also have the option to add documentation in the form keyboarded text and also uploaded documentation attachments. You can also specify a license for the documentation.
• Copy and Paste. This selection contains many of the same fields as importing a file with an additional textarea to copy and paste data to.
• Empty dataset. Start with a blank dataset and manually add data with the Picostat dataset editor.
• Excel file - Choose this option if you would like to convert your Excel spreadsheet to a Picostat dataset. With this selection, you will have the option to specify whether to use the first row in the Excel file as column names. If you would like to choose a specific sheet to use, you can also specify with entering its name in the text input.
• sas7bdat file - SAS is a powerful statistical software package that has its own proprietary file format. Choose this option if you are importing a SAS file.
• SPSS sav file - SPSS is a statistical package owned by IBM. You can import SPSS files by choosing this option.
6. Choose whether or not the dataset contains a header. Some of the dataset input methods allow you to specify whether or not a Header exists on the file. Sometimes dataset files contain a Header as the first row which names the columns. If you choose "Yes" to this, the first row in the dataset will become column headers.
7. You can also add documentation and specify a documentation license. This can be used to help explain your dataset to those unfamiliar with it.
8. Choose whether or not to upload an supporting attachments.
9. Pass the captcha. To prevent spam submissions, Picostat has a captcha which is used to prevent automated submissions by bots.
10. Choose a privacy setting for the dataset. You can also specify whether or not the dataset is Public. If you uncheck this setting, only you and the Picostat administrator will be able to view the dataset.
11. Submit the form. Once the form is validated, you will be redirected to the dataset homepage where you can choose to edit or perform statistical operations on the dataset.
###### How To Perform Statistical Analysis with Picostat
1. Go to any dataset homepage. You can get a full list at the dashboard.
2. Near the top of the page there will be two drop downs. One for analysis and one for education. Here we will choose Analyis. Choose from one of the following:
• Numerical Summaries - Here you can get the:
1. Arithmetic mean
2. Median
3. Quartiles
4. Minimum and Maximum
5. Stem-and-leaf plot
6. Standard deviation and Variance
7. IQR
8. Cumulative frequencies
• Plot - a plot of two columns on the cartesian coordinate system
• Boxplot - a Boxplot (box-and-whisker plot) of a column.
• Correlation Coefficient - Compute the correlation coefficient between two columns.
• Cumulative Frequency Histogram - Display a cumulative frequency histogram
• Dotplot
• Hollow Histogram - Plot two columns on the same histogram with a different color for each column.
• Pie Chart
• Regression - Perform a simple linear regression and compute the p-value and regression line. Also plots the data with the regression line.
• Stem and Leaf Plots - Plot a one or two-sided stem-and-leaf plot from one or two columns respectively.
• Visual Summaries - plots the following:
1. Frequency Histogram
2. Relative Frequency Histogram
3. Cumulative Frequency Histogram
4. Boxplot (Box-and-whisker plot)
5. Dotplot
3. PDF - Check this box if you want to download a PDF of the output from the statistical analysis application.
4. HTML - Check this box if you want a link to the HTML page created with the application.
5. R File - Check this box if you want a link to download the R commands used to generate output from this statistical application. You should be able to copy the commands verbatim into R Studio or R to recreate the analysis.
###### How To Use Educational Applications with Picostat
1. Go to any dataset homepage. You can get a full list at the dashboard.
2. Near the top of the page there will be two drop downs. One for analysis and one for education. Here we will choose Education. Choose from one of the following:
• How To Create a Barplot - This will show you how to create a bar chart after selecting a column with the mouse.
• How To Create a Stacked Barplot - This application will show you how to create a stacked bar plot from a column vector.
• How To Create a Pie Chart - This application will show you how to create a pie chart from a column of data
• How To Compute the Mean - This application will show you how to compute the mean from a column vector
• How To Create a Plot - This app will show you how to plot two columns in the cartesian coordinate system.
• How To Compute the Media - This statisistical app will show you how to compute the median from a column vector.
3. PDF - Check this box if you want to download a PDF of the output from the education application.
4. HTML - Check this box if you want a link to the HTML page created with the application.
Recent Queries For This Dataset

No queries made on this dataset yet.

Title Authored on Content type
OpenIntro Statistics Dataset - husbands_wives August 9, 2020 - 12:25 PM Dataset
R Dataset / Package wooldridge / twoyear March 9, 2018 - 1:06 PM Dataset
R Dataset / Package Stat2Data / WordMemory March 9, 2018 - 1:06 PM Dataset
OpenIntro Statistics Dataset - ppp_201503 August 9, 2020 - 2:30 PM Dataset
R Dataset / Package car / Greene March 9, 2018 - 1:06 PM Dataset