Package 'remify'

a remify object.

Any R object.

is a remify object.

one or more numbers between 1 and 7. Plots described in order: (1) distribution of the inter-event times (histogram), (2) tile plot titled 'activity plot', with in-degree and out-degree activity line plots on the sides (or total-degree on the top side if the network is undirected). Tiles' color is scaled based on the count of the directed (or undirected) dyad, (3) for directed networks two plots of normalized out-degree and in-degree (values ranging in [0,1]) over a set of n_intervals (evenly spaced). For undirected networks one plot of normalized total-degree over the n_intervals (also here values ranging in [0,1]). The normalization is calculated in each interval as the (degree-min(degree))/(max(degree)-min(degree))) for each actor considering minimum and maximum degree (in-, out- or total-) observed in the interval (opacity and size of the points is proportional to the normalized measure), (4) four plots: (i) number of events (# events) per time interval, (ii) proportion of observed dyads (# dyads / x$D) per time interval, (iii) and (iv) (for directed network only) proportion of active senders and receivers per time interval (calculated as # senders/ x$N and # receiver/x$N per interval), (5) two networks: (i) network of events where edges are considered undirected (edges' opacity is proportional to the counts of the undirected events, vertices' opacity is proportional to the total-degree of the actors), (ii) visualization of directed network (edges' opacity is proportional to the counts of the directed events, vertices' opacity is proportional to the in-degree of the actors), (6) counting process N(t): cumulative number of events against real (non-ordinal) time, with a dashed constant-rate reference line. A roughly straight curve indicates a stationary baseline rate; systematic departures or kinks indicate non-stationarity. Skipped for ordinal event histories, (7) dyad event raster: one row per dyad with a tick per event, capped at the top 30 most active dyads. For directed histories the two directions of a pair are kept as adjacent rows, so reciprocity appears as neighbouring rows firing in sequence and inertia as horizontal clustering of ticks, (8) distribution of dyad recurrence times: pooled gaps between successive events on the same dyad, with a dashed same-mean exponential reference. Excess mass at short gaps (above the reference) indicates temporal clustering, suggesting window/decay memory effects rather than plain inertia. For directed histories recurrence is measured on directed dyads; for ordinal histories the gap is the event-index distance and the exponential reference is omitted, (9) dyad concentration (Lorenz curve): cumulative share of events against cumulative share of dyads, taken over all x$D possible dyads so that never-active dyads enter as zeros. The dashed diagonal is perfect equality; a curve bowing toward the bottom-right indicates events concentrated on few dyads, which motivates frailty and makes per-dyad effects harder to identify. A Gini coefficient is reported in the corner, (10) reciprocity scatter (directed histories only): for each unordered pair, the event count in the busier direction against the count in the quieter direction, with point area proportional to the pair's total volume. Points near the dashed diagonal indicate balanced exchange (reciprocity has little asymmetry to explain); points near the horizontal axis indicate one-directional flow (reciprocity is informative), (11) distribution of event durations (remify_durem objects only): histogram of spell lengths (end - time) with a same-mean exponential reference for interval histories; for ordinal histories the duration is an event-index gap and the reference is omitted, (12) concurrency over time (remify_durem only): number of simultaneously-active spells as a step function, with the peak marked; spells are counted as active over the closed interval [time, end], (13) event spells / Gantt (remify_durem only): one lane per dyad with a segment per spell from start to end, capped at the 30 most active dyads and coloured by event type when the history is typed. Note, plots 2–4 use multi-panel layouts and need a sufficiently large graphics device. If you see "figure margins too large", enlarge the plot pane or open a larger device with e.g. dev.new().

default is 15L and it describes the number of cells of the histogram plot for the inter-event times. It can be specified in the same way as the argument used by the function graphics::hist() (see ?graphics::hist for more details).

a palette from grDevices::hcl.pals() (default is the "Purples" palette).

number of time intervals for time plots (default is 10).

default is TRUE (reverse order of the color specified in palette)

default is the set of actors in the network (see .dict[["actors"]]). The user can specify a subset of actors on which to run the descriptive plots. If the set contains more than 50 actors, then the function will select the 50 most active actors from the set provided.

default is 20. Shape of the points for the degree plots (in-degree, out-degree, total-degree).

color of the edges in visualization of the network with vertices and nodes. The user can specify the hex value of a color, the color name or use the functiongrDevices::rgb() which returns the hex value.

color of the vertices in visualization of the network with vertices and nodes. The user can specify the hex value of a color, the color name or use the function grDevices::rgb() which returns the hex value.

other graphical parameters

a remify object.

further arguments.

A remify_durem object.

Passed to summary.remify_durem.

the relational event history. An object of class data.frame with first three columns corresponding to time, and actors forming the dyad. The first three columns will be re-named "time", "actor1", "actor2" (where, for directed networks, "actor1" corresponds to the sender and "actor2" to the receiver of the relational event). Optional columns that can be supplied are: 'type' and 'weight'. If one or both exist in edgelist, they have to be named accordingly.

logical value indicating whether events are directed (TRUE) or undirected (FALSE). (default value is TRUE)

logical value indicating whether only the order of events matters in the model (TRUE) or also the exact timing must be considered in the model (FALSE). (default value is FALSE). If TRUE, then the column "time" of edgelist is still used to extract the order.

either "tie" (default) or "actor" oriented modeling. For "tie", the riskset is at the dyad level. For "actor", the model has two sub-processes: a sender rate model (who sends next?) and a receiver choice model (who does the sender choose?). Actor-oriented modeling requires directed=TRUE. The returned object includes sender_riskset, receiver_riskset, and activeN (see @return).

