Brochopred OR
OR_V9 <- read.csv2("C:/Users/mallah.s/Desktop/Stats et Theses/Brochopred_Arvin/OR_V9.csv", stringsAsFactors=TRUE)
BP_OR<-OR_V9
dependent_1A="OMS_1_2"
dependent_1B="OMS_3"
dependent_2A="Histo_ADK_oui"
dependent_2B="Histo_ADK_Non"
dependent_3A="Histo_Epidermoide_oui"
dependent_3B="Histo_Epidermoide_Non"
dependent_4A="traitement_Oui"
dependent_4B="Traitement_Non"
dependent_5A="Stade_NYHA_II"
dependent_5B="Stade_NYHA_III"
dependent_5C="Stade_NYHA_IV"
dependent_6A="Site_Trachee_carene_Oui"
dependent_6B="Site_Trachee_carene_Non"
dependent_7A="BP_OR$Programmee"
dependent_7B="BP_OR$Urgente"
dependent_8A="BP_OR$Atelectasie_RT_Oui"
dependent_8B="BP_OR$Atelectasie_RT_Non"
explanatory_1 = c("Delta.Total.QLQ30.V3.V1","Delta.Total.QLQ30_amelioration","Delta.Total.LC13.V3.V1"," Delta.Total.LC13.V3.V2_amelioration","Delta.QOL","Delta_QOL_amelioration","Delta.DyspnŽe.QLQ30","Delta.Dyspnee.QLQ31_Amelioration","Delta.DyspnŽe.LC13", "Delta.Dyspnee.LC14_amelioration")
explanatory_2 = c("Delta.Total.QLQ30_amelioration","Delta.Total.LC13.V3.V2_amelioration","Delta_QOL_amelioration","Delta.Dyspnee.QLQ31_Amelioration", "Delta.Dyspnee.LC14_amelioration")
OR OMS 1-2
kable(res_glm_uni_1A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
1.11 (0.11-10.91, p=0.924) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.58 (0.06-5.76, p=0.626) |
| Delta_QOL_ameliorationOui |
1.11 (0.11-10.91, p=0.924) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.90 (0.09-8.82, p=0.924) |
| Delta.Dyspnee.LC14_ameliorationOui |
1.38 (0.14-13.58, p=0.773) |
OR OMS 3
kable(res_glm_uni_1B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.90 (0.09-8.82, p=0.924) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
1.71 (0.17-17.10, p=0.626) |
| Delta_QOL_ameliorationOui |
0.90 (0.09-8.82, p=0.924) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
1.11 (0.11-10.91, p=0.924) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.73 (0.07-7.14, p=0.773) |
OR Histo ADK: Oui
kable(res_glm_uni_2A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.42 (0.07-2.20, p=0.309) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.50 (0.08-2.75, p=0.434) |
| Delta_QOL_ameliorationOui |
0.08 (0.01-0.48, p=0.012) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.27 (0.04-1.46, p=0.141) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.13 (0.02-0.75, p=0.031) |
OR Histo ADK : Non
kable(res_glm_uni_2B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
2.40 (0.46-14.05, p=0.309) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
2.00 (0.36-12.71, p=0.434) |
| Delta_QOL_ameliorationOui |
13.33 (2.07-130.86, p=0.012) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
3.73 (0.69-24.59, p=0.141) |
| Delta.Dyspnee.LC14_ameliorationOui |
7.78 (1.33-60.08, p=0.031) |
OR Histo epidermoîde: Oui
kable(res_glm_uni_3A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
1.25 (0.23-7.06, p=0.795) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.67 (0.11-3.73, p=0.649) |
| Delta_QOL_ameliorationOui |
2.67 (0.49-17.18, p=0.270) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.80 (0.14-4.33, p=0.795) |
| Delta.Dyspnee.LC14_ameliorationOui |
4.67 (0.79-39.77, p=0.111) |
OR Epidermoide: Non
kable(res_glm_uni_3B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.80 (0.14-4.33, p=0.795) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
1.50 (0.27-9.51, p=0.649) |
| Delta_QOL_ameliorationOui |
0.38 (0.06-2.05, p=0.270) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
1.25 (0.23-7.06, p=0.795) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.21 (0.03-1.27, p=0.111) |
OR traitement : Oui
kable(res_glm_uni_4A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
1.25 (0.23-7.06, p=0.795) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
1.44 (0.25-8.24, p=0.676) |
| Delta_QOL_ameliorationOui |
1.25 (0.23-7.06, p=0.795) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.80 (0.14-4.33, p=0.795) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.94 (0.17-5.28, p=0.940) |
OR traitement: Non
kable(res_glm_uni_4B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.80 (0.14-4.33, p=0.795) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.69 (0.12-3.95, p=0.676) |
| Delta_QOL_ameliorationOui |
0.80 (0.14-4.33, p=0.795) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
1.25 (0.23-7.06, p=0.795) |
| Delta.Dyspnee.LC14_ameliorationOui |
1.07 (0.19-5.88, p=0.940) |
OR Stade NYHA II
kable(res_glm_uni_5A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.88 (0.15-5.03, p=0.879) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
2.00 (0.34-12.24, p=0.438) |
| Delta_QOL_ameliorationOui |
4.50 (0.75-38.59, p=0.122) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
6.00 (0.98-52.70, p=0.068) |
| Delta.Dyspnee.LC14_ameliorationOui |
1.46 (0.26-9.31, p=0.673) |
OR Stade NYHA III
kable(res_glm_uni_5B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
1.33 (0.22-8.71, p=0.753) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
1.25 (0.19-7.77, p=0.809) |
| Delta_QOL_ameliorationOui |
1.33 (0.22-8.71, p=0.753) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.75 (0.11-4.52, p=0.753) |
| Delta.Dyspnee.LC14_ameliorationOui |
2.50 (0.40-21.30, p=0.347) |
OR Stade NYHA IV
kable(res_glm_uni_5C,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.88 (0.15-5.03, p=0.879) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.38 (0.05-2.31, p=0.318) |
| Delta_QOL_ameliorationOui |
0.17 (0.02-1.02, p=0.068) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.22 (0.03-1.34, p=0.122) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.30 (0.04-1.72, p=0.187) |
OR : trachée ou caréne : Oui
kable(res_glm_uni_6A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.83 (0.16-4.35, p=0.827) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
1.25 (0.23-7.06, p=0.795) |
| Delta_QOL_ameliorationOui |
3.50 (0.66-21.61, p=0.153) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.60 (0.11-3.09, p=0.538) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.57 (0.10-3.00, p=0.511) |
OR : trachée ou caréne : Non
kable(res_glm_uni_6B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
1.20 (0.23-6.41, p=0.827) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.80 (0.14-4.33, p=0.795) |
| Delta_QOL_ameliorationOui |
0.29 (0.05-1.52, p=0.153) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
1.68 (0.32-9.23, p=0.538) |
| Delta.Dyspnee.LC14_ameliorationOui |
1.75 (0.33-9.93, p=0.511) |
OR Programmée
kable(res_glm_uni_7A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
1.12 (0.16-7.70, p=0.901) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.55 (0.08-3.76, p=0.528) |
| Delta_QOL_ameliorationOui |
2.86 (0.43-25.09, p=0.292) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
2.25 (0.34-19.60, p=0.414) |
| Delta.Dyspnee.LC14_ameliorationOui |
1.43 (0.21-9.89, p=0.708) |
OR Urgente
kable(res_glm_uni_7B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.89 (0.13-6.06, p=0.901) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
1.83 (0.27-12.94, p=0.528) |
| Delta_QOL_ameliorationOui |
0.35 (0.04-2.32, p=0.292) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.44 (0.05-2.94, p=0.414) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.70 (0.10-4.78, p=0.708) |
library(rmarkdown)
library(markdown)
library(tidyverse)
library(finalfit)
library(survival)
library(survminer)
library(survminer,quietly=True)
OR Atelectasie : Oui
kable(res_glm_uni_8A,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
0.57 (0.10-3.00, p=0.511) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
0.44 (0.07-2.43, p=0.353) |
| Delta_QOL_ameliorationOui |
0.57 (0.10-3.00, p=0.511) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
0.42 (0.07-2.20, p=0.309) |
| Delta.Dyspnee.LC14_ameliorationOui |
0.78 (0.14-4.14, p=0.768) |
OR atelectasie : Non
kable(res_glm_uni_8B,row.names=FALSE, align=c("l", "l", "r", "r", "r", "r"))
| Delta.Total.QLQ30_ameliorationOui |
1.75 (0.33-9.93, p=0.511) |
| Delta.Total.LC13.V3.V2_ameliorationOui |
2.25 (0.41-13.37, p=0.353) |
| Delta_QOL_ameliorationOui |
1.75 (0.33-9.93, p=0.511) |
| Delta.Dyspnee.QLQ31_AmeliorationOui |
2.40 (0.46-14.05, p=0.309) |
| Delta.Dyspnee.LC14_ameliorationOui |
1.29 (0.24-7.21, p=0.768) |
PFS tous les patients
PFS avec ou sans traitement:
PFS_ttt<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~traitement_Oui,data=OR_V9)
PFS_ttt
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## traitement_Oui, data = OR_V9)
##
## n events median 0.95LCL 0.95UCL
## traitement_Oui=0 14 13 18.90 12.43 35.7
## traitement_Oui=1 9 7 2.07 1.83 NA
Courbe de survie avec ou sans traitement
ggsurvplot(PFS_ttt,xlab="Time(month)")

PFS_total QLQ30
PFS_tot_QlQ30<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Total.QLQ30_amelioration,data=OR_V9)
PFS_tot_QlQ30
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Total.QLQ30_amelioration, data = OR_V9)
##
## n events median 0.95LCL 0.95UCL
## Delta.Total.QLQ30_amelioration=Non 11 10 3.07 1.07 NA
## Delta.Total.QLQ30_amelioration=Oui 12 10 12.90 4.73 NA
Courbe de survie avec ou sans amélioration total QLq30
ggsurvplot(PFS_tot_QlQ30,xlab="Time(month)")

PFS_ Total.LC13.V3.V2
PFS_tot_LC13<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Total.LC13.V3.V2_amelioration,data=OR_V9)
PFS_tot_LC13
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Total.LC13.V3.V2_amelioration, data = OR_V9)
##
## n events median 0.95LCL 0.95UCL
## Delta.Total.LC13.V3.V2_amelioration=Non 14 12 13.30 2.33 NA
## Delta.Total.LC13.V3.V2_amelioration=Oui 9 8 7.07 2.07 NA
Courbe de survie avec ou sans amélioration Total.LC13.V3.V2
ggsurvplot(PFS_tot_LC13,xlab="Time(month)")

PFS_ Delta.Dyspnee.QLQ31_Amelioration
PFS_dysp_Qlq30<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Dyspnee.QLQ31_Amelioration,data=OR_V9)
PFS_dysp_Qlq30
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Dyspnee.QLQ31_Amelioration, data = OR_V9)
##
## n events median 0.95LCL 0.95UCL
## Delta.Dyspnee.QLQ31_Amelioration=Non 12 10 13.3 1.07 NA
## Delta.Dyspnee.QLQ31_Amelioration=Oui 11 10 12.4 4.73 NA
Courbe de survie avec ou sans amélioration Delta.Dyspnee.QLQ31
ggsurvplot(PFS_dysp_Qlq30,xlab="Time(month)")

PFS_ Dyspnee .LC13_amelioration
PFS_dys_LC13<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Dyspnee.LC14_amelioration,data=OR_V9)
PFS_dys_LC13
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Dyspnee.LC14_amelioration, data = OR_V9)
##
## n events median 0.95LCL 0.95UCL
## Delta.Dyspnee.LC14_amelioration=Non 10 9 2.47 1.07 NA
## Delta.Dyspnee.LC14_amelioration=Oui 13 11 14.17 4.73 NA
Courbe de survie avec ou sans amélioration dyspnée LC13
ggsurvplot(PFS_dys_LC13,xlab="Time(month)")

PFS_ QOL
PFS_QOL<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta_QOL_amelioration,data=OR_V9)
PFS_QOL
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta_QOL_amelioration, data = OR_V9)
##
## n events median 0.95LCL 0.95UCL
## Delta_QOL_amelioration=Non 11 8 12.4 1.07 NA
## Delta_QOL_amelioration=Oui 12 12 10.6 2.33 NA
Courbe de survie avec ou sans amélioration dyspnée LC13
ggsurvplot(PFS_QOL,xlab="Time(month)")

PFS patients traités
OR_V9_TTT <- read.csv2("C:/Users/mallah.s/Desktop/Stats et Theses/Brochopred_Arvin/OR_V9_TTT.csv", stringsAsFactors=TRUE)
PFS patients traités
PFS_ttt1<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~1,data=OR_V9_TTT)
PFS_ttt1
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## 1, data = OR_V9_TTT)
##
## n events median 0.95LCL 0.95UCL
## [1,] 9 7 2.07 1.83 NA
Courbe de survie avec ou sans traitement
ggsurvplot(PFS_ttt1,xlab="Time(month)")

PFS_total QLQ30_ttt
PFS_tot_QlQ30_t<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Total.QLQ30_amelioration,data=OR_V9_TTT)
PFS_tot_QlQ30_t
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Total.QLQ30_amelioration, data = OR_V9_TTT)
##
## n events median 0.95LCL 0.95UCL
## Delta.Total.QLQ30_amelioration=Non 4 3 1.63 0.90 NA
## Delta.Total.QLQ30_amelioration=Oui 5 4 2.07 1.87 NA
Courbe de survie avec ou sans amélioration total QLq30_ttt
ggsurvplot(PFS_tot_QlQ30_t,xlab="Time(month)")

PFS_ Total.LC13.V3.V2_ttt
PFS_tot_LC13_t<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Total.LC13.V3.V2_amelioration,data=OR_V9_TTT)
PFS_tot_LC13_t
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Total.LC13.V3.V2_amelioration, data = OR_V9_TTT)
##
## n events median 0.95LCL 0.95UCL
## Delta.Total.LC13.V3.V2_amelioration=Non 5 3 2.33 0.93 NA
## Delta.Total.LC13.V3.V2_amelioration=Oui 4 4 1.97 1.83 NA
Courbe de survie avec ou sans amélioration Total.LC13.V3.V2
ggsurvplot(PFS_tot_LC13_t,xlab="Time(month)")

