Evaluate how reliable the distinction is between two archaeological spatial units (layers or any other type of spatial units).
frag.layers.admixture(graph, layer.attr, morphometry, x, y, z, verbose=TRUE)A numeric vector with the admixture of each pair of spatial units.
An undirected igraph object.
Character. The name of the vertex attribute giving the spatial unit of each fragment (e.g. layer).
Character. Optional, to pass to the `frag.edges.weighting` function: name of the vertex attribute with the morphometric value of the fragments.
Character. Optional, to pass to the `frag.edges.weighting` function: name of the vertex attribute with the "x" coordinate of the fragments.
Character. Optional, to pass to the `frag.edges.weighting` function: name of the vertex attribute with the "y" coordinate of the fragments.
Character. Optional, to pass to the `frag.edges.weighting` function: name of the vertex attribute with the "z" coordinate of the fragments.
Logical. Whether to print or not warning messages.
Sebastien Plutniak <sebastien.plutniak at posteo.net>
This function returns a value reflecting the robustness of the distinction between two spatial units (layers or any other relevant archaeological spatial unit). The admixture value is computed as:
1 - cohesion(layer 1) - cohesion(layer 2)
The admixture of two spatial units is equal to the cohesion of a virtual third unit, which is defined by the fragments and the connection relationships intersecting the two spatial units.
Results range in [0;1] with 0 for two completely independent spatial units and values towards 1 as the robustness of the boundary between the two spatial units is lower.
As it appears, this function calls the frag.layers.cohesion function.
The basic use of this function is intended for a graph with two layers, whose edges have been previously weighted using the frag.edges.weighting function (an error message is displayed if the vertice attribute "layer" contains less than two spatial units, and if the graph does not have an edge attribute "weight").
However, the function can also be used for a graph with more than two spatial units. In this case, a subgraph is generated for each pair of layers (using the frag.get.layers.pair function), the frag.edges.weighting function is applied to weight their edges (a warning message is displayed), and the admixture is computed.
An error message is displayed if the vertex "layer" attribute has less than two values.
frag.edges.weighting,
frag.layers.cohesion
g <- frag.simul.process(n.components=20, vertices=50, disturbance = .15)
g <- frag.edges.weighting(g, layer.attr="layer")
frag.layers.admixture(g, "layer")
Run the code above in your browser using DataLab