Comments (3)
Hi,
-Inf should be log(0). These must be segments where FREEC predicts copy number of zero.
By visualization, I meant to visualize the results as .png to visually evaluate the amount of noise after the normalization and the quality of FREEC's CNA calls.
from freec.
Hi,
I have a rough idea of your advice.I should visualize ratio.txt to remove the noise and outlier.Later I'll take care of it.
Our research is drawing to a close.We reviewed all the steps and collated them.I am currently working on the first draft of my article.As I review the FREEC, I have a few questions.
- 0 means Zero copies of DNA in this region predited in _CNVs file. And 0 always corresponds to a loss according to my output.But how to explain it? 0 fragment is lossing from this region? Or too much 0 means noise,and I should filter them, as you
adviced before?
Part of my _CNVs file:
NC_048218.1 68935104 69748872 1 loss
NC_048218.1 72301368 72327936 0 loss
NC_048218.1 89715216 89736864 0 loss
NC_048218.1 94680480 94689336 7 gain
NC_048218.1 98317344 98352768 1 loss
NC_048218.1 98685360 98714880 0 loss
NC_048218.1 100452624 100478208 0 loss
NC_048218.1 109823256 109882296 0 loss
NC_048218.1 115849272 116669928 3 gain
NC_048218.1 130496112 130513824 0 loss
- The CNVs FEEEC predicted show two types: loss and gain. I wonder if FREEC writes the normal region into the output -- the copy number does not change compared to the reference genome. If 0 means no change ,why does "loss" appears?
Wish your reply!
from freec.
Hi, I have a rough idea of your advice.I should visualize ratio.txt to remove the noise and outlier.Later I'll take care of it.
Our research is drawing to a close.We reviewed all the steps and collated them.I am currently working on the first draft of my article.As I review the FREEC, I have a few questions.
- 0 means Zero copies of DNA in this region predited in _CNVs file. And 0 always corresponds to a loss according to my output.But how to explain it? 0 fragment is lossing from this region? Or too much 0 means noise,and I should filter them, as you
adviced before?Part of my _CNVs file: NC_048218.1 68935104 69748872 1 loss NC_048218.1 72301368 72327936 0 loss NC_048218.1 89715216 89736864 0 loss NC_048218.1 94680480 94689336 7 gain NC_048218.1 98317344 98352768 1 loss NC_048218.1 98685360 98714880 0 loss NC_048218.1 100452624 100478208 0 loss NC_048218.1 109823256 109882296 0 loss NC_048218.1 115849272 116669928 3 gain NC_048218.1 130496112 130513824 0 loss
- The CNVs FEEEC predicted show two types: loss and gain. I wonder if FREEC writes the normal region into the output -- the copy number does not change compared to the reference genome. If 0 means no change ,why does "loss" appears?
Wish your reply!
Hi,
I 'm sorry to trouble you. I found that I forgot a truth -- there are 2 chromosomes.
And I missed the meaning of CN. I misunderstood it to the number of copy number region. I always tried to understand CN as how many losses or gains there are. But in reality the two are not equal.
So the anwser is 0 and 1 both represent loss. And CN=2 represents the normal fragments ,which FREEC writes in the ratio.txt.
Thanks to this review, I have to make some adjustments to my research.Best wishes!
from freec.
Related Issues (20)
- Questions about the output of CNVs and the use of assess_significance HOT 3
- run with multiple control files HOT 2
- failed to run segmentation on chr*
- Error: to calculate BAF values HOT 2
- New release including makeGraph 2.0? HOT 1
- Make graph error HOT 1
- Problem making plots HOT 2
- Somatic cnv issue
- All zero reads with segfault HOT 1
- freec can't extract reads from bam file HOT 1
- questions about the output HOT 7
- missing p values in assess_significance.R
- Some settings for WGS explaind in manual are not included in provided config file for WGS.
- error in running gccount HOT 1
- No normalization although it is requested HOT 2
- Understanding _ratio.txt and _BAF.txt for counting informative markers HOT 1
- What is the meaning of suclonal information in the _ratio.txt file HOT 1
- get the float copy number values HOT 2
- Segmentation fault after samtools or sambamba read the BAM file HOT 2
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 freec.