Re: [R-sig-Geo] how is variance defined in krige {gstats}

Hi all, from what I could devine (with the help of collegues) out of the gstat documentation (should have looked at www.gstat.org earlier...) the kriging variance (var1.var) is the kriging standard error ^2. I am still trying to get Chiles (1999) as it will be able to enlighten me on other things

[R-sig-Geo] Help with the Kmeasure function from the Spatstat package

Debarchana Ghosh writes: >I'm doing exploratory point pattern data analysis with 285 points in a >irregular polygon. I'm really exploring and hence trying out all possible >clustering functions on a single data set. In the 'Kmeasure' function, one >of the arguments that the user needs to be suppl

Re: [R-sig-Geo] Can you krige with a nominal response variable

I might consider using gam() in R to smooth the 0/1 indicator for each class individually based on a Bernoulli model - doing each class individually. A downside of course is that the resulting smoothed maps of the proportions for each class would not add to one at an individual location, but th

[R-sig-Geo] local K function

Hello all, is the local K function available in any R package? Thanks, Colin Colin Robertson MSc Candidate Dept of Geography University of Victoria www.geog.uvic.ca/spar ph. (250)885-8065 [[alternative HTML version deleted]] _

Re: [R-sig-Geo] map to shapefile with projection info

On Wed, 6 Dec 2006, Sebastian P. Luque wrote: > On Wed, 06 Dec 2006 14:00:38 -0600, > "Sebastian P. Luque" <[EMAIL PROTECTED]> wrote: > > > On Wed, 6 Dec 2006 20:28:42 +0100 (CET), > > Roger Bivand <[EMAIL PROTECTED]> wrote: > > [...] > > >> names <- sapply(slot(canada_sp, "polygons"), function

Re: [R-sig-Geo] map to shapefile with projection info

On Wed, 06 Dec 2006 14:00:38 -0600, "Sebastian P. Luque" <[EMAIL PROTECTED]> wrote: > On Wed, 6 Dec 2006 20:28:42 +0100 (CET), > Roger Bivand <[EMAIL PROTECTED]> wrote: [...] >> names <- sapply(slot(canada_sp, "polygons"), function(i) slot(i, "ID")) >> canada_spdf <- SpatialPolygonsDataFrame(can

Re: [R-sig-Geo] map to shapefile with projection info

On Wed, 6 Dec 2006 20:28:42 +0100 (CET), Roger Bivand <[EMAIL PROTECTED]> wrote: [...] > canada <- map("world2Hires", "Canada", fill=TRUE, plot=FALSE) > fill=TRUE retrives polygons, otherwise you get line segments. For this > example it now works. You can also make several retrieved island belon

[R-sig-Geo] Can you krige with a nominal response variable

Is there a method for kriging a response with multiple states, for instance vegetation class in 10 classes? I know indicator kriging would work for a boolean variable. Is there an extension for multiple responses? Kenneth B. Pierce Jr. Research Ecologist Landscape Ecology, Modeling, Mapping and

Re: [R-sig-Geo] map to shapefile with projection info

On Wed, 6 Dec 2006, Sebastian P. Luque wrote: > Hi, > > Trying to extract a map (package maps) to a shapefile, IIUC the procedure > is to get it into a SpatialPolygonsDataFrame and then use > writePolyShape(). This seems to require an intermediate step through > SpatialPolygons. However, I'm ru

[R-sig-Geo] map to shapefile with projection info

Hi, Trying to extract a map (package maps) to a shapefile, IIUC the procedure is to get it into a SpatialPolygonsDataFrame and then use writePolyShape(). This seems to require an intermediate step through SpatialPolygons. However, I'm running into some problems: R> canada <- map("world2Hires",

[R-sig-Geo] Live e-seminar

Edzer Pebesma's e-seminar is just beginning (1700 UTC) on: www.giscience.info If you miss it, you should be able to access a recording from the same site or from www.wun.ac.uk/ggisa/seminars.html later. -- Roger Bivand Economic Geography Section, Department of Economics, Norwegian School of Ec

Re: [R-sig-Geo] project the lat./long. to x/y

On Wed, 6 Dec 2006, zhijie zhang wrote: > Dear friends, > I had a point pattern dataset, whose locations was defined by lat./long.. > In the past, I used ARCGIS to specify its coordinate system and the > projection method, now i want to do it in R. I don't know how to do it, > because it needs

[R-sig-Geo] project the lat./long. to x/y

Dear friends, I had a point pattern dataset, whose locations was defined by lat./long.. In the past, I used ARCGIS to specify its coordinate system and the projection method, now i want to do it in R. I don't know how to do it, because it needs me to specify the *Geographic Coordinate System and

Re: [R-sig-Geo] krig.image beyond data points

Dear Micheal, I think you can also use the spplot command. It is not present in the gstat package but in the sp package. The most recent version of gstat (0.9-34) needs sp present, so I think you already have this package. For more info on spplot use: ?spplot I think that the ylim and xlim o

Re: [R-sig-Geo] Help with the Kmeasure function from the Spatstat package

Dear Debarchana, I think that you could check any of the good books about point patterns that there are currently out there. Diggle (2003), for example, describes what you ask and most of the examples in the book can be reproduced using package Splancs. Check manual page of function mse2d and ref

Re: [R-sig-Geo] how is variance defined in krige {gstats}

Dear Herry, Details on kriging variance are given in several books on geostatistics. For example the following book: (Chiles1999) Chiles, J. & Delfiner, P. Geostatistics: Modeling Spatial Uncertainty John Wiley & Sons, New York, 1999 A paragraph on the kriging variance is on page 176-177. How