It is difficult to work just with the str(MM).
No libraries are needed. You could try something like:
# input: df - output: subset by plot and some other criteria
f <- function(x){
x.splitted <- split(x, x$plot)
lapply(x.splitted, function(y) y[y$Especie == 24 &
y$Especie == 28,])
}
#
MM.subset <- lapply(MM, f)
Optionally you may like to have all data consolidated in a single data.frame
# In case you prefer a data.frame instead of a list of data.frames
do.call(rbind, MM.subset)