varix sensitivity analysis
市田 親正
2023/4/20
必要なlibraryの読み込み
library(tidyverse)## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.2 ✔ readr 2.1.4
## ✔ forcats 1.0.0 ✔ stringr 1.5.0
## ✔ ggplot2 3.4.1 ✔ tibble 3.2.1
## ✔ lubridate 1.9.2 ✔ tidyr 1.3.0
## ✔ purrr 1.0.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errorslibrary(readxl)
library(tableone)
library(lubridate)
library(skimr)
library(svglite)
library(haven)
library(summarytools)## Warning in system2("/usr/bin/otool", c("-L", shQuote(DSO)), stdout = TRUE):
## 命令 ''/usr/bin/otool' -L
## '/Library/Frameworks/R.framework/Resources/library/tcltk/libs//tcltk.so''
## の実行は状態 1 を持ちました##
## 次のパッケージを付け加えます: 'summarytools'
##
## 以下のオブジェクトは 'package:tibble' からマスクされています:
##
## viewlibrary(naniar)##
## 次のパッケージを付け加えます: 'naniar'
##
## 以下のオブジェクトは 'package:skimr' からマスクされています:
##
## n_completelibrary(devtools)## 要求されたパッケージ usethis をロード中ですlibrary(reader)## 要求されたパッケージ NCmisc をロード中です
##
## 次のパッケージを付け加えます: 'reader'
##
## 以下のオブジェクトは 'package:NCmisc' からマスクされています:
##
## cat.path, get.ext, rmv.extlibrary(stringr)
library(missForest)
library(mice)##
## 次のパッケージを付け加えます: 'mice'
##
## 以下のオブジェクトは 'package:stats' からマスクされています:
##
## filter
##
## 以下のオブジェクトは 'package:base' からマスクされています:
##
## cbind, rbindlibrary(geepack)
library(ggplot2)
library(cobalt)## cobalt (Version 4.4.1, Build Date: 2022-11-03)library(tableone)
library(sandwich)
library(survey)## 要求されたパッケージ grid をロード中です
## 要求されたパッケージ Matrix をロード中です
##
## 次のパッケージを付け加えます: 'Matrix'
##
## 以下のオブジェクトは 'package:tidyr' からマスクされています:
##
## expand, pack, unpack
##
## 要求されたパッケージ survival をロード中です
##
## 次のパッケージを付け加えます: 'survey'
##
## 以下のオブジェクトは 'package:graphics' からマスクされています:
##
## dotchartlibrary(WeightIt)
library(twang)## To reproduce results from prior versions of the twang package, please see the version="legacy" option described in the documentation.library(naniar)
pacman::p_load(norm2, miceadds, lmtest, car, ROCR, pROC)データ読み込み
data <- read_excel("varix_icu_ok.xlsx")データの整理
df <-
data|> #miceパッケージが順序ありの因子型でないと読めないためこれを指定。
mutate(
hosp_id=as.integer(hosp_id),
year=as.integer(year),
pt_id=as.integer(pt_id),
ad_num=as.integer(ad_num),
antibio_exposure=as.integer(antibio_exposure), #geeglmを回すにはfactorでは回らない
antibio_type_a= factor(antibio_type_a, levels = c("a", "b","c","d","e","f","g","h","i","99")),
antibio_type_b= factor(antibio_type_b, levels = c("c","d","f","h","i","99")),
antibio_term= as.integer(antibio_term ),
age_num=as.integer(age_num),
age_cate=factor(age_cate,levels = c("0", "1", "2", "3")),
sex= factor(sex, levels = c("M", "F")),
smoking= as.integer(smoking),
brthel_cate= factor(brthel_cate, levels = c("0", "1", "2")),
child_numl= as.integer(child_numl),
child_score=factor(child_score, levels = c("a", "b", "c")),
jcs_all=factor(jcs_all, levels = c("0", "1", "2","3","10","20","30","100","200","300")),
cci_num=as.integer(cci_num),
malignancy=factor(malignancy),
hd=factor(hd),
hcc=factor(hcc),
alocohol=factor(alocohol),
past_varix_rup=factor(past_varix_rup),
antiplate_user_01=factor(antiplate_user_01),
anticoag_user01=factor(anticoag_user01),
nsaids_user01=factor(nsaids_user01),
steroid_user01=factor(steroid_user01),
ppi_user01=factor(ppi_user01),
vaso=factor(vaso),
map_day0= as.integer(map_day0),
shock=factor(shock),
eGFR30=factor(eGFR30),
outcome=factor(outcome),
death_day2=factor(death_day2),
rebleed_end_hemo=factor(rebleed_end_hemo),
sbp=factor(sbp),
infection=factor(infection),
infection_all=factor(infection_all),
los=as.integer(los)
)感度分析 抗生剤投与期間2日以上
# Define the conditions for the datasets
conditions <- list(
"Unrestricted" = df,
"2 days or more" = df %>% filter(!(antibio_term == 1)),
"3 days or more" = df %>% filter(!(antibio_term == 1 | antibio_term == 2)),
"4 days or more" = df %>% filter(!(antibio_term == 1 | antibio_term == 2 | antibio_term == 3))
)
# Perform IPTW analysis for each dataset
results <- list()
for (condition_name in names(conditions)) {
imp <- mice(conditions[[condition_name]], m = 5, maxit = 5, method = "pmm", seed = 10)
gee_results <- list()
summary_tables <- list()
conf_ints <- list()
weighted_data <- list()
ps_data <- list()
for (i in 1:5) {
imputed_data <- complete(imp, i)
# Drop unused levels in factors
imputed_data <- droplevels(imputed_data)
# Fit GEE propensity score model
propensity_model <- geeglm(antibio_exposure ~ age_cate + sex + brthel_cate + child_numl + cci_num + malignancy + hd + hcc + alocohol + past_varix_rup + antiplate_user_01 + anticoag_user01 + nsaids_user01 + steroid_user01 + ppi_user01 + vaso + map_day0 + shock + t_bil + ast + alt + wbc + hb + plt + alb + eGFR30 + crp + pt + aptt,
data = imputed_data,
id = hosp_id,
family = binomial("logit"),
corstr = "exchangeable")
# Calculate propensity scores
ps <- predict(propensity_model, type = "response", newdata = imputed_data)
imputed_data$ps <- ps
# Calculate IPTW and fit the logistic regression model
stabilized_weights <- with(imputed_data, antibio_exposure/(ps) + (1 - antibio_exposure)/(1 - ps))
iptw_model <- glm(outcome ~ antibio_exposure, data = imputed_data, weights = stabilized_weights, family = binomial("logit"))
# Store the results
gee_results[[i]] <- iptw_model
summary_tables[[i]] <- summary(iptw_model)
conf_ints[[i]] <- confint(iptw_model, level = 0.95, method = "Wald")
# Store weighted data
imputed_data$weight <- stabilized_weights
weighted_data[[i]] <- imputed_data
ps_data[[i]] <- imputed_data[, c("ps", "antibio_exposure")]
}
# Save the results from each imputed dataset
params <- list()
for (i in 1:5) {
params[[i]] <- data.frame(Estimate = coef(gee_results[[i]])[2],
Std.Error = coef(summary_tables[[i]])[2, "Std. Error"],
df = Inf)
}
# Calculate the pooled estimates
Qbar <- mean(sapply(params, function(x) x$Estimate))
se_Qbar <- sqrt(sum(sapply(params, function(x) x$Std.Error^2))/length(params))
# Calculate the pooled odds ratio and confidence intervals
pooled_or <- exp(Qbar)
pooled_lower <- exp(Qbar - 1.96 * se_Qbar)
pooled_upper <- exp(Qbar + 1.96 * se_Qbar)
# Calculate the pooled p-value (two-tailed)
pooled_p <- 2 * (1 - pnorm(abs(Qbar / se_Qbar)))
# Store the results in a data frame
results[[condition_name]] <- data.frame(
Condition = condition_name,
OR = pooled_or,
Lower_CI = pooled_lower,
Upper_CI = pooled_upper,
P_value = pooled_p
)
}##
## iter imp variable
## 1 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt## Warning: Number of logged events: 450## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!##
## iter imp variable
## 1 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt## Warning: Number of logged events: 450
## Warning: non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!##
## iter imp variable
## 1 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 1 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 2 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 3 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 4 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 1 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 2 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 3 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 4 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt
## 5 5 bmi_num smoking brthel_cate child_numl child_score shock t_bil ast alt wbc hb plt alb eGFR eGFR30 crp pt aptt## Warning: Number of logged events: 450
## Warning: non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!##
## iter imp variable
## 1 1 bmi_num* smoking* brthel_cate* child_numl child_score* shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 1 2 bmi_num smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 1 3 bmi_num smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 1 4 bmi_num* smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 1 5 bmi_num smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 2 1 bmi_num* smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 2 2 bmi_num smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 2 3 bmi_num* smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 2 4 bmi_num* smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 2 5 bmi_num smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 3 1 bmi_num smoking* brthel_cate* child_numl child_score* shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 3 2 bmi_num* smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 3 3 bmi_num smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 3 4 bmi_num* smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 3 5 bmi_num smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 4 1 bmi_num* smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 4 2 bmi_num smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 4 3 bmi_num* smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 4 4 bmi_num* smoking* brthel_cate* child_numl child_score* shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 4 5 bmi_num* smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 5 1 bmi_num* smoking* brthel_cate* child_numl child_score* shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 5 2 bmi_num smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 5 3 bmi_num smoking* brthel_cate* child_numl child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt
## 5 4 bmi_num smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp pt* aptt
## 5 5 bmi_num smoking* brthel_cate* child_numl* child_score shock* t_bil* ast* alt* wbc hb plt alb* eGFR eGFR30 crp* pt* aptt## Warning: Number of logged events: 687
## Warning: non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!## Waiting for profiling to be done...## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!
## Warning in eval(family$initialize): non-integer #successes in a binomial glm!# Combine the results into a single data frame
results_table <- do.call(rbind, results)
# Display results as a table
print(results_table)## Condition OR Lower_CI Upper_CI P_value
## Unrestricted Unrestricted 1.230067 0.8779559 1.723394 0.2287818
## 2 days or more 2 days or more 1.317591 0.9308002 1.865112 0.1198151
## 3 days or more 3 days or more 1.328498 0.9271164 1.903652 0.1217017
## 4 days or more 4 days or more 1.312365 0.9089354 1.894856 0.1469175# Create forest plot data
forest_plot_data <- data.frame(
Outcome = results_table$Condition,
OR = results_table$OR,
Lower_CI = results_table$Lower_CI,
Upper_CI = results_table$Upper_CI,
P_value = results_table$P_value
)
# Create forest plot
forest_plot_ordered <- ggplot(forest_plot_data, aes(x = Outcome, y = OR, ymin = Lower_CI, ymax = Upper_CI)) +
geom_pointrange(size = 0.5, fatten = 0.5) +
geom_hline(yintercept = 1, linetype = "dashed") +
scale_y_log10(limits = c(0.5, 4)) + # オッズ比の幅を指定
coord_flip() +
theme_bw() +
labs(
x = "Condition",
y = "Odds Ratio (95% Confidence Interval)"
) +
geom_text(aes(label = sprintf("p = %.3f", P_value)), vjust = -1.0, size = 3) +
scale_x_discrete(limits = c("4 days or more", "3 days or more", "2 days or more", "Unrestricted"))
# Display the forest plot
print(forest_plot_ordered)