sim.snp.expsurv.sctest {survSNP}

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

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