the local model to be fitted using the same syntax used in the zeroinfl function of the R package pscl. This is a sting (a symbolic description of the model) that is passed to the sub-models' zeroinfl function. For more details look at the details of the zeroinfl function.
family
a specification of the count model family to be used in the local model as in the zeroinfl function. Currently the only option tested is "poisson".
dframe
a numeric data frame of at least two suitable variables (one dependent and one independent)
coords
a numeric matrix or data frame of two columns giving the X,Y coordinates of the observations
kernel
the kernel to be used in the regression. Options are "adaptive" or "fixed". The weighting scheme used here is defined by the bi-square function (weight = (1-(ndist/H)^2)^2 for distances less than or equal to H, 0 otherwise)
algorithm
a character argument that specifies whether the function will use an exhaustive or a heuristic algorithm. In the first case all possible bandwidths within a range are being tested. In the second case the optim function is being used allowing for the choice of various optimisation methods (such as Brent or BFGS) that may find a global or local optimum. The default algorithm is "exhaustive"
optim.method
the optimisation method to be used. A detailed discussion is available at the 'Details' section of the function optim (stats). Example methods are "Nelder-Mead", "Brent", "BFGS", "CG" and "L-BFGS-B". The default method is "Nelder-Mead".
b.min
the minimum bandwidth. This is important for both algorithms. In the case of the exhaustive algorithm it sets the lower boundary for the range in which the function will compute the CV score for each possible bandwidth. In the case of the heuristic algorithm it provides the initial value for the bandwidth to be optimised which is very important. In the latter case b.min and b.max should be provided if the optimisation method "L-BFGS-B" or "Brent" has been selected.
b.max
the maximum bandwidth. This is important for both algorithms. In the case of the exhaustive algorithm it sets the upper boundary for the range in which the function will compute the CV score for each possible bandwidth. In the case of the heuristic algorithm b.max and b.min should be provided if the optimisation method "L-BFGS-B" or "Brent" has been selected.
step
this numeric argument is used only in the case of a fixed kernel indicating the increment of the sequence of bandwidths in between the b.min and the b.max. In the case of the adaptive kernel the increment is 1 neighbour.
Details
Please carefully read the function optim(stats) when using a heuristic algorithm.
Value
bw
The optimal bandwidth (fixed or adaptive)
CV
The corresponding Cross Validation score for the optimal bandwidth
CVs
Available only in the case of the exhaustive algorithm. This is a numeric matrix in which the first column refers to the bandwidth in test and the second to the corresponding CV score.
Warning
Large datasets increase the processing time.
Note
Please select the optimisation algorithm carefully. This function needs further testing. Please report any bugs!
Author(s)
Stamatis Kalogirou <skalo@hua.gr>
References
Kalogirou, S. (2015) Destination Choice of Athenians: an application of geographically weighted versions of standard and zero inflated Poisson spatial interaction models, Geographical Analysis, On line first. DOI: 10.1111/gean.12092