SETSe_bicomp: SETSe embedding on each bi-connected component using...

Description Usage Arguments Details Value See Also Examples

View source: R/SETSe_bicomp.R

Description

Performs the SETS embedding on a network using the bi-connected component/block-tree method. This is the most reliable method to perform SETSe embeddings and can be substantially quicker on certain network topologies.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
SETSe_bicomp(
  g,
  force = "force",
  distance = "distance",
  edge_name = "edge_name",
  k = "k",
  tstep = 0.02,
  tol,
  max_iter = 20000,
  mass = NULL,
  sparse = FALSE,
  sample = 100,
  static_limit = NULL,
  hyper_iters = 100,
  hyper_tol = 0.1,
  hyper_max = 30000,
  drag_min = 0.01,
  drag_max = 100,
  tstep_change = 0.2,
  verbose = FALSE,
  noisy_termination = TRUE
)

Arguments

g

An igraph object

force

A character string. This is the node attribute that contains the force the nodes exert on the network.

distance

A character string. The edge attribute that contains the original/horizontal distance between nodes.

edge_name

A character string. This is the edge attribute that contains the edge_name of the edges.

k

A character string. This is k for the moment don't change it.

tstep

A numeric. The time interval used to iterate through the network dynamics.

tol

A numeric. The tolerance factor for early stopping.

max_iter

An integer. The maximum number of iterations before stopping. Larger networks usually need more iterations.

mass

A numeric. This is the mass constant of the nodes in normalised networks. Default is set to NULL and call mass_adjuster to set the mass for each biconnected component

sparse

Logical. Whether sparse matrices will be used. This becomes valuable for larger networks

sample

Integer. The dynamics will be stored only if the iteration number is a multiple of the sample. This can greatly reduce the size of the results file for large numbers of iterations. Must be a multiple of the max_iter

static_limit

Numeric. The maximum value the static force can reach before the algorithm terminates early. This prevents calculation in a diverging system. The value should be set to some multiple greater than one of the force in the system. If left blank the static limit is the system absolute mean force.

hyper_iters

integer. The hyper parameter that determines the number of iterations allowed to find an acceptable convergence value.

hyper_tol

numeric. The convergence tolerance when trying to find the minimum value

hyper_max

integer. The maximum number of iterations that SETSe will go through whilst searching for the minimum.

drag_min

integer. A power of ten. The lowest drag value to be used in the search

drag_max

integer. A power of ten. if the drag exceeds this value the tstep is reduced

tstep_change

numeric. A value between 0 and 1 that determines how much the time step will be reduced by default value is 0.5

verbose

Logical. This value sets whether messages generated during the process are suppressed or not.

noisy_termination

Stop the process if the static force does not monotonically decrease.

Details

This approach can be faster for larger graphs or graphs with many nodes of degree 2, or networks with a low clustering coefficient. This is because although the algorithm is very efficient the topology of larger graphs make them more difficult to converge. Large graph tend to be made of 1 very large biconnected component and many very small biconnected components. As the mass of the system is concentrated in the major biconnected component the small ones can be knocked around by minor movements of the major. This leads to long convergence times. By solving all biconnected components separately and then reassembling the block tree at the end, the system can be converged considerably faster. In addition the smaller biconnected components iterate faster than a single large one.

Setting mass to the absolute system force divided by the total nodes, often leads to faster convergence. As such When mass is left to the default of NULL, the mean force value is used.

Value

A list containing 5 dataframes.

  1. The node embeddings. Includes all data on the nodes the forces exerted on them position and dynamics at simulation termination

  2. The network dynamics describing several key figures of the network during the convergence process, this includes the static_force

  3. memory_df A dataframe recording the iteration history of the convergence of each component.

  4. A data frame giving the time taken for the simulation as well as the number of nodes and edges. Node and edge data is given as this may differ from the total number of nodes and edges in the network depending on the method used for convergence. For example if SETSe_bicomp is used then some simulations may contain as little as two nodes and 1 edge

  5. time taken. the amount of time taken per component, includes the edge and nodes of each component

See Also

SETSe_auto SETSe

Examples

1
2
3
4
5
6
7
8
9
set.seed(234) #set the random see for generating the network
g <- generate_peels_network(type = "E")
embeddings <- g %>%
#prepare the network for a binary embedding
prepare_SETSe_binary(., node_names = "name", k = 1000, 
                     force_var = "class", 
                     positive_value = "A") %>%
#embed the network using  SETSe_bicomp
  SETSe_bicomp(tol = 0.02)

rsetse documentation built on Nov. 12, 2020, 5:08 p.m.