#为 ATAC 和 RNA seq 接入号对制作 1 个 library.csvdef read_rna_accession(): with open('rna_accessions.txt') as f: samples = [sample for sample in f.read().split('\n')...

#Make 1 library.csv for a pair of ATAC and RNA seq accession numbers
def read_rna_accession():
    with open('rna_accessions.txt') as f:
        samples = [sample for sample in f.read().split('\n') if len(sample) > 0]  # Remove empty lines
        return samples

def read_atac_accession():
    with open('atac_accessions.txt') as f:
        samples = [sample for sample in f.read().split('\n') if len(sample) > 0]  # Remove empty lines
        return samples

# Read ATAC and RNA accession IDs
atac_SRRs = read_atac_accession()
rna_SRRs = read_rna_accession()


# Define all rule for generating libraries.csv
rule all:
    input:
        expand("{atac_srr}_{rna_srr}_libraries.csv", atac_srr=atac_SRRs, rna_srr=rna_SRRs)

# Rule to create libraries.csv
rule create_libraries_csv:
    output:
        "{atac_srr}_{rna_srr}_libraries.csv"
    run:
        atac_srr = wildcards.atac_srr
        rna_srr = wildcards.rna_srr
        with open(output[0], "w") as f:
            f.write("fastqs,sample,library_type\n")
            f.write(f"atac_seq/{atac_srr},{atac_srr},Chromatin Accessibility\n")
            f.write(f"rna_seq/{rna_srr},{rna_srr},Gene Expression\n")

您好，我正在尝试为 atac_seq 和 rna_seq 编号对 {SRRXXXXXX} 创建一个 libraries.csv 文件。目前，我有 2 个 atac_seq 的 SRR 和 2 个 rna_seq 的 SRR；我想为每一对生成一个 libraries.csv 文件，因此最终应该会得到 2 个 {atac_srr}_{rna_srr}libraries.csv 文件。

使用我拥有的代码，它会为 atac_seq 和 rna_seq 接入号之间的所有可能组合生成它；我想按顺序生成它，这意味着 atac 接入号列表中的第一个 SRR 对应于 rna 接入号列表中 SRR 中的第一个元素；有没有可能做到这一点？

这是文档相关部分的链接，其中提到了这一点：snakemake.readthedocs.io/en/stable/snakefiles/…

您可以按照 Swifty 的建议将其作为 for 循环来执行。或者如果您想保留扩展函数：

expand("{atac_srr}_{rna_srr}_libraries.csv",zip, atac_srr=atac_SRRs, rna_srr=rna_SRRs)

您是否在 IDEA 中配置了 SDK，以便它知道您将 Coretto 放在哪里？