运行代码的时候出错,我确认文件都存在,报错代码如下所示

 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

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

输入文件格式不标准吧,你截图贴一下 输入文件,看看分组名字在不在table和metadata里面

请先 登录 后评论