楼主 | 收藏 | 举报 2018-07-01 00:00   浏览:158   回复:3

cnvnator安装及使用方法简介

1、安装gcc

推荐yum安装,自己寻找依存关系,版本当在4.8以上。

2、安装root

a、下载root_v6.03.04.source.tar.gz

b、cd root-6.03.04

3、安装CNVnator

4、使用方法

打赏
沙发 | 回复 | 举报 2017-09-30 01:26
bin size的选择不是经验值,bin size的选择有两个标准:1. mean to sigma的比率大概在4-5之间, 计算mean to sigma的命令为:cnvnator -root my.root -eval2. 跟测序深度有关系:given the same data quality and read length, we observed that the optimal bin size, and thus breakpoint resolution accuracy, scales roughly inversely with the coverage, resulting in ~100-bp bins for 20–30X coverage; ~500-bp bins for 4–6X coverage, and ~30-bp bins for ~100X coverage.至于具体的bin size的选择,就要根据这两条去权衡。
藤椅 | 回复 | 举报 2017-12-14 15:01
非常感谢您的解答! 我再请教您一个问题啊,就是我发现好像CNVnator不能处理case-control的sample, 那么对于这种有normal-tumor的成对样本,您觉得用什么做CNV分析比较好呢?
板凳 | 回复 | 举报 2017-10-06 02:40
**How to select the bin size?1. the ratio of mean to sigam is roughly 4-5; -eval returns values of mean and sigma once histograms are calculated eg: cnvnator -root my.root -his 1000 -d dir_with_chrom_fas cnvnator -root my.root -eval 10002. Specifically, given the same data quality and read length, we observed that the optimal bin size, and thus breakpoint resolution accuracy, scales roughly inversely with the coverage, resulting in ~100-bp bins for 20–30X coverage; ~500-bp bins for 4–6X coverage, and ~30-bp bins for ~100X coverage.**How to filter the results from cnvnator? 1. q0<0.5 2. the first e-value(it is like the blast e-value) is small**批量获得genotype 和 Rd value awk '{print $2}END{print "exit"}' out.cnv | cnvnator -root out.root -genotype 100**-1 for q0 means it couldn’t be calculated. The likely reason is that no reads were mapped in that region.
网站首页 | 关于我们 | 联系方式 | 使用协议 | 版权隐私 | 网站地图  |  排名推广  |  广告服务  |  积分换礼  |  网站留言  |  RSS订阅  |  违规举报
 
免责声明:本站有部分内容来自互联网,如无意中侵犯了某个媒体 、公司 、企业或个人等的知识产权,请来电或致函告之,本网站将在规定时间内给予删除等相关处理。