Living Conditions and Dropout (step 3)


Multivariate Parametric Survival Regression (Jul 2023)

Show code
ss.gomp.rmst_noboot_diff_tab<-
ss.gomp.rmst_noboot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.gomp.rmst_boot_diff_tab<-
ss.gomp.rmst_boot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.gomp.rmst_noboot_tab<-
ss.gomp.rmst_noboot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.gomp.rmst_boot_tab<-
ss.gomp.rmst_boot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))


ss.gomp.surv_noboot_diff_tab<-
ss.gomp.surv_noboot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.gomp.surv_boot_diff_tab<-
ss.gomp.surv_boot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.gomp.surv_noboot_tab<-
ss.gomp.surv_noboot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.gomp.surv_boot_tab<-
ss.gomp.surv_boot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

#_________________________
#rp04

ss.rp04.rmst_noboot_diff_tab<-
  ss.rp04.rmst_noboot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.rp04.rmst_boot_diff_tab<-
  ss.rp04.rmst_boot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.rp04.rmst_noboot_tab<-
  ss.rp04.rmst_noboot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.rp04.rmst_boot_tab<-
  ss.rp04.rmst_boot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))


ss.rp04.surv_noboot_diff_tab<-
  ss.rp04.surv_noboot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.rp04.surv_boot_diff_tab<-
  ss.rp04.surv_boot_diff %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))



ss.rp04.surv_noboot_tab<-
  ss.rp04.surv_noboot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

ss.rp04.surv_boot_tab<-
  ss.rp04.surv_boot %>% 
  dplyr::rename_with(~ paste0(.x, "_est"), .cols =  names(.)[!str_ends(names(.), "time|_lci|_uci")]) %>% 
  dplyr::rename_with(~ sub("2_1","21",.x), .cols =  names(.)[grepl("2_1",names(.))]) %>% 
  dplyr::rename_with(~ sub("3_1","31",.x), .cols =  names(.)[grepl("3_1",names(.))]) %>% 
  pivot_longer(!time, names_to=c("at","est"), names_sep="_", values_to="rmst") %>% 
  pivot_wider(id_cols=c("time", "at"), names_from="est", values_from="rmst") %>% 
  dplyr::filter(time %in% c(.5,1,3))

Time for this code chunk to run: 0 minutes

Show code
combined_df <- bind_rows(
  cbind(mod="Gompertz, RMST, No bootstrap, Difference", ss.gomp.rmst_noboot_diff_tab),
  cbind(mod="Gompertz, RMST, Bootstrap, Difference", ss.gomp.rmst_boot_diff_tab),
  cbind(mod="Gompertz, RMST, No bootstrap, Ratio", ss.gomp.rmst_noboot_tab),
  cbind(mod="Gompertz, RMST, Bootstrap, Ratio", ss.gomp.rmst_boot_tab),
  cbind(mod="Gompertz, Survival, No bootstrap, Difference", ss.gomp.surv_noboot_diff_tab),
  cbind(mod="Gompertz, Survival, Bootstrap, Difference", ss.gomp.surv_boot_diff_tab),
  cbind(mod="Gompertz, Survival, No bootstrap, Ratio", ss.gomp.surv_noboot_tab),
  cbind(mod="Gompertz, Survival, Bootstrap, Ratio", ss.gomp.surv_boot_tab),
  cbind(mod="Royston-Parmar 04, RMST, No bootstrap, Difference", ss.rp04.rmst_noboot_diff_tab),
  cbind(mod="Royston-Parmar 04, RMST, Bootstrap, Difference", ss.rp04.rmst_boot_diff_tab),
  cbind(mod="Royston-Parmar 04, RMST, No bootstrap, Ratio", ss.rp04.rmst_noboot_tab),
  cbind(mod="Royston-Parmar 04, RMST, Bootstrap, Ratio", ss.rp04.rmst_boot_tab),
  cbind(mod="Royston-Parmar 04, Survival, No bootstrap, Difference", ss.rp04.surv_noboot_diff_tab),
  cbind(mod="Royston-Parmar 04, Survival, Bootstrap, Difference", ss.rp04.surv_boot_diff_tab),
  cbind(mod="Royston-Parmar 04, Survival, No bootstrap, Ratio", ss.rp04.surv_noboot_tab),
  cbind(mod="Royston-Parmar 04, Survival, Bootstrap, Ratio", ss.rp04.surv_boot_tab)
) %>% 
  dplyr::mutate_at(c("est","lci","uci"),~sprintf("%.3f",.)) %>% 
  dplyr::mutate(est_lci_uci= paste0(est, " (",lci,", ", uci,")")) %>%
  tidyr::pivot_wider(id_cols=c(mod,time), names_from = at, values_from = est_lci_uci)

