On Thu, 26 Aug 2010, Dorothea Lemke wrote:

Dear list,

I've a problem with plotting a numeric object in a SpatialPixelsDataFrame.
I've done 2 Kerneldensity estimations with the package splancs and I've
calculated the relative risk by taking the ratio of this two intensities. When
I use the spplot() function I get one window with 3 graphics on the same
scale. But I want to plot the graphics with seperate scales. How can I do
this?

Save the three output objects of spplot, then see:

library(lattice)
?plot.trellis

especially the more= and split= arguments.

Hope this helps,

Roger


My Rcode is:

 > kcasesC67m<-spkernel2d(cases, pRBMS, h0=3500, gt)
 > kcontrolsC67m<-spkernel2d (controls, pRBMS, h0=3500, gt)
 > df0<-data.frame (kcasesC67m = kcasesC67m, kcontrolsC67m = kcontrolsC67m)
 > spkratio0<-SpatialGridDataFrame (gt,  data = df0)
 > spkratio<-as(spkratio0, "SpatialPixelsDataFrame")
 > spkratio$kratio<-spkratio$kcasesC67m/spkratio$kcontrolsC67m
 > is.na(spkratio$kratio)<-!is.finite(spkratio$kratio)

Many thanks for your help

Dorothea

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo


--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: roger.biv...@nhh.no

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to