Roleren / ORFik

MIT License
33 stars 9 forks source link

QCreport(df.rfp) throws an error, when working with ribo-seq data #135

Closed yeroslaviz closed 1 year ago

yeroslaviz commented 1 year ago

before i can p-shift my samples, I try to run the QCreport(df.rfp) command, but it ends up throwing an error when plotting the meta coverage plots:

> QCreport(df.rfp)
Started ORFik QC report for experiment: ORFiK_Testing_Ribo-seq
Saving, ofst files to: /fs/pool/pool-cox-projects-bioinformatics/Dept_Hartl/Rongqin/ORFik_testing/processed_data/Ribo-seq/ORFiK_Testing/aligned/ofst/
...
- Meta coverage plots
  - Metacoverage of CDS region only
RFP_TRiC-total_r1

RFP_Ssb-total_r2

RFP_TRiC-bound-pyromycin_r1

RFP_Ssb-total_r1

Error: BiocParallel errors
  4 remote errors, element index: 1, 2, 3, 4
  0 unevaluated and other errors
  first remote error:
Error: BiocParallel errors
  1 remote errors, element index: 1
  2 unevaluated and other errors
  first remote error:
Error in S4Vectors:::normarg_names(value, class(x), length(x)): attempt to set too many names (2) on PartitioningByEnd object of length
  0

In addition: Warning messages:
1: In dir.create(stats_folder, recursive = TRUE) :
  '/fs/pool/pool-cox-projects-bioinformatics/Dept_Hartl/Rongqin/ORFik_testing/processed_data/Ribo-seq/ORFiK_Testing/aligned/QC_STATS' already exists
2: In correlation.plots(df, stats_folder, region, plot.ext = plot.ext,  :
  Complex correlation plots are for now deprecated!
3: In QCplots(df, "mrna", stats_folder, plot.ext = plot.ext, complex.correlation.plots = complex.correlation.plots,  :
  No 5' UTRs or 3' of significant length defined, UTR metacoverage plots can not be made, check your annotation file. In case no UTRs exist in your annotation,
    you can add pseudo UTRs, to also see coverage profiles over those areas.

Not sure if this is because I'm working with Ribo-Seq or not, but I would appreciate the help fixing it.

thanks Assa

Roleren commented 1 year ago

Did you manage to get this to work ?

Roleren commented 1 year ago

This issue btw, is because of you leader annotation being lacking, I will close this issue in a few days if you don't need anything else :)

yeroslaviz commented 1 year ago

I will close it now and if needed will reopen, when I need further help.

Thanks for everything.