Calculate_Row_Weights {bgafun} | R Documentation |
This will calculate the sequence weights for each group using the Heinkoff and Heinkoff method. Each residue in the sequence is assigned a weight depending on how unique it is in the column. The sequence weight is then the sum of these weights, and the total weight is the number of groups
Calculate_Row_Weights(my_amino,label)
my_amino |
Matrix representation of alignment generated by convert_aln_amino |
label |
Vector or factor that shows the group representation for each sequence in the alignment |
Henikoff, S. and J. G. Henikoff (1994). "Position-based sequence weights." J Mol Biol 243(4): 574-8.
library("bgafun") data(LDH.amino.gapless) data(LDH.groups) LDH.weights=Calculate_Row_Weights(LDH.amino.gapless,LDH.groups) sum(LDH.weights)