geozoning (version 1.0.0)

SigmaL2: SigmaL2

Description

SigmaL2

Usage

SigmaL2(zlab, listZonePoint, tabVal, surfVoronoi)

Arguments

zlab

list with zone numbers for each zone label

listZonePoint

list of indices of data points within zones, result of call to calNei

tabVal

SpatialPointsDataFrame containing data values

surfVoronoi

Surfaces of the Voronoi polygons corresponding to data pts

Value

a list with components

cL

weighted (with Voronoi surfaces) average of per label variances

SigmaL2

vector of per label variances

SL

vector of per label Voronoi surfaces

mL

vector of weighted (with Voronoi surfaces) per label average values

voroLab

vector of per label data

Details

compute overall mean and variance of all zones for each label plus sum of them for all labels

Examples

# NOT RUN {
data(mapTest)
# run zoning with 2 quantiles corresponding to probability values 0.4 and 0.7
# save initial zoning and last level zonings
criti=correctionTree(c(0.4,0.7),mapTest,SAVE=TRUE)
K=criti$zk[[2]][[1]]
uni=unique(K$lab)
zlab=sapply(uni,function(x){(1:length(K$lab))[K$lab==x]})
sig=SigmaL2(zlab,K$listZonePoint,mapTest$krigData,mapTest$krigSurfVoronoi)
# }