setwd("C:/Users/jmorel/RCharlson/test")
library(testthat)
library(RCharlson)
abstract1_expected <- read.csv("abstract1_expected.csv",header=TRUE,stringsAsFactors=FALSE)
abstract2_expected <- read.csv("abstract2_expected.csv",header=TRUE,stringsAsFactors=FALSE)
abstract3_expected <- read.csv("abstract3_expected.csv",header=TRUE,stringsAsFactors=FALSE)
expect_equal(charlsonindex(abstract1_expected[,1:51],diag_types=TRUE)$charlson_index,abstract1_expected$charlson_index)
expect_equal(charlsonindex(abstract2_expected[,1:23],diag_types=TRUE)$charlson_index,abstract2_expected$charlson_index)
expect_equal(charlsonindex(abstract3_expected[,1:26])$charlson_index,abstract3_expected$charlson_index)
# if(charlson_index_total == 0) {
# comorbidity_index = 0
# } else if (charlson_index_total %in% c(1,2)) {
# comorbidity_index = 1
# } else {comorbidity_index = 2}
#Column index of first diagnosis code
diag_index = grep("diag_code_01",colnames(visits))
actual_columns = colnames(visits)[diag_index:(diag_index+49)]
if (mean(actual_columns == actual_columns) < 1 | is.na(mean(actual_columns == actual_columns))) {
stop("diag_code_01 to diag_code_25 and diag_type_01 to diag_type_25 are required in that order")}
setwd("C:/Users/Josh.Josh-PC/Documents/healthcarerisk/test")
visits = read.csv("sourcedata_for_charlson.csv",header=TRUE,stringsAsFactors = FALSE)
#Testing Readmissions
setwd("~/../RCharlson/test")
visits_tests <- read.csv("abstract_readmit1_expected.csv",header=TRUE,stringsAsFactors = FALSE)
visits_tests <- predict_readmissions(visits_tests,
prov_ref_rate = 0.164359254690341,
readm_model_params_path = "~/../RCharlson/data/hig_readmit_input_model_parameters.csv")
cbind(visits_tests$cohort,round(visits_tests$readm_prob_expected - visits_tests$readm_prob,4))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.