goodStrata: Find the strata that have the smallest distance at the end

In RJauslin/SystematicSampling: Systematic Sampling

Description

Find the strata that have the smallest distance at the end

Usage

goodStrata(u_Strata)

Arguments

X	sptaial coordinates
pik_tmp	temporary pik (unit k has a cutted pik)
toreBound

Examples

############################# 1D

rm(list = ls())
set.seed(1)
eps <- 1e-13
N <- 10
n <- 3
X <- as.matrix(seq(1,N,1))
tb = 0.001
pik <- rep(0.3,N)
ref <- 6
pik[ref] <- runif(1,min = eps,max = pik[ref] - eps)
tore = TRUE
toreBound = 10
l_Strata <- allStrata(X,pik,ref,tb)
u_Strata <- unionStrata(l_Strata,ref)
s <- goodStrata(u_Strata)
plot(s,X)

################################## 2D

rm(list = ls())
set.seed(1)
eps <- 1e-13
N <- 36
n <- 8
pik <- rep(n/N,N)


tb <- runif(2)/100
X <- as.matrix(cbind(runif(N),runif(N)))
tore = FALSE
toreBound = FALSE
ref = 1
pik[ref] <- runif(1,min = eps,max = pik[ref] - eps)

l_Strata <- allStrata(X,pik,ref,tb)
u_Strata <- unionStrata(l_Strata)
s <- goodStrata(u_Strata)
plot(s,X)

RJauslin/SystematicSampling documentation built on Feb. 10, 2021, 9:37 a.m.