然后有小伙伴留言，想知道每一個(gè)個(gè)體對(duì)應(yīng)的單倍型是哪一個(gè)？發(fā)現(xiàn)，10號(hào)染色體的，這個(gè)區(qū)間，有幾個(gè)位點(diǎn)位于一個(gè)block里面。## 導(dǎo)入基因型vcf數(shù)據(jù)library(geneHapR)
vcf = import_vcf("df2.vcf")
# 單倍型分型hapResult <- vcf2hap(vcf)
write.csv(hapResult,"df2-1-hapresult.csv")
上面結(jié)果中，共有11個(gè)單倍型，每個(gè)個(gè)體都會(huì)給出具體的單倍型分型。分割線
1，快來領(lǐng)取 | 飛哥的GWAS分析教程 2，飛哥匯總 | 入門數(shù)據(jù)分析資源推薦
3，數(shù)量遺傳學(xué)，分享幾本書的電子版 4，R語言學(xué)習(xí)看最新版的電子書不香嘛？ 5，書籍及配套代碼領(lǐng)取--統(tǒng)計(jì)遺傳分析導(dǎo)論 6，十一在家把GWAS分析學(xué)會(huì)吧！
|
