hu.mackey.thomas.linear.power: Random coefficient regression models (RCRM) sample size...

In longpower: Sample Size Calculations for Longitudinal Data

Random coefficient regression models (RCRM) sample size calculations

Description

This function computes sample size and power needed for the random coefficient regression models (RCRM) based on the formula from Hu, Mackey, and Thomas (2021). The RCRM assumes that the experimental and control arms have the same population baseline value.

Usage

hu.mackey.thomas.linear.power(
  n = NULL,
  delta = NULL,
  power = NULL,
  t = NULL,
  lambda = 1,
  sig2.i = 0,
  cor.s.i = NULL,
  sig2.s = 0,
  sig2.e = NULL,
  p = NULL,
  sig.level = 0.05,
  alternative = c("two.sided", "one.sided"),
  tol = .Machine$double.eps^2
)

Arguments

n	sample size, group 1. This formula can accommodate unbalanced group allocation via lambda.
delta	Effect size (absolute difference in rate of decline between tx and placebo)
power	power
t	Vector of visit time points (including time 0)
lambda	allocation ratio (sample size group 1 divided by sample size group 2)
sig2.i	Variance of random intercept
cor.s.i	Correlation between random intercept & slope
sig2.s	Variance of random slope
sig2.e	Variance of pure error
p	proportion vector for both groups; if i indexes visits, p[i] = the proportion whose last visit was at visit i (p sums to 1)
sig.level	type one error
alternative	one- or two-sided test
tol	numerical tolerance used in root finding

Details

See Hu. Mackey, and Thomas (2021) for parameter details.

See Equations (7) and (8) in Hu, Mackey, and Thomas (2021)

Value

One of the number of subject required per arm, the power, or detectable effect size given sig.level and the other parameter estimates.

Author(s)

Monarch Shah

References

Hu, N., Mackey, H., & Thomas, R. (2021). Power and sample size for random coefficient regression models in randomized experiments with monotone missing data. Biometrical Journal, 63(4), 806-824.

See Also

lmmpower, diggle.linear.power, liu.liang.linear.power, edland.linear.power

Examples

## Not run: 
browseVignettes(package = "longpower")

## End(Not run)
# An Alzheimer's Disease example using ADAS-cog pilot estimates
t <- seq(0,1.5,0.25)
p <- c(rep(0, 6),1)

hu.mackey.thomas.linear.power(delta=1.5, t=t, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p, power=0.80)
hu.mackey.thomas.linear.power(n=180, t=t, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p, power=0.80)
hu.mackey.thomas.linear.power(n=180, delta=1.5, t=t, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p)

hu.mackey.thomas.linear.power(delta=1.5, t=t, lambda=2, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p, power=0.80)
hu.mackey.thomas.linear.power(n=270, t=t, lambda=2, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p, power=0.80)
hu.mackey.thomas.linear.power(n=270, delta=1.5, t=t, lambda=2, 
  sig2.s=24, sig2.e=10, p=p, cor.s.i=0.5)

hu.mackey.thomas.linear.power(delta=1.5, t=t, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p, power=0.80, alternative='one.sided')
hu.mackey.thomas.linear.power(n=142, t=t, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p, power=0.80, alternative='one.sided')
hu.mackey.thomas.linear.power(n=142, delta=1.5, t=t, 
  sig2.s=24, sig2.e=10, cor.s.i=0.5, p=p, sig.level=0.05, alternative='one.sided')

longpower documentation built on Jan. 11, 2023, 5:12 p.m.