Dear Bruno, Can you provide a worked example? By the way, when looking at the code it seems that the simulated BM values are constrained to be positive (l. 30). Is it normal?
All the best, Julien Note: the sigma value computed is the REML estimate. If you need the ML use sig2 <- sum(pic(mean.perf, tree)^2)/length(mean.perf) instead. ---------------------------------------- > From: bruno.freder...@ulg.ac.be > Date: Tue, 1 Sep 2015 23:15:21 +0200 > To: r-sig-phylo@r-project.org > Subject: [R-sig-phylo] phylo.ANOVA.intra > > Dear all, > > I am using the function phylo.ANOVA.intra developed by Hernandez-Lopez et al. > 2014. PLoS ONE. 9(3): e89832, allowing the use of individual data. I have > compared groups using various variables and the tests are running well with > most variables. > However, for two of them (both are Log-transformed data), I got the following > error message: > > Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : > NA/NaN/Inf in ‘y' > > It looks like the problem is coming from the data... probably > log-transformation. > > Does someone have experience with this test? > Do you have tips to solve this problem? > > Thank you in advance for your help! > All the best, > > Bruno > > ---------------------------------------- > Dr. Bruno Frédérich > Postdoctoral Researcher F.R.S.-FNRS > Lab. of Functional and Evolutionary Morphology > University of Liège, Institut de Chimie (B6c) > 4000 Liège, BELGIUM > > Phone: 003243665133 > Fax: 003243663715 > E-mail: bruno.freder...@ulg.ac.be <mailto:bruno.freder...@ulg.ac.be> > Website: http://www2.ulg.ac.be/morfonct/ <http://www2.ulg.ac.be/morfonct/> > ---------------------------------------- > Open Repository and Bibliography (ORBi) http://orbi.ulg.ac.be/ > <http://orbi.ulg.ac.be/> > Personal Publications (Open Access) http://tinyurl.com/Bruno-Reef > <http://tinyurl.com/Bruno-Reef> > > > > > > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-phylo mailing list - R-sig-phylo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-phylo > Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/ _______________________________________________ R-sig-phylo mailing list - R-sig-phylo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
