Tracing the likelihood calculation of a Gaussian model2 years ago
Example tree and data | A mixed Gaussian model | Likelihood values for the three example variants | Tracing the likelihood calculation using the function PCMLikTrace | Variant 1 | Variant 2 | Variant 3 | A step by step description of the log-likelihood calculation | Step 1: Calculating $\vec | Calculating $\vec{\omega}$, $\mathbf{\Phi}$ and $\mathbf{V}$ for a node in an OU regime | Calculating $\vec{\omega}$, $\mathbf{\Phi}$ and $\mathbf{V}$ for a node in a BM regime | Step 2: Calculating $\mathbf{A}$, $\vec{b}$, $\mathbf{C}$, $\vec{d}$, $\mathbf{E}$ and $f$ for tips and internal nodes | Step 3: Calculating $\mathbf{L}$, $\vec{m}$, $r$ for the internal nodes and the root | Step 4: Calculating the log-likelihood value using $\mathbf{L}{0}$, $\vec{m}{0}$, and $r_{0}$. | References
PCMBase 1.2.15Venelin MitovPCMTracePruning.Rmd