Re: [R-sig-phylo] ancestral trait reconstruction and OU model

2009-10-26 Thread Marguerite Butler
-- the data which you are trying to fit). Hope this is not too confusing. Take care, Marguerite Butler On Oct 26, 2009, at 8:53 AM, Andrew Hipp wrote: Dear Fabio, In ouch v2 and up, the root state is not treated as a free parameter and is not returned. One way to think about your question might

Re: [R-sig-phylo] Identifying nodes in ouch tree

2011-01-24 Thread Marguerite Butler
Hi Alejandro, The info you want is given in the example script of the dataset bimac, which is included for ouch. To see the help page: require(ouch) data(bimac) ?bimac Then what you need is: tree - with(bimac,ouchtree(node,ancestor,time/max(time),species)) plot(tree,node.names=TRUE) The

Re: [R-sig-phylo] Dealing with Bounded Trait Measures

2011-03-05 Thread Marguerite Butler
Hi David, Liam and everyone, Reflecting traits at boundaries or absorbing them is something that can be done, but I guess I'd like to encourage everyone to think carefully about the interpretation of such simulations. What are you trying to model and what does it mean at the end? Doing these

Re: [R-sig-phylo] Dealing with Bounded Trait Measures

2011-03-05 Thread Marguerite Butler
, absorbing boundaries don't seem like a good model to me. I have never been able to think up a good biological mechanism for a reflecting boundary. Sorry for the extra email from the incomplete thought. M On Mar 5, 2011, at 8:50 PM, Marguerite Butler wrote: Hi David, Liam and everyone

Re: [R-sig-phylo] Dealing with Bounded Trait Measures

2011-03-06 Thread Marguerite Butler
different from BM expectation. Perhaps, the threshold model might be of a lot of use here; I'm going to have to look into that much more. Thanks for you advice, all! -Dave On Sun, Mar 6, 2011 at 12:50 AM, Marguerite Butler mbut...@hawaii.edu wrote: Hi David, Liam and everyone, Reflecting

Re: [R-sig-phylo] multi-state categorical predictor variables in PGLS

2011-03-07 Thread Marguerite Butler
Hi Andrew, Does this sidestep the degrees of freedom problem discussed by Garland et al.? Can anybody point me to references discussing the mechanics of this process and why this is an appropriate thing to do? Others on this list will disagree with me, but it's not a degrees of

Re: [R-sig-phylo] How to detect phylogenetic signal (lambda) in one unscaled trait?

2011-04-01 Thread Marguerite Butler
On Apr 1, 2011, at 7:57 AM, Marguerite Butler wrote: Hi Alberto and Ted, and others, Maybe this is too late to jump into the discussion, but I wanted to add a few comments regarding analysis of the evolution of size and shape. I agree with Ted that there are a number of important

Re: [R-sig-phylo] multiple traits measured within species

2012-03-10 Thread Marguerite Butler
Dear Kaspar, Are your traits the same? Basically, are the seven color patches on the same individual all one trait, or are they seven traits? Do all species have the seven color patches? Or are the number of color patches variable? (this could also be a character). I think before you attempt

Re: [R-sig-phylo] Normality requirement for assessment of lambda with phylosig (phytools) and fitContinuous (geiger)

2012-04-25 Thread Marguerite Butler
Hi Nina and everyone, One thing to consider is that not all zero data are the same. Zeros under a model of continuous trait evolution with a gaussian process as assumed under Brownian motion and OU processes would occasionally cross zero, maybe go negative, etc. For example if you were

Re: [R-sig-phylo] Ouch v.2.8-2 on the same tree gives different results (maybe)

2012-09-14 Thread Marguerite Butler
Dear Elena, If I understand you correctly, you are trying to fit a unique regime to each branch of the tree? This is not advisable, as there are more branches than there are more branches than there are terminal taxa (data points). So you are trying to fit many more parameters than you have

Re: [R-sig-phylo] on evolutionary models

2013-01-24 Thread Marguerite Butler
Dear Pas, I am wondering what the general aim of your study is? I agree with Ted, that you should be careful about your assumptions and the tools you are using. Ted's point that K is a descriptive statistic is an important one. K describes the expected divergence among tip traits assuming a

Re: [R-sig-phylo] Question on ace ML reconstruction of discrete binary character

2013-07-30 Thread Marguerite Butler
Oops. Sorry the citation is Schluter, Price, Mooers, Ludwig 1997. Likelihood of ancestor states in adaptive radiation. Evolution 51:1699-1711. This issue has been known for a long time. On Jul 30, 2013, at 6:51 AM, Marguerite Butler mbutler...@gmail.com wrote: Hi Tom, One thing to keep

Re: [R-sig-phylo] compar.ou

