tag:blogger.com,1999:blog-8499895524521663926.post713891082425502465..comments2024-03-27T07:13:39.236-04:00Comments on Phylogenetic Tools for Comparative Biology: Negative estimated variances from the likelihood surfaceLiam Revellhttp://www.blogger.com/profile/04314686830842384151noreply@blogger.comBlogger4125tag:blogger.com,1999:blog-8499895524521663926.post-15757919316760468912011-09-06T15:33:43.035-04:002011-09-06T15:33:43.035-04:00Thanks Rafael.Thanks Rafael.Liam Revellhttps://www.blogger.com/profile/04314686830842384151noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-25194759061756310312011-09-04T11:58:35.622-04:002011-09-04T11:58:35.622-04:00in case anyone runs into the same problem, in my c...in case anyone runs into the same problem, in my case it was because I was using the hessian() function from the maxLik package, and not from the numDeriv package as Liam suggested in his post (I ran evol.vcv and got an error that the function wasn't found, so I searched ??hessian and the hit I got was for the maxLik package). Once I loaded the numDeriv package, everything ran smoothly. Thanks!Rafael Maiahttps://www.blogger.com/profile/02395303011366399867noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-46241085042551836372011-09-03T09:39:14.040-04:002011-09-03T09:39:14.040-04:00Hi Rafael.
I have not been able to reproduce your...Hi Rafael.<br /><br />I have not been able to reproduce your error. If you send me your data and tree I will try with that (liam.revell@umb.edu).<br /><br />- LiamLiam Revellhttps://www.blogger.com/profile/04314686830842384151noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-76428179003316182282011-09-02T13:42:45.227-04:002011-09-02T13:42:45.227-04:00hey Liam,
sorry I cant help, and instead I'm...hey Liam, <br /><br />sorry I cant help, and instead I'm bringing you another issue... I am trying to use the evol.vcv function and I keep getting the error:<br /><br />Error in x$hessian : object of type 'closure' is not subsettable<br /><br />if I use v0.1 it runs ok... any ideas? Thanks!Rafael Maiahttps://www.blogger.com/profile/02395303011366399867noreply@blogger.com