接上一篇我进行修改完成了一个批量处理过程

1.将文件导入到新的路径并重新命名

cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR3.txt 522chr3.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR3.txt 680chr3.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR3.txt forecastchr3.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR4.txt 467chr4.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR4.txt 522chr4.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR4.txt 680chr4.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR4.txt forecastchr4.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR5.txt 467chr5.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR5.txt 522chr5.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR5.txt 680chr5.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR5.txt forecastchr5.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR7.txt 467chr7.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR7.txt 522chr7.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR7.txt 680chr7.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR7.txt forecastchr7.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR8.txt 467chr8.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR8.txt 522chr8.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR8.txt 680chr8.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR8.txt forecastchr8.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR9.txt 467chr9.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR9.txt 522chr9.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR9.txt 680chr9.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR9.txt forecastchr9.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR10.txt 467chr10.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR10.txt 522chr10.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR10.txt 680chr10.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR10.txt forecastchr10.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR11.txt 467chr11.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR11.txt 522chr11.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR11.txt 680chr11.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR11.txt forecastchr11.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR12.txt 467chr12.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR12.txt 522chr12.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR12.txt 680chr12.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR12.txt forecastchr12.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR13.txt 467chr13.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR13.txt 522chr13.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR13.txt 680chr13.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR13.txt forecastchr13.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR14.txt 467chr14.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR14.txt 522chr14.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR14.txt 680chr14.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR14.txt forecastchr14.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR15.txt 467chr15.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR15.txt 522chr15.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR15.txt 680chr15.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR15.txt forecastchr15.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR17.txt 467chr17.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR17.txt 522chr17txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR17.txt 680chr17.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR17.txt forecastchr17.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR18.txt 467chr18.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR18.txt 522chr18.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR18.txt 680chr18.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR18.txt forecastchr18.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR19.txt 467chr19.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR19.txt 522chr19.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR19.txt 680chr19.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR19.txt forecastchr19.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/467_ind/CHR20.txt 467chr20.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/522_ind/CHR20.txt 522chr20.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/04.try/680_ind/CHR20.txt 680chr20.txt
cp /home/DL_2/CUTech/comput_2/Tech_No1/group4/yjc/workspace/SCT_data/var/01.combine_gvcf/04.deep_allele_maf/02.ROH/02.genetic_load/01.ROH_num_len/test/forecast/CHR20.txt forecastchr20.txt

2.执行批量化绘图脚本

import pandas as pd
import matplotlib.pyplot as plt

def plot_roh_distribution(chromosome_number):
    # 根据染色体编号构造文件名
    file_paths = [
        f'467chr{chromosome_number}.txt',
        f'522chr{chromosome_number}.txt',
        f'680chr{chromosome_number}.txt',
        f'forecastchr{chromosome_number}.txt'
    ]

    # 加载数据
    data_list = []
    for file_path in file_paths:
        data = pd.read_csv(file_path, delim_whitespace=True)
        individual_id = file_path.split('chr')[0]
        data_list.append((data, f"Individual {individual_id}"))

    # 初始化图形
    fig, axes = plt.subplots(nrows=4, ncols=1, figsize=(10, 8))

    # 在每个子图上绘制每个个体的ROH分布
    for ax, (data, title) in zip(axes, data_list):
        for _, row in data.iterrows():
            ax.bar((row['POS1'] + row['POS2']) / 2, 1, width=(row['POS2'] - row['POS1']),
                   align='center', color='#4c924c', edgecolor='black', linewidth=0.5)
        ax.set_title(title)
        ax.set_xlim(0, data['POS2'].max() + 100000)
        ax.set_xlabel('Position along chr{} (bp)'.format(chromosome_number))
        ax.set_yticks([])

    plt.tight_layout()
    plt.savefig(f'ROH_Distribution_chr{chromosome_number}.pdf')
    plt.close()  # 关闭图形，避免内存泄漏

# 定义要处理的染色体列表
chromosomes = [1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 12, 13, 14, 15, 18, 19, 20]

# 批量处理所有染色体
for chromosome in chromosomes:
    plot_roh_distribution(chromosome)

最后编辑于：2024.04.19 13:33:37

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 199,271评论 5赞 466
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 83,725评论 2赞 376
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 146,252评论 0赞 328
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 53,634评论 1赞 270
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 62,549评论 5赞 359
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 47,985评论 1赞 275
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 37,471评论 3赞 390
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 36,128评论 0赞 254
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 40,257评论 1赞 294
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 35,233评论 2赞 317
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 37,235评论 1赞 328
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 32,940评论 3赞 316
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 38,528评论 3赞 302
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 29,623评论 0赞 19
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 30,858评论 1赞 255
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 42,245评论 2赞 344
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 41,790评论 2赞 339

接上一篇我进行修改完成了一个批量处理过程

1.将文件导入到新的路径并重新命名

2.执行批量化绘图脚本

推荐阅读更多精彩内容