[R-sig-Geo] hscattergram between two different variables

[Edzer Pebesma]

facrimas at libero.it 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: Roger.Bivand at nhh.no
>> Data: 03/04/2009 21.40
>> A: "milton ruser"<milton.ruser at gmail.com>
>> Cc: <r-sig-geo at stat.math.ethz.ch>
>> 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
>>> 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
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>>     
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> 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 e.pebesma at wwu.de