4

我正在尝试生成包含 10 个具有随机质量分数的随机序列的FASTQ 文件。我最初使用以下代码，效果很好：

my @seq  = (rand_fa_seq() for ^10);
my @qual = (rand_qual()   for ^10);

@seq.perl.say;
@qual.perl.say;

sub rand_fa_seq
{
    return join("", roll(20,"ACGT".comb));

}

sub rand_qual
{
    return join("", roll(20,"EFGHIJ".comb))
}

但是，为了进一步简化它，我想也许我可以从右边的语句中删除括号。@seq当我这样做时，我在and中只有一个元素@qual。

my @seq  = rand_fa_seq() for ^10;
my @qual = rand_qual()   for ^10;

@seq.perl.say;
@qual.perl.say;

sub rand_fa_seq
{
    return join("", roll(20,"ACGT".comb));

}

sub rand_qual
{
    return join("", roll(20,"EFGHIJ".comb))
}

这是一个错误还是它应该表现的方式？没有括号，这是一个标量上下文吗？Great List Refactor会改变这种行为吗？

我见过这种行为的 Perl 6 解释器版本：

MoarVM：

perl6 版本 2015.03-204-g8578022 基于 MoarVM 版本 2015.03-60-g36d56f7

虚拟机：

perl6 版本 2015.03-305-ga95107d 基于 JVM java 版本“1.7.0_79”
OpenJDK 运行时环境 (rhel-2.5.5.1.el7_1-x86_64 u79-b14)
OpenJDK 64 位服务器 VM（构建 24.79-b02，混合模式）

4

2 回答 2

5

语句修饰符是 Perl 语法的标志之一（参见 Perl 5 的perldoc），它们在向版本 6 的过渡中幸存下来。

这意味着您的陈述等同于

my @seq;
for ^10 { @seq = rand_fa_seq() }

my @qual;
for ^10 { @qual = rand_qual() }

即，您连续 10 次分配新值（并且只有最后一个存在）。

请注意，这也可以更简洁地写成

my @seq  = rand_fa_seq() xx 10;
my @qual = rand_qual() xx 10;

另请注意，默认情况下，sub 中的最后一条语句提供返回值。

假设您可能需要长度不是 20 的序列，您可以对该值进行参数化，最终得到以下结果：

sub rand-fa-seq($n = 20) { <A C G T>.roll($n).join }
sub rand-qual($n = 20)   { <E F G H I J>.roll($n).join }

my @seq  = rand-fa-seq() xx 10;
my @qual = rand-qual() xx 10;

我使用引号而不是拆分字符串。

于 2015-05-16T15:41:56.627 回答
2

第二个示例类似于行尾的经典 perl5 控制语句。

像

say "true" if something();
say "yeah" for ^10;

如果您将打印语句放入其中，rand_qual()您会注意到它仍在执行 10 次。

于 2015-05-16T13:02:30.350 回答