## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----eval=FALSE---------------------------------------------------------------
# install.packages('NetworkInference')

## ----eval=FALSE---------------------------------------------------------------
# install.packages(c('dplyr', 'igraph'))

## ----message=FALSE------------------------------------------------------------
library(NetworkInference)
# Load the `policies` dataset (?policies for details).
data('policies')

## ----eval=FALSE---------------------------------------------------------------
# head(policies)

## ----results="asis", echo=FALSE-----------------------------------------------
pander::pandoc.table(head(policies))

## -----------------------------------------------------------------------------
length(unique(policies$policy))

## -----------------------------------------------------------------------------
nrow(policies)

## -----------------------------------------------------------------------------
unique(policies$policy)[1:5]

## ----eval=FALSE---------------------------------------------------------------
# library(dplyr)
# filter(policies_metadata, policy %in% unique(policy)[100:104]) %>%
#     select(-source)

## ----results="asis", echo=FALSE, message=FALSE--------------------------------
library(dplyr)
pander::pandoc.table(filter(policies_metadata, policy %in% unique(policy)[100:104]) %>%
                         select(-source))

## -----------------------------------------------------------------------------
policy_cascades <- as_cascade_long(policies, cascade_node_name = 'statenam',
                                   event_time = 'adopt_year', 
                                   cascade_id = 'policy')

## -----------------------------------------------------------------------------
class(policy_cascades)
length(policy_cascades)
names(policy_cascades)

## -----------------------------------------------------------------------------
policy_cascades$cascade_nodes[2:3]

## -----------------------------------------------------------------------------
policy_cascades$cascade_times[2:3]

## -----------------------------------------------------------------------------
selected_policies <- subset_cascade(cascade = policy_cascades, 
                                    selection = c('clinic_access', 'cogrowman'))
selected_policies[1:2]

## -----------------------------------------------------------------------------
time_constrained <- subset_cascade_time(cascade = selected_policies, 
                                        start_time = 1990, end_time = 2000)
time_constrained[1:2]

## -----------------------------------------------------------------------------
less_nodes <- drop_nodes(cascades = time_constrained, 
                         nodes = c('Maryland', 'Washington'))
less_nodes[1:2]

## -----------------------------------------------------------------------------
summary(policy_cascades)

## ----fig.align='center', fig.width=7, fig.height=4----------------------------
selection <- c('guncontrol_assaultweapon_ba', 'guncontrol_licenses_dealer')
plot(policy_cascades, label_nodes = TRUE, selection = selection)

## ----fig.align='center', fig.width=7, fig.height=4----------------------------
selection <- c('waiting', 'threestrikes', 'unionlimits', 'smokeban', 
               'paperterror', 'miglab', 'methpre', 'lott', 'lemon', 'idtheft',
               'harass', 'hatecrime', 'equalpay')
plot(policy_cascades, label_nodes = FALSE, selection = selection)

## -----------------------------------------------------------------------------
results <- netinf(policy_cascades, trans_mod = "exponential", n_edges = 100, 
                  params = 0.5, quiet = TRUE)

## -----------------------------------------------------------------------------
npe <- count_possible_edges(policy_cascades)
npe

## -----------------------------------------------------------------------------
results <- netinf(policy_cascades, trans_mod = "exponential", 
                  p_value_cutoff = 0.1, params = 0.5, quiet = TRUE)
nrow(results)

## -----------------------------------------------------------------------------
results <- netinf(policy_cascades, trans_mod = "exponential", 
                  p_value_cutoff = 0.1, quiet = TRUE)
nrow(results)

## ----eval=FALSE, echo=TRUE----------------------------------------------------
# head(results)

## ----results = "asis", echo=FALSE---------------------------------------------
pander::pandoc.table(head(results))

## ----fig.align='center', fig.width=7, fig.height=4----------------------------
plot(results, type = "improvement")

## ----fig.align='center', fig.width=7, fig.height=4----------------------------
plot(results, type = 'p-value')

## ----fig.width=7, fig.height=5.5----------------------------------------------
#install.packages('igraph')
# For this functionality the igraph package has to be installed
# This code is only executed if the package is found:
if(requireNamespace("igraph", quietly = TRUE)) {
    plot(results, type = "network")
}

## ----message=FALSE, eval=FALSE------------------------------------------------
# if(requireNamespace("igraph", quietly = TRUE)) {
#     library(igraph)
#     g <- graph_from_data_frame(d = results[, 1:2])
#     plot(g, edge.arrow.size=.3, vertex.color = "grey70")
# }