Entering edit mode
3.1 years ago
MEJIA.MAZA.ALAN
•
0
Hi All,
I am new in pacbio demultiplexing strategy but I was able to use Lima (split named). I obtained individual .bam files but I can not generate plots using Rscript_report.
The command I used in the RStudio terminal was: Rscript --vanilla report_detail.R [prefix].demux.lima.report
I tried with different samples and it keeps giving me the same error which I don't completely understand where it comes.
The error output I get is this:
WARNING: ignoring environment value of R_HOME
Error: Problem with filter() input ..2.
x Input ..2 must be a logical vector, not a integer.
i Input ..2 is PassedFilters.
Backtrace:
█
├─%>%(...)
├─dplyr::mutate(., cs = cumsum(n))
├─dplyr::count(., BarcodePair, ScoreCombined)
├─dplyr::arrange(., BarcodePair, desc(ScoreCombined))
├─dplyr::group_by(., BarcodePair)
├─dplyr::select(., BarcodePair, ScoreCombined)
├─dplyr::filter(., BarcodePair %in% unique_bps$BarcodePair)
├─dplyr::filter(., Barcoded, PassedFilters)
├─dplyr:::filter.data.frame(., Barcoded, PassedFilters)
│ └─dplyr:::filter_rows(.data, ...)
│ ├─base::withCallingHandlers(...)
│ └─mask$eval_all_filter(dots, env_filter)
├─dplyr:::abort_glue(...)
│ ├─rlang::exec(abort, message = message, class = class, !!!data)
│ └─(function (message = NULL, class = NULL, ..., trace = NULL, parent = NULL, ...
│ └─rlang:::signal_abort(cnd)
│ └─base::signalCondition(cnd)
└─(function (e) ...
Execution halted
I'd highly appreciate any help. Thanks in advance