##Bipedal dinosaurs
data(dinosbip)
#Estimates for Tyrannosaurus (FMNH PR 2081 "Sue")
sue<-which(dinosbip$Taxon=="TyrannosaurusFMNH2081")
cQE(dinosbip$FC[sue], equation = "raw", cor = 2) #default correction factor
cQE(dinosbip$FC[sue], equation = "raw", cor = 1.815) #based on eccentricity of the femur
##Estimates for the sample of bipedal dinosaurs using phylogeneteically corrected formula
cQE(dinosbip$FC, equation = "phylocor", cor = dinosbip$cor, data = dinosbip)Run the code above in your browser using DataLab