步驟1：數(shù)據(jù)質(zhì)控： fastp代碼如下殿托，保存為fastp.sh 執(zhí)行fastp命令伞广，批量處理fastq文件 步驟2：call SNP戳吝，使用GATK的流程genome文件構(gòu)建索...
發(fā)簡(jiǎn)信
IP屬地：香港
-
2024-06-20 QTL-seq從原始數(shù)據(jù)開始的流程
-
2024-03-02從gff文件中獲取多區(qū)域基因ID
注意，region.xls和 final_sorted_with_promoter.change.gene文件均為bed格式急鳄。bed格式前3列為必須有的，分別為染色體乐导，起始位...
-
2023-10-20尋找文件并壓縮
將當(dāng)前目錄下大于200Mb典阵，且沒有壓縮的文件均壓縮 這個(gè)命令使用find命令來查找當(dāng)前目錄下大于200MB且不以".gz"結(jié)尾的文件奋渔。然后使用-exec選項(xiàng)將找到的文件傳遞給...
-
2023-02-20更改vcf.gz 格式文件名稱
如果你想替換特定的樣本名稱，你可以創(chuàng)建一個(gè)包含舊名稱和新名稱的文件壮啊，然后使用bcftools reheader命令進(jìn)行替換嫉鲸。這個(gè)文件的格式應(yīng)該是每行一個(gè)舊名稱和一個(gè)新名稱，中...
-
二代群體遺傳與重測(cè)序（中）二.變異結(jié)果注釋與統(tǒng)計(jì) 書接上回http://www.reibang.com/p/ab6a35502786[http://www.reibang.com/p/ab6a35...
-
Cytoscape繪制相關(guān)網(wǎng)絡(luò)圖導(dǎo)讀 用R軟件corr.test函數(shù)進(jìn)行兩組數(shù)組的相關(guān)性分析也榄，用Cytoscape繪制相關(guān)網(wǎng)絡(luò)圖。 一司志、輸入數(shù)據(jù) 代謝數(shù)據(jù)：metabo.txt 16S數(shù)據(jù)：genus.tx...
-
annovar注釋除人類以外SNP
快速注釋已經(jīng)得到的vcf文件 所需文件 ref.fa gff3文件或者gtf vcf文件 簡(jiǎn)單流程 1.建立一文件夾 2. gff3文件轉(zhuǎn)變格式 gff3ToGenePred...
-
-
-
2023-02-21 安裝agricolae報(bào)錯(cuò) non-zero exit status
安裝 R 包agricolae報(bào)錯(cuò)：installation of package ‘a(chǎn)gricolae’ had non-zero exit status 往上查看時(shí) 發(fā)現(xiàn)...
-
2023-02-20 SnpSift 安裝及使用
SnpSift 可以分割染色體，截取部分vcf文件激才。具體使用方法見：https://blog.csdn.net/weixin_43569478/article/details...
-
您好 PCA.R 程序哪里有呢瘸恼？謝謝劣挫。是自己寫的嗎？
重測(cè)序分析（10）群體進(jìn)化分析之PCA分析
PCA PCA(Principal Components Analysis)即主成分分析东帅，也稱主分量分析或主成分回歸分析法压固，是一種無監(jiān)督的數(shù)據(jù)降維方法。首先利用線性變換靠闭，將數(shù)...
-
-
理解OPLS-DA分析的原理和結(jié)果圖OPLS-DA：正交篇最小二乘法判別分析 原理： OPLS-DA不同于PCA，它是一種有監(jiān)督的判別分析統(tǒng)計(jì)方法媚朦。運(yùn)用偏最小二乘回歸建立代謝物表達(dá)量與樣本類別之間的關(guān)系模型捡絮，來...
-
GATK流程生成多態(tài)性的 SNP/InDel
目的：提取群體中多態(tài)性較好的SNP和InDel。將GVCF文件進(jìn)行合并莲镣，并提取SNP和InDel，然后使用vcftools獲得在群體中多態(tài)性較好的SNP和InDel涎拉。 第一步...
-
2022-06-21EXCEL批量轉(zhuǎn)pdf
import osimport xlwings as xw def excel2pdf(file):app = xw.App(visible=True,add_book=Fa...
