Generating plots from lima report: demultiplex pacbio .bam files
0
0
Entering edit mode
3.0 years ago

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

demultiplex Barcode PacBio • 527 views
ADD COMMENT

Login before adding your answer.

Traffic: 2545 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6