mgwrsar (version 1.1)

find_TP: Search of a suitable set of target points. find_TP is a wrapper function that identifies a set of target points based on spatial smoothed OLS residuals.

Description

Search of a suitable set of target points. find_TP is a wrapper function that identifies a set of target points based on spatial smoothed OLS residuals.

Usage

find_TP(formula, data,coords,kt,ks=16,Wtp=NULL,type='residuals',
model_residuals=NULL,verbose=0,prev_TP=NULL,nTP=NULL)

Value

find_TP returns an index vector of Target Points set.

Arguments

formula

a formula

data

a dataframe or a spatial dataframe (SP package)

coords

a dataframe or a matrix with coordinates, not required if data is a spatial dataframe

kt

the minimum number of first neighbors with lower (resp.higer) absolute value of the smoothed residuals.

ks

the number of first neighbors for computing the smoothed residuals, default 16.

Wtp

a precomputed matrix of weights, default NULL.

type

method for choosing TP, could be 'residuals', ' kdtree','random', default 'residuals'

model_residuals

(optional) a vector of residuals.

verbose

verbose mode, default FALSE.

prev_TP

index of already used TP (version length(kt)>1), default NULL.

nTP

numbeer of target points for random choice of target points, default NULL.

Details

find_TP is a wrapper function that identifies a set of target points, based on spatial smoothed residuals by default. If no vector of residuals are provided, OLS residuals are computed. The function first computes the smooth of model residuals using a Shepard's kernel with ks neighbors (default 16). Then it identifies local maxima (resp. minima) that fits the requirement of having at least kt neighbors with lower (resp.higer) absolute value of the smoothed residuals. As kt increases the number of target points decreases.

Examples

# \donttest{
 library(mgwrsar)
 ## loading data example
 data(mydata)
 coords=as.matrix(mydata[,c("x","y")])
 TP=find_TP(formula = 'Y_gwr~X1+X2+X3', data =mydata,coords=coords,kt=6,
 type='residuals')
 # only 60 targets points are used
 length(TP)

 model_GWR_tp<-MGWRSAR(formula = 'Y_gwr~X1+X2+X3', data = mydata,
 coords=coords, fixed_vars=NULL,kernels=c('gauss'),  H=0.03, Model = 'GWR',
 control=list(SE=TRUE,TP=TP,ks=12))
 summary(model_GWR_tp@Betav)
 # }