Comments (2)
Hi @andreykoch
I have included GC and map wig files for 50kb and 10kb. For the normal panel, unfortunately, the ultra-low coverage (0.1x) is not sufficient to generate a panel for 50kb. I am not sure the read coverage of your data but if it is also ~0.1x I would not recommend using 50kb. If you are using much deeper coverage, then using a matched normal would be a much better suited for normalization.
Hope this helps,
Gavin
from ichorcna.
Hi Gavin,
Thanks for your input. The data I have has x0.4 coverage, so hopefully it should be sufficient to generate PoN with 50Kb bins.
On the other hand, I'm experiencing hickups with generating PoN with the error message similar to the issue24.
I have previously generated 4 normal wig files using
$readCounter --window 50000 --quality 20 --chromosome "chr1,chr2,chr3,chr4,chr5,chr6,chr7,chr8,chr9,chr10,chr11,chr12,chr13,chr14,chr15,chr16,chr17,chr18,chr19,chr20,chr21,chr22,chrX,chrY" $normal_sample.bam > $normal_sample.wig
And they look OK, e.g.
fixedStep chrom=chr1 start=1 step=50000 span=50000
5
47
31
The actual command line for the PoN generation looks like this:
Rscript /rdseqdata/bin/ichorCNA-master/scripts/createPanelOfNormals.R
--filelist /rndhome/akoch/SamBam_file_example/TMP/ref/wigPoNfiles.txt
--gcWig /rdseqdata/bin/ichorCNA-master/inst/extdata/gc_hg19_50kb.wig
--mapWig /rdseqdata/bin/ichorCNA-master/inst/extdata/map_hg19_50kb.wig
--centromere /rdseqdata/bin/ichorCNA-master/inst/extdata/GRCh37.p13_centromere_UCSC-gapTable.txt
--outfile /rdseqdata/bin/ichorCNA-master/inst/extdata/PoN_hg19_50Kb_median_normAutosome_median.rds
And the output before blowup looks like this:
Loading normal file:/rndhome/akoch/SamBam_file_example/TMP/ref/1801100551-NIK-3-26_TCATCTCC_L003_R1_export.wig
Slurping: /rndhome/akoch/SamBam_file_example/TMP/ref/1801100551-NIK-3-26_TCATCTCC_L003_R1_export.wig
Parsing: fixedStep chrom=chr1 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr2 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr3 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr4 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr5 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr6 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr7 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr8 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr9 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr10 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr11 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr12 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr13 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr14 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr15 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr16 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr17 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr18 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr19 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr20 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr21 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chr22 start=1 step=50000 span=50000
Parsing: fixedStep chrom=chrX start=1 step=50000 span=50000
Parsing: fixedStep chrom=chrY start=1 step=50000 span=50000
Sorting by decreasing chromosome size
Reading GC and mappability files
Slurping: /rdseqdata/bin/ichorCNA-master/inst/extdata/gc_hg19_50kb.wig
Parsing: fixedStep chrom=1 start=1 step=50000 span=50000
Error in seq.default(from = as.integer(tokens[2]), by = as.integer(tokens[3]), :
'from' must be a finite number
Calls: wigToRangedData -> seq -> seq.default
Execution halted
I would appreciate any hints on what is going on.
Andrey
from ichorcna.
Related Issues (20)
- Slurping
- Duplicate levels error after loading GenomeInfoD using mouse reference HOT 1
- question about coverage computation
- logR and copy.number in cna.seg. files HOT 1
- New version release HOT 3
- Panel of normals for 10kb and 50kb
- Intuitive understanding of `txnE` and `txnStrength` HOT 1
- map_hg38_1k.wig HOT 2
- PON - unable to find an inherited method for function ‘seqinfo’ for signature ‘"data.table"’ HOT 2
- Is cutoff 0.03 suitable for cfDNA 5hmC sequencing?
- Failed to connect to ftp.ncbi.nlm.nih.gov HOT 1
- correctedDepth.txt File interval missing
- PON- missing .rds metadata columns
- Final segmentation from alternate ichorCNA solutions HOT 1
- Error while using a BED file
- Rumtime error (Error in lambdas[, , i] <- na.omit(param$lambda) : incorrect number of subscripts) HOT 9
- Are duplicated records in BAM file counted or not? HOT 1
- Add config option for people with color vision defficiency HOT 1
- Error in .order_seqlevels(chrom_sizes[, "chrom"]) :
- The subclonal fraction (fraction of tumor-derived DNA that is subclonal) is 0.446 but tumor fraction is 0
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from ichorcna.