censbnd: Calculate bound for censoring distribution.

Description Usage Arguments Details Value Author(s) References Examples

Description

This function computes the bound

Usage

1
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 λ_0, λ_1 and λ_2 with mixing proportion p_0,p_1 and p_2. 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. http://onlinelibrary.wiley.com/doi/10.1002/gepi.21645/abstract

Examples

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

Example output

Loading required package: survival
Loading required package: Rcpp
Loading required package: lattice
Loading required package: foreach
Loading required package: xtable
[1] 2.145579

survSNP documentation built on May 2, 2019, 2:48 p.m.