DNAseqtest (version 1.0)

Distance: Paralinear Distances

Description

This function calculates the paralinear distance between K matched DNA sequences.

Usage

Distance(F4)

Arguments

F4

a $4^K$ array containing the joint distribution array $F(t)$ or the observed array N

Value

Details

This function calculates the paralinear distances between K matched DNA sequences, depending on the joint distribution array for these K sequences or on the observed divergence array N.

References

Faisal Ababneh, Lars S Jermiin, Chunsheng Ma, John Robinson (2006). Matched-pairs tests of homogeneity with applications to homologous nucleotide sequences. Bioinformatics, 22(10), 1225-1231.

See Also

gn2, gn, Fmatrix, Ntml

Examples

merge2<-matrix(c(-1,-4,-3,2,-2,-5,1,3), 4, 2)
theta<-c(rep(.25,3), rep(.25,3), rep(.25,3), c(.2,.35,.79,.01,.93,.47), 3,.1,.5,.8)
F1<-gn(theta,merge2)
dn<-Distance(F1)
dn