## Fatty Acids
data(FAset)
fa.set = as.vector(unlist(FAset))
## prey.cluster requires full prey database.
data(preyFAs)
prey.sub=(preyFAs[,4:(ncol(preyFAs))])[fa.set]
prey.sub=prey.sub/apply(prey.sub,1,sum)
group=as.vector(preyFAs$Species)
prey.matrix=cbind(group,prey.sub)
prey.cluster(prey.matrix,method="average",dist.meas=3)
Run the code above in your browser using DataLab