[Mattice-commits] r121 - pkg/inst/doc

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Tue Dec 30 08:09:19 CET 2008


Author: andrew_hipp
Date: 2008-12-30 08:09:19 +0100 (Tue, 30 Dec 2008)
New Revision: 121

Modified:
   pkg/inst/doc/maticce.Rnw
Log:
adding analysis details to vignette

Modified: pkg/inst/doc/maticce.Rnw
===================================================================
--- pkg/inst/doc/maticce.Rnw	2008-12-30 06:45:38 UTC (rev 120)
+++ pkg/inst/doc/maticce.Rnw	2008-12-30 07:09:19 UTC (rev 121)
@@ -118,5 +118,23 @@
 \caption{Simulated character on \code{ovales.tree} at model-averaged theta values, with coloring according to \code{ou2}}
 \end{figure}
 
+What is the relative support for the brownian motion vs. OU-2 model? We find which model has a change only at node 2 by inspecting the regime matrix
 
+<<ha, fig = FALSE, echo = TRUE>>=
+ha.4.2[['regMatrix']][['overall']]
+@
+
+Then we can find the likelihood of these models on a given tree:
+
+<<haLnl, fig = FALSE, echo = TRUE>>=
+ha.4.2[['hansens']][[1]] 
+@
+
+or the information criterion weights:
+
+<<haWeights, fig=FALSE, echo=TRUE>>=
+summary(ha.4.2)[['modelsMatrix']][[1]]
+@
+
+
 \end{document}



More information about the Mattice-commits mailing list