data(landsnail, vhabitat, vsample, vtable)
v1 <- as.xcount(vtable, FALSE, segment="mixed")
spec <- xorder(v1, which="species", landsnail, 2)
sampl <- xorder(v1, which="samples", vsample, 1)
vmf <- mefa(v1, sampl, spec)
v2 <- strify(vmf, strata = "site.descr", which = "samples")
v2 ## output is 'mefa'
v3 <- strify(v2, strata = "familia", which = "species")
v3 ## output is 'xcount'
check.attrib(as.xcount(v2), which="samples", vhabitat, 2)
habi <- xorder(as.xcount(v2), which="samples", vhabitat, 2)
vmf2 <- mefa(as.xcount(v2), habi, spec)
vmf2 ## 'mefa' again with attributes
Run the code above in your browser using DataLab