DESeq2如何使用

這篇文章主要介紹“DESeq2如何使用”的相關(guān)知識，小編通過實際案例向大家展示操作過程，操作方法簡單快捷，實用性強，希望這篇“DESeq2如何使用”文章能幫助大家解決問題。

deseq_analysis.r 差異基因分析DESeq2

使用方法:

$Rscript $scriptdir/deseq_analysis.r -h
usage: /work/my_stad_immu/scripts/deseq_analysis.r [-h] -i filepath -m
                                                   filepath -t treatname
                                                   --control CONTROL --case
                                                   CASE [-f fdr] [-c fc]
                                                   [-s size] [-a alpha]
                                                   [-X x.lab] [-Y y.lab]
                                                   [-T title] [-H height]
                                                   [-W width] [-o path]
                                                   [-p prefix]

DESeq2 analysis : https://www.億速云.com/article/1500

optional arguments:
  -h, --help            show this help message and exit
  -i filepath, --input filepath
                        input read count file [required]
  -m filepath, --metadata filepath
                        metadata file , required
  -t treatname, --treatname treatname
                        treat colname in group file, required
  --control CONTROL     set control group name required
  --case CASE           set case group name required
  -f fdr, --fdr fdr     set fdr threshold [default 0.05]
  -c fc, --fc fc        set fold change threshold [default 2]
  -s size, --size size  point size [optional, default: 0.7]
  -a alpha, --alpha alpha
                        point transparency [0-1] [optional, default: 1]
  -X x.lab, --x.lab x.lab
                        the label for x axis [optional, default: log2FC]
  -Y y.lab, --y.lab y.lab
                        the label for y axis [optional, default: -log10(FDR)]
  -T title, --title title
                        the label for main title [optional, default: Volcano]
  -H height, --height height
                        the height of pic inches [default 5]
  -W width, --width width
                        the width of pic inches [default 5]
  -o path, --outdir path
                        output file directory [default
                        /work/my_stad_immu/05.enrich]
  -p prefix, --prefix prefix
                        out file name prefix [default Volcano]

參數(shù)說明：

-i 輸入基因表達矩陣文件，必須為count表達文件：

-k 輸入基因表達矩陣文件，為fpkm 或者 tpm文件 用于相關(guān)性分析：

-m metadata文件路徑,樣本的分組信息，第一列必須和表達文件的樣本名稱對應(yīng)：

-t subtype.hclust   --case S1 --control  S2  ： 指定metadata 分組列名，分組里面的比較組名字 ，如果分組名字有空格，應(yīng)該用引號引起來：  “Stage IA”

--fdr 0.01 --fc 2  設(shè)置差異基因的篩選條件： 顯著性和差異倍數(shù)

使用舉例：

Rscript $scriptdir/deseq_analysis.r  -i ../01.TCGA_download/TCGA-STAD_gene_expression_Counts.tsv \
  -k ../01.TCGA_download/TCGA-STAD_gene_expression_TPM.tsv --fdr 0.01 --fc 2 \
  -m ../03.TIME/metadata.group.tsv -t subtype.hclust   --case S1 --control  S2 -p S1_vs_S2

關(guān)于“DESeq2如何使用”的內(nèi)容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業(yè)相關(guān)的知識，可以關(guān)注億速云行業(yè)資訊頻道，小編每天都會為大家更新不同的知識點。