embxyz.c


Function embXyzSeqsetNR

Reads a list of AjPSeq's and writes an array describing the redundancy in the list. Elements in the array correspond to sequences in the list, i.e. the array[0] corresponds to the first sequence in the list.

Sequences are classed as redundant (0 in the array, i.e. they are possibly to be discarded later) if they exceed a threshold (%) level of sequence similarity to any other in the set (the shortest sequence of the current pair will be discarded).

The set output will always contain at least 1 sequence.

Synopsis

Prototype
AjBool embXyzSeqsetNR (
      const AjPList input,
      AjPInt* keep,
      ajint* nset,
      const AjPMatrixf matrix,
      float gapopen,
      float gapextend,
      float thresh
);

TypeNameRead/WriteDescription
const AjPListinputInputList of ajPSeq's
AjPInt*keepOutput0: rejected (redundant) sequence, 1: the sequence was retained
ajint*nsetOutputNumber of sequences in nr set (no. of 1's in the keep array)
const AjPMatrixfmatrixInputResidue substitution matrix
floatgapopenInputGap insertion penalty
floatgapextendInputGap extension penalty
floatthreshInputThreshold residue id. for "redundancy"
AjBool RETURNajTrue on success

Input
input:(Input)List of ajPSeq's
matrix:(Input)Residue substitution matrix
gapopen:(Input)Gap insertion penalty
gapextend:(Input)Gap extension penalty
thresh:(Input)Threshold residue id. for "redundancy"
Output
keep:(Output)0: rejected (redundant) sequence, 1: the sequence was retained
nset:(Output)Number of sequences in nr set (no. of 1's in the keep array)
Returns
AjBool:ajTrue on success

Description

Reads a list of AjPSeq's and writes an array describing the redundancy in the list. Elements in the array correspond to sequences in the list, i.e. the array[0] corresponds to the first sequence in the list.

Sequences are classed as redundant (0 in the array, i.e. they are possibly to be discarded later) if they exceed a threshold (%) level of sequence similarity to any other in the set (the shortest sequence of the current pair will be discarded).

The set output will always contain at least 1 sequence.

Usage

See source code

Example

In preparation

Errors

See source code

Dependencies

See source code

See Also

See other functions in this section

Availability

In release 5.0.0

Function embXyzSeqsetNRRange

Reads a list of AjPSeq's and writes an array describing the redundancy in the list. Elements in the array correspond to sequences in the list, i.e. the array[0] corresponds to the first sequence in the list.

Sequences are classed as redundant (0 in the array, i.e. they are possibly to be discarded later) if they lie outside a range of threshold (%) sequence similarity to others in the set (the shortest sequence of the current pair will be discarded).

Synopsis

Prototype
AjBool embXyzSeqsetNRRange (
      const AjPList input,
      AjPInt* keep,
      ajint* nset,
      const AjPMatrixf matrix,
      float gapopen,
      float gapextend,
      float thresh1,
      float thresh2
);

TypeNameRead/WriteDescription
const AjPListinputInputList of ajPSeq's
AjPInt*keepOutput0: rejected (redundant) sequence, 1: the sequence was retained
ajint*nsetOutputNumber of sequences in nr set (no. of 1's in the keep array)
const AjPMatrixfmatrixInputResidue substitution matrix
floatgapopenInputGap insertion penalty
floatgapextendInputGap extension penalty
floatthresh1InputThreshold lower limit
floatthresh2InputThreshold upper limit
AjBool RETURNajTrue on success

Input
input:(Input)List of ajPSeq's
matrix:(Input)Residue substitution matrix
gapopen:(Input)Gap insertion penalty
gapextend:(Input)Gap extension penalty
thresh1:(Input)Threshold lower limit
thresh2:(Input)Threshold upper limit
Output
keep:(Output)0: rejected (redundant) sequence, 1: the sequence was retained
nset:(Output)Number of sequences in nr set (no. of 1's in the keep array)
Returns
AjBool:ajTrue on success

Description

Reads a list of AjPSeq's and writes an array describing the redundancy in the list. Elements in the array correspond to sequences in the list, i.e. the array[0] corresponds to the first sequence in the list.

Sequences are classed as redundant (0 in the array, i.e. they are possibly to be discarded later) if they lie outside a range of threshold (%) sequence similarity to others in the set (the shortest sequence of the current pair will be discarded).

Usage

See source code

Example

In preparation

Errors

See source code

Dependencies

See source code

See Also

See other functions in this section

Availability

In release 5.0.0