PFS_ Delta.Dyspnee.QLQ31_Amelioration_ttt
PFS_dysp_Qlq30_t<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Dyspnee.QLQ31_Amelioration,data=OR_V9_TTT)
PFS_dysp_Qlq30_t
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Dyspnee.QLQ31_Amelioration, data = OR_V9_TTT)
##
## n events median 0.95LCL 0.95UCL
## Delta.Dyspnee.QLQ31_Amelioration=Non 5 3 1.83 0.93 NA
## Delta.Dyspnee.QLQ31_Amelioration=Oui 4 4 2.20 1.87 NA
Courbe de survie avec ou sans amélioration Delta.Dyspnee.QLQ31
ggsurvplot(PFS_dysp_Qlq30_t,xlab="Time(month)")

PFS_ Dyspnee .LC13_amelioration_ttt
PFS_dys_LC13_t<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta.Dyspnee.LC14_amelioration,data=OR_V9_TTT)
PFS_dys_LC13_t
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta.Dyspnee.LC14_amelioration, data = OR_V9_TTT)
##
## n events median 0.95LCL 0.95UCL
## Delta.Dyspnee.LC14_amelioration=Non 4 3 1.40 0.90 NA
## Delta.Dyspnee.LC14_amelioration=Oui 5 4 2.33 2.07 NA
Courbe de survie avec ou sans amélioration dyspnée LC13
ggsurvplot(PFS_dys_LC13_t,xlab="Time(month)")

PFS_ QOL ttt
PFS_QOL_t<-survfit(Surv(delais_DC_dernieres_nouvelles_mois,DCD)~Delta_QOL_amelioration,data=OR_V9_TTT)
PFS_QOL_t
## Call: survfit(formula = Surv(delais_DC_dernieres_nouvelles_mois, DCD) ~
## Delta_QOL_amelioration, data = OR_V9_TTT)
##
## n events median 0.95LCL 0.95UCL
## Delta_QOL_amelioration=Non 4 2 0.93 0.90 NA
## Delta_QOL_amelioration=Oui 5 5 2.07 1.87 NA
Courbe de survie avec ou sans amélioration dyspnée LC13
ggsurvplot(PFS_QOL_t,xlab="Time(month)")