combined_df %>% 
  knitr::kable(format="markdown",caption="Parametric survival regression, RMSTs and Survival Marginal Estimations", col.names=c("model", "time (years)", "With couple/children", "Alone", "Family of origin", "Alone vs. Couple/children","Fam of origin vs. Alone"))
Table 1: Parametric survival regression, RMSTs and Survival Marginal Estimations
model time (years) With couple/children Alone Family of origin Alone vs. Couple/children Fam of origin vs. Alone
Gompertz, RMST, No bootstrap, Difference 0.5 0.361 (0.356, 0.367) 0.365 (0.362, 0.368) 0.350 (0.342, 0.357) 0.004 (0.000, 0.007) -0.012 (-0.024, 0.001)
Gompertz, RMST, No bootstrap, Difference 1.0 0.593 (0.579, 0.608) 0.603 (0.595, 0.611) 0.563 (0.545, 0.583) 0.009 (0.001, 0.018) -0.030 (-0.062, 0.002)
Gompertz, RMST, No bootstrap, Difference 3.0 1.267 (1.217, 1.320) 1.300 (1.272, 1.329) 1.164 (1.101, 1.230) 0.033 (0.002, 0.063) -0.104 (-0.214, 0.007)
Gompertz, RMST, Bootstrap, Difference 0.5 0.361 (0.355, 0.367) 0.365 (0.361, 0.368) 0.350 (0.341, 0.357) 0.004 (-0.003, 0.010) -0.012 (-0.023, -0.001)
Gompertz, RMST, Bootstrap, Difference 1.0 0.593 (0.578, 0.608) 0.603 (0.594, 0.611) 0.563 (0.543, 0.581) 0.009 (-0.007, 0.025) -0.030 (-0.056, -0.004)
Gompertz, RMST, Bootstrap, Difference 3.0 1.267 (1.217, 1.316) 1.300 (1.272, 1.328) 1.164 (1.094, 1.223) 0.033 (-0.025, 0.087) -0.104 (-0.197, -0.015)
Gompertz, RMST, No bootstrap, Ratio 0.5 0.361 (0.356, 0.367) 0.365 (0.362, 0.368) 0.350 (0.342, 0.357) 1.010 (1.001, 1.020) 0.967 (0.934, 1.002)
Gompertz, RMST, No bootstrap, Ratio 1.0 0.593 (0.579, 0.608) 0.603 (0.595, 0.611) 0.563 (0.545, 0.583) 1.016 (1.001, 1.031) 0.949 (0.898, 1.004)
Gompertz, RMST, No bootstrap, Ratio 3.0 1.267 (1.217, 1.320) 1.300 (1.272, 1.329) 1.164 (1.101, 1.230) 1.026 (1.001, 1.051) 0.918 (0.838, 1.006)
Gompertz, RMST, Bootstrap, Ratio 0.5 0.361 (0.355, 0.367) 0.365 (0.361, 0.368) 0.350 (0.341, 0.357) 1.010 (0.993, 1.028) 0.967 (0.938, 0.996)
Gompertz, RMST, Bootstrap, Ratio 1.0 0.593 (0.578, 0.608) 0.603 (0.594, 0.611) 0.563 (0.543, 0.581) 1.016 (0.989, 1.043) 0.949 (0.906, 0.993)
Gompertz, RMST, Bootstrap, Ratio 3.0 1.267 (1.217, 1.316) 1.300 (1.272, 1.328) 1.164 (1.094, 1.223) 1.026 (0.981, 1.071) 0.918 (0.850, 0.987)
Gompertz, Survival, No bootstrap, Difference 0.5 0.545 (0.528, 0.562) 0.556 (0.546, 0.565) 0.511 (0.489, 0.533) 0.011 (0.001, 0.021) -0.034 (-0.071, 0.002)
Gompertz, Survival, No bootstrap, Difference 1.0 0.406 (0.388, 0.425) 0.418 (0.408, 0.428) 0.369 (0.346, 0.393) 0.012 (0.001, 0.023) -0.037 (-0.077, 0.002)
Gompertz, Survival, No bootstrap, Difference 3.0 0.311 (0.293, 0.330) 0.322 (0.312, 0.333) 0.275 (0.253, 0.298) 0.011 (0.001, 0.022) -0.037 (-0.075, 0.002)
Gompertz, Survival, Bootstrap, Difference 0.5 0.545 (0.527, 0.561) 0.556 (0.546, 0.565) 0.511 (0.487, 0.531) 0.011 (-0.008, 0.029) -0.034 (-0.065, -0.005)
Gompertz, Survival, Bootstrap, Difference 1.0 0.406 (0.388, 0.424) 0.418 (0.408, 0.428) 0.369 (0.344, 0.390) 0.012 (-0.009, 0.031) -0.037 (-0.070, -0.006)
Gompertz, Survival, Bootstrap, Difference 3.0 0.311 (0.292, 0.329) 0.322 (0.312, 0.334) 0.275 (0.251, 0.296) 0.011 (-0.009, 0.031) -0.037 (-0.070, -0.007)
Gompertz, Survival, No bootstrap, Ratio 0.5 0.545 (0.528, 0.562) 0.556 (0.546, 0.565) 0.511 (0.489, 0.533) 1.020 (1.001, 1.039) 0.937 (0.874, 1.005)
Gompertz, Survival, No bootstrap, Ratio 1.0 0.406 (0.388, 0.425) 0.418 (0.408, 0.428) 0.369 (0.346, 0.393) 1.029 (1.001, 1.058) 0.908 (0.819, 1.007)
Gompertz, Survival, No bootstrap, Ratio 3.0 0.311 (0.293, 0.330) 0.322 (0.312, 0.333) 0.275 (0.253, 0.298) 1.037 (1.001, 1.075) 0.883 (0.772, 1.009)
Gompertz, Survival, Bootstrap, Ratio 0.5 0.545 (0.527, 0.561) 0.556 (0.546, 0.565) 0.511 (0.487, 0.531) 1.020 (0.986, 1.054) 0.937 (0.882, 0.991)
Gompertz, Survival, Bootstrap, Ratio 1.0 0.406 (0.388, 0.424) 0.418 (0.408, 0.428) 0.369 (0.344, 0.390) 1.029 (0.979, 1.079) 0.908 (0.833, 0.985)
Gompertz, Survival, Bootstrap, Ratio 3.0 0.311 (0.292, 0.329) 0.322 (0.312, 0.334) 0.275 (0.251, 0.296) 1.037 (0.973, 1.106) 0.883 (0.788, 0.976)
Royston-Parmar 04, RMST, No bootstrap, Difference 0.5 0.366 (0.363, 0.369) 0.362 (0.360, 0.364) 0.337 (0.333, 0.341) -0.004 (-0.005, -0.003) -0.029 (-0.033, -0.024)
Royston-Parmar 04, RMST, No bootstrap, Difference 1.0 0.488 (0.481, 0.495) 0.478 (0.474, 0.483) 0.418 (0.411, 0.426) -0.010 (-0.013, -0.006) -0.070 (-0.081, -0.058)
Royston-Parmar 04, RMST, No bootstrap, Difference 3.0 0.556 (0.545, 0.567) 0.541 (0.533, 0.548) 0.448 (0.438, 0.458) -0.015 (-0.021, -0.010) -0.108 (-0.126, -0.090)
Royston-Parmar 04, RMST, Bootstrap, Difference 0.5 0.366 (0.363, 0.369) 0.362 (0.360, 0.364) 0.337 (0.334, 0.341) -0.004 (-0.007, 0.000) -0.029 (-0.034, -0.024)
Royston-Parmar 04, RMST, Bootstrap, Difference 1.0 0.488 (0.481, 0.495) 0.478 (0.473, 0.483) 0.418 (0.411, 0.425) -0.010 (-0.017, -0.002) -0.070 (-0.079, -0.060)
Royston-Parmar 04, RMST, Bootstrap, Difference 3.0 0.556 (0.545, 0.566) 0.541 (0.533, 0.548) 0.448 (0.438, 0.458) -0.015 (-0.027, -0.002) -0.108 (-0.121, -0.094)
Royston-Parmar 04, RMST, No bootstrap, Ratio 0.5 0.366 (0.363, 0.369) 0.362 (0.360, 0.364) 0.337 (0.333, 0.341) 0.989 (0.986, 0.993) 0.921 (0.909, 0.934)
Royston-Parmar 04, RMST, No bootstrap, Ratio 1.0 0.488 (0.481, 0.495) 0.478 (0.474, 0.483) 0.418 (0.411, 0.426) 0.980 (0.973, 0.987) 0.857 (0.836, 0.879)
Royston-Parmar 04, RMST, No bootstrap, Ratio 3.0 0.556 (0.545, 0.567) 0.541 (0.533, 0.548) 0.448 (0.438, 0.458) 0.972 (0.963, 0.982) 0.806 (0.778, 0.835)
Royston-Parmar 04, RMST, Bootstrap, Ratio 0.5 0.366 (0.363, 0.369) 0.362 (0.360, 0.364) 0.337 (0.334, 0.341) 0.989 (0.981, 1.000) 0.921 (0.909, 0.935)
Royston-Parmar 04, RMST, Bootstrap, Ratio 1.0 0.488 (0.481, 0.495) 0.478 (0.473, 0.483) 0.418 (0.411, 0.425) 0.980 (0.966, 0.997) 0.857 (0.840, 0.876)
Royston-Parmar 04, RMST, Bootstrap, Ratio 3.0 0.556 (0.545, 0.566) 0.541 (0.533, 0.548) 0.448 (0.438, 0.458) 0.972 (0.951, 0.997) 0.806 (0.786, 0.830)
Royston-Parmar 04, Survival, No bootstrap, Difference 0.5 0.405 (0.396, 0.415) 0.392 (0.385, 0.399) 0.307 (0.297, 0.318) -0.014 (-0.018, -0.009) -0.098 (-0.114, -0.082)
Royston-Parmar 04, Survival, No bootstrap, Difference 1.0 0.135 (0.129, 0.142) 0.127 (0.122, 0.132) 0.075 (0.069, 0.080) -0.009 (-0.012, -0.006) -0.061 (-0.071, -0.051)
Royston-Parmar 04, Survival, No bootstrap, Difference 3.0 0.004 (0.003, 0.004) 0.003 (0.003, 0.004) 0.001 (0.001, 0.001) -0.000 (-0.001, -0.000) -0.003 (-0.004, -0.002)
Royston-Parmar 04, Survival, Bootstrap, Difference 0.5 0.405 (0.395, 0.414) 0.392 (0.385, 0.398) 0.307 (0.296, 0.317) -0.014 (-0.024, -0.002) -0.098 (-0.111, -0.083)
Royston-Parmar 04, Survival, Bootstrap, Difference 1.0 0.135 (0.129, 0.142) 0.127 (0.122, 0.131) 0.075 (0.069, 0.080) -0.009 (-0.017, -0.001) -0.061 (-0.068, -0.053)
Royston-Parmar 04, Survival, Bootstrap, Difference 3.0 0.004 (0.003, 0.004) 0.003 (0.003, 0.004) 0.001 (0.001, 0.001) -0.000 (-0.001, 0.000) -0.003 (-0.003, -0.002)
Royston-Parmar 04, Survival, No bootstrap, Ratio 0.5 0.405 (0.396, 0.415) 0.392 (0.385, 0.399) 0.307 (0.297, 0.318) 0.966 (0.955, 0.977) 0.758 (0.724, 0.794)
Royston-Parmar 04, Survival, No bootstrap, Ratio 1.0 0.135 (0.129, 0.142) 0.127 (0.122, 0.132) 0.075 (0.069, 0.080) 0.935 (0.914, 0.957) 0.550 (0.498, 0.608)
Royston-Parmar 04, Survival, No bootstrap, Ratio 3.0 0.004 (0.003, 0.004) 0.003 (0.003, 0.004) 0.001 (0.001, 0.001) 0.888 (0.848, 0.930) 0.213 (0.164, 0.278)
Royston-Parmar 04, Survival, Bootstrap, Ratio 0.5 0.405 (0.395, 0.414) 0.392 (0.385, 0.398) 0.307 (0.296, 0.317) 0.966 (0.941, 0.996) 0.758 (0.728, 0.791)
Royston-Parmar 04, Survival, Bootstrap, Ratio 1.0 0.135 (0.129, 0.142) 0.127 (0.122, 0.131) 0.075 (0.069, 0.080) 0.935 (0.880, 0.993) 0.550 (0.508, 0.599)
Royston-Parmar 04, Survival, Bootstrap, Ratio 3.0 0.004 (0.003, 0.004) 0.003 (0.003, 0.004) 0.001 (0.001, 0.001) 0.888 (0.760, 1.064) 0.213 (0.175, 0.265)

