Hi Lara, You can use the make.era.map function in phytools to slice your tree in simmap like format. Then you multiply the mapped time slices by the required rate (eg, tree$mapped.edge[,1]*sigma) and then use tree$edge.length<-rowSum(tree$mapped.edge) or something like that.
HTH Julien > From: [email protected] > To: [email protected] > Date: Mon, 3 Nov 2014 18:36:34 +0100 > Subject: [R-sig-phylo] Rescaling trees > > Hello R Phylo list > > I have a question regarding rescaling trees. > > "rescale" function in geiger has an option to rescale the > tree in which different time slices have a different rate > of evolution. > > I would like to rescale each time slice with a different > rate but also with a specific alpha parameter (for a > time-dependent OU model). > > Do you know if that would be possible somehow? > > I will greatly appreciate any thoughts and suggestions. > > > > Lara Budic > PhD student at Biometry and Environmental System Analysis > Faculty of Environment and Natural Resources > University of Freiburg > Tennenbacherstr. 4 > 79106 Freiburg i. Br. > > _______________________________________________ > R-sig-phylo mailing list - [email protected] > https://stat.ethz.ch/mailman/listinfo/r-sig-phylo > Searchable archive at http://www.mail-archive.com/[email protected]/ [[alternative HTML version deleted]] _______________________________________________ R-sig-phylo mailing list - [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/[email protected]/
