R/rectangle_adapt_sorted.R

In mgwrsar: GWR and MGWR with Spatial Autocorrelation

#' rectangle_adapt_sorted
#' to be documented
#' @usage rectangle_adapt_sorted(d,h)
#' @param d a vector of distances
#' @param h a scalar or vector of number of neighbours as bandwidths (lenght(h)=1 or length(h)=length(d))
#' @noRd
#' @return a vector of weights.
rectangle_adapt_sorted<-function (d, h) {
  if(is.unsorted(d[1,])) {
    dd=t(apply(d,1,sort))
  } else dd=d
if(any(h!=h[1])) {h <- dd[cbind(1:nrow(dd),h+1)]} else h=dd[,h[1]+1]
x<-d/h
x[abs(x)> 1]<- 0
x[abs(x)>0]<- 1
x[,1]<-1
x<-x/rowSums(x)
x
}