data(elbaue)
clust <- vector('integer', nrow(elbaue.env))
clust[elbaue.env$MGL < -50 & elbaue.env$SDGL < 50] <- 1
clust[elbaue.env$MGL < -50 & elbaue.env$SDGL >= 50] <- 2
clust[elbaue.env$MGL >= -50 & elbaue.env$SDGL >= 50] <- 3
clust[elbaue.env$MGL >= -50 & elbaue.env$SDGL < 50] <- 4
# syntab(elbaue, clust, limit=30)
levels(clust) <- c('dry.ld','dry.hd', 'wet.hd','wet.ld')
syntab(elbaue, clust, limit=30, mupa=TRUE)Run the code above in your browser using DataLab