Hi, We have a recent paper on bias caused by intraspecific variation/measurement error in phylogenetic regression setups and whether one should correct for it and if so how :
T. F. Hansen and K. Bartoszek. Interpreting the evolutionary regression: the interplay between observational and biological errors in phylogenetic comparative studies. Syst Biol., 61(3): 413-425, 2012. with R code to do it : GLSME.R http://www.math.chalmers.se/~krzbar/GLSME/GLSME.R Hope this will be useful for you! Cheers Krzysztof Bartoszek PhD Student Mathematical Statistics Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg ________________________________________ From: r-sig-phylo-boun...@r-project.org [r-sig-phylo-boun...@r-project.org] on behalf of r-sig-phylo-requ...@r-project.org [r-sig-phylo-requ...@r-project.org] Sent: Thursday, September 27, 2012 12:00 PM To: r-sig-phylo@r-project.org Subject: R-sig-phylo Digest, Vol 56, Issue 20 Send R-sig-phylo mailing list submissions to r-sig-phylo@r-project.org To subscribe or unsubscribe via the World Wide Web, visit https://stat.ethz.ch/mailman/listinfo/r-sig-phylo or, via email, send a message with subject or body 'help' to r-sig-phylo-requ...@r-project.org You can reach the person managing the list at r-sig-phylo-ow...@r-project.org When replying, please edit your Subject line so it is more specific than "Re: Contents of R-sig-phylo digest..." Today's Topics: 1. Lack of fit between regression lines and data && phylogenetic regression with covariate (Fabio de A. Machado) ---------------------------------------------------------------------- Message: 1 Date: Wed, 26 Sep 2012 13:34:00 -0300 From: "Fabio de A. Machado" <macfa...@gmail.com> To: r-sig-phylo@r-project.org Subject: [R-sig-phylo] Lack of fit between regression lines and data && phylogenetic regression with covariate Message-ID: <157d92f5-f6ef-440f-8cbd-2836b26be...@gmail.com> Content-Type: text/plain Dear list, There is something bothering me in a recent analysis that I am performing. I am currently investigating the association between a morphological variable with an ecological one through gls( ), using a size covariate. I am performing a separate analysis for both sexes. I have 17 species in my sample. I've found a significant interaction term between the ecological variable and size, and this result is consistent between sexes. However, when I visually inspect the interaction term (using 'effects' package) I find that, for one sex, the regression lines passes through the middle of my data cloud, but for the other, the lines seems to miss the majority of the points. This might be a visualization problem, but I would appreciate if anyone can pinpoint any obvious issue with this analysis. Aditionally, I would like to know if there is any way to perform a phylogenetic regression with size as a covariate using intraspecific variation. I have over 1000 specimens measured and maybe this could solve my problems. Best regards, ----------------------------------------------------- Fabio Andrade Machado Laborat?rio de Evolu??o de Mam?feros Departamento de Gen?tica e Biologia Evolutiva Universidade de S?o Paulo Rua do Mat?o, trav.14, n?277, Edificio Minas Gerais 05422-970 f.mach...@usp.br ; macfa...@gmail.com +55 11 82631029 Lattes: http://lattes.cnpq.br/3673327633303737 Google Scholar: http://scholar.google.com/citations?hl=en&user=2l6-VrQAAAAJ [[alternative HTML version deleted]] ------------------------------ _______________________________________________ R-sig-phylo mailing list R-sig-phylo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-phylo End of R-sig-phylo Digest, Vol 56, Issue 20 _______________________________________________ R-sig-phylo mailing list R-sig-phylo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-phylo