File: epi.pooled.R

package info (click to toggle)
r-cran-epir 2.0.80%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, trixie
  • size: 2,332 kB
  • sloc: makefile: 5
file content (16 lines) | stat: -rw-r--r-- 332 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
epi.pooled <- function(se, sp, P, m, r){
   # Herd specificity:
   PlSp <- sp^m
   HSp <- (PlSp)^r

   # Herd sensitivity:
   HSe <- 1 - ((1 - (1 - P)^m) * (1 - se) + (1 - P)^m * PlSp)^r
   
   # Herd level apparent prevalence:
   HAPneg <- 1 - HSp

   rval <- list(HAPneg = HAPneg, HSe = HSe, HSp = HSp)
   rval
}