if (FALSE) {
comm = rbind(c(1,3,0,5,3), c(3,2,5,0,0))
colnames(comm) = c("SpA", "SpB", "SpC", "SpD", "SpE")
rownames(comm) = c("Site 1", "Site 2")
trait = data.frame(body = c(1,2,3,4,4), beak = c(1,5,4,1,2))
rownames(trait) = colnames(comm)
hv = kernel.build(comm[1,], trait)
plot(hv)
kernel.alpha(hv)
hot = kernel.hotspots(hv, 0.5)
plot(hot)
kernel.alpha(hot)
hvlist = kernel.build(comm, trait)
hot = kernel.hotspots(hvlist, 0.1)
kernel.alpha(hot)
}
Run the code above in your browser using DataLab