Method
consensus_z
Classpackage
G::Seq::Consensus;
Function
Analyzes a consensus through Z-values.
Description
Consensus with Z value greater than int high_z (default is 1) is
capitalized, and consensus with Z value less than int low_z
(default is 0.2) is shown as '-'.
pointer array_of_sequences is a reference to an array of sequences. The sequences must be of same lengths.
eg. my @seq = ("atgc","ctaa","tttt","cttg"); &consensus_z(\@seq);
returns 'cTt-'
Requirements
applications
SutOpt.pm
Statistics::Descriptive;
applications
none
data
none
Parameters
essential
first pointer array_of_sequences Reference to an array of sequences.
option
high int z_value greater than which is significant default=1
low int z_value less than which is insignificant default=0.2
filename string out filename default='stdout'
output string [g or f or show] default='g'
Returns
first string consensus
Usage
(string consensus) = &consensus_z(pointer array_of_sequences, high=>1, low=>0.2, output=>'show',
filename=>'consensus_z.gif');
Input file
none
Output file
none
References
none
Comments
none
Authors
This method was written by Kazuharu Arakawa(gaou@g-language.org) at Institute for Advanced Biosciences Keio Univ.
History
Written 13 July 2001 - Kazuharu Arakawa