2013-10-29 Thread Marguerite Butler
. Is that problem resolved in these functions? Best, Sandra. 2013/10/26 Marguerite Butler mbutler...@gmail.com Dear Sandra, You might also want to look at the papers that go along with slouch, ouch, and ouwie. Here are some pdfs, along with some tutorials. On my Rcompstart, the ouch stuff

Re: [R-sig-phylo] estimate ancestral state with OUwie models

2018-06-12 Thread Marguerite Butler
FAX: 808-956-4745 http://butlerlab.org http://manoa.hawaii.edu/biology/people/marguerite-butler http://www2.hawaii.edu/~mbutler > On Jun 11, 2018, at 7:33 PM, Simone Blomberg wrote: > > This sounded wrong to me, as the OU process should be agnostic to the > dataset: There are no

Re: [R-sig-phylo] Interpretation of standard errors of parameter estimates in OUwie models

2018-04-04 Thread Marguerite Butler
gt; > > > > -- > > William Gearty > > PhD Candidate, Paleobiology > > Department of Geological Sciences > > Stanford School of Earth, Energy & Environmental Sciences > > williamgearty.com <http://williamgearty.com/> > > > >

Re: [R-sig-phylo] OU for non-ultrametric trees

2018-12-16 Thread Marguerite Butler
hz.ch/mailman/listinfo/r-sig-phylo > Searchable archive at > http://www.mail-archive.com/r-sig-phylo@r-project.org/ > -- Marguerite A. Butler Professor Department of Biology 2538 McCarthy Mall, Edmondson Hall 216 Honolulu, HI 96822 Office:

Re: [R-sig-phylo] model averaging for discrete character evolution

2019-08-09 Thread Marguerite Butler
7803sdata=ogCaep2DN92bk6%2FZwOOcpmHaCPQwX7BzNgKo2UMcWtw%3Dreserved=0 > > > >> > > > [[alternative HTML version deleted]] > > _______ > R-sig-phylo mailing list - R-sig-phylo@r-project.org >

Re: [R-sig-phylo] How to sort trait data according to tree

2021-05-17 Thread Marguerite Butler
is called. > > So you make it work by removing "dat$" in the call to corBrownian(): > > cbm <- corBrownian(1, tree, form = ~Species) > > Best, > > Emmanuel > > - Le 17 Mai 21, à 6:22, Marguerite Butler mbutler...@gmail.com a > écrit : >

Re: [R-sig-phylo] How to sort trait data according to tree

2021-05-17 Thread Marguerite Butler
. Thanks again for your kind replies, Marguerite On Sun, May 16, 2021 at 8:58 PM Emmanuel Paradis wrote: > - Le 17 Mai 21, à 13:41, Marguerite Butler a > écrit : > > Thank you very much for the reply Emmanuel! > > OK, yes I just tried and Iʻm surprised that this

Re: [R-sig-phylo] How to sort trait data according to tree

2021-05-15 Thread Marguerite Butler
archive at > http://www.mail-archive.com/r-sig-phylo@r-project.org/ > -- ________ Marguerite A. Butler Professor Department of Biology 2538 McCarthy Mall, Edmondson Hall 216 Honolulu, HI 96822 Office: 808-956-4713 Dept: 808-956-8617 Lab: 808

Re: [R-sig-phylo] How to sort trait data according to tree

2021-05-16 Thread Marguerite Butler
lculate correlation separately gls(Y ~ X, dat, correlation = cbm) ## works! Thanks, Marguerite On Fri, May 14, 2021 at 10:21 PM Marguerite Butler wrote: > Aloha Oliver, > > From the cor.Brownian help page, the explanation for the form argument is > this: > > a one si

Re: [R-sig-phylo] If my trait X cannot be regressed by body size, how can I rescue residuals corrected by the phylogeny and SE?

2021-05-28 Thread Marguerite Butler
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:/

Re: [R-sig-phylo] units of sigsq

2021-03-21 Thread Marguerite Butler
Biology 2538 McCarthy Mall, Edmondson Hall 216 Honolulu, HI 96822 Office: 808-956-4713 Dept: 808-956-8617 Lab: 808-956-5867 FAX: 808-956-4745 http://butlerlab.org http://manoa.hawaii.edu/biology/people/marguerite-butler http://www2.hawaii.edu/~mbutler [[alternative HTML version deleted]]

Re: [R-sig-phylo] units of sigsq

2021-03-19 Thread Marguerite Butler
partment of Biology 2538 McCarthy Mall, Edmondson Hall 216 Honolulu, HI 96822 Office: 808-956-4713 Dept: 808-956-8617 Lab: 808-956-5867 FAX: 808-956-4745 http://butlerlab.org http://manoa.hawaii.edu/biology/people/marguerite-butler http://www2.hawaii.edu/~mbutler [[alternative HTML version

Re: [R-sig-phylo] issue with reading matrix in R package DispRity

2021-12-06 Thread Marguerite Butler
> 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/ > -- Marguerite