Estimate the A.C.E. model from a phyloseq object

Estimates the A.C.E. model of heritabily (with confidence intervals) from a phyloseq object using the twinlm function.

calculate_ACE(ps, zygosity, dz, family_id, verbose = TRUE)

Arguments

ps	A phyloseq object that contains sample_data.
zygosity	The name of the column in sample_data() that contains the zygosity variable (either MZ or DZ).
dz	The level in the zygosity variable corresponding to dyzogitic twins.
family_id	The name of the column in sample_data() that contains the family_id variable. In other words, both members of a twin pair should have the same value in this column.
verbose	Logical. If FALSE, messages will be supressed.

Value

This function returns a data.frame with A.C.E. values (and confidence intervals) for each ASV.

Examples

df <- calculate_ACE(ps, zygosity = "Zygosity", dz = "DZ", family_id = "Family_ID", verbose = TRUE)
#> Error in calculate_ACE(ps, zygosity = "Zygosity", dz = "DZ", family_id = "Family_ID",     verbose = TRUE): object 'ps' not found