- 博客(64)
- 收藏
- 关注
原创 SV merge 整理
vg的github问答区也有类似现象https://github.com/vgteam/vg/issues/4318 怀疑这个人也用的SURVIVOR。构建泛基因组时候,将每个样本的ins 和del提取出来,就可以不添加–ignore-types了。可能其他类型的变异这个软件有的类型不识别吧。我的变异是用minimap2比对,使用svim-asm haploid 执行的变异检测 不知道啥情况。这里的sample_files里面是vcf文件路径,一行一个vcf文件,不使用压缩格式。
2024-12-03 15:57:20
361
原创 gwas数据根据eaf Z 和N 求beta和se
https://www.nature.com/articles/s41590-023-01588-w#Sec10
2023-10-01 13:25:40
1054
原创 根据GWAS数据估算样本量N和使用千人基因组填充maf的参考文献
https://github.com/GenomicSEM/GenomicSEM/wiki/2.1-Calculating-Sum-of-Effective-Sample-Size-and-Preparing-GWAS-Summary-Statistics
2023-10-01 10:24:42
575
原创 R语言学习——线性拟合、检查正态分布、秩和检验、百分比检验、卡方检验、方差检验
【代码】R语言学习——线性拟合、检查正态分布、秩和检验、百分比检验、卡方检验、方差检验。
2023-07-05 20:55:17
810
原创 linux删除带有NA的行
在此命令中,-v选项用于反转匹配,即只输出不包含匹配项的行。'NA’是要匹配的字符串,该命令将从input_file中读取数据,并将不包含’NA’的行写入output_file中。
2023-06-26 19:39:51
303
原创 linux将制表符分隔的文件将outukb-b-6134.txt文件中第一行的ID替换为SNP,ALT替换为A1,REF替换为A2,并将结果保存回原文件。
sed -i ‘1s/variant_id/SNP/’
2023-06-24 19:35:02
248
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人