# NOT RUN {
df <- data.frame(id = seq(1:20),
n_tested = rpois(20, 6),
N = rpois(20, 50),
test_Se = 0.3,
dp = 0.05)
## Calculate the herd level sensitivity for each of these herds. If
## the ratio of the number tested to number of animals in the herd
## exceeds the threshold then the finite method is used, otherwise the
## infinite method is used.
herd_Se <- hse(df$id,
df$n_tested,
df$N,
df$test_Se,
df$dp,
threshold = 0.1)
## Calculate the system sensitivity given the testing and sensitivity
## in these herds:
sysse(dp = rep(0.10, nrow(herd_Se)),
hse = herd_Se$HSe)
# }
Run the code above in your browser using DataLab