符合過濾條件的snp應(yīng)該符合以下參數(shù)：
QD > 2.0
FS < 60.0
MQ > 40.0
MQRankSum > -12.5
ReadPosRankSum > -8.0
SOR < 3.0
而不是FS > 60.0,SOR > 3.0。對應(yīng)硬過濾中視為需要過濾的位點(diǎn)的表達(dá)式為"QD < 2.0 || MQ < 40.0 || FS > 60.0 || SOR > 3.0 || MQRankSum < -12.5 || ReadPosRankSum < -8.0"。沒考慮過為什么你在table的時(shí)候結(jié)果剛好相反嗎？
「簡化基因組」如何過濾用GATK分析得到的SNPGATK官方提供了一個(gè)SNP過濾的標(biāo)準(zhǔn)，howto-apply-hard-filters-to-a-call-set,如果你按照它的要求來過濾簡化基因組中的SNP數(shù)據(jù)，也就是...
plink的pca有人用過嗎，效果如何？感覺plink計(jì)算的很快，是算法不一樣？
PCA的分析常用的工具 PCA分析中常用的工具有GCTA中的PCA模塊，老牌的軟件EIGENSOFT中的smartpca，還有很多最近推出的R包都能夠做PCA分析。然后分析完后，可視化的...
大家好，本周給大家分享的是2021年12月份發(fā)表在Current biology上與中國家馬體型相關(guān)主效基因突變的一篇文章。 文章題目：A single-nucleotide...
pixy用于從vcf變異文件中計(jì)算pi（π）,Dxy和Fst的值。 1.使用conda安裝pixy github[https://github.com/ksamuk/pixy...
終于超過1k粉絲啦～之前說要小小慶祝一下，害怕直接放文章里會被封文章，所以福利大家可以見評論。前幾周忙著答辯畢業(yè)事宜沒有時(shí)間更新，雖然這幾天也挺忙哈哈哈哈，不過還是要定期更新...
許多R包都可以畫?；鶊D和輪狀圖，如 circlize, networkD3, and ggalluvial。這些包都非常好，非常有用。問題來了，吃飽了撐的再弄一個(gè)同樣的畫圖包...
aplot(https://github.com/YuLab-SMU/aplot[https://github.com/YuLab-SMU/aplot])是Y叔開發(fā)的一個(gè)基于...
MOVICS包整合多組學(xué)數(shù)據(jù)進(jìn)行分析，其安裝的過程太讓人心酸，大家安裝的過程慢慢體會，下面就來學(xué)習(xí)下這個(gè)由“大魚海棠”老師開發(fā)的R包吧。 1.示例數(shù)據(jù) brca.tcga.R...
論文是 Global burden of lung cancer attributable to ambient fine particulate matter pollut...
每次有正事的時(shí)候就想逃避寫個(gè)簡書，假裝自己沒有干不正經(jīng)的事。（明明我還有一些數(shù)據(jù)需要分析呢，但是毫無分析的動(dòng)力。）好的，今天我們來介紹一個(gè)很酷炫的圖：弧形圖（Arc Diag...
@斬毛毛 請問如何隨機(jī)選取20ksnp呢？
admixture 群體結(jié)構(gòu)分析tructure是與PCA、進(jìn)化樹相似的方法，就是利用分子標(biāo)記的基因型信息對一組樣本進(jìn)行分類，分子標(biāo)記可以是SNP、indel、SSR。相比于PCA，進(jìn)化樹，群體結(jié)構(gòu)分析可明...
請問小鼠的全基因組重測序需要基因型填充嗎？
群體遺傳學(xué)習(xí)筆記-基因型缺失數(shù)據(jù)的填充Genotype Imputation是在高通量測序中常出現(xiàn)的定義，按照義譯就是基因型填充。要真正理解imputation這個(gè)概念，我們就需要先理解基因型缺失（genotyp...
請問圖一熱圖的配色是什么呢？感覺很好看
10X單細(xì)胞WGCNA分析聯(lián)合10X空間轉(zhuǎn)錄組WGCNA分析各位同學(xué)，大家好，今天我們繼續(xù)WGCNA的分析點(diǎn)，這一次，我們需要聯(lián)合兩種技術(shù)，10X單細(xì)胞和10X空間轉(zhuǎn)錄組，也就是在多組學(xué)層面上分析WGCNA得到的模塊，關(guān)于WGCNA，...
