sim.snp.expsurv.sctest: Simulation of Cox Score Statistic

Description

This function simulates the asymptotic P-value for the Cox score statistic under a specified genetic risk model

Usage

sim.snp.expsurv.sctest(n, gtprev, lam, a, b, ztest, diag = FALSE)

Value

A vector of length 2 containing the observed event rate and the asymptotic P-value for the Cox score test

Arguments

n: Sample size
gtprev: Relative genotypic frequency
lam: Exponential hazard rates for conditional time to event survival functions
a: Lower bound for the uniform censoring bound
b: Upper bound for the uniform censoring bound
ztest: Assumed genetic model
diag: Set to TRUE if print out of diagnostics is desired

Author

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

Details

This function simulates the genotypes 0, 1 or 2 from relative genotypic frequencies. Then conditional on genotype, it simulates the time to event from an exponential distribution with parameter \(lambdag\). The censoring time is drawn from a uniform law on the interval \([a,b]\). 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