rn_phases: Extract Phases from weighted RN

In FredHasselman/casnet: A Toolbox for Studying Complex Adaptive Systems and Networks

Extract Phases from weighted RN

Description

This function will extract phases (regions of attraction in phase space) based on a weighted RN object created with function rn. The assumption is that coordinates in state space that are either close in terms of distance, or are re-visited with short recurrence times, or with high-frequency, are regions of attraction for the system.

Usage

rn_phases(
  RN,
  maxPhases = NA,
  minStatesinPhase = 2,
  maxStatesinPhase = NROW(RN),
  useDegree = FALSE,
  inverseWeight = TRUE,
  cleanUp = TRUE,
  returnCentroid = c("no", "mean.sd", "median.mad", "centroid")[1],
  removeSingularities = FALSE,
  standardise = c("none", "mean.sd", "median.mad", "unit")[4],
  returnGraph = FALSE,
  doPhaseProfilePlot = FALSE,
  plotCentroid = FALSE,
  dimNames = NULL,
  colOrder = FALSE,
  phaseColours = NULL,
  doSpiralPlot = FALSE,
  doPhaseSeriesPlot = FALSE,
  doPhaseSpacePojectionPlot = FALSE,
  showEpochLegend = TRUE,
  epochColours = NULL,
  epochLabel = "Phase",
  excludeOther = FALSE,
  excludeNorec = TRUE
)

Arguments

RN	A matrix produced by the function rn
maxPhases	The maximum number of phases to extract. These will be the phases associated with the highest node degree or node strength. All other recurrent points will be labelled with "Other". If NA, the value will be set to NROW(RN), this will return all the potential phases in the data irrespective of their frequency/strength of recurrence (default = 10)
minStatesinPhase	A parameter applied after the extraction of phases (limited by maxPhases). If any extracted phases do not have a minimum number of minStatesinPhase + 1 (= the state that was selected based on node strength), the phase will be removed from the result (default = 1)
maxStatesinPhase	A parameter applied after the extraction of phases (limited by maxPhases). If any extracted phases exceeds a maximum number of maxStatesinPhase + 1 (= the state that was selected based on node strength), the phase will be removed from the result (default = NROW(RN))
useDegree	By default, node strength will be used to determine the phases. Set to TRUE to use the node degree (default = FALSE)
inverseWeight	Whether to perform the operation 1/weight on the edge weights. The default is TRUE, if the matrix was weighted by a distance metric (weightedBy = "si") edges with smaller distances (recurring coordinates closer to the current coordinate) have greater impact on the node strength calculation used to select the phases. If the matrix was weighted by recurrence time (weightedBy = "rt") and inverseWeight = TRUE, recurrent points with shorter recurrence times will have greater impact on the strength calculation. If weightedBy = "rf", lower frequencies will end up having more impact. (default = TRUE)
cleanUp	Try to assign states to phases that were not assigned by the algorithm. If FALSE, these states will be added to category "Other" (default = TRUE)
returnCentroid	Values can be "no", "mean.sd", "median.mad", "centroid". Any other value than "no" will return a data frame with the central tendency and deviation measures for each phase (default = "no")
removeSingularities	Will remove states that recur only once (nodes with degree(g) == 1) (default = FALSE)
standardise	Standardise the series using ts_standardise() with adjustN = FALSE (default = "mean.sd")
returnGraph	Returns all the graph object objects of the plots that have been produced (default = FALSE)
doPhaseProfilePlot	Produce a profile plot of the extracted phases (default = TRUE)
plotCentroid	Plot the centroid requested in returnCentroid? (default = FALSE)
dimNames	A vector of titles to use for the dimensions. If NULL the values will be read from the attribute of RN.
colOrder	Should the order of the dimensions reflect the order of the columns in the dataset? If FALSE the order will be based on the values of the dimensions observed in the first extracted phase (default = FALSE)
phaseColours	Colours for the different phases in the phase plot. If epochColours also has a value, phaseColours will be used instead (default = NULL)
doSpiralPlot	Produce a plot of the recurrence network with the nodes coloured by phases (default = FALSE)
doPhaseSeriesPlot	Produce a time series of the phases as they occur with a marginal histogram of their frequency (default = FALSE)
showEpochLegend	Should a legend be shown for the epoch colours? (default = TRUE)
epochColours	A vector of length vcount(g) with colour codes (default = NULL)
epochLabel	A title for the epoch legend (default = "Epoch")
excludeOther	Should the phase "Other" be excluded from plots? (default = FALSE)
excludeNorec	Should the category "No recurrence" be excluded from plots? (default = TRUE)
doPhaseSpaceProjectionPlot	produce a 2D umap projection of the phases (default = FALSE)

Details

The method used for the identification of phases is on the properties of the RN object:

• If weighted by distance "si", the inverse distance will be used, which means higher weights correspond to closer states.

• If weighted by recurrence time "rt", the inverse time will be used, which means higher weights correspond to faster recurrence times.

The procedure is as follows:

1. Identify the node with the highest strength

2. Identify the nodes that connect to this node

3. Identify the node with highest strength that does not connect to the node identified in step 1.

4. Repeat until criteria set in maxPhases, minStatesinPhase and maxStatesinPhase are triggered.

5. Clean up. If cleanUp = TRUE the remaining states that connect to nodes of one of the identified phases, but not to the node with highest strength identified in step 1 will be added to that phase. Otherwise these nodes will end up in category "Other".

Value

A data frame with information about the phases or a list object with data and graph objects (if requested).

See Also

Other Distance matrix operations (recurrence network): mat_di2bi(), mat_di2ch(), mat_di2we(), rn(), rn_plot(), rn_recSpec()

Examples

# Use the ManyAnalysts dataset to create a phase plot with default settings
data("manyAnalystsESM")
df <- manyAnalystsESM[4:10]
RN <- rn(y1 = df, doEmbed = FALSE, weighted = TRUE, weightedBy = "si", emRad = NA)

# This returns 6 phases which have minimally 2 states
rn_phases(RN, maxPhases = 10, doPhaseProfilePlot = TRUE)

# Use min. number of states as the extraction criterion
rn_phases(RN, maxPhases = NA, minStatesinPhase = 7, doPhaseProfilePlot = TRUE)

FredHasselman/casnet documentation built on Jan. 6, 2025, 3:18 a.m.