The Perl Toolchain Summit needs more sponsors. If your company depends on Perl, please support this very important event.

NAME

bc_measure_distance - Measure beta-diversity between communities

SYNOPSIS

  bc_measure_distance -input_files   communities.generic      \
                      -dist_type     'hellinger'              \
                      -output_prefix community_distance

DESCRIPTION

This script reads files containing biological communities and calculate how dissimilar they are (beta diversity). The output is a tab-delimited file containing the beta diversity between all pairs of communities. The columns of this file contain the name of the first community, of the second community, and their beta diversity, respectively. See Bio::Community::Meta::Beta for more details. Note that beta diversity metrics are based on relative abundances. Hence, any weight you provide will affect the results.

REQUIRED ARGUMENTS

-if <input_files>... | -input_files <input_files>...

Input file containing the communities to manipulate. When providing communities in a format that supports only one community per file (e.g. gaas), you can provide multiple input files.

OPTIONAL ARGUMENTS

-wf <weight_files>... | -weight_files <weight_files>...

Tab-delimited files containing weights to assign to the community members.

-wa <weight_assign> | -weight_assign <weight_assign>

When using a files of weights, define what to do for community members whose weight is not specified in the weight file (default: weight_assign.default):

* $num : assign to the member the arbitrary weight $num provided

* average : assign to the member the average weight in this file.

* ancestor: go up the taxonomic lineage of the member and assign to it the weight of the first ancestor that has a weight in the weights file. Fall back to the 'average' method if no taxonomic information is available for this member (for example a member with no BLAST hit).

-op <output_prefix> | -output_prefix <output_prefix>

Path and prefix for the output files. Default: output_prefix.default

-dt <dist_type> | -dist_type <dist_type>

The type of distance or beta-diversity metric to calculate: 1-norm, euclidean (2-norm), hellinger, infinity-norm, bray-curtis... Default: dist_type.default

-pf <pair_files>... | -pair_files <pair_files>...

Input file specifying the pairs of communities for which to calculate the beta-diversity. Each line of a file of pairs should have the name of two communities, separated by a tab. For each file of pairs, there will be a corresponding output file.

FEEDBACK

Mailing Lists

User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to one of the Bioperl mailing lists.

Your participation is much appreciated.

  bioperl-l@bioperl.org                  - General discussion
  http://bioperl.org/wiki/Mailing_lists  - About the mailing lists

Support

Please direct usage questions or support issues to the mailing list:

bioperl-l@bioperl.org

rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible.

Reporting Bugs

Report bugs to the Bioperl bug tracking system to help us keep track the bugs and their resolution. Bug reports can be submitted via the web:

  http://bugzilla.open-bio.org/

AUTHOR - Florent Angly

Email florent.angly@gmail.com