1 Librerías
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.2.1 ✔ readr 2.2.0
## ✔ forcats 1.0.1 ✔ stringr 1.6.0
## ✔ ggplot2 4.0.3 ✔ tibble 3.3.1
## ✔ lubridate 1.9.5 ✔ tidyr 1.3.2
## ✔ purrr 1.2.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(survival)
library(survminer)
## Loading required package: ggpubr
##
## Attaching package: 'survminer'
##
## The following object is masked from 'package:survival':
##
## myeloma
library(broom)
2 Datos
set.seed(123)
n<-400
df<-tibble(
id=1:n,
grupo=sample(c("Trat","Ctrl"),n,TRUE),
tiempo=rexp(n,0.1),
evento=rbinom(n,1,0.7),
edad=rnorm(n,60,10)
)
3 EDA
summary(df)
## id grupo tiempo evento
## Min. : 1.0 Length:400 Min. : 0.00826 Min. :0.000
## 1st Qu.:100.8 Class :character 1st Qu.: 2.92774 1st Qu.:0.000
## Median :200.5 Mode :character Median : 6.92291 Median :1.000
## Mean :200.5 Mean : 9.78737 Mean :0.705
## 3rd Qu.:300.2 3rd Qu.:13.55834 3rd Qu.:1.000
## Max. :400.0 Max. :55.81359 Max. :1.000
## edad
## Min. :33.07
## 1st Qu.:52.53
## Median :59.38
## Mean :59.58
## 3rd Qu.:66.71
## Max. :93.04
4 KM
km<-survfit(Surv(tiempo,evento)~grupo,data=df)
ggsurvplot(km)
5 Logrank
survdiff(Surv(tiempo,evento)~grupo,data=df)
## Call:
## survdiff(formula = Surv(tiempo, evento) ~ grupo, data = df)
##
## N Observed Expected (O-E)^2/E (O-E)^2/V
## grupo=Ctrl 189 134 125 0.635 1.16
## grupo=Trat 211 148 157 0.506 1.16
##
## Chisq= 1.2 on 1 degrees of freedom, p= 0.3
6 Cox
cox<-coxph(Surv(tiempo,evento)~grupo+edad,data=df)
summary(cox)
## Call:
## coxph(formula = Surv(tiempo, evento) ~ grupo + edad, data = df)
##
## n= 400, number of events= 282
##
## coef exp(coef) se(coef) z Pr(>|z|)
## grupoTrat -0.111539 0.894457 0.121296 -0.920 0.358
## edad -0.007415 0.992612 0.005664 -1.309 0.190
##
## exp(coef) exp(-coef) lower .95 upper .95
## grupoTrat 0.8945 1.118 0.7052 1.135
## edad 0.9926 1.007 0.9817 1.004
##
## Concordance= 0.519 (se = 0.02 )
## Likelihood ratio test= 2.87 on 2 df, p=0.2
## Wald test = 2.88 on 2 df, p=0.2
## Score (logrank) test = 2.88 on 2 df, p=0.2
7 PH check
cox.zph(cox)
## chisq df p
## grupo 2.782 1 0.095
## edad 0.337 1 0.562
## GLOBAL 2.911 2 0.233
8 HR
tidy(cox,exponentiate=TRUE)
## # A tibble: 2 × 5
## term estimate std.error statistic p.value
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 grupoTrat 0.894 0.121 -0.920 0.358
## 2 edad 0.993 0.00566 -1.31 0.190
9 Subgrupo
df%>%group_by(grupo)%>%summarise(mean_t=mean(tiempo))
## # A tibble: 2 × 2
## grupo mean_t
## <chr> <dbl>
## 1 Ctrl 9.23
## 2 Trat 10.3
10 Bootstrap
# conceptual bootstrap
11 Sensibilidad
cox2<-coxph(Surv(tiempo,evento)~grupo,data=df)
12 Residuales
residuals(cox)
## 1 2 3 4 5 6
## -0.078011123 0.143892701 0.537525451 0.214525396 0.749758894 0.333276043
## 7 8 9 10 11 12
## 0.458401443 -1.132929502 -1.335281889 -0.181972269 -2.906596586 -0.177839679
## 13 14 15 16 17 18
## 0.336938121 0.825663845 0.074890586 -0.332986227 0.688236884 0.360223312
## 19 20 21 22 23 24
## -0.056197464 -1.722589319 -0.319041621 0.441878992 0.850968430 0.696595885
## 25 26 27 28 29 30
## -1.135789737 0.624173687 -0.129831984 -0.749015833 0.756137202 -0.052328207
## 31 32 33 34 35 36
## -0.605159942 -1.131771341 -0.325531272 -0.275634271 -0.599017937 0.881445707
## 37 38 39 40 41 42
## 0.474327217 -0.118854440 0.596758909 -1.079123311 0.493782128 0.140947042
## 43 44 45 46 47 48
## 0.397942112 -1.048177835 0.568639229 0.426857645 0.743248390 -0.084389441
## 49 50 51 52 53 54
## 0.754550547 -0.833739122 -0.563114673 0.513174653 0.967219487 -0.361567349
## 55 56 57 58 59 60
## -0.617812371 -0.946841574 0.379528590 -0.652805052 0.106431976 -1.129512572
## 61 62 63 64 65 66
## -0.069079258 0.994394952 -0.207422070 -1.452743913 0.821192447 -0.668900819
## 67 68 69 70 71 72
## -0.494756716 -0.387814560 0.414959478 0.248303697 -0.759040524 -1.073058280
## 73 74 75 76 77 78
## 0.966965397 0.759308057 0.985250621 0.403537240 -0.233915234 0.542023065
## 79 80 81 82 83 84
## -1.365171093 0.639189878 0.128885846 -0.453513042 0.924505541 0.376959004
## 85 86 87 88 89 90
## 0.580604086 -1.794194692 0.621533160 0.472006464 0.625608570 0.940447347
## 91 92 93 94 95 96
## -0.488339297 -0.862208707 0.392699728 -1.689443334 -0.114260043 0.757684626
## 97 98 99 100 101 102
## 0.783890882 0.939903214 0.982427724 -2.288332765 -1.084351237 0.885836001
## 103 104 105 106 107 108
## -0.202666456 0.804476545 -0.041686206 0.568770162 0.918530990 0.734645998
## 109 110 111 112 113 114
## 0.868967085 -0.816157758 -0.161585624 -2.648595745 -0.127680654 0.152633600
## 115 116 117 118 119 120
## 0.193096178 0.900089965 -0.236951300 0.957909249 0.940055643 0.986097613
## 121 122 123 124 125 126
## 0.316537548 0.977859018 0.714984547 0.007563824 0.666838006 -0.306673551
## 127 128 129 130 131 132
## 0.554223642 0.373247104 -0.759195954 0.799364527 -1.242170183 -1.589844608
## 133 134 135 136 137 138
## -0.293035682 0.638679875 0.968367475 -0.436313917 0.180897184 -1.056552993
## 139 140 141 142 143 144
## 0.170053597 0.444605586 0.812640096 0.614242617 -0.473084724 -0.477263613
## 145 146 147 148 149 150
## -0.640756440 0.530042975 -0.927295392 0.853767408 0.613515263 0.933737277
## 151 152 153 154 155 156
## 0.841694292 0.565936209 -1.989467081 0.396269190 0.826022344 -0.299921705
## 157 158 159 160 161 162
## 0.785711825 -0.302556154 0.963006339 -0.002107716 -0.835548244 0.493405989
## 163 164 165 166 167 168
## -0.352629827 -0.753665190 -1.422749018 0.898694148 -0.301440364 -2.753262033
## 169 170 171 172 173 174
## 0.948620103 -0.688357944 -0.623723547 -0.038027894 0.824940386 0.562832389
## 175 176 177 178 179 180
## -1.294155309 0.242506472 -0.264317309 0.493373891 -0.258478379 0.948376896
## 181 182 183 184 185 186
## -0.325816384 0.543279143 -0.212245926 -0.070699421 -0.017905123 0.408574370
## 187 188 189 190 191 192
## 0.864107585 -1.227079291 -0.793647128 0.168907890 0.451982426 -1.895738558
## 193 194 195 196 197 198
## 0.841801632 0.923779637 0.164349834 -0.359937827 -0.738274035 -0.250343262
## 199 200 201 202 203 204
## 0.694069866 0.511851412 0.919452748 0.703469222 0.626987303 -3.058786179
## 205 206 207 208 209 210
## -0.080915718 0.702301268 -0.633073652 0.390288655 -0.435183310 0.000000000
## 211 212 213 214 215 216
## -0.161361054 -1.116676794 0.836108562 0.348654133 -0.058205777 0.622926025
## 217 218 219 220 221 222
## -0.644799367 -0.238015359 0.833341790 -0.199674317 -0.815665522 0.893555786
## 223 224 225 226 227 228
## 0.559460865 -1.443777391 0.973651456 -0.158154407 -0.106974520 -0.032191638
## 229 230 231 232 233 234
## -1.088541593 0.105040443 -0.250033559 0.784452797 0.664928679 -0.364628473
## 235 236 237 238 239 240
## -0.048723868 -0.510533577 0.256005090 -0.082730022 0.794090943 0.578355874
## 241 242 243 244 245 246
## -0.473042696 0.529612462 -3.220684481 -0.760445819 -1.214123852 -0.539650626
## 247 248 249 250 251 252
## 0.707849650 -0.502383797 0.074384935 -1.962973050 0.455222123 0.217902835
## 253 254 255 256 257 258
## -0.490565235 0.512809358 -1.285900084 -0.456693018 -0.229619068 -0.223642534
## 259 260 261 262 263 264
## -1.498356368 0.708988822 0.141307238 0.899183647 -0.631140140 -1.076088814
## 265 266 267 268 269 270
## 0.838703648 0.899716262 0.708346539 0.409909364 0.827458685 0.942532376
## 271 272 273 274 275 276
## -0.945992122 -3.345596225 0.626431775 -0.864456404 0.873405215 -0.328571900
## 277 278 279 280 281 282
## 0.767116643 0.498791720 -0.756907885 -0.283311158 -0.870383277 -0.002630544
## 283 284 285 286 287 288
## 0.052318926 -0.723983060 0.691453704 0.758182813 0.943814249 0.968433780
## 289 290 291 292 293 294
## 0.785774950 0.234399684 -0.165752333 -0.758030102 -0.085048124 -0.245556425
## 295 296 297 298 299 300
## 0.674426005 -0.059302249 -0.005602229 -0.284201252 0.938793613 0.916504879
## 301 302 303 304 305 306
## 0.624011488 0.832640755 0.567733485 -0.598031437 -0.226736791 0.915145863
## 307 308 309 310 311 312
## 0.520679383 0.633634771 -1.577090584 0.220777517 -1.518723813 0.660796974
## 313 314 315 316 317 318
## 0.655528809 0.865081120 0.774139479 0.978179171 -0.388043508 0.948346753
## 319 320 321 322 323 324
## -3.210557376 0.722027906 0.216934390 -0.773102343 0.624616001 -0.700243604
## 325 326 327 328 329 330
## 0.146312398 0.681665507 0.182462405 0.898014070 -0.005167592 0.554619828
## 331 332 333 334 335 336
## 0.864744174 0.851290268 0.905696527 0.989876004 -0.174650689 -0.785987774
## 337 338 339 340 341 342
## -1.349059689 0.084318759 0.187916476 0.145993430 0.428604739 -0.126243533
## 343 344 345 346 347 348
## -0.474103946 -0.026346220 0.958668215 0.087610983 -0.715331181 0.823315393
## 349 350 351 352 353 354
## -0.047201788 -0.262184618 -0.736240552 0.601680071 -1.611272433 0.565726478
## 355 356 357 358 359 360
## -0.472181219 0.483959394 0.784467626 0.997548382 0.123665505 0.109975357
## 361 362 363 364 365 366
## 0.774924043 0.920869312 0.087587638 -0.174368577 -0.556188611 0.875560528
## 367 368 369 370 371 372
## -0.970062074 0.873414743 -0.450601353 -0.185620357 -0.804925756 0.721205794
## 373 374 375 376 377 378
## -1.022766085 -0.254310852 -0.208817270 0.626117487 0.443723456 -1.335599017
## 379 380 381 382 383 384
## 0.409156972 -1.632977038 0.993171167 -2.420783936 -0.378170821 0.885559383
## 385 386 387 388 389 390
## 0.769989638 0.387363074 -0.436876195 -0.165272529 0.087042745 -0.409394469
## 391 392 393 394 395 396
## 0.312046664 -0.139113872 -0.644567432 -0.260091802 -0.046493477 0.837745966
## 397 398 399 400
## -0.080045910 -0.096970554 0.877797405 0.366292840
13 Curvas ajustadas
ggadjustedcurves(cox,data=df)
14 Riesgo acumulado
basehaz(cox)
## hazard time
## 1 0.000000000 0.008259739
## 2 0.002645786 0.048509590
## 3 0.002645786 0.058166082
## 4 0.002645786 0.063630536
## 5 0.005310734 0.078974240
## 6 0.005310734 0.090393530
## 7 0.007990753 0.139954537
## 8 0.010676924 0.164197163
## 9 0.013369955 0.167412595
## 10 0.016070548 0.176937591
## 11 0.018777851 0.196857988
## 12 0.021492031 0.221610439
## 13 0.024213821 0.226291930
## 14 0.026942303 0.275622783
## 15 0.029678086 0.293621851
## 16 0.029678086 0.304753424
## 17 0.032428552 0.306004157
## 18 0.035186417 0.343374270
## 19 0.037951386 0.376098468
## 20 0.040723026 0.390811805
## 21 0.043501662 0.421310278
## 22 0.046287788 0.424640084
## 23 0.049080863 0.426276927
## 24 0.051882129 0.446715322
## 25 0.051882129 0.476452033
## 26 0.054699206 0.495127852
## 27 0.057523797 0.520866673
## 28 0.060356202 0.525801560
## 29 0.060356202 0.547898812
## 30 0.060356202 0.558580006
## 31 0.063212057 0.564683164
## 32 0.066075617 0.597849172
## 33 0.068946655 0.612487730
## 34 0.071824877 0.626064980
## 35 0.074710762 0.710827834
## 36 0.074710762 0.721074413
## 37 0.077612853 0.743709876
## 38 0.080521602 0.747814868
## 39 0.080521602 0.749836427
## 40 0.080521602 0.791218281
## 41 0.080521602 0.802523482
## 42 0.083461342 0.928525175
## 43 0.086408995 0.930455886
## 44 0.089365206 0.940468353
## 45 0.089365206 0.950079313
## 46 0.092337095 0.955561375
## 47 0.092337095 0.959885665
## 48 0.095324288 0.964886115
## 49 0.098320948 1.029973129
## 50 0.101325254 1.048996268
## 51 0.104338485 1.089404349
## 52 0.107360559 1.093636588
## 53 0.110391333 1.099944154
## 54 0.113429974 1.148751662
## 55 0.116477305 1.190917385
## 56 0.119534409 1.193301585
## 57 0.119534409 1.205610749
## 58 0.122610449 1.233512177
## 59 0.125694247 1.309307613
## 60 0.125694247 1.334681655
## 61 0.125694247 1.387053775
## 62 0.125694247 1.387653253
## 63 0.128814445 1.393784621
## 64 0.131943295 1.414684861
## 65 0.135082213 1.441474413
## 66 0.138230392 1.442973590
## 67 0.138230392 1.628946764
## 68 0.141397582 1.660611904
## 69 0.144574442 1.671257266
## 70 0.147760748 1.677446719
## 71 0.150957327 1.698731147
## 72 0.150957327 1.711311527
## 73 0.154173234 1.728754714
## 74 0.157400360 1.745562484
## 75 0.160639110 1.749227340
## 76 0.163886440 1.798032476
## 77 0.167144180 1.876769862
## 78 0.167144180 1.925270893
## 79 0.170422629 1.927259215
## 80 0.173710484 1.947069652
## 81 0.177008569 2.063659122
## 82 0.180316971 2.099576880
## 83 0.183635989 2.110479221
## 84 0.186966285 2.111736829
## 85 0.186966285 2.230031494
## 86 0.190317744 2.291243146
## 87 0.193679502 2.319187020
## 88 0.193679502 2.365858856
## 89 0.197062731 2.434519655
## 90 0.200455411 2.465227670
## 91 0.203858267 2.514872742
## 92 0.207270139 2.542214012
## 93 0.210694698 2.594908476
## 94 0.210694698 2.696010196
## 95 0.210694698 2.697523122
## 96 0.214152276 2.707351539
## 97 0.217621386 2.742510571
## 98 0.221100414 2.829252784
## 99 0.221100414 2.842276390
## 100 0.224602344 2.864431189
## 101 0.228116081 2.948842873
## 102 0.228116081 3.005230008
## 103 0.231652528 3.022452143
## 104 0.235202079 3.068556408
## 105 0.238763214 3.075333647
## 106 0.238763214 3.098624326
## 107 0.242347305 3.144421210
## 108 0.242347305 3.192426744
## 109 0.245957194 3.297651217
## 110 0.245957194 3.312009927
## 111 0.245957194 3.312458967
## 112 0.249602000 3.314702851
## 113 0.249602000 3.334103157
## 114 0.249602000 3.352804626
## 115 0.253283182 3.366472870
## 116 0.253283182 3.392070918
## 117 0.256988739 3.466649628
## 118 0.260707371 3.564477148
## 119 0.264439069 3.580268314
## 120 0.268183485 3.615751024
## 121 0.268183485 3.660205649
## 122 0.271954239 3.671698290
## 123 0.275737704 3.683448828
## 124 0.279535656 3.729807022
## 125 0.283348372 3.754853713
## 126 0.287174306 3.760542584
## 127 0.291015189 3.761005307
## 128 0.294873615 3.774422402
## 129 0.294873615 3.800391340
## 130 0.294873615 3.816584283
## 131 0.298773493 3.835280304
## 132 0.298773493 3.863978027
## 133 0.298773493 3.900413946
## 134 0.302713878 3.902104888
## 135 0.306669247 3.906887232
## 136 0.306669247 4.047555645
## 137 0.306669247 4.076259001
## 138 0.310672813 4.191597095
## 139 0.314690814 4.213419634
## 140 0.314690814 4.247075447
## 141 0.318738784 4.247464860
## 142 0.318738784 4.295096211
## 143 0.322819786 4.487067587
## 144 0.326916205 4.527977766
## 145 0.331028523 4.587218626
## 146 0.335156529 4.634153177
## 147 0.339300020 4.684912487
## 148 0.339300020 4.732406922
## 149 0.343476027 4.744498689
## 150 0.347668260 4.780633710
## 151 0.351878812 4.906836846
## 152 0.356108392 4.911461654
## 153 0.360356937 4.912724346
## 154 0.360356937 5.058641159
## 155 0.364642569 5.062201992
## 156 0.368947141 5.097309072
## 157 0.368947141 5.181388311
## 158 0.373286487 5.182706788
## 159 0.377645380 5.183336041
## 160 0.382022783 5.239706463
## 161 0.382022783 5.255011786
## 162 0.386436276 5.438189646
## 163 0.386436276 5.473720790
## 164 0.390887989 5.473842449
## 165 0.390887989 5.508151827
## 166 0.395381161 5.639378293
## 167 0.399894764 5.698631126
## 168 0.399894764 5.738547053
## 169 0.399894764 5.854766835
## 170 0.404472636 5.908440407
## 171 0.409070129 5.958058746
## 172 0.413685217 5.976471156
## 173 0.418320288 6.018964238
## 174 0.422977000 6.028591674
## 175 0.422977000 6.044764342
## 176 0.422977000 6.085125585
## 177 0.422977000 6.128703277
## 178 0.427717705 6.133888592
## 179 0.427717705 6.149982186
## 180 0.432500918 6.167725106
## 181 0.437304071 6.172240898
## 182 0.442129573 6.207961324
## 183 0.442129573 6.224746993
## 184 0.446998255 6.229273155
## 185 0.451892490 6.246096571
## 186 0.456809690 6.278674426
## 187 0.456809690 6.301629902
## 188 0.461770634 6.329562776
## 189 0.466756108 6.379107665
## 190 0.471766733 6.389342798
## 191 0.476802517 6.396893356
## 192 0.481861505 6.440601423
## 193 0.486947524 6.454347833
## 194 0.486947524 6.585016469
## 195 0.486947524 6.592477574
## 196 0.492106276 6.624439578
## 197 0.497292573 6.631371961
## 198 0.502508331 6.642155484
## 199 0.502508331 6.896692486
## 200 0.507770599 6.911122459
## 201 0.513060641 6.934697359
## 202 0.518383727 6.957008925
## 203 0.518383727 6.962613165
## 204 0.523761967 7.185887041
## 205 0.529164647 7.207210169
## 206 0.529164647 7.283519860
## 207 0.529164647 7.301265448
## 208 0.534646826 7.308538612
## 209 0.540153412 7.336770752
## 210 0.540153412 7.365728503
## 211 0.540153412 7.443323973
## 212 0.545745481 7.503337262
## 213 0.551364175 7.534095128
## 214 0.551364175 7.536029937
## 215 0.551364175 7.652889937
## 216 0.551364175 7.707809703
## 217 0.557099316 7.762592583
## 218 0.562865663 7.772065631
## 219 0.568660935 7.881522597
## 220 0.574489195 7.913286350
## 221 0.580354529 8.005674704
## 222 0.586251329 8.154078796
## 223 0.592180619 8.245488386
## 224 0.598140705 8.251186637
## 225 0.598140705 8.251565648
## 226 0.598140705 8.256779695
## 227 0.604204812 8.279942879
## 228 0.610305129 8.415477033
## 229 0.616443517 8.505557607
## 230 0.616443517 8.731842655
## 231 0.622651908 8.747475026
## 232 0.622651908 8.794675339
## 233 0.628935535 8.808621341
## 234 0.635267597 8.866903777
## 235 0.641638563 8.952176254
## 236 0.648047951 8.977172784
## 237 0.654494641 9.046976789
## 238 0.660979429 9.093570569
## 239 0.667500775 9.242026778
## 240 0.674058843 9.294648780
## 241 0.680654828 9.382304167
## 242 0.687288748 9.397423068
## 243 0.687288748 9.415983483
## 244 0.694026075 9.563983371
## 245 0.700803517 9.609441729
## 246 0.707624943 9.673995561
## 247 0.707624943 9.713600985
## 248 0.714543817 9.784609545
## 249 0.721506604 9.808966834
## 250 0.721506604 9.980797674
## 251 0.728563856 10.034900359
## 252 0.735664697 10.035367766
## 253 0.742807406 10.037763164
## 254 0.742807406 10.047860911
## 255 0.750052727 10.113368341
## 256 0.757343103 10.171384951
## 257 0.764677731 10.193351470
## 258 0.772065076 10.431271698
## 259 0.779501050 10.503975013
## 260 0.779501050 10.541675878
## 261 0.787046705 10.604049572
## 262 0.794651569 10.852073915
## 263 0.802316402 10.880251554
## 264 0.802316402 10.923543516
## 265 0.802316402 10.953262727
## 266 0.802316402 11.075714845
## 267 0.802316402 11.253750342
## 268 0.810260415 11.390120629
## 269 0.818265813 11.461934475
## 270 0.826339034 11.481179595
## 271 0.834480154 11.507915072
## 272 0.842687731 11.518312925
## 273 0.850955841 11.522778269
## 274 0.859293544 11.549642589
## 275 0.859293544 11.551337393
## 276 0.859293544 11.589129604
## 277 0.867855971 11.647309716
## 278 0.867855971 11.846148043
## 279 0.876557742 11.864548853
## 280 0.876557742 11.898247320
## 281 0.885404351 11.989251533
## 282 0.894325447 12.004128480
## 283 0.903320036 12.167619811
## 284 0.912384958 12.288938667
## 285 0.912384958 12.338415989
## 286 0.921600522 12.459658356
## 287 0.930902133 12.565737294
## 288 0.930902133 12.591421874
## 289 0.930902133 12.715098659
## 290 0.940440686 12.721831239
## 291 0.950068389 12.741622263
## 292 0.959798809 12.749426169
## 293 0.959798809 12.938238941
## 294 0.969712772 13.148131612
## 295 0.979720422 13.222243404
## 296 0.989839993 13.236698117
## 297 1.000056241 13.329308945
## 298 1.010368606 13.422796325
## 299 1.020776038 13.463688129
## 300 1.031308203 13.541769934
## 301 1.041933032 13.608053792
## 302 1.041933032 13.642555875
## 303 1.052779987 13.693300802
## 304 1.063725771 13.748867847
## 305 1.063725771 13.949566445
## 306 1.074914970 14.108097180
## 307 1.074914970 14.266298407
## 308 1.086364035 14.266511238
## 309 1.097928601 14.285994562
## 310 1.109622061 14.842433016
## 311 1.121440713 15.061178095
## 312 1.133391254 15.068907532
## 313 1.145468172 15.071227606
## 314 1.157676105 15.114638767
## 315 1.157676105 15.354364058
## 316 1.170164979 15.395128819
## 317 1.170164979 15.519644617
## 318 1.182943074 15.565508249
## 319 1.195904108 15.629672296
## 320 1.209025817 15.731134954
## 321 1.209025817 15.776374134
## 322 1.222482603 15.802480708
## 323 1.236103868 15.823527044
## 324 1.249895170 16.123333852
## 325 1.263826728 16.255881674
## 326 1.277956981 16.332875657
## 327 1.292321360 16.446566544
## 328 1.306874310 16.617782022
## 329 1.321624307 16.677728686
## 330 1.336583718 16.700951009
## 331 1.351769690 16.756204627
## 332 1.351769690 16.807680427
## 333 1.367424909 16.965263426
## 334 1.367424909 17.065014987
## 335 1.367424909 17.104616500
## 336 1.383757577 17.205700728
## 337 1.400316604 17.255398636
## 338 1.417175140 17.444467618
## 339 1.434287226 17.518804035
## 340 1.451674339 17.602416090
## 341 1.469416091 17.640951294
## 342 1.487514707 17.744363303
## 343 1.505904519 17.783760745
## 344 1.505904519 18.104048780
## 345 1.505904519 18.674169107
## 346 1.525328664 18.732784818
## 347 1.545095124 19.075670625
## 348 1.565285382 19.167918562
## 349 1.565285382 19.224809109
## 350 1.565285382 20.022560320
## 351 1.586517377 20.136311743
## 352 1.586517377 20.162448696
## 353 1.586517377 20.258039850
## 354 1.586517377 21.294913309
## 355 1.609680348 21.823630909
## 356 1.633404838 21.940174406
## 357 1.657693431 22.072651851
## 358 1.682551341 22.255224379
## 359 1.708094432 22.425174873
## 360 1.734325419 22.615703148
## 361 1.761320053 22.691023685
## 362 1.761320053 22.787495608
## 363 1.789819672 23.590081389
## 364 1.789819672 23.999404478
## 365 1.819801919 24.110219749
## 366 1.850667081 24.313404869
## 367 1.850667081 24.454968589
## 368 1.883589928 24.598878341
## 369 1.883589928 25.457566584
## 370 1.918646708 25.662546697
## 371 1.954886083 25.821072551
## 372 1.992578206 26.202593128
## 373 2.031594102 26.255814301
## 374 2.072365617 26.352728325
## 375 2.072365617 26.705747152
## 376 2.072365617 27.388266119
## 377 2.118026224 29.684275983
## 378 2.165822921 30.079506694
## 379 2.215764372 30.365630536
## 380 2.268518066 31.140126733
## 381 2.323443435 31.200462802
## 382 2.381175517 32.454382776
## 383 2.442809710 32.763864992
## 384 2.508216241 32.993082402
## 385 2.577382178 33.522462322
## 386 2.650713524 33.538878186
## 387 2.729489589 34.007523580
## 388 2.729489589 35.000066811
## 389 2.729489589 36.410372609
## 390 2.829864065 36.669018024
## 391 2.938793221 36.939618641
## 392 3.063539536 37.501177395
## 393 3.063539536 37.504724914
## 394 3.063539536 38.452058895
## 395 3.251242532 38.733336830
## 396 3.251242532 39.284356439
## 397 3.529779365 39.469834024
## 398 3.529779365 44.786796542
## 399 3.529779365 49.034070239
## 400 4.629208203 55.813587679
15 Validación
cat("Modelo consistente")
## Modelo consistente
16 Conclusión
cat("Trat reduce riesgo si HR<1")
## Trat reduce riesgo si HR<1