Re: [R-sig-Geo] Normal score transformation

2008-07-08 Thread Laura Poggio
Thank you very much for your code. I am looking for a method to do normal score transform on raster data, and I will try your code. For the moment I am not interested in the back transformation , but it could be useful for the future. I am sorry if the topic was already considered in the pas, but

Re: [R-sig-Geo] why are idw and gstat different? [SEC=UNCLASSIFIED]

2008-07-08 Thread Edzer Pebesma
[EMAIL PROTECTED] wrote: > Edzer, > > Thanks for identifying the problem. The problem showed how big the impact of > the choice of nmax value on the estimation is. > > In both krige and gstat, the default value for nmax is Inf. I am doing a > simulation experiment now to assess the performance of

Re: [R-sig-Geo] why are idw and gstat different? [SEC=UNCLASSIFIED]

2008-07-08 Thread Jin.Li
Edzer, Thanks for identifying the problem. The problem showed how big the impact of the choice of nmax value on the estimation is. In both krige and gstat, the default value for nmax is Inf. I am doing a simulation experiment now to assess the performance of several spatial interpolators using a

[R-sig-Geo] Japanese translation of "Applied Spatial Data Analysis with R" planned?

2008-07-08 Thread Hisaji ONO
Hello. Some Springer's R related books has been published in Japan. I think some demand for a book about spatial statistics in my country. Any plan or contact of Japanese translation of this interesting book prepared? Regards. ___ R-sig-Geo mai

Re: [R-sig-Geo] Normal score transformation

2008-07-08 Thread Ashton Shortridge
Hi Laura, I've written code to do normal score and back-transforms. Edzer dropped this code from me last month on this list, but I'll repost. Normalizing data is extremely simple in R. The back-transform is the problematic part. Please read the documentation in the file to learn more about the

Re: [R-sig-Geo] why are idw and gstat different? [SEC=UNCLASSIFIED]

2008-07-08 Thread Edzer Pebesma
Jin, you specified nmax = 7 in the second call, but not in the first. Comparing idp values is easiest when other things remain equal. -- Edzer [EMAIL PROTECTED] wrote: Dear there, I tried to compare idw and gstat in library(gstat). I specified idp as 0.5 in both functions, I expected identical

[R-sig-Geo] Normal score transformation

2008-07-08 Thread Laura Poggio
Dear all, I have a small image with three bands and I would like to normalise the variance of the three bands in order to get a more homogeneous dataset for further elaborations. I am thinking to use the normal score transform, but I have some questions about how to implement in R and especially ho

[R-sig-Geo] why are idw and gstat different? [SEC=UNCLASSIFIED]

2008-07-08 Thread Jin.Li
Dear there, I tried to compare idw and gstat in library(gstat). I specified idp as 0.5 in both functions, I expected identical results, but what I got are different as shown for the first five samples. Please help. Thanks. > data(meuse) > coordinates(meuse) = ~x+y > data(meuse.grid) > grid