{
### Data preparation ###
data(testotu)
groupinformation <- data.frame(
group = c(rep("a", 10), rep("b", 10)),
factor1 = rnorm(10),
factor2 = rnorm(mean = 100, 10),
subject = factor(c(1:10, 1:10)),
group2 = c(rep("e", 5), rep("f", 5), rep("e", 5), rep("f", 5))
)
### Packaging metafile, community data, and taxonomy table ###
test_object <- tax_summary(
groupfile = groupinformation,
inputtable = testotu[, 2:21],
reads = TRUE,
taxonomytable = testotu[, c(1, 22)]
)
### Object configuration ###
test_object_plan1 <- object_config(
taxobj = test_object,
treat_location = 1,
rep_location = 4
)
### Facet configuration ###
test_object_plan2 <- object_config(
taxobj = test_object,
treat_location = 1,
rep_location = 4,
facet_location = 5
)
}
Run the code above in your browser using DataLab