## EXAMPLE 1 (from Machin, Campbell, Tan and Tan 2009 p. 113):
## It is anticipated that patients on a particular drug have a mean diastolic
## blood pressure of 96 mmHg, as against 94 mmHg on an alternative. It is also
## anticipated that the standard deviation of diastolic BP is approximately
## 8 mmHg. If one wishes to confirm that the difference is likely to be less
## than 5 mmHg, that is, one wishes to show equivalence, how many patients
## are needed to be enrolled in the trial? Assume 80% power and
## 95% significance.
epi.ssequc(treat = 94, control = 96, sigma = 8, delta = 5, n = NA,
power = 0.80, r = 1, type = "equivalence", nfractional = FALSE,
alpha = 0.05)
## A total of 244 subjects need to be enrolled in the trial, 122 in the
## treatment group and 122 in the control group.
## EXAMPLE 2 (from Chow S, Shao J, Wang H 2008, p. 64):
## A pharmaceutical company is interested in conducting a clinical trial
## to compare two cholesterol lowering agents for treatment of patients with
## congestive heart disease using a parallel design. The primary efficacy
## parameter is the LDL. In what follows, we will consider the situation
## where the intended trial is for testing equivalence of mean responses
## in LDL. Assume that 80% power is required at a 5% level of significance.
## In this example, we assume a 5 unit (i.e., delta = 5) change of LDL is
## considered of clinically meaningful difference. Assume the standard
## of LDL is 10 units and the LDL concentration in the treatment group is 20
## units and the LDL concentration in the control group is 21 units.
epi.ssequc(treat = 20, control = 21, sigma = 10, delta = 5, n = NA,
power = 0.80, r = 1, type = "equivalence", nfractional = FALSE,
alpha = 0.05)
## A total of 216 subjects need to be enrolled in the trial, 108 in the
## treatment group and 108 in the control group.
## EXAMPLE 2 (cont.):
## Suppose only 150 subjects were enrolled in the trial, 75 in the treatment
## group and 75 in the control group. What is the estimated study power?
epi.ssequc(treat = 20, control = 21, sigma = 10, delta = 5, n = 150,
power = NA, r = 1, type = "equivalence", nfractional = FALSE,
alpha = 0.05)
## With only 150 subjects enrolled the estimated study power is 0.58.
## EXAMPLE 3:
## A study is to be carried out to determine if the daily weight gain of calves
## (in kilograms) on two diets differ. The null hypothesis is that there is
## no difference in daily weight gain. The alternative hypothesis is that
## daily weight gains for the two groups differ. How many calves need to
## be enrolled into a study to detect a statistically significant difference
## in daily weight gain for the two groups? Assume daily weight gain for the
## treatment group is 0.70 kg per day and daily weight gain for the control
## group is 0.65 kg per day and the standard deviation of daily weight gain in
## calves of this age is 0.06. Assume power of 0.80 and alpha of 0.05.
epi.ssequc(treat = 0.70, control = 0.65, sigma = 0.06, delta = NA, n = NA,
power = 0.80, r = 1, type = "equality", nfractional = FALSE,
alpha = 0.05)
## A total of 48 calves need to be enrolled into the study: 24 in the treatment
## group and 24 in the control group.
Run the code above in your browser using DataLab