enforceDomainOnParameters {RCEIM}

Enforce domain boundaries

Description

A small function to assure that the domains are respected during the optimization process. If any of them not respected, the ofending parameters are replaced by the value of the nearest border.

Usage

enforceDomainOnParameters(paramsArray, domain)

Arguments

paramsArray

The array with the parameters to check.

domain

The domain boudaries.

Value

The parameter array, with ofending values replaced if necessary.

Author(s)

Alberto Krone-Martins

Examples

# Creates a random set of parameters in an interval larger than a certain domain
# and apply the enforceDomainOnParameters function and represent graphically
# the parameters before and after the function.
dev.new()
paramArr <- matrix((runif(100)-0.5)/0.5*13, nrow=50)
domain <- matrix(c(-10, -10, 10, 10), ncol=2)
newParamArr <- enforceDomainOnParameters(paramArr, domain)
plot(paramArr[,1], paramArr[,2], xlab="x", ylab="y", main="black: input\n red: output")
points(newParamArr[,1], newParamArr[,2], col="red", pch=19, cex=0.7)