inst/tinytest/test-interact.R

In remstats: Computes Statistics for Relational Event History Data

# Small edgelist
edgelist <- data.frame(
  time = 1:5,
  actor1 = c(1, 1, 2, 2, 3),
  actor2 = c(2, 3, 1, 3, 2)
)

# Actor info
info <- data.frame(
  name = 1:3,
  time = rep(0, 3),
  x1 = c(10, 20, 30),
  x2 = c(0, 1, 1)
)

# Tie model
reh <- remify(edgelist, model = "tie")
effects <- ~ send(variable = "x1"):inertia()
stats <- remstats(reh, tie_effects = effects, attr_actors = info)

interaction <- stats[, , "send_x1"] * stats[, , "inertia"]
expect_equal(stats[, , "send_x1:inertia"], interaction)

# # Actor model
# reh <- remify(edgelist, model = "actor")
# sender_effects <- ~ send(variable = "x1"):outdegreeSender()
# receiver_effects <- ~ receive(variable = "x1"):inertia()
# stats <- remstats(reh, sender_effects = sender_effects, 
#     receiver_effects = receiver_effects, attr_actors = info)
# sender_stats <- stats$sender_stats
# receiver_stats <- stats$receiver_stats
# 
# sender_interaction <- sender_stats[, , "send_x1"] * 
#   sender_stats[, , "outdegreeSender"]
# expect_equal(sender_stats[, , "send_x1:outdegreeSender"], 
#   sender_interaction)
# 
# receiver_interaction <- receiver_stats[, , "receive_x1"] * 
#   receiver_stats[, , "inertia"]
# expect_equal(receiver_stats[, , "receive_x1:inertia"], receiver_interaction)