r - Rで別のベクトルに基づいてデータをマップする方法

Question

次のようなデータセット (pts) があります。

x <- seq(-124.25,length=115,by=0.5)    
y <- seq(26.25,length=46,by=0.5)
z <- 1:5290

longlat <- expand.grid(x = x, y = y)
pts <- data.frame(longlat,z) 
names(pts) <- c( "lon","lat","data")

次のようにして、データフレーム（pts）をマップにマップできることを知っていました。

library(sp)
library(rgdal)
library(raster)
library(maps)
coordinates(pts)=~x+y

proj4string(pts)=CRS("+init=epsg:4326") # set it to long, lat

pts = spTransform(pts,CRS(" +init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0"))
pts <- as(pts, "SpatialPixelsDataFrame")
r = raster(pts)
projection(r) = CRS(" +init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0")

plot(r)
map("usa",add=T)

もう1つの要素「クラス」がある場合

pts$class=rep(1:5,1058)

同じマップ上に 5 セットのモノカラー スキームを配置できるように、「クラス」の情報を追加するにはどうすればよいですか?