Package {syrona}

Type: Package
Title: Stratified Prevalence Comparison Across OMOP CDM Datasets
Version: 0.2.1
Description: Derives stratified prevalence tables from the condition, procedure, and drug records in OMOP CDM (Observational Medical Outcomes Partnership Common Data Model) databases, computes log2 prevalence ratios between paired datasets, and synthesizes them via random-effects meta-analysis at multiple aggregation levels (year, age group, and sex). Between-study variance is estimated with the Paule-Mandel method, as described in Paule and Mandel (1982) <doi:10.6028/jres.087.022>.
License: MIT + file LICENSE
URL: https://github.com/HealthInformaticsUT/Syrona
BugReports: https://github.com/HealthInformaticsUT/Syrona/issues
Encoding: UTF-8
Language: en-US
Depends: R (≥ 4.1.0)
Imports: dplyr (≥ 1.1.0), dbplyr (≥ 2.3.0), DBI (≥ 1.2.0), CDMConnector (≥ 2.0.0), omopgenerics (≥ 1.3.0), tibble, tidyr, cli, rlang, shiny, ggplot2, ggiraph, ggtext, scales, grDevices, stats, readr (≥ 2.0.0), meta (≥ 6.0.0)
Suggests: testthat (≥ 3.0.0), duckdb (≥ 0.9.0), RPostgres, CohortConstructor (≥ 0.6.0), DT, shinycssloaders, knitr, rmarkdown
Config/testthat/edition: 3
VignetteBuilder: knitr
Config/roxygen2/version: 8.0.0
NeedsCompilation: no
Packaged: 2026-08-06 12:34:32 UTC; maarjapajusalu
Author: Maarja Pajusalu [aut, cre]
Maintainer: Maarja Pajusalu <maarja.pajusalu@ut.ee>
Repository: CRAN
Date/Publication: 2026-09-04 20:10:02 UTC

syrona: Stratified prevalence comparison across OMOP CDM datasets

Description

Derives stratified prevalence tables from the condition, procedure, and drug records in OMOP CDM databases, computes log2 prevalence ratios between paired datasets, and synthesizes them via random-effects meta-analysis at multiple aggregation levels (year, age group, sex).

Author(s)

Maintainer: Maarja Pajusalu maarja.pajusalu@ut.ee

Authors:

See Also

Useful links:

Build the SELECT SQL for care-site cohort generation.

Description

Build the SELECT SQL for care-site cohort generation.

Usage

.build_caresite_sql(
  care_site_id,
  cohort_id,
  cdm_schema,
  collapse_strategy,
  restrict_to_observation
)

Age decade clamping: all ages >= this value merge into a single "80+" group.

Description

Age decade clamping: all ages >= this value merge into a single "80+" group.

Usage

AGE_CLAMP_MAX

ATC vocabulary constants for chapter lookup.

Description

ATC vocabulary constants for chapter lookup.

Usage

ATC_CHAPTER_VOCAB

Anchor concept_ids for condition chapter assignment.

Description

Anchor concept_ids for condition chapter assignment.

Usage

CHAPTER_ROOTS

SNOMED relationship types to extract as condition attributes.

Description

SNOMED relationship types to extract as condition attributes.

Usage

CONDITION_RELATIONSHIPS

Drug attribute relationships (empty - RxNorm relationships are structural).

Description

Drug attribute relationships (empty - RxNorm relationships are structural).

Usage

DRUG_RELATIONSHIPS

Gender concept_id to label mapping.

Description

Gender concept_id to label mapping.

Usage

GENDER_LABELS

k-anonymity threshold. Any stratum cell with fewer patients is suppressed.

Description

k-anonymity threshold. Any stratum cell with fewer patients is suppressed.

Usage

K_ANONYMITY

Anchor concept_ids for procedure chapter assignment.

Description

Anchor concept_ids for procedure chapter assignment.

Usage

PROCEDURE_CHAPTER_ROOTS

SNOMED relationship types to extract as procedure attributes.

Description

SNOMED relationship types to extract as procedure attributes.

Usage

PROCEDURE_RELATIONSHIPS

Base directory for extracted source data.

Description

Base directory for extracted source data.

Usage

SOURCES_DIR

Add prevalence ratio columns (log2 + natural + inference) to a data frame.

Description

Expects columns: prevalence_d1, prevalence_d2, denominator_d1, denominator_d2.

Usage

add_pr_columns(df)

Arguments

df

Data frame with prevalence and denominator columns.

Details

Formulas: log2_pr = log2(p2 / p1) SE(ln) = sqrt((1-p1)/(p1*n1) + (1-p2)/(p2*n2)) SE(log2)= SE(ln) / ln(2) CI = log2_pr +/- 1.96 x SE

Value

Data frame with PR columns appended.

Apply OHDSI cohort filtering to CDM table references.

Description

Modifies db$cdm table references to filter by a standard OHDSI cohort table. Filters person, observation_period, death, and event tables to cohort members and their cohort windows.

Usage

apply_cohort_filter(db, cohort_id, cohort_schema = NULL)

Arguments

db

Connection list (from syrona_connect).

cohort_id

Integer cohort_definition_id to filter by.

cohort_schema

Schema containing the cohort table (NULL = default).

Value

Modified db list with filtered CDM table references.

Apply k-anonymity suppression to all extracted source tables.

Description

Apply k-anonymity suppression to all extracted source tables.

Usage

apply_k_anonymity(tables, k = K_ANONYMITY)

Arguments

tables

Named list of tibbles (before k-anonymity).

k

Minimum cell count (default K_ANONYMITY = 5).

Value

Named list of tibbles with k-anonymity applied.

Build a per-year forest plot for a single concept

Description

Shows yearly points + CIs, meta-analysis diamond, and colored heatmap tile row. Faceted by sex x age group.

Usage

build_forest_detail(
  yearly_data,
  meta_data,
  concept_name,
  concept_code,
  name1,
  name2,
  scale_mode = "log2"
)

Arguments

yearly_data

Yearly comparison data for one concept.

meta_data

Meta-analysis data (across years) for one concept.

concept_name

Concept name string.

concept_code

Concept code string.

name1

Reference dataset name.

name2

Comparison dataset name.

scale_mode

"fold" or "log2".

Value

A girafe object.

Build an overview forest plot paired with the heatmap

Description

Compact forest showing Female, Male, Both meta PR for every concept. Concept y-axis ordering matches the heatmap.

Usage

build_forest_overview(
  meta_sex_df,
  meta_summary_df,
  concept_ids,
  concept_info,
  name1,
  name2,
  scale_mode = "log2",
  concept_order = NULL
)

Arguments

meta_sex_df

Meta-analysis by sex data (F/M rows).

meta_summary_df

Meta-analysis summary data (Both rows).

concept_ids

Vector of concept_ids to include.

concept_info

Data frame with concept_id, concept_name, concept_code, pop_weight.

name1

Reference dataset name.

name2

Comparison dataset name.

scale_mode

"fold" or "log2".

concept_order

Vector of concept_ids in desired display order.

Value

A girafe object.

Build a compact F/M/Both forest plot for a single concept

Description

Build a compact F/M/Both forest plot for a single concept

Usage

build_forest_summary(
  meta_sex_data,
  meta_summary_data,
  concept_name,
  concept_code,
  name1,
  name2,
  scale_mode = "log2"
)

Arguments

meta_sex_data

Meta-analysis data by sex (F, M) for one concept.

meta_summary_data

Meta-analysis summary (Both) for one concept.

concept_name

Concept name string.

concept_code

Concept code string.

name1

Reference dataset name.

name2

Comparison dataset name.

scale_mode

"fold" or "log2".

Value

A girafe object.

Build an interactive heatmap of prevalence ratios

Description

Build an interactive heatmap of prevalence ratios

Usage

build_heatmap(
  data,
  name1,
  name2,
  scale_mode = "log2",
  concept_weights = NULL,
  concept_order = NULL
)

Arguments

data

Data frame with meta_agegroups data (concept_id, concept_name, concept_code, sex, age_group, log2_pr, ci_low, ci_high, p_value).

name1

Reference dataset name.

name2

Comparison dataset name.

scale_mode

"fold" or "log2" for axis labels.

concept_weights

Data frame with concept_id, concept_name, pop_weight.

concept_order

Vector of concept_ids in desired display order.

Value

A girafe object.

Build a one-row tibble from run_meta output + aggregated source columns.

Description

Build a one-row tibble from run_meta output + aggregated source columns.

