censbnd: Calculate bound for censoring distribution.

View source: R/censbnd.R

censbndR Documentation

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 λ_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. https://onlinelibrary.wiley.com/doi/full/10.1002/gepi.21645

Examples

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

survSNP documentation built on Feb. 16, 2023, 10:10 p.m.