Nothing
build.lut <-
function(LEVEL=6, REPSIM=5, RAJZ=FALSE, CIM="", ENV="data") {
#--------------------------------------------------------------
#
# TITLE: build.lut()
# AUTHOR: TARMO REMMEL (ORIGINAL BY FERKO CSILLAG)
# DATE: 23 January 2020
# CALLS: wtest.run()
# CALLED BY: NA
# NEEDS: NA
# NOTES: USES wtest.run() TO RUN A LOOP ALONG AUTOCORRELATION
# (rho) AND PROPORTION AND THE RESULTS ARE SUMMARIZED
# AS A LOOKUP TABLE FOR BIAS ADJUSTMENT ON RHO
# ESTIMATION ON A BINARY MAP.
#--------------------------------------------------------------
# SAVE GRAPHIC PARAMETERS AND RESTATE THEM ON EXIT
opar <- par(no.readonly =TRUE)
on.exit(par(opar))
DIFF <- rep(0, 110)
dim(DIFF) <- c(10, 11)
IX <- 0
IY <- 0
for(luprho in seq(0, 0.2499999, 0.0277777)) {
RHO <- luprho
IX <- IX + 1
IY <- 0
for(lupcprop in seq(0.1, 0.9, 0.1)) {
CPROP <- lupcprop
IY <- IY + 1
for(lup in 1:REPSIM) {
RESULTT <- wtest.run(REPSIM = REPSIM, LEVEL = LEVEL, RHO= RHO, CPROP = CPROP, RAJZ = RAJZ, CIM = CIM, ENV=ENV)
} # END FOR: lup
DIFF[IX, IY] <- median(RESULTT[1, ]) - median(RESULTT[2,])
} # END FOR: lupcprop
} # END FOR: luprho
# RETURN THE BIAS CORRECTION MATRIX TO THE USER
return(DIFF)
} # END FUNCTION: build.lut
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.