Question: Potential Problems with BAM file header, r package input, or RNASeq data pipeline
gravatar for ss650012
2.5 years ago by
ss6500120 wrote:

Hi, I am trying to make the GLMVC pipeline ( work and am incurring a few issues. For this pipeline, you simply need to download the package and run it in the command line with two BAM files, one normal and one tumor. I have done that and am receiving the following error that is very persistent:

20180719 12:23:23 - filter process started at 07/19/2018 12:23:23
20180719 12:23:24 - E:Loading required package: profileModel
20180719 12:23:24 - E:'brglm' will gradually be superseded by 'brglm2' (, which provides utilities for mean and median bias reduction for all GLMs and methods for the detection of infinite estimates in binomial-response models.
20180719 12:23:24 - E:Error in read.table(inputfile, header = T, stringsAsFactors = F) : 
20180719 12:23:24 - E:  no lines available in input
20180719 12:23:24 - E:Execution halted
Error: R command failed to genearte result as /Users/as/GLMVC/rep1_r1.unfiltered.tsv. Check the log file /Users/as/GLMVC/rep1_r1.tsv.r.log.

You can manully run the R script file /Users/as/GLMVC/rep1_r1.tsv.r to find out the problem.

Trace:   at RCPA.R.RProcessor.Process () [0x00219] in <2447b922059f4ad28c04638fa7048d73>:0 
  at CQS.Genome.SomaticMutation.FilterProcessor.Process () [0x00084] in <b7c50a2fd30542368a2470adaa3a0af4>:0 
  at CQS.Genome.SomaticMutation.PipelineProcessor.Process () [0x000e5] in <b7c50a2fd30542368a2470adaa3a0af4>:0 
  at RCPA.Commandline.AbstractCommandLineCommand`1[T].DoProcess (System.String[] args, System.Boolean result) [0x0005c] in <2447b922059f4ad28c04638fa7048d73>:0 
  at RCPA.Commandline.AbstractCommandLineCommand`1[T].Process (System.String[] args) [0x00015] in <2447b922059f4ad28c04638fa7048d73>:0

I have tried editing the BAM file header using Picard AddOrReplaceReadGroups, and I have ensured that all of my BAM files were processed using the same pipeline.

Any help would be appreciated.

Thank you so much in advance!

R rna-seq bam glmvc software error • 576 views
ADD COMMENTlink modified 2.5 years ago by _r_am32k • written 2.5 years ago by ss6500120

What does the log file ( /Users/as/GLMVC/rep1_r1.tsv.r.log ) says?

Did you try manually running the script "/Users/as/GLMVC/rep1_r1.tsv.r", as the error message suggested?

ADD REPLYlink written 2.5 years ago by h.mon32k

Followed your advise and ran R script manually and unfortunately I am getting error

No lines are available for input I believe this error is being produced because I think GLMVC suppose to make input for the Read Table R package internally. In order to produce this file GLMVC run samtools Mpileup method after running this function manually on the files I receive

[mpileup] 2 samples in 2 input files

is this correct I am unsure since the file which is outputed does not seem to be right. Please help.

Thanks in advance

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by KG0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2285 users visited in the last hour