frag.simul.compare {archeofrag}

From an observed fragmentation graph, simulates two series of graphs corresponding to two deposition hypotheses.

Description

Given an observed fragmentation graph, simulates two series of graphs corresponding to two deposition hypotheses, compares their properties and returns a summary table.

Usage

frag.simul.compare(graph, layer.attr, iter, summarise=TRUE, ...)

Arguments

graph

An undirected igraph object. The 'observed' graph to compare to simulated graphs.

layer.attr

Character. The name of the vertices attribute giving the layer of the fragments.

iter

Numerical. The number of simulated graphs to generate for each hypothesis (minimal value: 30).

summarise

Logical. Whether to report a comparative summary of the results.

...

Further arguments passed to the 'frag.simul.process' function.

Details

This function is a convenient wrapper integrating several functions of the archeofrag package to compare an observed fragmentation graph to similar simulated graphs. The frag.simul.process is used to generate two series of graphs from the properties of the observed graph: the first series is generated under the formation hypothesis H1 (one initial spatial unit) and the second series is generated under the hypothesis H2 (two initial spatial units). The edge count, edge weights sum, balance, disturbance, admixture, and cohesion values of the generated graphs are measured.

By default, the results are post-processed with the frag.simul.summarise function and a summary data frame is printed and included in the list of results which is silently returned. If the summarise parameter is set to FALSE, then the function returns a list of two data frames containing the numeric values measured for H1 and H2.

Value

A named list with three items: "h1.data", a data frame with the numerical values measured on the graphs generated for H1; "h2.data", a data frame with the numerical values measured on the graphs generated for H2; "summary", a data frame summarising the comparison between the results for the two hypotheses and the values measured on the empirical graph.

Author(s)

Sebastien Plutniak <sebastien.plutniak at posteo.net>

See Also

frag.simul.process, frag.simul.summarise

Examples

g <- frag.simul.process(n.components=20, vertices=50, disturbance=.15)
g <- frag.edges.weighting(g, layer.attr="layer")
## Not run: frag.simul.compare(g, layer.attr="layer", iter=30)