Hi All,
I've got another question for the group. I'm using the bubble function
in GSTAT to plot the residuals from a krige.cv object. I'd like to
place on top the outlines of my counties of interest (or even better
I'd like to plot them in the reverse order). What I'm finding is that
even if I eliminate the key and eliminate the "asp" call from the
bubble function, my two plots don't align as they do under more
traditional plotting circumstances
Here's the idea:
kcv28<-krige.cv(avg.pm25~1,pm.all28, vFitRaw, nmax=dim(pm.all28)[1],
nfold=dim(pm.all28)[1])
bubble(kcv28, "residual", do.sqrt=F)
map.nyc<-read.shape("D:/junk/counties/all_counties.shp", dbf.data="">
plot(map.nyc, add=T, fg="transparent", ol="gray50")
Any suggestions? Zev
--
Untitled Document
Zev
Ross
ZevRoss Spatial Analysis
303 Fairmount Ave
Ithaca, NY 14850
(607) 277-0004 (phone)
(866) 877-3690 (fax toll-free)
zev@xxxxxxxxxxx
www.zevross.com
|
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@xxxxxxxxxxxxxxxxx
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
| 