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

## ----setup, , echo = FALSE, results = "hide", message = FALSE-----------------
library(reappraised)

## -----------------------------------------------------------------------------
 pval_cont_data <- load_clean(import= "no", file.cont = "SI_pvals_cont", pval_cont= "yes", format.cont = "wide")$pval_cont_data

## ----eval = FALSE-------------------------------------------------------------
# # get path for example files
# path <- system.file("extdata", "reappraised_examples.xlsx",
#              package = "reappraised", mustWork = TRUE)
# # delete file name from path
# path <- sub("/[^/]+$", "", path)
# 
# # load data
# pval_cont_data <- load_clean(import= "yes", pval_cont = "yes", dir = path,
#         file.name.cont = "reappraised_examples.xlsx",
#         sheet.name.cont = "SI_pvals_cont",
#         range.name.cont = "A1:O51", format.cont = "wide")$pval_cont_data

## ----fig.width=4, fig.align='center', fig.asp = 0.7---------------------------
base_pval_cont <- pval_cont_fn(pval_cont_data, btsp= 500, verbose = FALSE)
base_pval_cont$all_results$pval_cont_calculated_pvalues

## ----fig.width=4, fig.align='center', fig.asp = 0.7---------------------------
base_pval_cont$all_results$pval_cont_auc_calculated_pvalues

## ----fig.width=4, fig.align='center', fig.asp = 0.7---------------------------
base_pval_cont$all_results$pval_cont_reported_pvalues

## ----fig.width=4, fig.align='center', fig.asp = 0.7---------------------------
base_pval_cont$all_results$pval_cont_auc_reported_pvalues

## -----------------------------------------------------------------------------
knitr::kable(head(base_pval_cont$pval_cont_check), format = "html",
             col.names = c("Study", "Variable", "Reported P", "Flag", "Min Possible P (Either)", "Max Possible P (Either)")) |>
kableExtra::kable_styling(bootstrap_options = c("striped", "condensed"), full_width = F, font_size = 12) |> kableExtra::column_spec(5:6, width = "5em")

## -----------------------------------------------------------------------------
match_data <- load_clean(import= "no", file.cont = "SI_pvals_cont", match= "yes", format.cont = "wide", verbose = FALSE)$match_data
match <- match_fn(match_data, verbose= FALSE)
match$match_ft_all

## -----------------------------------------------------------------------------
knitr::kable(match$match_probbystudy |> dplyr::filter (Study == "H11"), format = "html") |> kableExtra::kable_styling(bootstrap_options = c("striped", "condensed"), full_width = F, font_size = 12) |> kableExtra::column_spec(c(2,6), width = "5em")

## ----fig.width=9, fig.align='center', fig.height=10---------------------------
match$match_graphs$comp

## ----warning=FALSE------------------------------------------------------------
cohort_data <- load_clean(import= "no", file.cont = "SI_cohort", cohort= "yes", format.cont = "long", verbose= FALSE)$cohort_data
cohort <- cohort_fn(cohort_data, seed = 10, sims = 100, verbose= FALSE)
cohort$cohort_ft

## ----fig.width=9, fig.align='center', fig.height=10---------------------------
cohort$cohort_all_graphs$all_graph

## ----fig.width=4, fig.align='center', fig.height=4, warning = FALSE-----------

anova_data <- load_clean(import= "no", file.cont = "SI_pvals_cont",anova= "yes", format.cont = "wide")$anova_data
anova_fn(anova_data, seed = 10, sims = 100, btsp = 100, verbose = FALSE)$anova_ecdf

## ----fig.width=7, fig.align='center', fig.height=8----------------------------
cat_data <- load_clean(import= "no", file.cat = "SI_cat", cat= "yes", format.cat = "wide", cat.names = c("n", "w"))$cat_data
cat_fn(cat_data, x_title="withdrawals", prefix="w", del.disparate="yes", verbose = FALSE)$cat_graph

## ----fig.width=4, fig.align='center', fig.asp=.67, warning = FALSE------------
sr_data <- load_clean(import= "no", file.cat = "SI_cat", sr= "yes", format.cat = "wide")$sr_data
sr_fn(sr_data, verbose = FALSE)$sr_graph

## ----fig.width=4, fig.align='center', fig.asp=.67, warning = FALSE------------
#create data set
block = data.frame(study = c(1,2,3,4,5,6,7,8,9,10), #study numbers
                   fb_sz= c(2,4,6,8,10,12,8,8,6,14), #final block size
                   n_fb = c(1,1,4,5,7,8,4,6,2,10), #number in final block
                   df=c(1,1,0,1,3,4,2,2,0,0)) #difference between groups


sr_fn(br = "yes", block = block, verbose = FALSE)$sr_graph

## ----fig.width=7, fig.align='center', fig.height=7, warning = FALSE-----------

cat_all_data <- load_clean(import= "No", file.cat = "SI_cat_all", cat_all= "yes", format.cat = "wide")$cat_all_data

cat_all_fn (cat_all_data, two_levels = "y", del.disparate = "y", excl.level = "yes", seed = 10, verbose = FALSE)$cat_all_graph_pc

## -----------------------------------------------------------------------------
pval_cat_data <- load_clean(import= "no", file.cat = "SI_cat_all", pval_cat= "yes", format.cont = "wide", verbose = FALSE)$pval_cat_data
pval_cat <- pval_cat_fn (pval_cat_data, match_pvals = "yes", miss_data_explain = "yes", verbose = FALSE)
pval_cat$pval_cat_ft_diff_calc_rep_p

## -----------------------------------------------------------------------------
row <- as.data.frame(t(pval_cat$pval_cat_check [9,]))
knitr::kable(row, col.names = "Value", caption = "Detailed view of the first result")

## ----fig.width=4, fig.align='center', fig.height=6, warning = FALSE-----------

pval_cat_data <- load_clean(import= "no", file.cat = "SI_cat_all", pval_cat= "yes", format.cont = "wide", verbose = FALSE)$pval_cat_data
pval_cat_fn(pval_cat_data, sims = 100, seed=10, btsp = 100, verbose = FALSE)$pval_cat_calculated_pvalues

## ----fig.width=4, fig.align='center', fig.asp=.7, warning = FALSE-------------

generic_data <- load_clean(import= "No", file.cont = "SI_pvals_cont", generic= "yes", gen.vars.del = c("p"), format.cont = "wide")$generic_data
final_digit_fn(generic_data, vars = c("m","s"), dec.pl = "n", verbose = FALSE)$digit_graph 

## ----fig.width=4, fig.align='center', warning = FALSE-------------------------

pval_cat_data <- load_clean(import= "no", file.cat = "SI_cat_all", pval_cat= "yes", format.cont = "wide", verbose = FALSE)$pval_cat_data
pval_cat <- pval_cat_fn(pval_cat_data, sims = 100, seed=10, btsp = 100, verbose = FALSE)
graphs <- pval_cat$all_results$graph_data_calc
graph_maker_fn(graphs) [[2]]