view release on metacpan or  search on metacpan

t/01.bigwig.t  view on Meta::CPAN

# does dumping work? 
my $fh = tempfile();
ok($wig->bigWigIntervalDump('I',1,5000,0,$fh));
seek ($fh,0,0);

my $count = 0;
$count++ while <$fh>;
ok($count>1);
undef $fh;

# test chromosome fetching
my $chr_head = $wig->chromList;
ok($chr_head);
my $nodes = 0;
for (my $i=$chr_head->head;$i;$i=$i->next) {
    $nodes++;
}
ok($nodes,7);

my $size = $wig->chromSize('I');
ok($size);
ok($size,$chr_head->head->size);
undef $chr_head;

# testing interval fetching
my $interval_head = $wig->bigWigIntervalQuery('I',1,5000);
ok($interval_head);
$nodes = 0;
for (my $i=$interval_head->head;$i;$i=$i->next) {
    $nodes++;
}
ok($nodes,$count-1);
undef $interval_head;

my $summary = $wig->bigWigSummaryArray('I',1,5000000,bbiSumMean,500);
ok($summary);
ok(ref $summary,'ARRAY');
ok(scalar @$summary,500);

my $should_fail = $wig->bigWigSummaryArray('chromFoo',1,5000000,bbiSumMean,500);
ok(!$should_fail);

my $extended_summary = $wig->bigWigSummary('I',1,5000000,500);
ok ($extended_summary);
ok ($extended_summary->isa('Bio::DB::bbiExtendedSummary'));
ok ($extended_summary->validCount(0) > 0);
ok ($extended_summary->sumData(1)/$extended_summary->validCount(1),$summary->[1]);

my $bin_sum = $wig->bigWigBinStats('I',1,5000000,500);
ok ($bin_sum);
ok (scalar @$bin_sum,500);
ok($bin_sum->[300]->sumData/$bin_sum->[300]->validCount,$summary->[300]);

my $es_list = $wig->bigWigSummaryArrayExtended('I',1,5000000,500);
ok($es_list);
ok(ref $es_list,'ARRAY');
ok(scalar @$es_list,500);
ok(join(' ',sort keys %{$es_list->[0]}),"maxVal minVal sumData sumSquares validCount");
ok($es_list->[300]{sumData}/$es_list->[300]{validCount},$summary->[300]);

my $value = $wig->bigWigSingleSummary('I',1,5000000,bbiSumMean,'NaN');
ok ($value > 0);

1;