## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")
library(MineSDG)

## -----------------------------------------------------------------------------
head(demo_mine_sites[, 1:8])

## -----------------------------------------------------------------------------
site <- demo_mine_sites[demo_mine_sites$site_id == "CU-ATAC" &
                          demo_mine_sites$year == 2024, ]

## -----------------------------------------------------------------------------
calculate_ghg_intensity(
  scope1_t = site$ghg_scope1_t,
  scope2_t = site$ghg_scope2_t,
  ore_processed_kt = site$ore_processed_kt
)

## -----------------------------------------------------------------------------
calculate_energy_intensity(
  energy_gj = site$energy_gj,
  ore_processed_kt = site$ore_processed_kt,
  renewable_gj = site$energy_gj * site$renewable_energy_pct / 100
)

## -----------------------------------------------------------------------------
calculate_safety_performance(
  hours_worked = site$hours_worked,
  recordable_injuries = site$recordable_injuries,
  lost_time_injuries = site$lost_time_injuries,
  fatalities = site$fatalities
)

## -----------------------------------------------------------------------------
calculate_water_efficiency(site$water_withdrawal_m3,
                           site$water_discharge_m3,
                           site$water_recycled_m3)

calculate_land_restoration(site$land_disturbed_ha,
                           site$land_rehabilitated_ha)

calculate_community_investment(site$community_investment_musd,
                               site$revenue_musd)

calculate_waste_intensity(site$ore_processed_kt,
                          site$tailings_kt,
                          waste_rock_kt = site$waste_rock_kt)

## -----------------------------------------------------------------------------
my_site <- data.frame(
  site_id = "MY-MINE", year = 2025,
  ore_processed_kt = 12000,
  ghg_scope1_t = 420000, ghg_scope2_t = 180000,
  hours_worked = 5.2e6, recordable_injuries = 18,
  lost_time_injuries = 6, fatalities = 0
)

calculate_ghg_intensity(my_site$ghg_scope1_t, my_site$ghg_scope2_t,
                        my_site$ore_processed_kt)$ghg_intensity