## Example from Alvarez (2017)
Wetlands_veg <- make_cocktail(Wetlands, Wetlands_veg, cover = "percen")
summary(as.factor(Wetlands_veg@header$Syntax))
## Same but only for two vegetation units
Wetlands_veg <- make_cocktail(Wetlands, Wetlands_veg,
which = c("HY1", "HY2"), cover = "percen"
)
summary(as.factor(Wetlands_veg$Syntax))
## Construct the 'shaker' object anew
Wetlands <- new("shaker")
## Set a pseudo-species
Wetlands <- set_pseudo(Wetlands, Wetlands_veg, c(
"Cyperus latifolius",
"Cyperus exaltatus"
))
## Set a species group
Wetlands <- set_group(Wetlands, Wetlands_veg,
group_id = "Cyperus papyrus",
group = c(
"Cyperus papyrus",
"Cyclosorus interruptus",
"Lepistemon owariense"
)
)
## Set a fromula
Wetlands <- set_formula(Wetlands, Wetlands_veg,
formula_id = "HE1",
formula = "groups:'Cyperus papyrus' | species:'Cyperus papyrus > 50'"
)
## Summaries
summary(Wetlands)
summary(Wetlands, Wetlands_veg)
Run the code above in your browser using DataLab