inst/doc/MedxR_vignette.R

In MedxR: Access Drug Regulatory Data via FDA and Health Canada APIs

## ----include = FALSE--------------------------------------------------------------------------------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  echo = TRUE,
  message = FALSE,
  warning = FALSE
)

# Save original user options
old_options <- options(
  width = 150,
  tibble.width = Inf,
  pillar.width = Inf
)


## ----setup------------------------------------------------------------------------------------------------------------------------------------------
library(MedxR)
library(ggplot2)
library(dplyr)

## ----fda-adverse-events, eval=FALSE-----------------------------------------------------------------------------------------------------------------
# aspirin_adverse_event <- get_fda_adverse_events("aspirin")
# aspirin_adverse_event %>%
#   select(report_id, date_received, country, serious, patient_sex) %>%
#   head(6)

## ----fda-approved-drugs, eval=FALSE-----------------------------------------------------------------------------------------------------------------
# lipitor_drugs_approved <- get_fda_drugs_approved("lipitor")
# lipitor_drugs_approved %>%
#   select(application_number, brand, generic, approval_date, strength, form)

## ----health-canada-din, eval=FALSE------------------------------------------------------------------------------------------------------------------
# hc_drug_din <- get_hc_drug_by_din("00000213")
# hc_drug_din %>%
#   select(drug_code, class_name, din, brand_name, company_name)

## ----health-canada-products-companies, eval=FALSE, message=FALSE, warning=FALSE, fig.width=10, fig.height=5-----------------------------------------
# # Retrieve data from Health Canada APIs
# drug_products <- get_hc_drug_products()
# companies <- get_hc_companies()
# 
# # Combine both datasets by company name
# merged_data <- drug_products %>%
#   left_join(companies, by = "company_name") %>%
#   filter(!is.na(company_name)) %>%
#   group_by(company_name) %>%
#   summarise(total_products = n()) %>%
#   arrange(desc(total_products)) %>%
#   slice_head(n = 10)
# 
# # Plot: Top 10 companies by number of registered drug products
# ggplot(merged_data, aes(x = reorder(company_name, total_products), y = total_products)) +
#   geom_col(fill = "steelblue") +
#   coord_flip() +
#   labs(
#     title = "Top 10 Canadian Companies by Number of Registered Drug Products",
#     x = "Company Name",
#     y = "Number of Drug Products",
#     caption = "Source: Health Canada Drug Product Database (via MedxR)"
#   ) +
#   theme_minimal(base_size = 13)

## ----cleanup, include = FALSE-------------------------------------------------
# Restore original user options

options(old_options)