require(pavo)
require(dplyr)
require(phytools)
require(geiger)
require(phylolm)
require(MuMIn)
require(ggplot2)
require(viridis)

Abbreviated methods

Caveats

  • Not everyone measured all the same body parts. I think most of them are standard, but some species have more or less body parts measured
  • samples per species per plumage patch ranged from 1 to 5
  • individuals measured per species varied a lot too
load('output/processedata.rdata')
plot(colspace(visdat), col=spec2rgb(alldatp, alpha=0.2))
## Warning: Quantum catch are not relative, and have been transformed

plot(xyzdatr, col=spec2rgb(alldatp, alpha=0.2), theta=275, arrow.p=3.5)