inst/script_example/script_example1.R

In DIETCOST: Calculate the Cost and Environmental Impact of a Ideal Diet

#Global parameters--------------------------------------------------------------
redmeat_ids <- c(71003, 71008, 71041, 81005, 81021, 81022, 81026, 81027, 81029)
filepath <- 'home/dir/dataset.xlsx'
linked_bread_low <- c("79034", "79015", "79016", "79032", "65048")
linked_bread_high <- c("80001", "80019", "81029", "81027", "81021", "70044", "69009")
linked_milk_low <- c("69029", "79012", "70038", "69043")
linked_milk_high <- c("79020", "79001", "79036", "79003", "79039")

#Initiates analysis-------------------------------------------------------------
foods_df <- DIETCOST::foodData(filepath = filepath, redmeat_ids = redmeat_ids, diets = c('C', 'PF', 'H'), max_scale = 2)
food_groups_df <- DIETCOST::foodGroupData(filepath = filepath, df_foods = foods_df, diets = c('C', 'PF', 'H'))
nutrient_targets <- DIETCOST::createNutrientTargets(filepath = filepath, allow_takeaway = FALSE, alcohol_perc_max = 20, discretionary_perc_max = 100)

#Runs a Monte Carlo Simulation--------------------------------------------------
DIETCOST::monteCarloSimulation('C:/Users/username/directory', 100000, foods_df, nutrient_targets, food_groups_df, person = 'woman', 'PF', c(1,2,3), 0.5, allow_discretionary = TRUE, allow_alcohol = TRUE, allow_takeaway = TRUE, linked_low_1 = linked_bread_low, linked_high_1 = linked_bread_high, linked_low_2 = linked_milk_low, linked_high_2 = linked_milk_high)

#Exports results----------------------------------------------------------------
DIETCOST::calculateResults('folder/results', 'C:/Users/username/directory', 0.05)
DIETCOST::calculateGroupedResults('folder/results', 'C:/Users/username/directory', 0.05)