Hi everybody,
Many thanks for all your inputs. After some digging, I think that MCMCglmm
is a good way to go. However, I am an illiterate in Bayesian modelling and
will probably have to seek for advice elsewhere, like r-sig-mixed-model.
Cheers,
Nicolay
On Wed, Aug 1, 2018 at 10:28 AM Masahito
Hi list,
I second Jon's suggestion for using phylogenetic mixed model. This seems to be
the closest to what Nicolay is after.
Yet another option is "SLOUCH" package in R, which can handle within-species
variation while performing phylogenetic corrections. Here are some example
codes.
You will want to try Bayesian mixed modeling which can handle your species
# problem by using “species” as a group level (i.e. “random”) effect. I
would recommend r packages brms and MCMCglmm and their vignettes. Both have
detailed “phylo” examples.
Maybe check out the r-sig-mixed-model listserv
In R the PhyloPars package might be interesting. Haven’t used it... though.
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2703881/
Cheers Franz
Sent from my mobile
> On 31. Jul 2018, at 21:04, Nicolay Cunha wrote:
>
> Dear list members,
>
> I would like to test a linear model between a
Two ways to go about this are here:
Ives, A. R., P. E. Midford, and T. Garland, Jr. 2007. Within-species
variation and measurement error in phylogenetic comparative methods.
Systematic Biology 56:252–270.
Felsenstein, J. 2008. Comparative methods with sampling error and
within-species variation: