survSNP (version 0.26)

censbnd: Calculate bound for censoring distribution.

Description

This function computes the bound

Usage

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

Value

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

Arguments

lambda

Baseline exponential hazard rate

p

Relative genotype frequency

crate

Desired censoring rate

rootint

Interval to be searched for the root

Author

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

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.

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