File: epi.dgamma.R

package info (click to toggle)
r-cran-epir 2.0.80%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 2,332 kB
  • sloc: makefile: 5
file content (8 lines) | stat: -rw-r--r-- 307 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
 
"epi.dgamma" <- function(rr, quantiles = c(0.05, 0.95)){
   fold.variation <- rr[2]/rr[1]
   low.p <- abs(qnorm(quantiles[1], mean = 0, sd = 1))
   up.p <- abs(qnorm(quantiles[2], mean = 0, sd = 1))
   p <- low.p + up.p
   tau <- (p^2) / (log(fold.variation) * log(fold.variation))
   return(tau)
}