Working in the mines got me interested in this especially since we just had a fatality at one of the mines nearby (not the one I work in). Their investigating and rumors are everywhere currently. It's different from my other research but I think it's time to branch out. To analyze the data I am using a survival analysis on mining accidents in the United States. So first we need to scrape some data from the Federal Mine Safety and Health Administration. We will scrape it in R. The code is below:

​

install.packages("rvest")
install.packages("dplyr")

library(rvest)
library(dplyr)

html_data0<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=0")
html_data1<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=1")
html_data2<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=2")
html_data3<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=3")
html_data4<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=4")
html_data5<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=5")
html_data6<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=6")
html_data7<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=7")
html_data8<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=8")
html_data9<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=9")
html_data10<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=10")
html_data11<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=11")
html_data12<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=12")
html_data13<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=13")
html_data14<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=14")
html_data15<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=15")
html_data16<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=16")
html_data17<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=17")
html_data18<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=18")
html_data19<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=19")
html_data20<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=20")
html_data21<-read_html("https://www.msha.gov/data-and-reports/fatality-reports/search?page=21")

data0<-html_data0 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data1<-html_data1 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data2<-html_data2 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data3<-html_data3 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data4<-html_data4 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data5<-html_data5 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data6<-html_data6 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data7<-html_data7 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data8<-html_data8 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data9<-html_data9 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data10<-html_data10 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data11<-html_data11 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data12<-html_data12 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data13<-html_data13 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data14<-html_data14 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data15<-html_data15 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data16<-html_data16 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data17<-html_data17 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data18<-html_data18 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data19<-html_data19 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data20<-html_data20 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()
data21<-html_data21 |>
html_elements(".views-field-field-primary-sic") |>
html_text2()

​

Then to print the data to make sure you grabbed it just type:

print(data0)

​

You should get this:

​[1] "Mined Material: Lime" [2] "Mined Material: Copper Ore NEC" [3] "Mined Material: Coal (Bituminous)" [4] "Mined Material: Crushed Broken Stone NEC" [5] "Mined Material: Construction Sand and Gravel" [6] "Mined Material: Coal (Bituminous)" [7] "Mined Material: Coal (Bituminous)" [8] "Mined Material: Crushed Broken Stone NEC" [9] "Mined Material: Coal (Bituminous)" [10] "Mined Material: Coal (Bituminous)" [11] "Mined Material: Sand industrial NEC" [12] "Mined Material: Copper Ore NEC" [13] "Mined Material: Coal (Bituminous)" [14] "Mined Material: Sand industrial NEC" [15] "Mined Material: Dimension Marble" [16] "Mined Material: Cement" [17] "Mined Material: Cement" [18] "Mined Material: Crushed, Broken Limestone NEC" [19] "Mined Material: Coal (Bituminous)" [20] "Mined Material: Crushed Broken Limestone NEC" [21] "Mined Material: Crushed Broken Limestone NEC" [22] "Mined Material: Salt" [23] "Mined Material: Coal (Bituminous)" [24] "Mined Material: Coal (Bituminous)" [25] "Mined Material: Crushed Broken Limestone NEC" [26] "Mined Material: Lime" [27] "Mined Material: Gold Ore" [28] "Mined Material: Crushed Broken Limestone NEC" [29] "Mined Material: Construction Sand and Gravel" [30] "Mined Material: Crushed Broken Limestone NEC"

​