dij {LHD}

Calculate the Inter-site Distance

Description

dij returns the inter-site distance of two design points of an LHD

Usage

dij(X, i, j, q = 1)

Arguments

X

A matrix object. In general, X stands for the design matrix.

i

A positive integer, which stands for the i^th row of X.

j

A positive integer, which stands for the j^th row of X. Both i and j should be in [1,nrow(X)] and they should not be equal to each other.

q

The default is set to be 1, and it could be either 1 or 2. If q is 1, dij is the Manhattan (rectangular) distance. If q is 2, dij is the Euclidean distance.

Value

If all inputs are logical, then the output will be a positive number indicating the distance. dij = \left\{ \sum_{k=1}^{m} \vert x_{ik}-x_{jk}\vert ^q \right\}^{1/q}

Examples

#create a toy LHD with 5 rows and 3 columns
toy=rLHD(n=5,k=3);toy

#Calculate the inter-site distance of the 2nd and the 4th row of toy (with default q)
dij(X=toy,i=2,j=4)

#Calculate the inter-site distance of the 2nd and the 4th row of toy (with q=2)
dij(X=toy,i=2,j=4,q=2)