data(fipslookup)
data(dplookup)
## Not run:
#
# censusData <- readCensus(state="Kentucky",fips=fipslookup,codes=dplookup)
# seerData <- readSEER(year=2005,state="Kentucky",cancer="Larynx",fips=fipslookup)
#
# library(diseasemapping)
# larynxRates <- cancerRates("USA", year=1998:2002,site="Larynx")
# kentucky <- getSMR(censusData,larynxRates,seerData,regoinCode="County")
#
# breaks = c(0,1,seq(2, ceiling(max(kentucky$SMR)),by=2))
# thecol = terrain.colors(length(breaks)-1)
# plot(kentucky, col = thecol[cut(kentucky$SMR, breaks,include.lowest=TRUE)] )
# ## End(Not run)
Run the code above in your browser using DataLab