hier | R Documentation |
The following methods allow the user to quickly change the hierarchy or population of a genind or genlight object.
hier(x, formula = NULL, combine = TRUE, value) hier(x) <- value
x |
a genind or genlight object |
formula |
a nested formula indicating the order of the population hierarchy to be returned. |
combine |
if |
value |
a formula specifying the full hierarchy of columns in the strata slot. (See Details below) |
You must first specify your strata before you can specify your hierarchies. Hierarchies are special cases of strata in that the levels must be nested within each other. An error will occur if you specify a hierarchy that is not truly hierarchical.
The preferred use of these functions is with a formula
object. Specifically, a hierarchical formula argument is used to name which
strata are hierarchical. An example of a hierarchical formula would
be:
~Country/City/Neighborhood |
This convention was
chosen as it becomes easier to type and makes intuitive sense when defining
a hierarchy. Note: it is important to use hierarchical formulas when
specifying hierarchies as other types of formulas (eg.
~Country*City*Neighborhood
) will give incorrect results.
Zhian N. Kamvar
strata
genind
as.genind
# let's look at the microbov data set: data(microbov) microbov # We see that we have three vectors of different names in the 'other' slot. ?microbov # These are Country, Breed, and Species names(other(microbov)) # Let's set the hierarchy strata(microbov) <- data.frame(other(microbov)) microbov # And change the names so we know what they are nameStrata(microbov) <- ~Country/Breed/Species # let's see what the hierarchy looks like by Species and Breed: hier(microbov) <- ~Species/Breed head(hier(microbov, ~Species/Breed))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.