Re: [R-sig-Geo] Holdout Sampling Adaptive Bandwidth SPGWR

2013-08-29 Thread Paul Bidanset
Thank you Luis. When calibrating the adaptive model, using adapt=t in the bandwidth selection created the proportion you speak of, which then allowed me to create a bandwidth matrix using gwr.adapt. However, this has not worked for me with holdout samples. Have you had success in this regard? I do

Re: [R-sig-Geo] Holdout Sampling Adaptive Bandwidth SPGWR

2013-08-29 Thread Luis Guerra
Dear Paul, I am dealing with this kind of problems right now, and if I am not wrong, when you want to apply an adaptative bandwidth, you should introduce a value for the "adapt" parameter instead of for the "bandwidth" parameter. This value will be between 0 and 1 and indicates the proportion of c

[R-sig-Geo] Holdout Sampling Adaptive Bandwidth SPGWR

2013-08-29 Thread Paul Bidanset
Hi Folks, I was curious if anyone has had experience applying an SPGWR model with an adaptive bandwidth matrix to a holdout or validation sample. I am using the "fit.points" command, which does not seem to allow for a new bandwidth calibrated around the holdout samples XY coordinates. Any directio

Re: [R-sig-Geo] Import custom projection in R

2013-08-29 Thread Pan
Dear all, Thanks for your help! CRS() is working fine now: > library(rgdal) > My.Proj<-"+proj=tmerc +datum=WGS84 +lat_0=0 +lon_0=17 +x_0=60 +y_0=1000" > CRS(My.Proj) CRS arguments: +proj=tmerc +datum=WGS84 +lat_0=0 +lon_0=17 +x_0=60 +y_0=1000 +ellps=WGS84 +towgs84=0,0,0 an

Re: [R-sig-Geo] Import custom projection in R

2013-08-29 Thread Roger Bivand
On Thu, 29 Aug 2013, Pan wrote: Dear Roger, Thank you for your valuable suggestions. First, I tried to enter my.proj as PROJ.4 string: > MyProj.PROJ.4<-"+proj=UTM +datum=WGS84 +lon_0=17 +x_0=60 +y_0=1000" > > library(rgdal) > > CRS("+proj=UTM +lon_0=17 +x_0=60 +y_0=1000") E

Re: [R-sig-Geo] Import custom projection in R

2013-08-29 Thread Barry Rowlingson
On Thu, Aug 29, 2013 at 11:46 AM, Pan wrote: > Dear Roger, > > Thank you for your valuable suggestions. First, I tried to enter my.proj as > PROJ.4 string: > > > MyProj.PROJ.4<-"+proj=UTM +datum=WGS84 +lon_0=17 +x_0=60 > +y_0=1000" > > > > library(rgdal) > > > > CRS("+proj=UTM

Re: [R-sig-Geo] Import custom projection in R

2013-08-29 Thread Pan
Dear Roger, Thank you for your valuable suggestions. First, I tried to enter my.proj as PROJ.4 string: > MyProj.PROJ.4<-"+proj=UTM +datum=WGS84 +lon_0=17 +x_0=60 +y_0=1000" > > library(rgdal) > > CRS("+proj=UTM +lon_0=17 +x_0=60 +y_0=1000") Error in CRS("+proj=UTM +lon_