alteval.rc {ALTopt}

Design evaluation with right censoring.

Description

alteval.rc calculates the objective function value (D, U or I) for a given design with right censoring plan.

Usage

alteval.rc(
  designTable,
  optType,
  tc,
  nf,
  alpha,
  formula,
  coef,
  useCond,
  useLower,
  useUpper
)

Arguments

designTable

a data frame containing the coordinates and the number of allocation of each design point. The design created by either altopt.rc or altopt.ic or any design matrix with the same form as those can be provided for this argument.

optType

the choice of "D", "U" and "I" optimality.

tc

the censoring time.

nf

the number of stress factors.

alpha

the value of the shape parameter of Weibull distribution.

formula

the object of class formula which is the linear predictor model.

coef

the numeric vector containing the coefficients of each term in formula.

useCond

the numeric vector of use condition. It should be provided when optType is "U". The length of the vector should be same as the number of stress factors.

useLower

the numeric vector of lower bound of use region. It should be provided when optType is "I". The length of the vector should be same as the number of stress factors.

useUpper

the numeric vector of upper bound of use region. It should be provided when optType is "I". The length of the vector should be same as the number of stress factors.

Value

The objective function value corresponded by optType for a given design with right censoring plan.

See Also

altopt.rc

Examples

# Evaluation of factorial design for right censoring.
x1 <- c(0, 1, 0, 1)
x2 <- c(0, 0, 1, 1)
allocation <- c(25, 25, 25, 25)
facDes <- data.frame(x1, x2, allocation)

alteval.rc(facDes, "D", 100, 2, 1, formula = ~ x1 + x2 + x1:x2,
coef = c(0, -4.086, -1.476, 0.01))

alteval.rc(facDes, "U", 100, 2, 1, formula = ~ x1 + x2 + x1:x2,
coef = c(0, -4.086, -1.476, 0.01), useCond = c(1.758, 3.159))

alteval.rc(facDes, "I", 100, 2, 1, formula = ~ x1 + x2 + x1:x2,
coef = c(0, -4.086, -1.476, 0.01), useLower = c(1.458, 2.859), useUpper = c(2.058, 3.459))