@[toc]
# Description
This is data derived as shapefiles.
Drought data for 1900-2015.
# Example R code
```{r}
library(rgdal)
flist <- dir("data_derived_shapefiles/", pattern = "shp$")
flist[1:3]
# first 6-mo are empty due to rolling 6-mo integration
flist <- flist[6:length(flist)]
for(fi in flist){
indat <- readOGR(path.expand(file.path("data_derived_shapefiles", fi)))
if(fi == flist[1]){
outdat <- indat@data[1,c(8,15)]
} else {
outdat <- rbind(outdat, indat@data[1,c(8,15)])
}
}
head(outdat)
with(outdat[as.numeric(outdat$count) >= 5 & as.Date(outdat$date) >= as.Date("1980-01-01"),], plot(as.Date(date), count, type = "h"))
```
# Citation
Please cite this project as:
Hanigan IC, Hutchinson MF, Porfirio L. 2013. The Hutchinson Drought Index Database and Algorithm. The Open Science Framework; doi: 10.17605/OSF.IO/PYTS3.