Displaying 1 result from an estimated 1 matches for "kmeansout_150".
2004 May 03
1
Speed up graphics output?
...------------------------
#Input file is comma-delimited and has GeogID, LON, LAT and CLUSTER
#LON should be in column 2 and LAT in column 3
library(maptools)
#Number of clusters in total
ClusCount <- 150
#Read in Cluster Assignments
datain <- as.data.frame( read.table("c:\\data\\Run1\\kmeansout_150.txt", header=TRUE ))
colnames(datain) <- c("geogid","long","lat","cluster")
#Set up screen device
#split 2x2 - permits 1 map per section
par(bg="white")
split.screen(c(2,2))
#Load maps once only
w <- read.shape("c:\\data\\region1...