Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error while using a BED file #116

Open
mokhouja opened this issue Dec 6, 2022 · 0 comments
Open

Error while using a BED file #116

mokhouja opened this issue Dec 6, 2022 · 0 comments

Comments

@mokhouja
Copy link

mokhouja commented Dec 6, 2022

Hi,

I am a complete newbie in data science and R but I have been trying to use the script to analyze CNA for a targeted capture sequencing in cfDNA as following:

/scratch/scripts/runIchorCNA.R --id tumor_sample   --WIG test.wig --ploidy "c(2,3)" --normal "c(0.5,0.6,0.7,0.8,0.9)" --maxCN 5   
--gcWig /scratch/inst/extdata/gc_hg38_1000kb.wig   --mapWig /scratch/inst/extdata/map_hg38_1000kb.wig   
--centromere /scratch/inst/extdata/GRCh38.GCA_000001405.2_centromere_acen.txt   
--includeHOMD False --chrs "c(1:22, \"X\")" --chrTrain "c(1:22)"   --estimateNormal True --estimatePloidy True 
--es[AtimateScPrevalence True   --scStates "c(1,3)" --txnE 0.9999 --txnStrength 10000 --outDir ./ 
--exons.bed  CNV.target.bed --genomeBuild hg38

and I am keep getting the following error:

Correcting Tumour
Removed 121 bins near centromeres.
Applying filter on data...
Correcting for GC bias...
Correcting for mappability bias...
Filtering low uniqueness regions with mappability score < 0.9
Removed 1 bins near centromeres.
Determining gender...Gender female
Outputting to: .//tumor_sample.correctedDepth.txt
Warning messages:
1: In simpleLoess(y, x, w, span, degree = degree, parametric = parametric,  :
  span too small.   fewer data values than degrees of freedom.
2: In simpleLoess(y, x, w, span, degree = degree, parametric = parametric,  :
  pseudoinverse used at 0.36574
3: In simpleLoess(y, x, w, span, degree = degree, parametric = parametric,  :
  neighborhood radius 0.0032647
4: In simpleLoess(y, x, w, span, degree = degree, parametric = parametric,  :
  reciprocal condition number  0
5: In simpleLoess(y, x, w, span, degree = degree, parametric = parametric,  :
  There are other near singularities as well. 0.0033244
6: In .Seqinfo.mergexy(x, y) :
  The 2 combined objects have no sequence levels in common. (Use
  suppressWarnings() to suppress this warning.)
runEM: Initialization
runEM iter1: Expectation
Error in runEM(dataMat, chr, chrInd, param, maxiter, verbose, estimateNormal = estimateNormal,  : 
  INTEGER() can only be applied to a 'integer', not a 'NULL'
Calls: HMMsegment -> runEM
Execution halted

and when running without using a bed file it works but gives the following warnings:

1: In KernSmooth::bkde2D(x, bandwidth = bandwidth, gridsize = nbin,  :
  Binning grid too coarse for current (small) bandwidth: consider increasing 'gridsize'
2: In KernSmooth::bkde2D(x, bandwidth = bandwidth, gridsize = nbin,  :
  Binning grid too coarse for current (small) bandwidth: consider increasing 'gridsize'
3: In KernSmooth::bkde2D(x, bandwidth = bandwidth, gridsize = nbin,  :
  Binning grid too coarse for current (small) bandwidth: consider increasing 'gridsize'
4: In KernSmooth::bkde2D(x, bandwidth = bandwidth, gridsize = nbin,  :
  Binning grid too coarse for current (small) bandwidth: consider increasing 'gridsize'

my .BED file looks like this

Chr	start	end
1	23558903	23559184
1	23559184	23559465
1	26696365	26696631

I'll appreciate your advice.
Bests,
Mouhamad

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant