tests/testthat/test_model1.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_table1_parameters(),
    get_evaluation_parameters(checking_parameters = T),
    get_predicition_parameters()
)

# Fit Model

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

# Check and evaluate model

ppc <- check_model(model_fit, data_list)

eval <- evaluate_model(model_fit, data_list)

# Plot models

p1 <- make_figure_three(model_fit)

p2 <- make_figure_four(model_fit)

p3 <- make_figure_five(model_fit)

# Make Table

tab1 <- make_table_one(model_fit)
head(tab1)