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

In survSNP: Power Calculations for SNP Studies with Censored Outcomes

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)

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

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 λ_g. 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.

Value

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

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

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