censbnd {survSNP}

Calculate bound for censoring distribution.

Description

This function computes the bound

Usage

censbnd(lambda, p, crate, rootint = c(0.1, 1000))

Arguments

lambda

Baseline exponential hazard rate

p

Relative genotype frequency

crate

Desired censoring rate

rootint

Interval to be searched for the root

Details

The time to event distribution T is assumed to be a mixture of exponentials with parameter lambda0, lambda1 and lambda2 with mixing proportion p0,p1 and p2. Suppose that the censoring distribution is uniform on the interval [0,b]. This function calculates b for a desired censoring rate. It is not intended to be called directly by the user. To conduct power calculations, use sim.snp.expsurv.power or the convenience wrapper function survSNP.power.table.

Value

This function returns a list from the uniroot function. The root component of this list in the

Author(s)

Kouros Owzar, Zhiguo Li, Nancy Cox, Sin-Ho Jung and Chanhee Yi

References

Kouros Owzar, Zhiguo Li, Nancy Cox and Sin-Ho Jung. Power and Sample Size Calculations for SNP Association Studies with Censored Time-to-Event Outcomes. https://onlinelibrary.wiley.com/doi/full/10.1002/gepi.21645

Examples

censbnd(0.1,hwe(0.1),0.9)$root