tests/testthat/test_model2.R

In philpatton/cowbird: Recreate analysis in Patton et al (TBD)

library(cowbird)

# Load Data

data("cowbird_data")

data_list <- make_data_list(cowbird_data)

# parameters for model 1
params <- c(
    get_table2_parameters(),
    get_evaluation_parameters(checking_parameters = T),
    get_predicition_parameters()
)

# Fit Model

model_fit <- fit_model('model_2', data_list, params)

# Check and evaluate model

ppc <- check_model(model_fit, data_list)

eval <- evaluate_model(model_fit, data_list)

# Plot models

# p1 <- plot_cooccur_probs(model_fit)

p2 <- make_figure_four(model_fit)

p3 <- make_figure_five(model_fit)

# Make Table

tab2 <- make_table_two(model_fit)