Re: [R-sig-phylo] obtaining reasonable values from OUCH (Krzysztof Bartoszek)

2013-07-03 Thread Krzysztof Bartoszek
Hi, For multiple co-evolving traits you can also look at the mvSLOUCH package, http://www.math.chalmers.se/~krzbar/mvSLOUCH/mvSLOUCH.html Bartoszek et. al., A phylogenetic comparative method for studying multivariate adaptation, Journal of Theoretical Biology, 314:204-215, 2012. However you can a

Re: [R-sig-phylo] obtaining reasonable values from OUCH

2013-07-01 Thread Pascal Title
Thanks Marguerite for your reply! To answer your questions, as for where I heard about trait variance approaching sigma2/(2*alpha), I had posted to r-sig-phylo about similar problems a little over a year ago, and Carl Boettiger had mentioned this as a way to see if the parameter estimates were in

Re: [R-sig-phylo] obtaining reasonable values from OUCH

2013-07-01 Thread Marguerite Butler
Good morning Pascal, We (Cressler, King, and myself) have a paper soon to be submitted that shows that model selection is very robust, but parameter estimates, esp. alpha and sigma are very difficult to estimate, because they have high variance. Although the bias is manageable, it is very har

Re: [R-sig-phylo] obtaining reasonable values from OUCH

2013-06-29 Thread Brian O'Meara
No, it's just univariate (I was responding to your example about doing a single character). For multivariate OU, I think OUCH is the only game in town. Best, Brian ___ Brian O'Meara Assistant Professor Dept. of Ecology & Evolutionary Biology U. of Tennessee, Kn

Re: [R-sig-phylo] obtaining reasonable values from OUCH

2013-06-29 Thread Pascal Title
Thanks for the reply, Brian. Does OUwie implement the multivariate model? >From the documentation, I don't see any mention of applying it to multiple traits at the same time. cheers, -Pascal On Sat, Jun 29, 2013 at 11:32 PM, Brian O'Meara wrote: > You're at best on the border of how many param

Re: [R-sig-phylo] obtaining reasonable values from OUCH

2013-06-29 Thread Brian O'Meara
You're at best on the border of how many parameters you can estimate given the size of your dataset. One thing you might try is running OUwie, which also implements that model but uses a different optimizer and starting values, though I wouldn't be surprised if you get similar answers. Trying a gri

[R-sig-phylo] obtaining reasonable values from OUCH

2013-06-29 Thread Pascal Title
Hello, Given a phylogeny with 29 tips, and 5 traits that each fit a OU model better than a BM model of trait evolution (according to fitContinuous in GEIGER), I would like to fit a multivariate OU model to these 5 traits using the hansen() function in the OUCH package in R. The goal is to get mult