Incidence and incidence difference statistics from trinomial prevalences of HIV and recency

1 2 3 |

`PrevH` |
Prevalence of HIV (vector/integer). |

`RSE_PrevH` |
Relative Standard Error (RSE) of estimate for population prevalence of HIV (vector/integer). |

`PrevR` |
Proportion of persons found to be 'recent' by biomarker assay among total persons found positive for HIV (vector/integer). |

`RSE_PrevR` |
Relative Standard Error (RSE) of estimate for population proportion of those testing positive for HIV who have been infected recently (vector/integer). |

`Boot` |
True/False variable indicating whether variance of point estimates is to be calculated by Empirical Bootstrapping (TRUE) or Delta Method (FALSE), the default setting. |

`BS_Count` |
Specifies number of bootstrap samples for bootstrapped confidence intervals of incidence. |

`alpha` |
test rejection threshold. |

`BMest` |
Biomarker estimation by one the 3 options 'same.test'(=default), 'FRR.indep', 'MDRI.FRR.indep' (string). |

`MDRI` |
mean duration of recent infection [days] (vector/integer). |

`RSE_MDRI` |
Relative standard error of MDRI [days] (vector/integer). |

`FRR` |
False recent rate (vector/integer). |

`RSE_FRR` |
Relative standard error of FRR (vector/integer). |

`BigT` |
post-infection time cut-off true vs false recent [days] default 730 days (integer). |

`Covar_HR` |
Covariance of probability of being positive and being categorized recent from survey (vector/integer). Note that as the variances of PrevH and PrevR are often quite small, only a suitably commensurate covariance will enable the inversion of the bootstrap covariance matrix for random number generation to proceed without error. |

Implements assay-based incidence estimation through cross-sectional prevalence and recency of infection tests as described by Kassanjee, et al. 'A new general biomarker-based incidence estimator,' *Epidemiology* (2012). Function parameters must be specified to include assay test characteristics and survey results as proportions. Confidence intervals are computed via Delta method approximation, except when Boot=TRUE is specified, in which case confidence intervals are generated by empirical bootstrap resampling. Inputs must be in appropriate ranges for appropriate units. Extreme input values may make calculation impossible, and if entered will elicit error notices.

The package contains long form documentation in the form of vignettes that cover the use of the main fucntions. Use browseVignettes(package="inctools") to access them.

Returns incidence and incidence difference statistics, annualised risk and implied CIs on test characteristics.

1 2 3 4 5 6 7 8 9 10 11 | ```
incprops(PrevH = 0.20, RSE_PrevH = 0.028, PrevR = 0.10, RSE_PrevR = 0.09,
BS_Count = 10000, Boot = TRUE, MDRI = 200, RSE_MDRI = 0.05,
FRR = 0.01,
RSE_FRR = 0.2, BigT = 730)
incprops(PrevH = c(0.20,0.21,0.18), RSE_PrevH = c(0.028,0.03,0.022),
PrevR = c(0.10,0.13,0.12), RSE_PrevR = c(0.094,0.095,0.05),
BS_Count = 10000, Boot = FALSE, BMest = 'MDRI.FRR.indep', MDRI = c(200,180,180),
RSE_MDRI = c(0.05,0.07,0.06), FRR = c(0.01,0.009,0.02), RSE_FRR = c(0.2,0.2,0.1),
BigT = 730)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.