国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 編程 > Perl > 正文

perl生成特定堿基比例的隨機序列的代碼

2020-06-04 20:29:56
字體:
來源:轉載
供稿:網友

方法一(不使用模塊,by agonyr)

代碼如下:
#!/usr/bin/perl -w

use strict;

my @seq = ( "A", "T", "C", "G" );
my $length = 10000;

undef my %hash;
$hash{"A"} = int( $length * 0.3 );
$hash{"C"} = int( $length * 0.3 );
$hash{"G"} = int( $length * 0.2 );
$hash{"T"} = int( $length * 0.2 );

my $i = 0;
while ( $i 《 $length ) {
    my $word = $seq[ rand(@seq) ];
    if ( $hash{$word} ) {
        print "$word";
        $i++;
    }
    $hash{$word}--;
}
print "n";

方法二(使用模塊,by yixf)

代碼如下:
#!/usr/bin/perl

use strict;
use warnings;

use BioX::SeqUtils::RandomSequence;

my $randomizer = BioX::SeqUtils::RandomSequence-》new(
    {
        l =》 10000,
        s =》 1,
        y =》 "dna",
        a =》 3,
        c =》 3,
        g =》 2,
        t =》 2
    }
);
print $randomizer-》rand_seq(), "n";

兩種方法比較

設定長度為10000,ACGT的比例為3:3:2:2。
代碼如下:
withoutModule Length=10000 GC=49.42% A=2558,C=2503,G=2439,T=2500,Others=0
withModule Length=10000 GC=50.00% A=3000,C=3000,G=2000,T=2000,Others=0

發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 苍山县| 桃园市| 三河市| 北川| 陆良县| 达尔| 颍上县| 肇东市| 奉化市| 大理市| 西宁市| 潮安县| 郓城县| 萍乡市| 阜宁县| 荃湾区| 花莲县| 西华县| 东港市| 大埔区| 中阳县| 安仁县| 栖霞市| 和静县| 海林市| 陆河县| 鄂温| 镇巴县| 沁水县| 桦南县| 福安市| 勃利县| 扬州市| 筠连县| 齐齐哈尔市| 肇州县| 石渠县| 尚义县| 新田县| 社旗县| 牟定县|