Time for this code chunk to run: 0 minutes


Session Info

June 29, 2023

Show code
Sys.getenv("R_LIBS_USER")
[1] "C:/Users/CISS Fondecyt/OneDrive/Documentos/R/win-library/4.1"
Show code
rstudioapi::getSourceEditorContext()
Document Context: 
- id:        '2728DA37'
- path:      'C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/joel/analisis_joel23.Rmd'
- contents:  <467 rows>
Document Selection:
- [427, 27] -- [427, 27]: ''
Show code
save.image("__analisis_joel23.RData")

unlink("*_cache", recursive = T, force = T, expand = TRUE)

sesion_info <- devtools::session_info()
dplyr::select(
  tibble::as_tibble(sesion_info$packages),
  c(package, loadedversion, source)
) %>% 
  DT::datatable(filter = 'top', colnames = c('Row number' =1,'Variable' = 2, 'Percentage'= 3),
              caption = htmltools::tags$caption(
        style = 'caption-side: top; text-align: left;',
        '', htmltools::em('Paquetes estadísticos utilizados')),
      options=list(
initComplete = htmlwidgets::JS(
      "function(settings, json) {",
      "$(this.api().tables().body()).css({'font-size': '80%'});",
      "}")))

Time for this code chunk to run: 0.1 minutes