输入文件格式不标准吧,你截图贴一下 输入文件,看看分组名字在不在table和metadata里面
Rscript $scriptdir/compare_errorbar.r -i $table \
> -m $metadata -t groupb \
> -a low -b high \
> --daa_method $m \
> --p_values_threshold 1 \
> --top 20 \
> -o compare_errorbar --prefix low_vs_high_L6
--> Q&A for bioinformatics, please visit the website: https://www.omicsclass.com/
--> R beginners ? I suggest your learning R language: https://study.omicsclass.com/index
Failed to create bus connection: No such file or directory
── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr 1.1.3 ✔ purrr 1.0.2
✔ forcats 1.0.0 ✔ stringr 1.5.0
✔ ggplot2 3.4.4 ✔ tidyr 1.3.0
✔ lubridate 1.9.3
── 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
Warning message:
In system("timedatectl", intern = TRUE) :
running command 'timedatectl' had status 1
New names:
• `Sepsis` -> `Sepsis...21`
• `Sepsis` -> `Sepsis...26`
Rows: 47 Columns: 33
── Column specification ────────────────────────────────────────────────────────
Delimiter: "\t"
chr (20): Wheel, BarcodeSequence, Genotype, Microbiota, Cohoused, Gender, Tr...
dbl (12): #SampleID, Donor_No, Exitus, Geschlecht, SofaAnstieg, Hauotdiagnos...
lgl (1): Groupf
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
Error in `select()`:
! Can't subset columns with `metaData[, 1, drop = T]`.
✖ Can't convert from `metaData[, 1, drop = T]` <double> to <integer> due to loss of precision.
Backtrace:
▆
1. ├─dplyr::select(abundance_data, metaData[, 1, drop = T])
2. ├─dplyr:::select.data.frame(abundance_data, metaData[, 1, drop = T])
3. │ └─tidyselect::eval_select(expr(c(...)), data = .data, error_call = error_call)
4. │ └─tidyselect:::eval_select_impl(...)
5. │ ├─tidyselect:::with_subscript_errors(...)
6. │ │ └─rlang::try_fetch(...)
7. │ │ └─base::withCallingHandlers(...)
8. │ └─tidyselect:::vars_select_eval(...)
9. │ └─tidyselect:::walk_data_tree(expr, data_mask, context_mask)
10. │ └─tidyselect:::eval_c(expr, data_mask, context_mask)
11. │ └─tidyselect:::reduce_sels(node, data_mask, context_mask, init = init)
12. │ └─tidyselect:::walk_data_tree(new, data_mask, context_mask)
13. │ └─tidyselect:::as_indices_sel_impl(...)
14. │ └─tidyselect:::as_indices_impl(...)
15. │ └─vctrs::vec_as_subscript(x, logical = "error", call = call, arg = arg)
16. └─rlang::cnd_signal(x)
Execution halted