missing value where TRUE/FALSE needed in finemap_enriched_windows.R

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
YeoLab / skipper

missing value where TRUE/FALSE needed in finemap_enriched_windows.R #12
