Prepare provider-level summary statistics for medication indicators
prepare_provider_summary.RdGenerates detailed provider-specific statistics for each medication indicator, including comparisons with overall averages, facility-level variations, and outlier analyses. For most indicators, statistics are calculated separately for all residents and by dementia status.
Arguments
- indicator
Character string specifying the medication indicator. Must be one of: "psychotropics", "antipsychotics", "antidepressants", "antiepileptics", "hypnotics", "polypharmacy", "antibiotics", "antivirals", or "opioids".
- data_summary
List containing:
provider_mean: Provider-specific statistics
overall_mean: Population-level statistics for comparison
Each element must contain columns for indicator, variant, and relevant statistics (mean_cmar, min_cmar, max_cmar).
- config
List containing provider configuration settings including:
provider_name: Full name of the provider
provider_abbr: Provider abbreviation (e.g., "AC", "BC", "SC")
- dta
data.table containing the full dataset needed for outlier analysis
Value
A nested list structure containing provider-specific statistics. The exact structure varies by indicator type:
For most medications: overall statistics plus dementia/non-dementia breakdowns
For polypharmacy: prescribed and administered statistics
For antibiotics: DOT rates and usage statistics
For antivirals: overall statistics plus COVID-specific information
Each section typically includes:
provider: Provider-specific rate
average: Overall population rate for comparison
comparison: Text describing how provider compares to average
facmin/facmax: Minimum and maximum facility rates
outliertxt: Description of facilities exceeding warning/action limits
Note
This function requires dta to be in the environment for the outlier calculations. This global dependency should be refactored in future versions.
Examples
if (FALSE) { # \dontrun{
# Prepare summary for antipsychotics
antipsych_summary <- prepare_provider_summary(
indicator = "antipsychotics",
data_summary = list(
provider_mean = provider_stats,
overall_mean = overall_stats
),
config = list(
provider_name = "Anglicare",
provider_abbr = "AC"
),
dta = facility_data
)
} # }