Usage

build_meta_row(meta_result, id_cols, g)

Arguments

meta_result

Named list from run_meta.

id_cols

One-row tibble with group ID columns.

g

Group data frame (for aggregated counts).

Value

One-row tibble.

Build a dumbbell chart showing absolute prevalence values

Description

Complements the forest detail (relative -> absolute). Same facet structure: sex x age_group.

Usage

build_point_diff(yearly_data, concept_name, concept_code, name1, name2)

Arguments

yearly_data

Yearly comparison data for one concept.

concept_name

Concept name string.

concept_code

Concept code string.

name1

Reference dataset name.

name2

Comparison dataset name.

Value

A girafe object.

Build a PR distribution histogram from meta_summary data

Description

Build a PR distribution histogram from meta_summary data

Usage

build_pr_distribution(
  data,
  d1_name,
  d2_name,
  scale_mode = "fold",
  fold_thresh = FOLD_THRESHOLD,
  domain_label = "conditions"
)

Arguments

data

Filtered meta_summary rows (one per concept). Expected columns: log2_pr, fold_diff, sig, concept_id.

d1_name

Name of reference dataset (shown on left side).

d2_name

Name of comparison dataset (shown on right side).

scale_mode

"fold" or "log2" - controls x-axis labels only.

fold_thresh

Fold difference threshold for zone coloring.

domain_label

Label for domain (e.g. "conditions", "drugs").

Value

A ggplot object.

Build a faceted PR distribution - one small histogram per chapter

Description

Interactive (ggiraph): clicking a chapter facet returns its chapter_id.

Usage

build_pr_distribution_chapters(
  data,
  d1_name,
  d2_name,
  scale_mode = "fold",
  fold_thresh = FOLD_THRESHOLD,
  domain_label = "conditions"
)

Arguments

data

Filtered meta_summary joined with chapters (one row per concept x chapter). Expected columns: log2_pr, fold_diff, sig, concept_id, chapter_name, chapter_id.

d1_name

Name of reference dataset.

d2_name

Name of comparison dataset.

scale_mode

"fold" or "log2".

fold_thresh

Fold difference threshold.

domain_label

Label for domain.

Value

A girafe object. Selection data_id = chapter_id (character).

Build an interactive population pyramid

Description

Build an interactive population pyramid

Usage

build_pyramid(demo_df, dataset_name)

Arguments

demo_df

Demographics data frame with birth_year, sex, patient_count.

dataset_name

Name of the dataset (used as title).

Value

A girafe object.

Get summary statistics for a cohort.

Description

Get summary statistics for a cohort.

Usage

cohort_summary(con, cohort_id, cohort_schema = NULL)

Arguments

con

DBI connection.

cohort_id

Integer cohort_definition_id.

cohort_schema

Schema containing the cohort table (NULL = default).

Value

A tibble with n_entries, n_persons, min_start, max_end.

Run the full comparison pipeline for two datasets across available domains.

Description

Loads both datasets, runs the 4-step comparison for each domain present in both, and optionally saves to CSV.

Usage

compare_all(
  d1_name,
  d2_name,
  domains = c("conditions", "procedures", "drugs"),
  save = TRUE
)

Arguments

d1_name

Name of dataset 1 (reference).

d2_name

Name of dataset 2 (comparison).

domains

Character vector of domains to compare.

save

If TRUE (default), writes CSV output.

Value

Named list of domain results (invisible).

Examples


# Run the bundled demo (a curated 32-concept subset) end to end.
# Copy it to a writable folder first, since the installed copy is read-only.
base <- tempdir()
file.copy(system.file("extdata", "demo", package = "syrona"),
          base, recursive = TRUE)
dir <- file.path(base, "demo")
old <- options(syrona.data_dir = dir)
# One domain keeps the example quick; drop `domains` to compare all three.
res <- compare_all("demo_population", "demo_selected", domains = "conditions")
res$condition_meta_summary
options(old)
# Then explore interactively: run_app(data_dir = dir)

Run the full comparison pipeline for one domain.

Description

Run the full comparison pipeline for one domain.

Usage

compare_domain(d1, d2, prev_table, domain_label)

Arguments

d1

Loaded dataset 1.

d2

Loaded dataset 2.

prev_table

Name of the prevalence table.