[optional] character vector of actors' names that may be observed interacting in the network. If NULL (default), actors' names will be taken from the input edgelist.

[optional] character value indicating the type of risk set to process: riskset = "full" (default) consists of all the possible dyadic events given the number of actors (and the number of event types) and it mantains the same structure over time. riskset = "active" considers at risk only the observed dyads and it mantains the same structure over time. riskset = "manual", allows the risk set to have a structure that is user-defined, and it is based on the instructions supplied via the argument manual_riskset. riskset = "active_saturated" extends the active riskset by adding the reverse direction for each observed dyad (if A->B is observed, B->A is also at risk) and includes all event types for each observed actor pair (type column is ignored). This reflects the assumption that observing any interaction between two actors implies both directions and all types are possible.

[optional] When riskset = "manual", this argument of class data.frame specifies which dyadic riskset to consider through the entire sequence. If observed dyads from the edgelist are missing, they will be automatically be added.

logical. FALSE (default). When event types are present (via event_type), controls whether the risk set is expanded over types. If TRUE (default when types are present), each actor pair is duplicated for each event type, so the risk set has size $D = N(N-1) \times C$ (directed) or $D = N(N-1)/2 \times C$ (undirected), and the type column appears in the decoded risk set. If FALSE, event type is treated as a mark on events only and does not expand the risk set: $D = N(N-1)$ (directed) or $D = N(N-1)/2$ (undirected), and no type column appears in the decoded risk set. This argument is ignored when no event types are present.

Optional. Either NULL (default) or a single character string giving the name of the column in edgelist that contains event types (marks).

If event_type is NULL, remify() uses edgelist$type if it exists; otherwise events are treated as untyped.

If event_type is a column name, that column is used as the event-type mark. If a column named type already exists and event_type != "type", the existing edgelist$type is overridden (with a warning).

When event types are present (via edgelist$type or event_type), the dyadic risk set is extended over types, i.e., each dyad is duplicated for each event type (dyad $\times$ type).

Optional. Either NULL (default) or a single character string giving the name of the column in edgelist that contains event weights.

If event_weight is NULL, remify() uses edgelist$weight if it exists; otherwise events are treated as unweighted (uniform weight = 1).

If event_weight is a column name, that column is used as the event weight. If a column named weight already exists and event_weight != "weight", the existing edgelist$weight is overridden (with a warning). This argument behaves analogously to event_type and applies to both standard and duration (duration = TRUE) models.

[optional] starting time point of the observation period (default is NULL). If it is supplied, it must have the same class of the 'time' column in the input edgelist. If unsupplied, the origin is set to the average waiting time in the sequence subtracted from the time of the first event.

Character string specifying the time unit for converting time values when 'edgelist$time' is of class Date or POSIXct; ignored for numeric or integer time. Default is "secs".

Integer >= 1. Event-time aggregated based on unique time points. Keeps every aggregate_time-th unique event time (after time translation) and maps each event time to the next kept time point (i.e., ceiling to the kept grid). This reduces the number of unique time points (and thus memory/computation in later steps).

Logical. If TRUE, attaches a list riskset_info to the returned remify object. The list contains the effective risk set representation used for estimation (e.g., riskset_idx, dyadIDactive, dictionaries, and basic risk set metadata). This is intended to make the returned object self-describing and easier to inspect/debug.

Character. Controls how (and whether) the included risk set dyads are decoded and attached in riskset_info$included.

"labels"

Attach a decoded dyad table including actor (and type) labels (e.g., actor1, actor2, and optional type).

"ids"

Attach a decoded dyad table with integer IDs only (e.g., actor1ID, actor2ID, optional typeID, and dyadID).

"none"

Do not attach a decoded dyad table.

Integer. Maximum number of included dyads (i.e., length(riskset_idx) / D_active) for which riskset_decode="labels" is allowed. If the included risk set exceeds this threshold, decoding to labels is skipped (typically falling back to "ids" with a warning) to avoid large memory usage.

Optional character vector of column names in edgelist to retain as additional event-level variables in the returned reh object.

These columns are stored as reh$event_attributes together with the corresponding time, actor1, and actor2 columns (and an internal .event_id). This is useful when downstream functions (e.g., in remstats) need access to event-level marks/attributes/covariates that are not part of the core reh$edgelist produced by remify().

Note: event_attributes does not affect risk set construction or type handling in remify(); it only preserves additional columns for later use. Currently there is no further support yet when event_attributes have been added.

[optional] number of cores used in the parallelization of the processing functions. (default is 1).

Logical. If TRUE, the edgelist is treated as a duration edgelist (each event has both a start time and an end time) and a remify_durem object is returned instead of a standard remify object. The edgelist must contain an end column (or end_time / duration as alternatives; see .durem_normalize_edgelist). Default FALSE.

Logical. Only used when duration = TRUE. If FALSE (default), the end process is undirected: either actor can terminate the event and only a combined dyad-level end rate is modelled. If TRUE, the end process is directed: a who_ended column in the edgelist ("actor1" / "actor2" / NA) records which actor terminated each event, enabling actor-level end-rate models. When who_ended is absent and dur_directed_end = TRUE, actor1 is assumed to terminate all events and a message is issued.

Logical. Only used when duration = TRUE, when typed events are present and when extend_riskset_by_type = TRUE. If TRUE, an active event of any type is a hard block on starting events of all other types for the same dyad (e.g. two actors cannot start a text conversation while already in a face-to-face conversation). If FALSE (default), types are treated as independent processes and a dyad can be simultaneously active in events of different types. Has no effect when extend_riskset_by_type = FALSE.

a remify object.

other arguments.

A remify_durem object.

Ignored.