[Mesquitelist] S-value

Monica Guerrini mguerrini at biologia.unipi.it
Mon Feb 2 00:14:56 PST 2009


Hi!
Thank you very much for your answers to my 
question some days ago! Now I have resolved my 
first problem, that is importing and using a gene 
tree reconstructed with other programs! It is very easy!!!!
But now I have another two questions.
First of all if I measure the distributions of 
test statistics (Slatkin & Maddison's S) on a 
simulated population tree, when I reject or 
accept this model drawn? I know that if the 
observed S.-value of my gene tree falls outside 
of the 95% confidence interval, that model are 
not supported by the data. On the contrary, if I 
have a case similar to the example of Mesquite's 
manual I don't reject the model (S = 9, observed 
values within the 95% confidence interval, that is 4 and 10).
But now I suppose this example:
I measure the observed S-value for my first gene 
tree and that is 9 (for example...!).
Now I draw some population trees corresponding to 
the evolution models I want to test. I measure 
the distribution of S-values for each model.
model n° 1: the S-value distribution has only one 
value, for example 2, so I reject the model (outside the boundary)
model n°2: the S-value distribution has two 
limits: 5 and 15, so I accept the model (inside the boundary)
model n°3: the S-value distribution has only one 
value, that is 9 (exactly my observed S-value)      ??????
model n°4: the S-value distribution has two 
limits: 9 and 15 (lower tail exactly my observed S-value)      ??????
What happen for models n° 3 and 4? I'm not sure 
to accept or reject these hypothesis.......

The second question is:
For use a gene tree (reconstructed with other 
software) and to measure S-value it is better to 
use a sequences tree or an haplotypes tree?

Thank you very very much again for your help!!!!!


Monica








More information about the Mesquitelist mailing list