Hi, I am trying to run ANCOM-BC but the following error (below) keeps popping up. The error says that I cannot load "'/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/gsl/libs/gsl.dylib':". I'm confused because when I am able to manually find this file in my folders following the same path. I have tried downloading and reinstalling miniconda3 and my qiime environment.
Update 6/16: just out of curiosity I tried running it in a qiime2 v2024.2 instead of v2025.2 and the same error showed up.
/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/q2_composition/_ancombc.py:77: FutureWarning: Series.__getitem__ treating keys as positions is deprecated. In a future version, integer keys will always be treated as labels (consistent with DataFrame behavior). To access a value by position, use `ser.iloc[pos]`
term_alpha_value = (metadata.get_column(term)
Running external command line application(s). This may print messages to stdout and/or stderr.
The command(s) being run are below. These commands cannot be manually re-run as they will depend on temporary files that no longer exist.
Command: run_ancombc.R --inp_abundances_path /var/folders/rt/6vtjnzqd64qbxyg4_fpkknjm0000gn/T/tmpc6d3msb7/input.biom.tsv --inp_metadata_path /var/folders/rt/6vtjnzqd64qbxyg4_fpkknjm0000gn/T/tmpc6d3msb7/input.map.txt --md_column_types {"Plate": "categorical", "Name": "categorical", "Infected": "categorical", "Oil": "categorical", "Trauma": "categorical", "CollectionDate": "categorical", "MultipleTimePoints": "categorical", "RemoveDuplicates": "categorical", "CollectionMonth": "numeric", "CollectionDay": "numeric", "CaptureDate": "categorical", "CaptureMonth": "numeric", "CaptureDay": "numeric", "County": "categorical", "City": "categorical", "Latitude": "numeric", "Longitude": "numeric", "Sex": "categorical", "WeightKg": "numeric", "Length": "numeric", "Status": "categorical", "Blank": "categorical", "Duplicate": "categorical", "SampleType": "categorical", "Reads": "numeric"} --formula Sex --p_adj_method holm --prv_cut 0.1 --lib_cut 0 --reference_levels ['Sex::Female'] --tol 1e-05 --max_iter 100 --conserve False --alpha 0.05 --output_loaf /var/folders/rt/6vtjnzqd64qbxyg4_fpkknjm0000gn/T/q2-DataLoafPackageDirFmt-ga3vswl8
── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr 1.1.4 ✔ readr 2.1.5
✔ forcats 1.0.0 ✔ stringr 1.5.1
✔ ggplot2 3.5.1 ✔ tibble 3.2.1
✔ lubridate 1.9.3 ✔ tidyr 1.3.1
✔ purrr 1.0.2
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag() masks stats::lag()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
Error: package or namespace load failed for ‘ANCOMBC’ in dyn.load(file, DLLpath = DLLpath, ...):
unable to load shared object '/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/gsl/libs/gsl.dylib':
dlopen(/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/gsl/libs/gsl.dylib, 0x0006): Library not loaded: @rpath/libgsl.25.dylib
Referenced from: <5451EAF0-839D-3B5C-80E8-426E99B788AD> /Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/gsl/libs/gsl.dylib
Reason: tried: '/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/gsl/libs/../../../lib/libgsl.25.dylib' (no such file), '/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/gsl/libs/../../../../libgsl.25.dylib' (no such file), '/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/gsl/libs/../../../lib/libgsl.25.dylib' (no such file), '/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/R/library/
8: stop(msg, call. = FALSE, domain = NA)
7: value[[3L]](cond)
6: tryCatchOne(expr, names, parentenv, handlers[[1L]])
5: tryCatchList(expr, classes, parentenv, handlers)
4: tryCatch({
attr(package, "LibPath") <- which.lib.loc
ns <- loadNamespace(package, lib.loc)
env <- attachNamespace(ns, pos = pos, deps, exclude, include.only)
}, error = function(e) {
P <- if (!is.null(cc <- conditionCall(e)))
paste(" in", deparse(cc)[1L])
else ""
msg <- gettextf("package or namespace load failed for %s%s:\n %s",
sQuote(package), P, conditionMessage(e))
if (logical.return && !quietly)
message(paste("Error:", msg), domain = NA)
else stop(msg, call. = FALSE, domain = NA)
})
3: library(ANCOMBC)
2: withCallingHandlers(expr, warning = function(w) if (inherits(w,
classes)) tryInvokeRestart("muffleWarning"))
1: suppressWarnings(library(ANCOMBC))
Traceback (most recent call last):
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/q2_composition/_ancombc.py", line 255, in _ancombc
run_commands([cmd])
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/q2_composition/_ancombc.py", line 32, in run_commands
subprocess.run(cmd, check=True)
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/subprocess.py", line 528, in run
raise CalledProcessError(retcode, process.args,
subprocess.CalledProcessError: Command '['run_ancombc.R', '--inp_abundances_path', '/var/folders/rt/6vtjnzqd64qbxyg4_fpkknjm0000gn/T/tmpc6d3msb7/input.biom.tsv', '--inp_metadata_path', '/var/folders/rt/6vtjnzqd64qbxyg4_fpkknjm0000gn/T/tmpc6d3msb7/input.map.txt', '--md_column_types', '{"Plate": "categorical", "Name": "categorical", "Infected": "categorical", "Oil": "categorical", "Trauma": "categorical", "CollectionDate": "categorical", "MultipleTimePoints": "categorical", "RemoveDuplicates": "categorical", "CollectionMonth": "numeric", "CollectionDay": "numeric", "CaptureDate": "categorical", "CaptureMonth": "numeric", "CaptureDay": "numeric", "County": "categorical", "City": "categorical", "Latitude": "numeric", "Longitude": "numeric", "Sex": "categorical", "WeightKg": "numeric", "Length": "numeric", "Status": "categorical", "Blank": "categorical", "Duplicate": "categorical", "SampleType": "categorical", "Reads": "numeric"}', '--formula', 'Sex', '--p_adj_method', 'holm', '--prv_cut', '0.1', '--lib_cut', '0', '--reference_levels', "['Sex::Female']", '--tol', '1e-05', '--max_iter', '100', '--conserve', 'False', '--alpha', '0.05', '--output_loaf', '/var/folders/rt/6vtjnzqd64qbxyg4_fpkknjm0000gn/T/q2-DataLoafPackageDirFmt-ga3vswl8']' returned non-zero exit status 1.
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/q2cli/commands.py", line 520, in __call__
results = self._execute_action(
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/q2cli/commands.py", line 581, in _execute_action
results = action(**arguments)
File "<decorator-gen-19>", line 2, in ancombc
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/qiime2/sdk/action.py", line 342, in bound_callable
outputs = self._callable_executor_(
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/qiime2/sdk/action.py", line 576, in _callable_executor_
output_views = self._callable(**view_args)
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/q2_composition/_ancombc.py", line 41, in ancombc
return _ancombc(
File "/Users/emilyyu/opt/miniconda3/envs/qiime2-amplicon-2024.5/lib/python3.9/site-packages/q2_composition/_ancombc.py", line 257, in _ancombc
raise Exception('An error was encountered while running ANCOM-BC'
Exception: An error was encountered while running ANCOM-BC in R (return code 1), please inspect stdout and stderr to learn more.
