From Code to Community: Sponsoring The Perl and Raku Conference 2025 Learn more

DBAURAIN / Bio-MUST-Core-0.250380 / bin / resample-ali.pl 

            

              
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
—
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143

              
#!/usr/bin/env perl
# PODNAME: resample-ali.pl
# ABSTRACT: Resample ALI files using (variable length) bootstrap or jackknife
use Modern::Perl '2011';
use autodie;
use File::Basename;
use Getopt::Euclid qw(:vars);
use Path::Class qw(dir file);
use POSIX;
use Smart::Comments;
use Bio::MUST::Core;
use Bio::MUST::Core::Utils qw(change_suffix append_suffix);
use aliased 'Bio::MUST::Core::Ali';
use aliased 'Bio::MUST::Core::SeqMask::Weights';
# setup resampling method
my $method = $ARGV_resampling . '_masks';
# setup replicate numbering format
my $field = ceil( log($ARGV_replicates) / log(10) );
for my $infile (@ARGV_infiles) {
    ### Processing: $infile
    my $ali = Ali->load($infile);
    # create output directory named after input file and resampling mode
    my $outdir = append_suffix(
        change_suffix($infile, q{}), "-$ARGV_resampling"
    );
    for my $width (@ARGV_width) {
        # create output subdirectory named after width
        my $subdir = dir( $outdir, 'width-' . $width )->relative;
        ### REPDIR: $subdir->stringify
        $subdir->mkpath();
        # generate resampling masks for Ali
        my @masks = Weights->$method(
            $ali, { 'rep' => $ARGV_replicates, 'width' => $width }
        );
        # apply masks to generate pseudo-replicates
        for my $rep (0..$ARGV_replicates-1) {
            my $new_ali = $masks[$rep]->filtered_ali($ali);
            my $filename = sprintf 'replicate-%0*d.ali', $field, $rep;
            ### REP file: $filename
            my $outfile = file($subdir, $filename);
            $new_ali->store($outfile);
        }
    }
}
__END__
=pod
=head1 NAME
resample-ali.pl - Resample ALI files using (variable length) bootstrap or jackknife
=head1 VERSION
version 0.250380
=head1 USAGE
    resample-ali.pl <infiles> [optional arguments]
=head1 REQUIRED ARGUMENTS
=over
=item <infiles>
Path to input ALI files [repeatable argument].
=for Euclid: infiles.type: readable
    repeatable
=back
=head1 OPTIONAL ARGUMENTS
=over
=item --resampling=<mode>
Mode of resampling to perform on the ALI [default: bootstrap]. The following
modes are available: bootstrap and jackknife.
=for Euclid: mode.type:       string, mode eq 'bootstrap' || mode eq 'jackknife'
    mode.type.error: <mode> must be one of bootstrap or jackknife (not mode)
    mode.default:    'bootstrap'
=item --rep[licates]=<n>
Number of pseudo-replicates [default: 100].
=for Euclid: n.type: int >= 0
        n.default: 100
=item --width=<n>...
Number of positions to be included in each pseudo-replicate [default: 1]. When
specified as a fraction between 0 and 1 (included), it is interpreted as
relative to the width of the ALI. Thus 10000 would mean 10,000 sites whereas
0.5 would mean half the width of the ALI. Hence, the default is to generate
pseudo-replicates of the same width as the ALI. Multiple values can be
provided if they are whitespace-separated.
=for Euclid: n.type: number
        n.default: [ 1 ]
=item --version
=item --usage
=item --help
=item --man
Print the usual program information
=back
=head1 AUTHOR
Denis BAURAIN <denis.baurain@uliege.be>
=head1 COPYRIGHT AND LICENSE
This software is copyright (c) 2013 by University of Liege / Unit of Eukaryotic Phylogenomics / Denis BAURAIN.
This is free software; you can redistribute it and/or modify it under
the same terms as the Perl 5 programming language system itself.
=cut