data(landsnail, vhabitat, vsample, vtable)
vt <- as.xcount(vtable, FALSE)
spec <- xorder(vt, which="species", landsnail, 2)
sampl <- xorder(vt, which="samples", vsample, 1)
vmf <- mefa(vt, sampl, spec)
v2 <- strify(vmf, strata = "site.descr", which = "samples")
habi <- xorder(as.xcount(v2), which="samples", vhabitat, 2)
vmf2 <- mefa(as.xcount(v2), habi, spec)
report.mefa("vspec-report.txt", vmf2, "species", c(3,2,4:7), "spec.name", "order", 0,
FALSE, FALSE, sep = c(", "," (",")","; "))
report.mefa("vsampl-report.txt", vmf2,"samples", c(3,2,4:7), "spec.name", "order", 0,
FALSE, FALSE, sep = c(", "," (",")",",\n"))
report.mefa("vspec-report.tex", vmf2,"species", c(3,2,4:7), "spec.name", "order", "author",
TRUE, FALSE, sep = c(", "," (",")","; "))
report.mefa("vsampl-report.tex", vmf2,"samples", c(3,2,4:7), "spec.name", "order", "author",
TRUE, TRUE, sep = c(", "," (",")",",\n"))
### Sweave example
mefadocs("SampleReport")
Run the code above in your browser using DataLab