Skip to content
Prev 424 / 29559 Next

kriging inside borders

My exact commands are :
library(maptools)
shapefile=read.shape('my.shp')
x=Map2poly(shapefile)

and the length is :
length(mappolys)
[1] 615
 But i'm not sure that it's the good polylist object i need because it's a polygon grid in fact. 
Thank you for your help
Chlo?
-----Message d'origine-----
De?: Roger Bivand [mailto:Roger.Bivand at nhh.no] 
Envoy??: jeudi 19 mai 2005 19:07
??: Chloe ARCHINARD
Cc?: r-sig-geo at stat.math.ethz.ch
Objet?: Re: [R-sig-Geo] kriging inside borders
On Thu, 19 May 2005, Chloe ARCHINARD wrote:

            
Firstly, it does help if the exact commands are given, rather than your 

choice of them. I think the main issue here is that a polylist object is a 

list of matrices with attributes, not a matrix. If your polylist object 

is:


x <- Map2poly(read.shape("my.shp"))


what is length(x)? If it is 1, you may be able to say x[[1]] to extract 

the first matrix in the list, but even then it may not work. If you can 

show us the commands you are giving, and the length of your polylist 

object, it should be possible to find something that works.
-- 

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



-- 
passerelle antivirus du campus CNRS de Montpellier
--