PEER ASSIGNMENT TWO

Loading Data…

The dataset shows events from January 1950 to November 2011.

getwd()

## [1] "C:/Users/user1/Downloads/Reproducible Research/Peer Assignment 2"

storm_data <- read.csv("./repdata-data-StormData.csv")
mycol <- c("EVTYPE", "FATALITIES", "INJURIES", "PROPDMG", "PROPDMGEXP", "CROPDMG", 
    "CROPDMGEXP")
mydata <- storm_data[mycol]

fatal <- aggregate(FATALITIES ~ EVTYPE, data = mydata, FUN = sum)
injury <- aggregate(INJURIES ~ EVTYPE, data = mydata, FUN = sum)
propdmg <- aggregate(PROPDMGEXP ~ EVTYPE, data = mydata, FUN = sum)

## Error in Summary.factor(structure(17L, .Label = c("", "-", "?", "+", "0", : sum not meaningful for factors

cropdmg <- aggregate(CROPDMGEXP ~ EVTYPE, data = mydata, FUN = sum)

## Error in Summary.factor(structure(1L, .Label = c("", "?", "0", "2", "B", : sum not meaningful for factors

Most harmful with respect to population health

fatal10 <- fatal[order(-fatal$FATALITIES), ][1:10, ]
# get top10 event with highest injuries
injury10 <- injury[order(-injury$INJURIES), ][1:10, ]
par(mfrow = c(1, 2), mar = c(12, 4, 3, 2), mgp = c(3, 1, 0), cex = 0.8)
barplot(fatal10$FATALITIES, las = 3, names.arg = fatal10$EVTYPE, main = "Weather Events With The Top 10 Highest Fatalities", 
    ylab = "number of fatalities", col = "red")
barplot(injury10$INJURIES, las = 3, names.arg = injury10$EVTYPE, main = "Weather Events With the Top 10 Highest Injuries", 
    ylab = "number of injuries", col = "red")

plot of chunk unnamed-chunk-3

The plots below show the Total and Per-Event impacts and the estimate of Human Economic Cost by type of event

propdmg10 <- propdmg[order(-propdmg$PROPDMGVAL), ][1:10, ]

## Error in eval(expr, envir, enclos): object 'propdmg' not found

# get top 10 events with highest crop damage
cropdmg10 <- cropdmg[order(-cropdmg$CROPDMGVAL), ][1:10, ]

## Error in eval(expr, envir, enclos): object 'cropdmg' not found

par(mfrow = c(1, 2), mar = c(12, 4, 3, 2), mgp = c(3, 1, 0), cex = 0.8)
barplot(propdmg10$PROPDMGVAL/(10^9), las = 3, names.arg = propdmg10$EVTYPE, 
    main = "Top 10 Events with Greatest Property Damages", ylab = "Cost of damages ($ billions)", 
    col = "red")

## Error in barplot(propdmg10$PROPDMGVAL/(10^9), las = 3, names.arg = propdmg10$EVTYPE, : object 'propdmg10' not found

barplot(cropdmg10$CROPDMGVAL/(10^9), las = 3, names.arg = cropdmg10$EVTYPE, 
    main = "Top 10 Events With Greatest Crop Damages", ylab = "Cost of damages ($ billions)", 
    col = "red")

## Error in barplot(cropdmg10$CROPDMGVAL/(10^9), las = 3, names.arg = cropdmg10$EVTYPE, : object 'cropdmg10' not found