[email protected] wrote: > Dear all, > > > I want to plot the cross hscattergram between two different variables. > > I used the following instruction: > > hscat(jura.pred$Cd~jura.pred$Zn, jura.pred, c(0, 0.1, 0.2, 0.3, 0.5, 1), > alpha=c(0,45)) > > where jura.pred$Cd and jura.pred$Ni are the two variables, but the result > obtained is the same of the univariate case when I put 1 instead of jura. > pred$Zn. > > Can you help me? > This cannot be done with hscat (package gstat), right now. It doesn't say so explicitly, but from the description:
formula: specifies the dependent variable you could have understood that a single dependent variable is considered. It doesn't seem terribly hard to program, either from scratch or building upon the sources in package gstat. -- Edzer > Thanks > > Cristiano > > >> ----Messaggio originale---- >> Da: [email protected] >> Data: 03/04/2009 21.40 >> A: "milton ruser"<[email protected]> >> Cc: <[email protected]> >> Ogg: Re: [R-sig-Geo] rasterize shape file and output as .img >> >> On Fri, 3 Apr 2009, milton ruser wrote: >> >> >>> Dear all, >>> >>> I have a set of shapefile and I need rasterized it using a collumn called >>> ClassCover, >>> and with the resolution of 30 meters. How can I do this on R? >>> >> Use bbox() of the SpatialPolygonDataFrame to find the extent. Create a >> GridTopology object, and overlay the polygons - possibly just with >> ClassCover - and the SpatialGrid built with the GridTopology object. This >> should yield a SpatialGridDataFrame with a single band, with NAs outside >> the polygons, and the polygon column values for raster cells where their >> centres fall inside the polygons. >> >> >>> Another thing is that I need to output the raster as IMG (erdas/arcgis) >>> format. >>> >> writeGDAL in rgdal. >> >> Hope this helps, >> >> Roger >> >> >>> Any help are welcome. >>> >>> Best wishes, >>> >>> milton >>> Brasil/Toronto >>> >>> >>>> require(maptools) >>>> require(sp) >>>> require(rgdal) >>>> todos<-readShapePoly("myshape.shp") >>>> sessionInfo() >>>> >>> R version 2.9.0 alpha (2009-03-26 r48224) >>> i386-pc-mingw32 >>> locale: >>> LC_COLLATE=English_Canada.1252;LC_CTYPE=English_Canada.1252; >>> > LC_MONETARY=English_Canada.1252;LC_NUMERIC=C;LC_TIME=English_Canada.1252 > >>> attached base packages: >>> [1] stats graphics grDevices utils datasets methods base >>> other attached packages: >>> [1] rgdal_0.6-7 maptools_0.7-20 sp_0.9-32 foreign_0.8-34 >>> loaded via a namespace (and not attached): >>> [1] grid_2.9.0 lattice_0.17-20 >>> >>> [[alternative HTML version deleted]] >>> >>> _______________________________________________ >>> R-sig-Geo mailing list >>> [email protected] >>> 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: [email protected] >> >> _______________________________________________ >> R-sig-Geo mailing list >> [email protected] >> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >> >> > > _______________________________________________ > R-sig-Geo mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de/ http://www.springer.com/978-0-387-78170-9 [email protected] _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
