怎么使用shapeit进行单倍型分析，相信很多没有经验的人对此束手无策，为此本文总结了问题出现的原因和解决方法，通过这篇文章希望你能解决这个问题。欢迎关注”生信修炼手册”!shapeit是一款单倍型分析工具，运算速度快，准确率高，是impute2官方推荐的pre-phasing工具。通过隐马可夫模型来分析单倍型，简化的模型示意如下从上到下依次有5个子图，用1到5来表示，需要分成3个部分来看。在1图中，表示的是8个位点构成的8种单倍型，每行表示一个单倍型，每一列代表一个位点，2图中用图状结构来表示上述的单倍型，每个节点表示一个SNP位点，依次用Z1到Z8表示，从1到8的完整路径代表一个单倍型。观察图1可以发现，前4个位点只有3种组成，后4个位点也是同样的，通过4号位点和5号位点的不同连线可以表示所有的单倍型，每条边上的数字代表对应的频数。图5表示的是某个样本的分型结果，分别用0,1,2表示不同的状态，0表示没有突变，1表示杂合突变，2表示纯合突变，根据分型结果拆分成单倍型的时候，杂合突变对应2个allel, 根据这个分型结果可以得到图4中对应的单倍型构成，图4中ref allel用空白方框表示，alt allel用黑色方框表示，对于前5个位点，存在了2个杂合突变，所以有4种路径，后3个位点也是4种。图3表示的是该软件的隐马可夫模型，将真实的单倍型当做是隐藏序列，将根据分型结果预测的单倍型当做是观察序列，通过建模之后来分析隐藏序列的组成，就得到了真实的单倍型分析结果。在文献中，将该软件与其他类似的工具进行了比较，结果示意如下采用了3个不同的数据集，比较了运行时间和错误率，shapeit错误率最低，运行速度最快。该软件的基本用法如下需要指定的参数分成了以下3个部分支持以下3种格式ped/mapbed/bim/famgen/samplevcf
前两种为plink软件的格式，是GWAS分析最常见的文件格式，第三种格式是WTCCC默认的文件格式 第四种是最常见的VCF格式。不同类型的输入文件对应的用法如下对于gen/sample文件格式，可以通过gtool这个软件来进行格式转换，参考基因组对应的连锁图谱，可以提高单倍型分析的准确性，官方提供了hapmap项目的连锁图谱供下载，链接如下http://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html#formats这个是一个可选参数，没有的情况下软件会根据线性模型来进行估算。默认用后缀为haps和sample的两个文件来描述单倍型， haps文件的内容如下所示每列之间用空格分隔，第一列为snp位点所在的染色体名称，第二列为snp id，第三列为染色体的位置，第四列为不同样本中该位点的分型结果，0代表ref allle, 1代表alt allel, 每两列对应一个样本。后缀为sample的文件内容如下所示用来描述样本的信息，同样的空格分隔，前两行内容固定，后续每一行代表一个样本，以上只是该文件最基本内容的展示，还可以有更多的列，用来描述样本的表型信息。在impute2中，phased reference panel会用hap/legend/sample3个文件来表示，通过下列代表可以进行格式转换不同格式的详细解释参考以下链接http://math 香港云主机gen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html#formats事先对需要填充的样本进行phasing, 可以有效提高填充的运行效率，如果后续使用impute2进行基因型填充，推荐使用shapeit对需要填充的样本进行单倍型分析。看完上述内容，你们掌握怎么使用shapeit进行单倍型分析的方法了吗？如果还想学到更多技能或想了解更多相关内容，欢迎关注开发云行业资讯频道，感谢各位的阅读！
这篇文章主要介绍“python中pdb调试的命令是什么”，在日常操作中，相信很多人在python中pdb调试的命令是什么问题上存在疑惑，小编查阅了各式资料，整理出简单好用的操作方法，希望对大家解答”python中pdb调试的命令是什么”的疑惑有所帮助！接下来，…
免责声明：本站发布的图片视频文字，以转载和分享为主，文章观点不代表本站立场，本站不承担相关法律责任；如果涉及侵权请联系邮箱：360163164@qq.com举报，并提供相关证据，经查实将立刻删除涉嫌侵权内容。
微信扫一扫 
