
Data set from NOAA's National Weather Service Indianapolis, IN Weather Forecast Office 6900 W. Hanna Ave.
data(tornados)
torn Tornado Touchdown points SpatialPointsDataFrame
- geographical projection
torn2 Tornado Touchdown points SpatialPointsDataFrame
- equal area projection
# Read in the data
data(tornados)
# Split the plot in two
opar <- par(mfrow=c(2,1))
# Plot US states
plot(us_states)
# Add Locations of observation stations
plot(torn,add=TRUE,pch=16,col='red')
# Plot a histogram of year of observation next to this
hist(torn$YEAR)
par(opar)
Run the code above in your browser using DataLab