# (1) load the data into the workspace
data(tie_data)
# (2) process event sequence with \code{remify}
tie_reh <- remify::remify(edgelist = tie_data$edgelist, model = "tie")
# (3) define linear predictor and claculate stastistcs with \code{remstats} package
## linear predictor
tie_model <- ~ 1 + remstats::indegreeSender() + remstats::inertia() + remstats::reciprocity()
## calculate statistics
tie_reh_stats <- remstats::remstats(reh = tie_reh, tie_effects = tie_model)
# (4) estimate model using method = "MLE" and print out summary
## estimate model
mle_tie <- remstimate::remstimate(reh = tie_reh, stats = tie_reh_stats, method = "MLE")
## print out a summary of the estimation
summary(mle_tie)
Run the code above in your browser using DataLab