Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS

2006-05-15 Thread Thomas Adams
Roger, That did it! The sessionInfo() command identified that I had older versions than what were available on CRAN, namely: lattice -- 0.13-8 gstat -- 0.9-29 (was older: 0.9-17) sp -- 0.8-14 (was older: 0.8-8) On my Mac, I have: lattice -- 0.11-8 gstat -- 0.9-21 sp -- 0.7-10 which works… BT

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS

2006-05-13 Thread Roger Bivand
Thomas: Please do: sessionInfo() on both machines, and see if any of the package versions are different. It often also helps to name function arguments explicitly, to be sure that they are being understood correctly inside the functions. Depending on the class of the data= arugment, as far as

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-05-01 Thread Thomas Adams
Roger, This is what the summary() command gives me before/after I use fullgrid(dem) <- FALSE dem<-readFLOAT6sp("ohrfc.dem") > summary(dem) Object of class SpatialGridDataFrame Coordinates: min max coords.x1 -351000 675000 coords.x2 -486000 486000 Is projected: TRUE proj4string : [+proj=lcc +la

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-04-28 Thread Roger Bivand
On Fri, 28 Apr 2006, Thomas Adams wrote: > Roger, > > Your suggestion: > > fullgrid(dem) <- FALSE > > did turn dem into class type SpatialGridDataFrame, but when I tried: > > z <- predict(UK_fit,newdata=dem) > > I got an error: > > Error in model.frame(... : > invalid variable type. > > I t

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-04-28 Thread Thomas Adams
Roger, Your suggestion: fullgrid(dem) <- FALSE did turn dem into class type SpatialGridDataFrame, but when I tried: z <- predict(UK_fit,newdata=dem) I got an error: Error in model.frame(... : invalid variable type. I think I should restate the problem: I have a file 'temps' which has class

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-04-28 Thread Roger Bivand
On Fri, 28 Apr 2006, Thomas Adams wrote: > Roger, > > This got me further along, but I am encountering a problem with: > > z <- predict(UK_fit, newdata=BMcD_SPx) > > The gstat step works for me, where I have: > > UK_fit<-gstat(formula=temps$temp~dem,data=temps,model=efitted) > > temps has cla

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-04-28 Thread Thomas Adams
Roger, This got me further along, but I am encountering a problem with: z <- predict(UK_fit, newdata=BMcD_SPx) The gstat step works for me, where I have: UK_fit<-gstat(formula=temps$temp~dem,data=temps,model=efitted) temps has class SpatialPointsDataFrame: coordinates cat x

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-04-27 Thread Thomas . Adams
age - From: Roger Bivand <[EMAIL PROTECTED]> Date: Thursday, April 27, 2006 4:21 pm Subject: Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed > On Thu, 27 Apr 2006, Thomas Adams wrote: > > > List: > > > > I can not seem to work out the

Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-04-27 Thread Roger Bivand
On Thu, 27 Apr 2006, Thomas Adams wrote: > List: > > I can not seem to work out the syntax for using R/gstat within a GRASS > 6.1 session to do universal kriging. I have a DEM (elevation data on a > grid) and point data for temperature; theoretically, the temperatures > should relate to elevat

[R-sig-Geo] Example of universal kriging with R/gstat in GRASS needed

2006-04-27 Thread Thomas Adams
List: I can not seem to work out the syntax for using R/gstat within a GRASS 6.1 session to do universal kriging. I have a DEM (elevation data on a grid) and point data for temperature; theoretically, the temperatures should relate to elevation. So, I am trying to spatially interpolate the tem