├── DIR1│   ├── smp1.fastq.gz│   ├── smp1_fastqc/│   ├── smp2.fastq.gz│   └── smp2_fastqc/└── DIR2    ├── smp3.fastq.gz    ├── smp3_fastqc/    ├── smp4.fastq.gz    └── smp4_fastqc/我想按樣本計(jì)算讀取次數(shù)，然后按目錄連接所有計(jì)數(shù)。我創(chuàng)建一個(gè)字典，將樣本 1 和 2 鏈接到目錄 1，將樣本 3 和 4 鏈接到目錄 2DIRS,SAMPLES = glob_wildcards(INDIR+'/{dir}/{smp}.fastq.gz')# Create samples missingdef filter_combinator(combinator, authlist):    def filtered_combinator(*args, **kwargs):        for wc_comb in combinator(*args, **kwargs):            if frozenset(wc_comb) in authlist:                yield wc_comb    return filtered_combinator# Authentificationcombine_dir_samples = []for dir in DIRS:    samples, = glob_wildcards(INDIR+'/'+dir+'/{smp}.fastq.gz')    for smp in samples:        combine_dir_samples.append( { "dir" : dir, "smp" : smp} )       combine_dir_samples = { frozenset( x.items() ) for x in combine_dir_samples }dir_samples = filter_combinator(product, combine_dir_samples)然后，我創(chuàng)建一個(gè)規(guī)則來按樣本計(jì)算我的讀取次數(shù)rule all:    input:        expand(INDIR+'/{dir}/{smp}_Nreads.txt', dir_samples, dir=DIRS, smp=SAMPLES)rule countReads:    input:        INDIR+'/{dir}/{smp}_fastqc/fastqc_data.txt'    output:        INDIR+'/{dir}/{smp}_Nreads.txt'    shell:        "grep 'Total\ Sequences' {input} | awk '{{print {wildcards.dir},$3}}' > {output}"---------------------------------------------------------------# result ok├── DIR1│   ├── smp1_Nreads.txt│   └── smp2_Nreads.txt└── DIR2    ├── smp3_Nreads.txt    └── smp4_Nreads.txt> cat smp1_Nreads.txtDIR1 15082186我嘗試為我的連接規(guī)則使用不同的輸入語(yǔ)法expand(OUTFastq+'/{dir}/FastQC/{{smp}}_Nreads.txt', dir_samples, dir=DIRS)lambda wildcards: expand(OUTFastq+'/{dir}/FastQC/{wildcards.smp}_Nreads.txt', dir_samples, dir=DIRS, smp=SAMPLES)expand(OUTFastq+'/{dir}/FastQC/{wildcards.smp}_Nreads.txt', dir_samples, dir=DIRS, smp=SAMPLES)我沒有找到任何解決方案，就像它不關(guān)心我的字典中的這條規(guī)則一樣。