uniform.crit {SFDesign}R Documentation

Uniform criterion

Description

This function calculates the wrap-around discrepancy of a design.

Usage

uniform.crit(design)

Arguments

design

a design matrix.

Details

uniform.crit calculates the wrap-around discrepancy of a design. The wrap-around discrepancy for a design D=[\bm x_1, \dots, \bm x_n]^T is defined as (Hickernell, 1998):

\phi_{wa} = -\left(\frac{4}{3}\right)^p + \frac{1}{n^2}\sum_{i,j=1}^n\prod_{k=1}^p\left[\frac{3}{2} - |x_{ik}-x_{jk}|(1-|x_{ik}-x_{jk}|)\right].

Value

wrap-around discrepancy of the design

References

Hickernell, F. (1998), “A generalized discrepancy and quadrature error bound,” Mathematics of computation, 67, 299–322.

Examples

n = 20
p = 3
D = randomLHD(n, p)
uniform.crit(D)


[Package SFDesign version 0.1.2 Index]