domain_label

Label for messages.

Value

Named list of 4 tibbles, or NULL if no overlapping concepts.

Meta-analysis across years: one row per concept_id x sex x age_group.

Description

Meta-analysis across years: one row per concept_id x sex x age_group.

Usage

compare_meta_agegroups(yearly_df)

Arguments

yearly_df

Output of compare_yearly.

Value

tibble with meta-analyzed rows.

Meta-analysis across age groups: one row per concept_id x sex.

Description

Meta-analysis across age groups: one row per concept_id x sex.

Usage

compare_meta_by_sex(meta_agegroups_df)

Arguments

meta_agegroups_df

Output of compare_meta_agegroups.

Value

tibble with meta-analyzed rows.

Meta-analysis across sexes: one row per concept_id (sex = "Both").

Description

Meta-analysis across sexes: one row per concept_id (sex = "Both").

Usage

compare_meta_summary(meta_by_sex_df)

Arguments

meta_by_sex_df

Output of compare_meta_by_sex.

Value

tibble with one summary row per concept.

Compare two datasets at the per-year stratum level.

Description

Matches on concept_id x year x sex x age_group (inner join). Only keeps strata where both datasets have prevalence > 0.

Usage

compare_yearly(d1, d2, prev_table = "condition_prevalence")

Arguments

d1

Named list from load_dataset (reference).

d2

Named list from load_dataset (comparison).

prev_table

Name of the prevalence table to compare.

Value

tibble with one row per matched stratum, full PR column set.

Create an OHDSI cohort from a care site.

Description

For each person who visited the specified care site, their cohort window runs from their first visit start date to their last visit end date at that site. Optionally clips to observation period overlap.

Usage

create_caresite_cohort(
  con,
  care_site_id,
  cohort_id,
  cohort_schema = NULL,
  cdm_schema,
  overwrite = FALSE,
  restrict_to_observation = TRUE,
  collapse_strategy = c("person_span", "visit_occurrence")
)

Arguments

con

DBI connection to the OMOP CDM database.

care_site_id

Integer care_site_id to define the cohort.

cohort_id

Integer cohort_definition_id to assign.

cohort_schema

Schema for the cohort table (NULL = default).

cdm_schema

Schema containing OMOP CDM tables.

overwrite

If TRUE, deletes existing entries for this cohort_id before inserting.

restrict_to_observation

If TRUE (default), clips cohort intervals to observation period overlap.

collapse_strategy

How to collapse visits per person: "person_span" (default) = first visit start to last visit end; "visit_occurrence" = one cohort entry per visit.

Details

This uses server-side INSERT...SELECT for performance on large databases (avoids downloading + re-uploading patient-level data).

Value

Number of rows inserted (invisible).

Create an OHDSI-standard cohort table.

Description

Creates the 4-column cohort table if it does not already exist. Works with both DuckDB and PostgreSQL.

Usage

create_cohort_table(con, cohort_schema = NULL)

Arguments

con

DBI connection.

cohort_schema

Schema for the cohort table. NULL for default.

Value

Invisible TRUE if created, FALSE if already existed.

Delete a cohort from the cohort table.

Description

Delete a cohort from the cohort table.

Usage

delete_cohort(con, cohort_id, cohort_schema = NULL)

Arguments

con

DBI connection.

cohort_id

Integer cohort_definition_id to delete.

cohort_schema

Schema containing the cohort table (NULL = default).

Value

Number of rows deleted (invisible).

Run the full source data extraction pipeline.

Description

Derives stratified prevalence tables from the records in an OMOP CDM database for one or more clinical domains. Applies k-anonymity suppression and optionally saves results to CSV.

Usage

extract_all(
  dataset_name,
  db,
  domains = c("conditions", "procedures", "drugs"),
  cohort_id = NULL,
  cohort_schema = NULL,
  save = TRUE
)

Arguments

dataset_name

Short label for the dataset (e.g. "EH30", "EstBB").

db

Either a DuckDB file path (string) or an existing connection list from syrona_connect / syrona_connect_pg.

domains

Character vector of domains to extract. Options: "conditions", "procedures", "drugs".

cohort_id

Integer cohort_definition_id to filter by. If NULL (default), extracts the full dataset.

cohort_schema

Schema containing the cohort table.

save

