run_meta: obtain input to rmeta meta.summaries, metaplot

In vjcitn/sars2app: COVID-19 data resources and analysis tools cut down for an interactive app

Description

obtain input to rmeta meta.summaries, metaplot

Usage

run_meta(nyd, opt_parms, ...)

Arguments

nyd	tibble, data source assumed to be nytimes_state_data()
opt_parms	list, such as that produced by min_bic_all_states()
...	passed to Arima_by_state

Examples

nyd = nytimes_state_data()
data(min_bic_2022_01_01)
m1 = run_meta(nyd, opt_parms=min_bic_2022_01_01) # must be relatively current
rmeta::meta.summaries(m1$drifts, m1$se.drifts)
names(m1$drifts) = gsub(".drift", "", names(m1$drifts))
nyind = which(names(m1$drifts) %in% c("New York", "New Jersey"))
rmeta::meta.summaries(m1$drifts[-nyind], m1$se.drifts[-nyind])
o = order(m1$drifts)
rmeta::metaplot(m1$drifts[o], m1$se.drifts[o], labels=names(m1$drifts)[o], cex=.7, 
  xlab="Infection velocity (CHANGE in number of confirmed cases/day)", ylab="State")
segments(rep(-350,46), seq(-49,-4), rep(-50,46), seq(-49,-4), lty=3, col="gray")

vjcitn/sars2app documentation built on Jan. 3, 2022, 12:19 a.m.