How to subset COMADRE Animal Matrix Database

subsetComadre.R

#How to subset the COMADRE Animal Matrix Database (works with COMPADRE too)

#Create a vector of indices to retain
subsetID <- which(comadre$metadata$Order %in% c("Monotremata","Didelphimorphia","Paucituberculata","Microbiotheria","Dasyurormorphia","Peramelemorphia","Diprotodontia"))

#Subset entire COMADRE object to JUST the above subset
#First make a copy

subset.comadre <- comadre

#Now subset the sub-parts of the database
#metadata
subset.comadre$metadata <- mammal.comadre$metadata[subsetID,]

#matrices
subset.comadre$mat <- mammal.comadre$mat[subsetID]

#matrix class information
subset.comadre$matmatrixClass <- mammal.comadre$matrixClass[subsetID]

This is now superceded by the subsetDB function in the Rcompadre package