[R-sig-Geo] expand the dist covariate of the prediction-grid meuse.grid (sp package)

2010-05-04 Thread Christoph Hofer
Dear list I would like to expand the covariate dist of the meuse.grid dataset, of the sp package, to the south-east direction. The value of the covariate dist is relative between [0,1]. Now, it is not clear for me how I shall calculate the distance between a new gird point and the river (is

Re: [R-sig-Geo] negative range from fit.variogram through autofitVariogram

2010-05-04 Thread Paul Hiemstra
Hi Mark, Thanks for the reproducible example. The problem is that when I look at the sample variogram, the semivariance values start high and end low. This is best illustrated by: plot(variogram(theta_percent~1, sparse)) You see that there are outliers in the data that cause high

[R-sig-Geo] sp2WB in Maptools - order of polygons

2010-05-04 Thread kath_o_reilly
Dear R users, I plan on using the command sp2WB to convert a .shp file to a .map file in order to carry out some analysis in WinBUGS. Thanks for the threads I've followed on how to do this. The command itself is working fine and produces the map I would expect. However to analyse the data for

Re: [R-sig-Geo] negative range from fit.variogram through autofitVariogram

2010-05-04 Thread Paul Hiemstra
I just uploaded a new version of automap (1.07) that fixes this problem. It deletes fitted variogram models with a negative sill/range/nugget. cheers, Paul Paul Hiemstra wrote: Hi Mark, Thanks for the reproducible example. The problem is that when I look at the sample variogram, the

Re: [R-sig-Geo] negative range from fit.variogram through autofitVariogram

2010-05-04 Thread Mark Connolly
I took the expedient of detecting the negative range and recording the offending time and space slice. Then I went back and checked each model type individually. This is perfectly acceptable but I did want to alert you to the possibility that autofitVariogram and therefore autoKrige could

[R-sig-Geo] enfa [adehabitat]

2010-05-04 Thread Consuelo Hermosilla
Hello guys, I have a couple of questions regarding enfa (adehabitat package). The first is related to this message: In predict.enfa(enfa1, octopus.hab$index, octopus.hab$attr) : the enfa *is not mathematically optimal for prediction*: please consider the madifa instead. I was wondering why

Re: [R-sig-Geo] enfa [adehabitat]

2010-05-04 Thread Clément Calenge
Hello Consuelo, Have a look at: Calenge, C., Darmon, G., Basille, M., Loison, A. and Jullien J.M. (2008) The factorial decomposition of the Mahalanobis distances in habitat selection studies. _Ecology_, *89*, 555-566. See especially the section The MADIFA and the

Re: [R-sig-Geo] enfa [adehabitat]

2010-05-04 Thread Consuelo Hermosilla
I see, thanks for the explanation and the reference Clément. I'll read it! Consuelo - Consuelo Hermosilla PhD student Departamento de Ecología y Biología Animal Departamento de Bioquímica, Genética e Inmunología, Área de Genética Facultad de Ciencias del Mar Campus de As

Re: [R-sig-Geo] negative range from fit.variogram through autofitVariogram

2010-05-04 Thread Mark Connolly
Thanks. I'll pick it up as it becomes available. Does the possibility exist that you would delete all models? In such a case, what does the function return? On 05/04/2010 05:49 AM, Paul Hiemstra wrote: I just uploaded a new version of automap (1.07) that fixes this problem. It deletes

Re: [R-sig-Geo] sp2WB in Maptools - order of polygons

2010-05-04 Thread Virgilio Gomez Rubio
Dear Kath, The command itself is working fine and produces the map I would expect. However to analyse the data for each polygon, I need to know which polygon is which. How do I determine the ID for each polygon generated? I'm hoping there is a better solution than figuring it out by hand...

[R-sig-Geo] how to get the values for corresponding neighbors

2010-05-04 Thread rusers.sh
Hi there, Take the data of columbus(spdep) as an example. Say “i” is one of the polygons, and “j” represents its neighbours. How to get the vlues of “sum(Wij*Xj)”? “Wij” means the weights, which can be get by “nb2listw$weights”. Say “Xj” to be “columbus$PLUMB”. The problem is how to get the

Re: [R-sig-Geo] expand the dist covariate of the prediction-grid meuse.grid (sp package)

2010-05-04 Thread Edzer Pebesma
sp does not do points-to-line distances, but can do point-to-point distances. In the second section I resampled the river boundary so we get a closer approximation, because more points. Maybe rgeos can be expanded do the real point-to-line things (R != GIS?). library(sp) data(meuse.riv) meuse.sr

Re: [R-sig-Geo] how to get the values for corresponding neighbors

2010-05-04 Thread Roger Bivand
On Tue, 4 May 2010, rusers.sh wrote: Hi there, Take the data of columbus(spdep) as an example. Say ?i? is one of the polygons, and ?j? represents its neighbours. How to get the vlues of ?sum(Wij*Xj)?? ?Wij? means the weights, which can be get by ?nb2listw$weights?. Say ?Xj? to be