Calculation of the utility index for case-control data

Share:

Description

Calculation of the utility index for case-control data

Usage

1
utility.cc(cases, contr, pD,pG)

Arguments

cases

a vector containing numbers of cases observed in the gene environment groups <G=0,E=0>, <G=1,E=0>, <G=0,E=1> and <G=1,E=1>

contr

a vector containing numbers of controls observed in the gene environment groups <G=0,E=0>, <G=1,E=0>, <G=0,E=1> and <G=1,E=1>

pD

disease prevalence should be either less than 1 (for a known value) or a vector of observed numbers of affected and unaffected individuals from an external data set

pG

frequency of the genetic risk variant should be either less than 1 (for a known value) or a vector of observed numbers of carriers and non-carriers observed from an external data set

Details

see Nguyen et al. 2013

Value

Utility.cc

a list containing the utility estimate, the corresponding variance and the boundaries of the 95% of the confidence interval

Author(s)

Thuy Trang Nguyen

References

Making Medical Decisions in Dependence of Genetic Background: Estimation of the Utility of DNA Testing in Clinical, Pharmaco-Epidemiological or Genetic Studies. Nguyen TT, Schaefer H, Timmesfeld N. Genet Epidemiol. 2013 Apr 4. doi: 10.1002/gepi.21701

See Also

gxe.coh, utility.coh, gxe.cc

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
## Parameter initialisation
cases <- c(36,10,84,25)
contr <- c(100,4,63,2)
pD <- c(22,9420)
pG <- c(52,618) 
## Calculation of the utility index
utility.cc(cases,contr,pD,pG)  # with pD and pG estimated from external data
utility.cc(cases,contr,pD,pG=0.078) # with known pG
utility.cc(cases,contr,pD=0.002,pG) #  with known pD
utility.cc(cases,contr,pD=0.002,pG=0.078) # with known pD and pG