Skip to content

Add howto on estimating ascertainment with EpiNow2 - #105

Open
jamesmbaazam wants to merge 7 commits into
mainfrom
add-frac-obs-example
Open

Add howto on estimating ascertainment with EpiNow2#105
jamesmbaazam wants to merge 7 commits into
mainfrom
add-frac-obs-example

Conversation

@jamesmbaazam

Copy link
Copy Markdown
Member

This PR adds a howto guide on how to estimate ascertainment or the fraction of observed cases from a time series of primary and secondary observations.

@jamesmbaazam
jamesmbaazam requested a review from avallecam March 27, 2025 14:56
@jamesmbaazam jamesmbaazam self-assigned this Apr 9, 2025

@avallecam avallecam left a comment

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks @jamesmbaazam for the contribution! I tagged you on comments that I need to have your input to proceed with the suggested changes.

I can take charge of some downstream pending actions:

  • create an estimate-ascertainment/ folder
  • move this entry to this new folder
  • add namespace to functions

Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd
# (scale) with a secondary case
res <- estimate_secondary(
df,
obs = obs_opts(

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@jamesmbaazam a question about the obs argument:

  • when used in estimate_infections() it refers to the % observed reported cases among infections
  • when used in estimate_secondary() it refers to the % observed secondary cases among primary cases.

Is this right? Slightly different if you confirm my interpretation, which could be more explicit in the documentation?

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Or does it always refer to the % of observed primary observations?

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Sorry for lagging this entry. I just filled in a discussion with a question I would like to clarify. epiforecasts/EpiNow2#1102

After it, my question would be:

  • To estimate the ascertainment rate as done by {cfr}, defined as ascertainment = defined / estimated (as in cfr::estimate_ascertainment()), should we add an extra step to the reprex code? For example, assuming that the prior of Normal(0.2,0.2) was the defined CFR.

Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Co-authored-by: James Azam <james.m.azam@gmail.com>

@avallecam avallecam left a comment

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Add some final edits

Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd Outdated
Comment thread internal/sandbox/estimate_ascertainment.qmd
@avallecam
avallecam self-requested a review June 18, 2025 17:45
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants