Question

我遇到了一个较大脚本的一小部分问题。我正在尝试将文件名存储在数组中作为哈希值，以便以后能够列出它们。下面是我构建散列％extramapped的循环：

for (my $i=0; $i<scalar(@refgenomes); $i++) {
            my $extraRefPath = $refgenomes[$i];
            my @extraFieldRef = split("/", $extraRefPath);
            my @extraFullRefName = split("\\.", $extraFieldRef[-1]);
            my $extraRefName = $extraFullRefName[0];
            my $extraTempOutFile = $gzFile.$extraRefName.".sort.bam";
            my $extraCmd = "$bwa bam2bam -n 0.01 -o 2 -l 16500 -g $extraRefPath <($fixpair -o /dev/stdout $tempOutFile) | samtools sort -o /dev/stdin ".$gzFile.".sort"." >$extraTempOutFile";
            push(@listTargets, $extraTempOutFile);
            push(@{$extramapped{$extraRefName}}, $extraTempOutFile);
            $stringMakeToPrint = $stringMakeToPrint.$extraTempOutFile.": ".$tempOutFile."\n\t$extraCmd\n\n";

当我使用Data :: Dumper检查哈希的内容时，它看起来像我期望的那样：

＆GT;

 $VAR1 = {
>           'rCRS' => [
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L001_R1_001.fastq.gzrCRS.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L001_R1_002.fastq.gzrCRS.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L001_R1_003.fastq.gzrCRS.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L002_R1_001.fastq.gzrCRS.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L002_R1_002.fastq.gzrCRS.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L002_R1_003.fastq.gzrCRS.sort.bam'
>                     ],
>           'Copy_NC' => [
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L001_R1_001.fastq.gzCopy_NC.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L001_R1_002.fastq.gzCopy_NC.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L001_R1_003.fastq.gzCopy_NC.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L002_R1_001.fastq.gzCopy_NC.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L002_R1_002.fastq.gzCopy_NC.sort.bam',
> '..//Sample_LP30/LP30_CCTAGGT-ACCTAGG_L002_R1_003.fastq.gzCopy_NC.sort.bam'
>                        ]
>         };

但是，我很难让这些文件名重新出现。我虽然可以通过join（“”，@ {extramapped {$ extraRefName}}）来做到这一点，但是它的行为并不像我期望的那样。当我只是要求它打印列表时：

my $extraFileList = join(" ", @{$extramapped{$extraRefName}});
print "This is the extra file list $extraFileList\n";

我得到了预期的清单：

> This is the extra file list ..//Sample_LP32/LP32_GCTCGAA-ACTTCAA_L001_R1_001.fastq.gzCopy_NC.sort.bam ..//Sample_LP32/LP32_GCTCGAA-ACTTCAA_L001_R1_002.fastq.gzCopy_NC.sort.bam ..//Sample_LP32/LP32_GCTCGAA-ACTTCAA_L002_R1_001.fastq.gzCopy_NC.sort.bam ..//Sample_LP32/LP32_GCTCGAA-ACTTCAA_L002_R1_002.fastq.gzCopy_NC.sort.bam

但是当我尝试将它合并到这样的循环中时：

for (my $i=0; $i<scalar(@refgenomes); $i++) {
    my $extraFileList = join(" ", @{$extramapped{$extraRefName}});
    my $extraMmfCmd = "$samtools merge $extraFileList | $retrieve /dev/stdin /dev/stdout | $rmdup -c -o $extraMapFileName /dev/stdin";

我收到以下错误，循环死亡。

> Can't use an undefined value as an ARRAY reference at line 152.

我不明白为什么在这种情况下输出是不同的。有谁可以向我解释一下？建议如何修复它？

干杯，安娜

Answer 1

@{...} 取消引用数组引用。所以，如果它在这里爆炸：

my $extraFileList = join(" ", @{$extramapped{$extraRefName}});

然后显然是$extramapped{$extraRefName}未将数组引用存储到取消引用的结果。

要么，

$extraRefName不再引用有效密钥
或，$extramapped{$extraRefName}被其他地方删除了。

有些人采用了使用空数组引用或来避免此错误的做法：@{ $hash{ $keyNotThere } || [] }。虽然，这在阵列创建中往往有点浪费。你真正想要的是一个更麻烦的事情：

ref( $aref ) eq 'ARRAY' ? @$aref : ();

（当然可能是Perl向导已经优化了空数组解除引用的情况。）

Perl - 尝试访问存储为哈希值的数组时出错

1 个答案: