DATA MINING Desktop Survival Guide by Graham Williams
Desktop Survival Project Home Rattle Introduction Getting Started Interacting with Rattle Data Exploring Data Test Transforming Data Cluster Analysis Association Analysis Classification and Regression Models Evaluation and Deployment Issues Moving into R Troubleshooting Beyond Rattle R Data Graphics in R Understanding Data Preparing Data Descriptive and Predictive Analytics Issues Evaluating Models Reporting Text Mining Text Mining Algorithms Bagging Bayes Classifier Cluster Analysis Conditional Trees Hierarchical Clustering K-Nearest Neighbours Linear Models Neural Networks Support Vector Machines Open Products AlphaMiner Borgelt Data Mining Suite KNime R Rattle Weka Closed Products C4.5 Clementine Equbits Foresight GhostMiner InductionEngine ODM Enterprise Miner Statistica Data Miner TreeNet Virtual Predict Installing Rattle Bibliography Index

Map Data

R provides a set of tools for reading geographic (map) data, particularly ESRI shapefiles, and plotting and manipulating such data. Maps for many countries are available, particularly the US, Europe, and New Zealand. Limited Australian map data is also freely available.

download.file("http://www.vdstech.com/mapdata/australia.zip", "australia.zip")
system("unzip australia.zip; rm australia.zip")

Rsourceget-australia.R

The data can be read in with Rfunction[]readShapePoly and displayed with Rfunction[]plot. Australia has a few outlying islands which we crop from the main focus of the map here using Roption[]xlim and Roption[]ylim.

library(maptools)
aus <- readShapePoly("australia.shp")
plot(aus, lwd=2, border="grey", xlim=c(115,155), ylim=c(-35,-20))
dev.off()

Rsourcemap-australia-plot.R

The class of the resulting object (aus) is SpatialPolygonsDataFrame. Such an object has a collection of slots. For example, the data slot includes meta information about the region recorded in the data frame.

> aus@data

  FIPS_ADMIN GMI_ADMIN                   ADMIN_NAME FIPS_CNTRY CNTRY_NAME
0       AS01   AUS-ACT Australian Capital Territory         AS  Australia
1       AS02   AUS-NSW              New South Wales         AS  Australia
2       AS03   AUS-NTR           Northern Territory         AS  Australia
3       AS04   AUS-QNS                   Queensland         AS  Australia
4       AS05   AUS-SAS              South Australia         AS  Australia
5       AS06   AUS-TSM                     Tasmania         AS  Australia
6       AS07   AUS-VCT                     Victoria         AS  Australia
7       AS08   AUS-WAS            Western Australia         AS  Australia

                 REGION CONTINENT POP_ADMIN  SQKM_ADMIN SQMI_ADMIN  TYPE_ENG
0 Australia/New Zealand Australia    292475    2342.295     904.36 Territory
1 Australia/New Zealand Australia   6338919  803110.812  310081.09     State
2 Australia/New Zealand Australia    161294 1352365.000  522148.09 Territory
3 Australia/New Zealand Australia   3107362 1733475.000  669294.69     State
4 Australia/New Zealand Australia   1445153  985308.500  380427.59     State
5 Australia/New Zealand Australia    472122   68131.477   26305.56     State
6 Australia/New Zealand Australia   4354611  227781.406   87946.40     State
7 Australia/New Zealand Australia   1655588 2533628.000  978233.81     State

   TYPE_LOC
0 Territory
1     State
2 Territory
3     State
4     State
5     State
6     State
7     State

The bounding box of the plot is available using Rfunction[]bbox.

> bbox(aus)
         min       max
r1 112.90721 159.10190
r2 -54.75389 -10.05139

Support further development through the purchase of the PDF version of the book.
PDF version is properly formatted and forms a comprehensive book (draft with over 700 pages).
Brought to you by Togaware. This page generated: Wednesday, 25 February 2009