require(pavo)
require(dplyr)
require(phytools)
require(geiger)
require(phylolm)
require(MuMIn)
require(ggplot2)
require(viridis)
Abbreviated methods
- Combined the data from Dale, Dunn, and myself
- Ran through a visual model assuming the average UV visual system and blue tit double cone
- Calculated \(\Delta\)S assuming a Weber fraction of 0.1 and receptor relative densities of 1:2:2:4
- Converted to Cartesian coordinates where Euclidean distances are in JND units
- considering only male color. In a few species, the sex of the only individual sampled was unknown; since these were monochromatic species I considered them males for the analyses.
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)