Rscript --vanilla skipper/bb63a25/bin/skipper/tools/finemap_enriched_windows.R output/finemapping/nt_coverage/Q0542b_CB.nt_coverage.bed output/finemapping/mapped_sites/ Q0542b_CB
── Attaching core tidyverse packages ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────── tidyverse 2.0.0 ──
✔ dplyr 1.1.2 ✔ readr 2.1.4
✔ forcats 1.0.0 ✔ stringr 1.5.0
✔ ggplot2 3.4.2 ✔ tibble 3.2.1
✔ lubridate 1.9.2 ✔ tidyr 1.3.0
✔ purrr 1.0.1
── 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 messages:
1: package ‘tidyverse’ was built under R version 4.2.3
2: package ‘ggplot2’ was built under R version 4.2.3
3: package ‘tibble’ was built under R version 4.2.3
4: package ‘tidyr’ was built under R version 4.2.3
5: package ‘readr’ was built under R version 4.2.3
6: package ‘purrr’ was built under R version 4.2.3
7: package ‘dplyr’ was built under R version 4.2.3
8: package ‘stringr’ was built under R version 4.2.3
9: package ‘forcats’ was built under R version 4.2.3
10: package ‘lubridate’ was built under R version 4.2.3
Error in `mutate()`:
ℹ In argument: `local_maximum = row_number() %in% localMaxima(enrichment_heuristic) | finemapped`.
ℹ In group 1: `window_n = 1`, `strand = "-"`, `chr = "chr1"`.
Caused by error in `if (x[[1]] == x[[2]]) ...`:
! missing value where TRUE/FALSE needed
Backtrace:
▆
1. ├─... %>% ...
2. ├─dplyr::filter(., enrichment_heuristic > enrichment_median, local_maximum)
3. ├─dplyr::mutate(...)
4. ├─dplyr:::mutate.data.frame(...)
5. │ └─dplyr:::mutate_cols(.data, dplyr_quosures(...), by)
6. │ ├─base::withCallingHandlers(...)
7. │ └─dplyr:::mutate_col(dots[[i]], data, mask, new_columns)
8. │ └─mask$eval_all_mutate(quo)
9. │ └─dplyr (local) eval()
10. ├─row_number() %in% localMaxima(enrichment_heuristic)
11. ├─global localMaxima(enrichment_heuristic)
12. └─base::.handleSimpleError(...)
13. └─dplyr (local) h(simpleError(msg, call))
14. └─rlang::abort(message, class = error_class, parent = parent, call = error_call)
Execution halted
