[R] mcmcsamp and variance ratios
Martin Henry H. Stevens
hstevens at muohio.edu
Wed Jan 3 19:44:09 CET 2007
I have assumed that ratios of variance components (Fst and Qst in
population genetics) could be estimated using the output of mcmcsamp
(the series on mcmc sample estimates of variance components).
What I have started to do is to use the matrix output that included
the log(variances), exponentiate, calculate the relevant ratio, and
apply either quantile or or HPDinterval to get confidence intervals.
This seems too simple but I can't think of what is wrong with it.
All thoughts appreciated.
Dr. Hank Stevens, Assistant Professor
338 Pearson Hall
Oxford, OH 45056
Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
"E Pluribus Unum"
More information about the R-help