# Example usage on toy metabolomics data:
data <- data.frame(
batch = rep(c("A","B","C"), each = 10),
injectionOrder = rep(1:30, times = 1),
metabolite1 = rnorm(30, mean = 100, sd = 10),
metabolite2 = rnorm(30, mean = 200, sd = 20)
)
result <- convex_analysis_of_variables(
data = data,
variable_columns = c("metabolite1", "metabolite2"),
batch_col = "batch",
sample_order_col = "injectionOrder",
impute_if_needed = "median",
mode = "global"
)
plot_all_convex_hulls(
target_file_path = file.path(tempdir(), "convex_hulls.pdf"),
convex_analysis_res = result,
show_points = TRUE,
mode = "global"
)
Run the code above in your browser using DataLab