R/mgwrsar_bootstrap_test_all.R

In mgwrsar: GWR and MGWR with Spatial Autocorrelation

#' A bootstrap test for testing nullity of all Betas for mgwrsar class model,
#' @usage mgwrsar_bootstrap_test_all(model,B=100,domc=NULL)
#' @param model A mgwrsar model
#' @param B number of bootstrap replications, default 100
#' @param domc  If TRUE, doMC parallelization
#' @seealso  mgwrsar_bootstrap_test
#' @return a matrix with statistical test values and p ratios
mgwrsar_bootstrap_test_all <-
function(model,B=100,domc=NULL){
x1=model
k=ncol(model$X)
names_x=colnames(model$X)
names_x<-names_x[-which(names_x=='Intercept')]
res<-matrix(nrow = length(names_x), ncol = 2)
for(i in 1:(k-1)){
name_x<-names_x[i]

x0<-MGWRSAR(formula=update(x1$formula,paste('~.-',name_x,sep='')),data=x1$data,coord=x1$coord,fixed_vars=x1$fixed_vars,kernels=x1$kernel,H=x1$H,Model=x1$Model,control=list(Method=x1$Method,W=x1$W,isgcv=FALSE,SE=TRUE))
cat(name_x,' ')
res[i,]<-unlist(mgwrsar_bootstrap_test(x0,x1,B=B,domc=domc,type='standard',eps='H1',df='H1',focal='median',D=NULL))
}
res
}