[Ouwie-commits] help with OUwie

[Jeremy Beaulieu]

Hi Ari~

Alter the following line:

model <- as.numeric(unlist(OUwie(sim10[[i]], data, model = c("OUM"), simmap.tree = TRUE)))

to:

model<-OUwie(sim10[[i]], data, model = c("OUM"), simmap.tree = TRUE))

Now to obtain the MLE of theta you would do:

model$theta[1,]

That should work.

All the best,
Jeremy

__________________________
Jeremy M. Beaulieu
Postdoctoral Fellow
National Institute for Mathematical
and Biological Synthesis (NIMBioS)
University of Tennessee
1122 Volunteer Blvd, Ste. 106
Knoxville, TN 37996
www.jeremybeaulieu.org