Description:

This function calculates the probability of seropositivity for a given age, based on the seroconversion rate (SCR) and the seroreversion rate (SRR). It uses a reversible catalytic model to calculate the probability that an individual will be seropositive at a specific age.

Parameters:

• SCR: The seroconversion rate.
• SRR: The seroreversion rate.
• t: The age of the individual.

Example:

prob_seropositive(0.03, 0.01, 45)

## [1] 0.6260258

Description:

This function calculates the seroprevalence (the proportion of individuals who are seropositive) in a population, considering an age distribution and using a reversible catalytic model. It aggregates the probability of seropositivity across different ages to provide a total estimate of seroprevalence for the population.

Parameters:

• ages: A vector representing the distribution of ages in the population.
• A_max: The maximum age in the population.
• SCR: The seroconversion rate.
• SRR: The seroreversion rate.

Example:

seroprevalence(rep(1 / 80, 80), 80, 0.03, 0.01)

## [1] 0.5296451

Description:

This function calculates the confidence interval for a seroprevalence estimate using a specified confidence level. It employs binomial confidence interval methods to compute the lower and upper bounds of the interval based on the given seroprevalence estimate and sample size. This function uses some of the methods available in the binom.confint function of the binom package (see the references for a link to the package).

Parameters:

• SP: The estimated seroprevalence value (proportion of seropositive individuals).
• n: The sample size used for the estimation.
• conf.level: The desired confidence level (default is 0.95, representing 95% confidence).
• method: The method to calculate the confidence interval (default is “asymptotic”). Available methods: c(“asymptotic”,“exact”,“ac”,“wilson”,“logit”,“cloglog”)

Example:

IC_SP(0.25, 100, conf.level = 0.95, method = "asymptotic")

## [1] 0.1651311 0.3348689

Description:

This function calculates the confidence interval for seroprevalence using the Wald method with a continuity correction. The continuity correction is applied by subtracting from the lower limit and adding to the upper limit 1/(2n).

Parameters:

• SP: The estimated seroprevalence value (proportion of seropositive individuals).
• n: The sample size used for the estimation.
• conf.level: The desired confidence level (default is 0.95, representing 95% confidence).

Example:

IC_SP_Waldcc(0.25, 100, conf.level = 0.95)

## [1] 0.1601311 0.3398689

Description:

This function calculates the confidence interval for the seroconversion rate (SCR) using the confidence interval of seroprevalence. The function uses numerical methods to solve for the SCR that results in the observed seroprevalence.

Parameters:

• SP_interval: A vector containing the lower and upper limits of the confidence interval for the seroprevalence.
• SRR: The seroreversion rate.
• ages: A vector with the distribution of ages in the population.
• A_max: The maximum age in the population.
• limits: A vector specifying the lower and upper bounds for the SCR (default is c(0, 1)).

Example:

IC_SCR(c(0.1, 0.2), 0.01, rep(1 / 80, 80), 80, limits = c(0, 1))

## [1] 0.003433113 0.007502935

Description:

This function estimates the required sample size so that the relative width of the confidence interval for the seroconversion rate (SCR) is equal to a specified value (RL). The function calculates the necessary sample size by iteratively adjusting the sample size until the confidence interval for SCR meets the desired width criteria. This calculation is based on the seroprevalence, confidence intervals for seroprevalence, and the seroconversion rate.

Parameters:

• SCR: The seroconversion rate.
• RL: The desired relative width for the confidence interval width for the seroconversion rate.
• SRR: The seroreversion rate.
• ages: A vector representing the distribution of ages in the population.
• A_max: The maximum age in the population.
• limits: The lower and upper limits for the SCR.
• max_iter: The maximum number of iterations for adjusting the sample size (default is 10000).
• conf.level: The confidence level for the confidence interval (default is 0.95).
• method: The method for calculating the confidence interval for seroprevalence, default is "asymptotic".

Example:

A_max <- 80
age_distribution <- rep(1 / A_max, A_max)
sample_s(0.03, 1, 0.01, age_distribution, A_max, limits = c(0, 1))

## $n
## [1] 57
## 
## $ci_sp
## [1] 0.4000717 0.6592185
## 
## $ci_scr
## [1] 0.0187055 0.0486105