arendsee / rhmmer

Simple R utilities for working with HMMER

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

failed to load tblout table

Valentin-Bio opened this issue · comments

Hello, I have previously used rhmmer to load hmmscan results using read_tblot() function. In my current case I want to load a hmmsearch table with nearly 26 thousand hits but I'm getting the following error:

Error in `$<-`:
! Assigned data `descriptions[!grepl(" *#", descriptions, perl = TRUE)]` must be compatible
  with existing data.
✖ Existing data has 25673 rows.
✖ Assigned data has 0 rows.
ℹ Only vectors of size 1 are recycled.
Caused by error in `vectbl_recycle_rhs_rows()`:
! Can't recycle input of size 0 to size 25673.

May this be due to the large amount of hits that report the table that I'm trying to load ?

bests,

Valentín.

Thanks for the report Valentin! rhmmer doesn't support the output of search. It should though. If your data isn't private, could you send me the first 20 lines or so? Then I can use that as a test case and write a new parser function.

Alternatively, if you like writing R code yourself, I would be happy to accept a pull request.