alteval.rc {ALTopt} R Documentation

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.

`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))
```

[Package ALTopt version 0.1.2 Index]