If TRUE (default), saves CSV to data/sources/<dataset_name>/.

Value

Named list of tibbles matching the Syrona schema (invisible).

Extract condition SNOMED attributes.

Description

Extract condition SNOMED attributes.

Usage

extract_condition_attributes(cdm, condition_ids)

Arguments

cdm

CDM reference object.

condition_ids

Integer vector of concept_ids.

Value

tibble with columns: concept_id, relationship, target_concept_id, target_concept_name

Extract condition chapter and sub-chapter assignments.

Description

Extract condition chapter and sub-chapter assignments.

Usage

extract_condition_chapters(cdm, condition_ids)

Arguments

cdm

CDM reference object.

condition_ids

Integer vector of concept_ids to assign chapters to.

Value

tibble with columns: concept_id, chapter_type, chapter_id, chapter_name, chapter_level, parent_chapter_id

Extract condition concept metadata.

Description

Extract condition concept metadata.

Usage

extract_condition_info(cdm)

Arguments

cdm

CDM reference object.

Value

tibble with columns: concept_id, concept_name, concept_code, vocabulary_id, concept_class_id, n_patients_total, n_records_total

Extract condition prevalence (ACHILLES-404 numerator + 116 denominator).

Description

Extract condition prevalence (ACHILLES-404 numerator + 116 denominator).

Usage

extract_condition_prevalence(cdm, denom_df)

Arguments

cdm

CDM reference object.

denom_df

Pre-computed denominator tibble.

Value

tibble with columns: concept_id, year, sex, age_group, patient_count, record_count, denominator, prevalence

Extract death counts by year x sex x age group (ACHILLES-504 pattern).

Description

Extract death counts by year x sex x age group (ACHILLES-504 pattern).

Usage

extract_death_counts(cdm, denom_df)

Arguments

cdm

CDM reference object.

denom_df

Pre-computed denominator tibble.

Value

tibble with columns: year, sex, age_group, death_count, denominator, mortality_rate

Extract demographics (birth year x sex).

Description

Extract demographics (birth year x sex).

Usage

extract_demographics(cdm)

Arguments

cdm

CDM reference object.

Value

tibble with columns: sex, birth_year, patient_count

Extract the ACHILLES-116 denominator: persons observed per year x sex x age group.

Description

Extract the ACHILLES-116 denominator: persons observed per year x sex x age group.

Usage

extract_denominators(cdm)

Arguments

cdm

CDM reference object.

Value

tibble with columns: year, sex, age_group, denominator

Extract drug chapter assignments (ATC 1st level).

Description

Extract drug chapter assignments (ATC 1st level).

Usage

extract_drug_chapters(cdm, drug_ids)

Arguments

cdm

CDM reference object.

drug_ids

Integer vector of Ingredient concept_ids.

Value

tibble with same columns as extract_condition_chapters.

Extract drug (ingredient) concept metadata.

Description

Extract drug (ingredient) concept metadata.

Usage

extract_drug_info(cdm)

Arguments

cdm

CDM reference object.

Value

tibble with same columns as extract_condition_info.

Extract drug prevalence at the Ingredient level.

Description

Extract drug prevalence at the Ingredient level.

Usage

extract_drug_prevalence(cdm, denom_df)

Arguments

cdm

CDM reference object.

denom_df

Pre-computed denominator tibble.

Value

tibble with same columns as extract_condition_prevalence.

Extract procedure SNOMED attributes.

Description

Extract procedure SNOMED attributes.

Usage

extract_procedure_attributes(cdm, procedure_ids)

Arguments

cdm

CDM reference object.

procedure_ids

Integer vector of concept_ids.

Value

tibble with same columns as extract_condition_attributes.

Extract procedure chapter and sub-chapter assignments.

Description

Extract procedure chapter and sub-chapter assignments.

Usage

extract_procedure_chapters(cdm, procedure_ids)

Arguments

cdm

CDM reference object.

procedure_ids

Integer vector of concept_ids.

Value

tibble with same columns as extract_condition_chapters.

Extract procedure concept metadata.

Description

Extract procedure concept metadata.

Usage

extract_procedure_info(cdm)

Arguments

cdm

CDM reference object.

Value

tibble with same columns as extract_condition_info.

Extract procedure prevalence.

Description

Extract procedure prevalence.

