[R] help with large gamm

Hi All, I was wondering if anyone has had any experience running gamm(mgcv) with a large spatial data set? I know the gamm help says this is a bad idea, but I wondered if there was a way to use something like filehash to store the matrices that gamm creates? -- David Depew PhD Candidate Depar

[R] gls and memory allocation

Hi all, A two part quick question regarding gls 1) I'd like to pass a formula to gnls from a gam(mgcv package). Is there a quick way to do this? I tried using >gnls(gam$call$formula, data=df, correlation=AC1) but it keeps outputting a message to the effect that I need x ~ y etc as the formula

Re: [R] spatialpixelsdataframe

Thanks Paul, Will keep that in mind for my next query :-) Paul Hiemstra wrote: > Dave Depew wrote: >> I have a spatialpixelsdataframe containing bathymetry data. Is it >> possible to save this as something similar to the meuse.grid so it >> can be loaded without havi

[R] spatialpixelsdataframe

I have a spatialpixelsdataframe containing bathymetry data. Is it possible to save this as something similar to the meuse.grid so it can be loaded without having to run the spatial interpolation each time the script is called? I don't see a straightforward method in the sp documentation. _

[R] memory issues

Hi all, I've read the R for windows FAQ and am a little confused re: memory.limit and memory.size to start using R 2.6.2 on WinXP, 2GB RAM, I have the command line "sdi --max-mem-size=2047M" Once the Rgui is open, memory.limit() returns 2047, memory.size() returns 11.315, and memory.size(max=T)

[R] memory issues

Hi all, I've read the R for windows FAQ and am a little confused re: memory.limit and memory.size to start using R 2.6.2 on WinXP, 2GB RAM, I have the command line "sdi --max-mem-size=2047M" Once the Rgui is open, memory.limit() returns 2047, memory.size() returns 11.315, and memory.size(max=T)

[R] minimum distances

Hi, I have a question about computing shortest Euclidean distances between two data frames of spatial points... I have 2 dataframes (not yet converted to spatial class) >Sewers<-data.frame(x=c(10,20,35,50),y=c(100,150,200,300)) >transect <- data.frame(x=seq(from=0, to=50, by=1),y=seq(from=100