ynegbinompowersim {NBDesign}R Documentation

Two-sample sample size calculation for negative binomial distribution with variable follow-up

Description

This will calculate the power for the negative binomial distribution for the 2-sample case under different follow-up scenarios: 1: fixed follow-up, 2: fixed follow-up with drop-out, 3: variable follow-up with a minimum fu and a maximum fu, 4: variable follow-up with a minimum fu and a maximum fu and drop-out.

Usage

ynegbinompowersim(nsize=200,r0=1.0,r1=0.5,shape0=1,shape1=shape0,pi1=0.5,
   alpha=0.05,twosided=1,fixedfu=1,type=1,u=c(0.5,0.5,1),ut=c(0.5,1.0,1.5),
   tfix=ut[length(ut)]+0.5,maxfu=10.0,tchange=c(0,0.5,1),
   ratec1=c(0.15,0.15,0.15),ratec0=ratec1,rn=10000)

Arguments

nsize

total number of subjects in two groups

r0

event rate for the control

r1

event rate for the treatment

shape0

dispersion parameter for the control

shape1

dispersion parameter for the treatment

pi1

allocation prob for the treatment

alpha

type-1 error

twosided

1: two-side, others: one-sided

fixedfu

fixed follow-up time for each patient

type

follow-up time type, type=1: fixed fu with fu time fixedfu; type=2: same as 1 but subject to censoring; type=3: depending on entry time, minimum fu is fixedfu and maximum fu is maxfu; type=4: same as 3 but subject to censoring

u

recruitment rate

ut

recruitment interval, must have the same length as u

tfix

fixed study duration, often equals to recruitment time plus minimum follow-up

maxfu

maximum follow-up time, should not be greater than tfix

tchange

a strictly increasing sequence of time points starting from zero at which the drop-out rate changes. The first element of tchange must be zero. The above rates and tchange must have the same length.

ratec1

piecewise constant drop-out rate for the treatment

ratec0

piecewise constant drop-out rate for the control

rn

Number of repetitions

Details

Let τmin\tau_{min} and τmax\tau_{max} correspond to the minimum follow-up time fixedfu and the maximum follow-up time maxfu. Let TfT_f, CC, EE and RR be the follow-up time, the drop-out time, the study entry time and the total recruitment period(RR is the last element of ut). For type 1 follow-up, Tf=τminT_f=\tau_{min}. For type 2 follow-up Tf=min(C,τmin)T_f=min(C,\tau_{min}). For type 3 follow-up, Tf=min(R+τminE,τmax)T_f=min(R+\tau_{min}-E,\tau_{max}). For type 4 follow-up, Tf=min(R+τminE,τmax,C)T_f=min(R+\tau_{min}-E,\tau_{max},C). Let ff be the density of TfT_f. Suppose that YiY_i is the number of event obsevred in follow-up time tit_i for patient ii with treatment assignment ZiZ_i, i=1,,ni=1,\ldots,n. Suppose that YiY_i follows a negative binomial distribution such that

P(Yi=yZi=j)=Γ(y+1/kj)Γ(y+1)Γ(1/kj)(kjui1+kjui)y(11+kjui)1/kj,P(Y_i=y\mid Z_i=j)=\frac{\Gamma(y+1/k_j)}{\Gamma(y+1)\Gamma(1/k_j)}\Bigg(\frac{k_ju_i}{1+k_ju_i}\Bigg)^y\Bigg(\frac{1}{1+k_ju_i}\Bigg)^{1/k_j},

where kj,j=0,1k_j, j=0,1 are the dispersion parameters for control j=0j=0 and treatment j=1j=1 and

log(ui)=log(ti)+β0+β1Zi.\log(u_i)=\log(t_i)+\beta_0+\beta_1 Z_i.

The data will be gnerated according to the above model. Note that the piecewise exponential distribution and the piecewise uniform distribution are genrated using R package PWEALL functions "rpwe" and "rpwu", respectively.

The parameters in the model are estimated by MLE using the R package MASS function "glm.nb".

Value

power

simulation power (in percentage)

Author(s)

Xiaodong Luo

Examples

##calculating the sample sizes
abc=ynegbinompowersim(nsize=200,r0=1.0,r1=0.5,shape0=1,
        pi1=0.5,alpha=0.05,twosided=1,fixedfu=1,
        type=4,u=c(0.5,0.5,1),ut=c(0.5,1.0,1.5),
        tchange=c(0,0.5,1),
        ratec1=c(0.15,0.15,0.15),rn=10)
###Power
abc$power

[Package NBDesign version 2.0.0 Index]