Usage

extract_procedure_prevalence(cdm, denom_df)

Arguments

cdm

CDM reference object.

denom_df

Pre-computed denominator tibble.

Value

tibble with columns: concept_id, year, sex, age_group, patient_count, record_count, denominator, prevalence

Format a log2 prevalence ratio as a fold difference string

Description

Vectorized. For values < 1, shows reciprocal in parentheses.

Usage

format_fold(log2_val)

Arguments

log2_val

Numeric vector of log2 PR values.

Value

Character vector of formatted fold strings.

Format p-values in Nature style

Description

Vectorized. Shows < 0.001 for very small values.

Usage

format_pval(p)

Arguments

p

Numeric vector of p-values.

Value

Character vector of formatted p-value strings.

Insert a cohort from a local data frame.

Description

Uploads a data frame to the database via omopgenerics::insertTable and marks it as a cohort table with omopgenerics::newCohortTable. Use this when cohort membership has already been computed locally (e.g. from a CSV or programmatic cohort definition).

Usage

insert_cohort(cdm, cohort_df, name = "cohort")

Arguments

cdm

A CDM reference (from syrona_connect()$cdm).

cohort_df

Data frame with columns: cohort_definition_id, subject_id, cohort_start_date, cohort_end_date.

name

Name for the cohort table in the database (default "cohort").

Details

Requires that the CDM connection was created with a writeSchema.

Value

Updated CDM reference with the cohort table attached.

List care sites with patient counts.

Description

Useful for identifying care_site_id values before creating cohorts.

Usage

list_care_sites(con, cdm_schema, min_patients = 100)

Arguments

con

DBI connection.

cdm_schema

Schema containing the OMOP CDM tables.

min_patients

Minimum number of patients to include (default 100).

Value

A tibble with care_site_id, care_site_name, n_patients.

List available comparisons.

Description

List available comparisons.

Usage

list_comparisons()

Value

Character vector of comparison directory names.

List all available datasets.

Description

List all available datasets.

Usage

list_datasets()

Value

Character vector of dataset names.

Load a previously saved comparison from CSV.

Description

Load a previously saved comparison from CSV.

Usage

load_comparison(d1_name, d2_name)

Arguments

d1_name

Dataset 1 name.

d2_name

Dataset 2 name.

Value

Named list of tibbles.

Load a previously saved dataset from CSV.

Description

Load a previously saved dataset from CSV.

Usage

load_dataset(dataset_name)

Arguments

dataset_name

Short label (subfolder name under data/sources/).

Value

Named list of tibbles.

Examples


base <- tempdir()
file.copy(system.file("extdata", "demo", package = "syrona"), base, recursive = TRUE)
old <- options(syrona.data_dir = file.path(base, "demo"))
d <- load_dataset("demo_population")
names(d)
options(old)

Load all Syrona design tokens into the calling environment

Description

Makes color constants, thresholds, and other design tokens available as variables. Used by the Shiny dashboard's global.R.

Usage

load_syrona_theme(envir = parent.frame())

Arguments

envir

