ONCOCNV 计算过程详解（1）

最新推荐文章于 2022-12-14 20:44:30 发布
最新推荐文章于 2022-12-14 20:44:30 发布
阅读量377 收藏
点赞数 1
分类专栏： 生信
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/dujidan/article/details/115675171
版权

ONCOCNV 计算过程详解（1）

计算原理

Control.stats.txt、Test.stats.txt 中 最后一列的，计算原理

readsCount / (sum_readslCount / sum_pos_len * pos_len)
基于指定的bed文件，bed
readsCount ： 指定区间内的 reads 数
sum_readslCount：所有区间内的 reads 总数
sum_pos_len：所有区间内的 碱基 总数；即，所有的 end - start 之和
pos_len：指定区间内的 碱基 数

背景的 Merge_control.stats.txt 合并到一个中的，计算原理

所有的最后一列取均值

ONCOCNV软件思路分析之control处理.
ONCOCNV软件思路分析之tumor处理.

毒鸡蛋
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
C语言编译过程总结详解
01-19
C语言的编译链接过程要把我们编写的一个c程序（源代码）转换成可以在硬件上运行的程序（可执行代码），需要进行编译和链接。编译就是把文本形式源代码翻译为机器语言形式的目标文件的过程。链接是把目标文件、操作...
ONCOCNV软件思路分析之control处理
weixin_33955681的博客
12-27 425
进行数据初步处理(perl) 统计amplicon的RC(read counts),并且相互overlap大于75%的amplicon合并起来 统计每个amplicon的GC含量,均值, 性别识别并校正，文库大小、长度、GC含量标准化 文库大小标准化 某个sample的文库大小（read counts * amplicon length） / 平均每个ampicon文库大小(sum(re...
ONCOCNV软件思路分析之tumor处理
weixin_33973609的博客
12-27 275
前期处理 perl脚本统计RC(RC(read counts)) 读入control baseline 和 sigma（最后baseline 预测的mad值） 将gc < 0.28或gc > 0.68，sigma乘上1.5,后来又乘以6，对于小于0.01或者大于0.99分位数，sigma取0.01和0.99分位点的sigma 将sigma转化为权重，SigmaForWeights =...
全基因组数据CNV分析简介
热门推荐
庐州月光的博客
08-23 1万+
欢迎关注”生信修炼手册”!除了利用aCGH和snp芯片来检测CNV之外，也可以通过NGS数据来分析CNV, 比如全基因组和全外显子测序。针对全基因组CNV的检测，还针对开发了一种称之为C...
CNV分析学习1（SNP芯片分型原理以及genomestudio部分结果展示）
sweet_yemen的博客
07-01 3436
CNV学习，持续更新中
nvchad自定义设置无法生效
01-21 2954
nvchad，在init.lua下修改或添加快捷键map，但是无法生效，可能是init.lua没有被正确加载，这种情况通常发生在第一次配置nvchad。 怎么看模块是否被正确加载？ 答… nvim -V 如果未被正确加载咋办？ 答…假如xxx是模块或者包名， lua通常定义模块的方式有两种，单个文件xxx.lua或者文件夹形式，xxx/init.lua，这和python有点类似。因此，最简单的方式就是把模块所在路径添加到LUA_PATH，注意通配符： export LUA_PATH="pat...
单片机执行指令过程详解
01-19
单片机执行指令过程详解 单片机执行程序的过程，实际上就是执行我们所编制程序的过程。即逐条指令的过程。计算机每执行一条指令都可分为三个阶段进行。即取指令-----分析指令-----执行指令。 取指令的...
视觉引导教程，计算过程详解
04-15
本教程将详细解释视觉引导中的计算过程，并以单相机和多相机引导为案例，提供具体的计算方法和步骤。 首先，视觉引导系统的核心在于座标系定位引导方法，它可以分为单相机引导和多相机引导。单相机引导通常包括两种...
软件开发过程详解-.pdf
11-17
软件开发过程详解-.pdf
Excel计算公式大全详解.pdf
08-21
Excel计算公式大全详解
nvcode:真正做到极致Neovim Config NVCode
03-17
NVCode 如果您要查找我的旧配置，请在此存储库中检出两个快照分支，其中一个用于CoC，一个用于本机LSP。 只需一条命令即可安装 bash <( curl -s https://raw.githubusercontent.com/ChristianChiarulli/nvcode/master/utils/installer/install-nv-code.sh ) 运行之后，您将可以访问nv命令，这不会覆盖您的nvim配置。 因此，您可以同时安装两者 获取最新版本的Neovim cd ~ sudo rm -r neovim git clone https://github.com/neovim/neovim cd neovim sudo make CMAKE_BUILD_TYPE=Release install cd ~ sudo rm -r neovim VSCode支持 在
nvchad添加插件
03-22 1674
首先创建 lua/custom/plugins/init.lua, 在下面添加插件的仓库名： return { { "max397574/better-escape.nvim", event = "InsertEnter"}, { "ggandor/lightspeed.nvim"}, } 修改 lua/custom/chadrc.lua，添加install，没有该文件的话从examples下面拷贝一份过去 local userPlugins = require "custom..
使用CNVkit进行CNV分析
庐州月光的博客
09-06 9413
欢迎关注”生信修炼手册”!CNVkit是一款CNV预测软件，适用于全外显子，目的区域靶向测序等数据的CNV检测，官网如下https://cnvkit.readthedocs.io/en/...
neovim安装配置
guohw
09-08 1843
neovim 简单自定义配置
【Python】利用滑动窗口计算全基因组每个窗口上CNV的拷贝数和Vst
每天都要学Python的博客
05-15 3153
目录Vst介绍计算每个窗口的绝对拷贝数1.文件准备2.编写脚本计算每个窗口的Vst Vst介绍 Vst是通过计算拷贝数的方差来衡量不同群体之间CNV的分化的一个指标，类似于Fst的概念，可以用来鉴定一些高分化的区域。 计算方法如下： Vpop1是指群体1的copy number的方差； Vpop2是指群体2的copy number的方差； Vtotal是全部个体的copy number的方差； Npop指的对应群体的个体数。 部分文献中，会按照10K的窗口，2K的步长，去计算全基因组上每个窗口的Vst。而
根据bed文件计算CG含量
somehappyness的博客
09-27 1276
看到一个临床小panel的测序数据的质控结果，CG含量 37%，（一般的结果都是45-55%的范围）那是实验和测序造成的，还是本来就是这个这个数值呢。 思路：根据bed文件提取基因组的的序列，然后再计算序列的CG含量。写脚本有点麻烦，下面这个方法稍微简单一点： 提取序列： for region in $( awk '{print $1":"$2"-"$3}' my_two_gene.bed ) ...
使用linux命令 根据序列ID快速提取fastq序列
dujidan的博客
04-02 6725
做生信的人经常会接触到一些fastq数据，有时想要从fastq文件中提取出某些序列来查看。 一般情况可以使用grep命令，就可以实现 grep -A 3 seq_id fastq.file # 匹配seq_id 并向下去3行 但是，如果需要的序列很多，在使用grep就会很慢了，所以这里给出了 awk 的命令，速度简直快的飞起。 awk -F ' ' 'NR==FNR {a[$1]=1; next} { if (a[substr($1,2)]) {print $0; getline b; print b.
R语言---生信分析---count转换成TPM、FPKM
最新发布
dujidan的博客
12-14 6478
R语言---生信分析---count转换成TPM、FPKM
如何下载NCBI的ftp数据
dujidan的博客
06-25 6146
如何下载NCBI的ftp数据 因为要从refseq中下载数据，知道 ftp 地址，浏览器打不开，用了好多下载工具都下不下来，所以有点难受。。。 一般的 ftp 下载链接是长这样的： ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/001/696/305/GCF_001696305.1_UCN72.1 但无意间找到NCBI的FQA，所以找到了下载的方法（手动狗头） 主要的方法就是：把下载链接开头的 ftp 替换成 rsync，然后用 rsync 进行传输 这个也是可以使用
DDR内存系统功率计算详解
"DDR功耗计算说明文档涵盖了DDR内存系统的日期表、系统配置、功耗计算详解以及DDR计算模块的各个方面。文档来源于Micron Technology，提醒读者产品和规格可能会不经通知而更改，仅供评估和参考，不提供任何担保。本...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值