Dear Nick, Although within-species error is not treated in the same way as in Rphylopars you can also impute missing values (and estimate ancestral states) for most multivariate models implemented in mvMORPH using the "estim" function (see ?estim).
Best wishes, Julien De : R-sig-phylo <r-sig-phylo-boun...@r-project.org> de la part de Liam J. Revell <liam.rev...@umb.edu> Envoyé : lundi 20 juillet 2020 18:47 À : Nick Crouch <ncro...@uchicago.edu>; r-sig-phylo@r-project.org <r-sig-phylo@r-project.org> Objet : Re: [R-sig-phylo] Rphylopars changing observed data values Dear Nick. I don't know what's going on with Rphylopars, but phytools also has a multivariate phylogenetic imputation function, phylo.impute. It's described here on my blog: http://blog.phytools.org/2019/03/phylogenetic-imputation-with.html (although, of course you should use the phytools version not the posted code). All the best, Liam Liam J. Revell Associate Professor, University of Massachusetts Boston Profesor Asistente, Universidad Católica de la Ssma Concepción web: http://faculty.umb.edu/liam.revell/, http://www.phytools.org Academic Director UMass Boston Chile Abroad: https://www.umb.edu/academics/caps/international/biology_chile On 7/16/2020 10:47 AM, Nick Crouch wrote: > [EXTERNAL SENDER] > > Hi all, > > I'm trying to use Rphylopars to impute some missing data, but am having some > strange results with my empirical data. > > The phylopars should return a `data.frame` of the original data with the > missing (NA) values calculated given the phylogenetic data and evolutionary > model. However, my known values are being changed in the output. For example: > >> head(empirical.data) > species Trait > 1 Species_1 14.0 > 2 Species_2 12.0 > >> head(p_BM$anc_recon) > species Trait > 1 Species_1 12.0 > 2 Species_2 8.9 > > The function works fine when I use different combinations of simulated data: > * simulated traits + simulated phylogeny works fine > * simulated traits + empirical phylogeny works fine > * empirical data + simulated phylogeny works fine > > The phylogeny has no 0 branch lengths, is binary, and ultrametric. I have > tried log-transforming the data. I have also tried sequentially dropping one > tip at a time from the tree and re-running to see if any single tip is > causing the problem. Different models also have the same effect. > > If anyone has any thoughts i'd appreciate it. > > Nick > > -- > > Nick Crouch > https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fnickcrouch.github.io%2F&data=02%7C01%7Cliam.revell%40umb.edu%7Cfb9cae70c8884b194ba308d82cc914ca%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637308589448752383&sdata=hXkX5aboOGy3cJaK6QZBYAnm3kNrs%2Fqj2hAocX9MTKU%3D&reserved=0 > Postdoctoral scholar, Jablonski lab > The University of Chicago > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-phylo mailing list - R-sig-phylo@r-project.org > https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-phylo&data=02%7C01%7Cliam.revell%40umb.edu%7Cfb9cae70c8884b194ba308d82cc914ca%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637308589448752383&sdata=VO6CC1WW3FuZg3geh0vZAL7NADKfHuxXxftuMW%2FTic4%3D&reserved=0 > Searchable archive at > https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.mail-archive.com%2Fr-sig-phylo%40r-project.org%2F&data=02%7C01%7Cliam.revell%40umb.edu%7Cfb9cae70c8884b194ba308d82cc914ca%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637308589448762381&sdata=skztD9PMN5AVLUr03a1TXg2iNjWL9Aw5NvA3aYCvoyg%3D&reserved=0 > _______________________________________________ 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/
