[R-sig-Geo] Plotting SpatialPixelDataFrame

[Roger Bivand]

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 at 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.Bivand at nhh.no