Environment to load into (default: caller's environment).

Value

No return value, called for its side effect of assigning the design tokens as variables in envir.

Run the Syrona dashboard

Description

Launches the Shiny dashboard for exploring prevalence comparisons. The app looks for data/sources/ and data/comparisons/ relative to your current working directory. Run this from the directory that contains your data/ folder.

Usage

run_app(data_dir = getwd(), port = NULL, launch.browser = TRUE, ...)

Arguments

data_dir

Path to the directory containing data/sources/ and data/comparisons/. Defaults to the current working directory.

port

Port to run the app on (default: auto-select).

launch.browser

Whether to open a browser window (default: TRUE).

...

Additional arguments passed to runApp.

Value

No return value, called for its side effect of launching the Shiny dashboard.

Examples

# Copy the bundled demo, generate the comparison, then launch the dashboard.
if (interactive()) {
  base <- tempdir()
  file.copy(system.file("extdata", "demo", package = "syrona"), base, recursive = TRUE)
  dir <- file.path(base, "demo")
  old <- options(syrona.data_dir = dir)
  compare_all("demo_population", "demo_selected")
  run_app(data_dir = dir)
  options(old)
}

Run meta-analysis on a vector of log2 prevalence ratios.

Description

Strategy: 1. Single study: Pass-Through (return as-is) 2. Try random-effects with Paule-Mandel tau 3. If RE fails: try fixed-effect 4. If both fail: return NULL

Usage

run_meta(te, se_te, studlab)

Arguments

te

Numeric vector of log2 prevalence ratios.

se_te

Numeric vector of standard errors (log2 scale).

studlab

Character/numeric vector of study labels.

Value

Named list with log2_pr, se, ci_low, ci_high, z, p_value, meta_model_type, tau2, I2, Q, pval_Q, lower_predict, upper_predict. NULL if meta-analysis fails.

Save comparison tables to CSV.

Description

Save comparison tables to CSV.

Usage

save_comparison(tables, d1_name, d2_name)

Arguments

tables

Named list of tibbles.

d1_name

Dataset 1 name.

d2_name

Dataset 2 name.

Save extracted tables to CSV.

Description

Save extracted tables to CSV.

Usage

save_dataset(tables, dataset_name, db_path = NA_character_)

Arguments

tables

Named list of tibbles.

dataset_name

Short label for the dataset.

db_path

Database path (stored in metadata).

Apply k-anonymity suppression to a single domain's tables.

Description

Apply k-anonymity suppression to a single domain's tables.

Usage

suppress_domain(
  info,
  prevalence,
  chapters,
  attributes,
  domain_label = "concept",
  k = K_ANONYMITY
)

Arguments

info

Info tibble (must have concept_id, n_patients_total, n_records_total).

prevalence

Prevalence tibble (must have concept_id, patient_count).

chapters

Chapters tibble (can be NULL).

attributes

Attributes tibble (can be NULL).

domain_label

Short label for messages.

k

Minimum cell count.

Value

Named list: prevalence, info, chapters, attributes, rare.

Connect to a DuckDB OMOP CDM.

Description

Opens a DuckDB file and creates a CDMConnector reference. By default the connection is read-only; set read_only = FALSE if you need to write cohort tables or temp tables into the database.

Usage

syrona_connect(
  db_path,
  cdm_schema = "main",
  write_schema = cdm_schema,
  read_only = TRUE
)

Arguments

db_path

Path to the DuckDB file.

cdm_schema

Schema containing OMOP CDM tables (default "main").

write_schema

Schema for writing temp/cohort tables. Defaults to cdm_schema. Set to a different schema if the CDM schema is read-only (common in production setups).

read_only

Logical. Open DuckDB in read-only mode? Default TRUE. Set to FALSE when you need to create cohort tables or use omopgenerics::insertTable.

Value

A list with components:

con

DBI connection object.

cdm

CDMConnector CDM reference (cdm_reference).

Connect to a PostgreSQL OMOP CDM.

Description

For remote databases, typically accessed via SSH tunnel: ssh -L 5432:localhost:5432 user@server.

Usage

syrona_connect_pg(
  host = "localhost",
  port = 5432,
  dbname,
  user,
  password = NULL,
  cdm_schema = "public",
  write_schema = cdm_schema
)

Arguments

host

Hostname (default "localhost" for SSH tunnel).

port

Port number (default 5432).

dbname

Database name.

user

PostgreSQL username.

password

Password. If NULL, uses .pgpass or PGPASSWORD env var.

cdm_schema

Schema with OMOP CDM tables (e.g. "cdm").

write_schema

Schema for cohort/temp tables (e.g. "results_maarja").

Value

Same structure as syrona_connect.

Disconnect from an OMOP CDM database.

Description

Disconnect from an OMOP CDM database.

Usage

syrona_disconnect(db)

Arguments

db

Connection list returned by syrona_connect or syrona_connect_pg.

Value

No return value, called for its side effect of closing the database connection.

Truncate a concept name to a maximum number of characters

Description

Adds "..." if truncated.

Usage

truncate_name(name, max = CONCEPT_NAME_MAX)

Arguments

name

Character vector of names.

max

Maximum character length (default CONCEPT_NAME_MAX).

Value

Character vector of (possibly truncated) names.

Map population weight (%) to a hex color on a log-scale gradient

Description

Light gray for rare concepts, dark charcoal for common ones.

Usage

weight_color(pct)

Arguments

pct

Numeric vector of population weight percentages (0-100).

Value

Character vector of hex color strings.