utility.cc {clinUtiDNA} | R Documentation |
Calculation of the utility index for case-control data
utility.cc(cases, contr, pD,pG)
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 |
see Nguyen et al. 2013
Utility.cc |
a list containing the utility estimate, the corresponding variance and the boundaries of the 95% of the confidence interval |
Thuy Trang Nguyen
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
## 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