twocon: Confidence Interval on the Response Rate from a Two Stage...

View source: R/phase2.R

twoconR Documentation

Confidence Interval on the Response Rate from a Two Stage Study

Description

Computes a confidence interval and several estimators of the response rate using data from a phase II study with two-stage sampling

Usage

twocon(n1, n2, r1, r, conf = 0.95, dp = 1)

Arguments

n1

Number of cases entered during the first stage

n2

Number of additional cases to be entered during the second stage

r1

max number of responses that can be observed in the first stage without continuing

r

total number responses observed

conf

two-sided confidence level (proportion) for the confidence interval

dp

Affects the ordering of outcomes within the sample space (see below)

Details

First n1 patients are entered on the study. If more than r1 responses are observed, then an additional n2 patients are entered. This function assumes that if the observed number of response r < r1, then only n1 patients were entered.

The estimators computed are the MLE (the observed proportion of responses), a bias corrected MLE, and an unbiased estimator, which is sometimes incorrectly described as the UMVUE.

The confidence interval is based on the exact sampling distribution. However, there is not a universally accepted ordering on the sample space in two-stage designs. The parameter dp can be used to modify the ordering by weighting points within the sample space differently. dp=0 will give the Atkinson and Brown procedure, and dp=1 will order outcomes base on the MLE. The Atkinson and Brown procedure orders outcomes based solely on the number of responses, regardless of the number cases sampled. The MLE ordering defines as more extreme those outcomes with a more extreme value of the MLE (the proportion of responses). Other powers of dp, such as dp=1/2, could also be used. Let R be the number of responses and N=n1 if R<=r1 and N=n1+n2 if R>r1. In general, the outcomes that are more extreme in the high response direction are those with R/(N^dp) >= r/(n^dp), where r and n are the observed values of R and N, and the outcomes that are more extreme in the low response direction are those with R/(N^dp) <= r/(n^dp).

Value

A vector with the lower confidence limit, the upper confidence limit, the bias corrected MLE, the MLE, and the unbiased estimator.

References

Atkinson and Brown (1985), BIOMETRICS 741-744.

See Also

binci

Examples

twocon(14, 18, 3, 4, dp = 0)
twocon(14, 18, 3, 4, dp = 1)


raredd/desmon documentation built on May 7, 2024, 3:46 p.m.