Smatrix {DNAseqtest} | R Documentation |
This function calculates the symmetric matrix S.
Smatrix(s, pix)
s |
a vector of variables containing the six free parameters in the S matrix |
pix |
a vector giving the stationary probabilities for the four nucleotides A, C, G and T |
This function calculates the matrix S, which we used to calculate the rate matrix R.
A 4 x 4 symmetric matrix
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.
Pt, Fmatrix, gn ,gn2
s<-c(.1,.2,.3,.4,.5,.6)
pi<-c(.1,.1,.1,.7)
sm<-Smatrix(s, pi)
sm