微生物与转录组关联分析

1. 相关性分析

1.1 相关性热图

test2.correlation.pdf

1.2 相关性网络图

network.pdf

1.3 相关性桑基图

sankey.pdf

2. CCA分析

3. O2PLS分析

流程脚本

python /TJPROJ6/RNA_SH/personal_dir/fengjie/Personal_analysis/Micro_vs_rna_asso/bin/prepareShell.py --help
usage: prepareShell.py [-h] -t TRAN_RESULT -m MICRO_RESULT -c COMPARES -s
                       SAMPLES -s2g SAMPLE2GROUP -g GROUPS -o OUTDIR

Prepare shells for micro vs tran associate analysis

optional arguments:
  -h, --help            show this help message and exit
  -t TRAN_RESULT, --tran_result TRAN_RESULT
                        转录组结果文件路径
  -m MICRO_RESULT, --micro_result MICRO_RESULT
                        微生物结果文件路径
  -c COMPARES, --compares COMPARES
                        比较组合
  -s SAMPLES, --samples SAMPLES
                        样本名称
  -s2g SAMPLE2GROUP, --sample2group SAMPLE2GROUP
                        分组方案,eg:sample1:sample2:sample3,sample4:samp
                        le5:sample6
  -g GROUPS, --groups GROUPS
                        分组名称,与上面分组方案对应
  -o OUTDIR, --outdir OUTDIR
                        shell脚本输出目录

示例:

python /TJPROJ6/RNA_SH/personal_dir/fengjie/Personal_analysis/Micro_vs_rna_asso/bin/prepareShell.py \
    --tran_result /TJPROJ6/RNA_SH/personal_dir/fengjie/Personal_analysis/Micro_vs_rna_asso/test/OSS_DOWN/Back_up/X101SC23105518-Z01-J001/03.Result_X101SC23105518-Z01-J001_bos_mutus/Result_X101SC23105518-Z01-J001_bos_mutus/ \
    --micro_result /TJPROJ6/RNA_SH/personal_dir/fengjie/Personal_analysis/Micro_vs_rna_asso/test/result_X101SC23011361-Z01-F002/ \
    --compares BJ10vsBJG,BJ30vsBJ10,BJ30vsBJG,BJ50vsBJ10,BJ50vsBJG,GZ10vsGZG,GZ30vsGZ10,GZ30vsGZG,GZ50vsGZ10,GZ50vsGZG \
    --samples BJ10_2,BJ30_1,BJ50_1,BJ10_3,BJ30_2,BJ50_2,BJ50_3,BJ30_3,BJG_1,BJG_2,BJG_3,BJ10_1,GZ10_2,GZ30_1,GZ50_1,GZ10_3,GZ30_2,GZ50_2,GZ50_3,GZ30_3,GZG_1,GZG_2,GZG_3,GZ10_1 \
    --sample2group BJ10_2:BJ10_3:BJ10_1,BJ30_1:BJ30_2:BJ30_3,BJ50_1:BJ50_2:BJ50_3,BJG_1:BJG_2:BJG_3,GZ10_2:GZ10_3:GZ10_1,GZ30_1:GZ30_2:GZ30_3,GZ50_1:GZ50_2:GZ50_3,GZG_1:GZG_2:GZG_3 \
    --groups BJ10,BJ30,BJ50,BJG,GZ10,GZ30,GZ50,GZG \
    --outdir /TJPROJ6/RNA_SH/personal_dir/fengjie/Personal_analysis/Micro_vs_rna_asso/bin/log \