@孟祥霄 應(yīng)該是R包和R版本的原因缤谎。如果你用4.1.0版本的R访娶，那這個R包就不會報錯
structure繪圖（Rscript）pophelper（R 包） Usage:Rscript ../../script/draw_admixture.R order admixture.fam root_tes...
發(fā)簡信
IP屬地：北京
-
-
群體結(jié)構(gòu)分析的繪圖腳本商叹，基于sns腳本 首先導(dǎo)入需要用到的包。--order和-n都是用于指定順序的蕴忆，有一個參數(shù)就可以了。其他都是必選參數(shù)。 核心的排序函數(shù)和繪圖函數(shù) 讀入數(shù)據(jù)颤绕，調(diào)顏色 畫圖 中間上色的部分和...
-
像是輸入文件的問題。你再檢查一下輸入
structure繪圖（Rscript）
pophelper（R 包） Usage:Rscript ../../script/draw_admixture.R order admixture.fam root_tes...
-
基因流分析——D統(tǒng)計基本概念 基因滲入(introgression)： 是在遺傳學(xué)中祟身，指兩個基因庫間的基因流動奥务，通常是經(jīng)過種間雜交產(chǎn)生⊥嗔颍基因滲入是一個長期的過程氯葬，它可能需要許多代雜交才能產(chǎn)生回交...
-
尋找SV.vcf文件中是否含有某個片段的插入現(xiàn)在已知這個片段在3號染色體上，基本思路是先將vcf文件變成fa文件婉陷，然后用blastn處理帚称。 第四列記錄的是INS的序列內(nèi)容 在實際操作過程中竟然出現(xiàn)了重復(fù)三次，并且在vc...
-
重測序-篩去部分樣品并統(tǒng)一編號
篩取vcf中的樣品 samplelist_vcf.txt中按列排樣本名即可 重命名秽澳，再轉(zhuǎn)為bed文件 rename_snp.txt是一個兩列文件闯睹，分別是舊ID和新ID 篩取b...
-
@小鹿不吃香菜 可能這個名字叫sub的其實已經(jīng)是ccs了
把HiFi reads變成轉(zhuǎn)錄本文件
這一次是主線任務(wù)，不過要簡單許多孩锡。官方幫助：IsoSeq/isoseq-clustering.md at master · PacificBiosciences/IsoSeq...
-
不愧是徐州更老師
說說R語言的轉(zhuǎn)置
在群里看到了一個問題酷宵，他用R語言的轉(zhuǎn)置函數(shù)t()對下面形式的數(shù)據(jù)進(jìn)行了轉(zhuǎn)置 轉(zhuǎn)置之后的數(shù)據(jù)結(jié)構(gòu)如下 他發(fā)現(xiàn)轉(zhuǎn)置之后多了一行，也就是有了geneid這一行躬窜，和他預(yù)期的不符合浇垦。按...
-
@今天能畢業(yè)了嗎 看看reads數(shù)量，然后問問你導(dǎo)師是不是不符合這個物種的規(guī)律
把HiFi reads變成轉(zhuǎn)錄本文件
這一次是主線任務(wù)荣挨，不過要簡單許多男韧。官方幫助：IsoSeq/isoseq-clustering.md at master · PacificBiosciences/IsoSeq...
-
2021-11-16 ggtree包繪制進(jìn)化樹余大神的力作之一：Guangchuang Yu (https://yulab-smu.top[https://yulab-smu.top/])。繪制出的進(jìn)化樹還是很專業(yè)美觀的...
-
兩年前會棒仍，現(xiàn)在不知道了
pasa安裝記錄和ORF提取
官方文檔：https://github.com/PASApipeline/PASApipeline/wiki/Pasa_installation_instructions 之...
-
@sky_info R繪圖我現(xiàn)在也只有草圖，沒有做得很精細(xì)臭胜。等以后調(diào)整好了會擇機發(fā)出來的莫其。
二代群體遺傳與重測序（中）
二.變異結(jié)果注釋與統(tǒng)計 書接上回http://www.reibang.com/p/ab6a35502786[http://www.reibang.com/p/ab6a35...
-
基于sort和uniq的集合操作（簡單明了）
sort:用于將文本文件內(nèi)容加以排序 uniq:刪除文件中重復(fù)的行，得到文件中唯一的行 交集 （用sort將a.txt b.txt文件進(jìn)行排序耸三，uniq使得兩個文件中的行唯一...
