data(bimac)
attach(bimac)
tree.plot(ancestor,time,names=species,regimes=OU.LP)
brown.fit(log(size),ancestor,time/max(time))
hansen.fit(log(size),ancestor,time/max(time),OU.1)
hansen.fit(log(size),ancestor,time/max(time),OU.3)
hansen.fit(log(size),ancestor,time/max(time),OU.4)
hansen.fit(log(size),ancestor,time/max(time),OU.LP)
Run the code above in your browser using DataLab