larssnip / micropan

R package for microbial pangenomics
21 stars 0 forks source link

Error in `vec_rbind()` occured while preparing the blast pairs pair.tbl #12

Open itsvetkova opened 1 year ago

itsvetkova commented 1 year ago

Hi, Lars, I encounter a problem with preparing the pair.tbl of blast pairs. This is related to vctrs package, but replacing the vctrs package version (0.4.1) with previous versions (to 0.3.8) does not help to solve this error. Here's the code:

self.tbl <- readBlastSelf(file.path("blast", list.files("blast", pattern = "txt$"))) readBlastSelf: ...received 50721 blast-files... ...found 318 self-alignment files... ...returns 14751398 alignment results pair.tbl <- readBlastPair(file.path("blast", list.files("blast", pattern = "txt$"))) readBlastPairs: ...received 50721 blast-files... ...found 50403 alignment files who are NOT self-alignments... Error in vec_rbind(): ! Negative n in compact_rep(). ℹ In file utils.c at line 897. ℹ Install the winch package to get additional debugging info the next time you get this error. ℹ This is an internal error that was detected in the vctrs package. Please report it at https://github.com/r-lib/vctrs/issues with a reprex and the full backtrace. Backtrace: ▆

  1. ├─micropan::readBlastPair(...)
  2. │ └─... %>% filter(.data$Evalue <= e.value)
  3. ├─dplyr::filter(., .data$Evalue <= e.value)
  4. ├─dplyr::rename(...)
  5. ├─dplyr::bind_rows(.)
  6. │ └─vctrs::vec_rbind(!!!dots, .names_to = .id, .error_call = current_env())
  7. └─rlang:::stop_internal_c_lib(...)
  8. └─rlang::abort(message, call = call, .internal = TRUE, .frame = frame)

Please could you help me with solution the problem. Best regards, Irina