## Example from Alvarez (2017)
Wetlands_veg@header <- make_cocktail(Wetlands, Wetlands_veg, cover="percen")
summary(as.factor(Wetlands_veg@header$Syntax))
## Same but only for two vegetation units
Wetlands_veg@header <- 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