DEM7223 Homework 1
Tobias Kuhn
10/5/2020
Libraries
library(haven);
library(sas7bdat);
library(plyr);
library(dplyr);##
## Attaching package: 'dplyr'## The following objects are masked from 'package:plyr':
##
## arrange, count, desc, failwith, id, mutate, rename, summarise,
## summarize## The following objects are masked from 'package:stats':
##
## filter, lag## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, unionlibrary(tidyr);
library(haven);
library(survival);
library(car);## Loading required package: carData##
## Attaching package: 'car'## The following object is masked from 'package:dplyr':
##
## recodelibrary(muhaz);
library(survminer);## Loading required package: ggplot2## Loading required package: ggpubr##
## Attaching package: 'ggpubr'## The following object is masked from 'package:plyr':
##
## mutatelibrary(ggplot2)Loading PM2.5 data, merging, and removing single years
PM2010 = read.csv("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/pm25TX2010.csv");
PM2011 = read.csv("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/pm25TX2011.csv");
PM2012 = read.csv("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/pm25TX2012.csv");
PM2013 = read.csv("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/pm25TX2013.csv");
PM2014 = read.csv("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/pm25TX2014.csv");
PM2015 = read.csv("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/pm25TX2015.csv");
PM25ALL <- rbind(PM2010,PM2011,PM2012,PM2013,PM2014,PM2015);
rm(PM2010,PM2011,PM2012,PM2013,PM2014,PM2015)Loading Mortality Data, merging, and removing single years
Mor2010 <- read_sas("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/death2010.sas7bdat");
Mor2011 <- read_sas("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/death2011.sas7bdat");
Mor2012 <- read_sas("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/death2012.sas7bdat");
Mor2013 <- read_sas("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/death2013.sas7bdat");
Mor2014 <- read_sas("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/death2014.sas7bdat");
Mor2015 <- read_sas("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/death2015.sas7bdat");
MorALL <- rbind(Mor2010,Mor2011,Mor2012,Mor2013,Mor2014,Mor2015);
rm(Mor2010,Mor2011,Mor2012,Mor2013,Mor2014,Mor2015)Loading a county / county number list in order to merge the data
County = read.csv("C:/Users/tobik/OneDrive/Documents/Air Pollution Research/data/TXCounties.csv")Adding to County to Fips codes in mortality data in order to merge the data sets (the county name may be redundant now, but that was a good way to check if it worked correctly)
MorALL$RESIDCOUNTY_CODE <- as.numeric(MorALL$RESIDCOUNTY_CODE);
County$ï..Number <- as.numeric(County$ï..Number);
MorALL <- merge(MorALL, County, by.x = "RESIDCOUNTY_CODE", by.y = "ï..Number", all.x = TRUE, all.y = FALSE)
rm(County)Recoding dates (strptime did not work on all dates, so, I had to take it apart as string and recombine it so I can merge the data)
PM25ALL$NewDateM <- substr(PM25ALL$Date, start = 1, stop = 2);
PM25ALL$NewDateD <- substr(PM25ALL$Date, start = 4, stop = 5);
PM25ALL$NewDateY <- substr(PM25ALL$Date, start = 7, stop = 10);
PM25ALL$Date <- paste(PM25ALL$NewDateY,PM25ALL$NewDateM,PM25ALL$NewDateD, sep = "")Creating averages per county for fine dust emissions per day
PM25ALL <- PM25ALL %>%
dplyr::group_by(Date, COUNTY_CODE) %>%
dplyr::mutate(
PM25Aver = mean(Daily.Mean.PM2.5.Concentration)) %>%
dplyr::ungroup()Merging mortality data with PM2.5 averages
PM25Sub <- PM25ALL[,c(1,17,18,24)]
FromDustToDust <- merge(x = MorALL, y = PM25Sub, by.x = c("DATEOFDEATH","Fips"), by.y = c("Date","COUNTY_CODE"), all.x = TRUE)
rm(MorALL,PM25ALL,PM25Sub)Removing some lines due to missing data
FromDustToDust <- subset(FromDustToDust, complete.cases (PM25Aver,AGE))
FromDustToDust <- subset(FromDustToDust, AGE < 999)
FromDustToDust <- subset(FromDustToDust, DEATHRESULTOFINJURY < 9)Creating Dummy Variable for premature death (below 65 = 1; 65Plus = 0)
FromDustToDust$AGE <- as.numeric(FromDustToDust$AGE)
#FromDustToDust$PrematureDeath <- as.numeric(FromDustToDust$AGE<65)
FromDustToDust$One <- 1 #dummy variable for survfit
names(FromDustToDust)[names(FromDustToDust) == "DEATHRESULTOFINJURY"]<-"Injury"This is for testing purposes if database reduction is needed #{r} FDTD<- FromDustToDust %>% select(AGE,SEX,PM25Aver,PrematureDeath) #
fit1 <- survfit(formula = Surv(AGE,One) ~ SEX, data = FromDustToDust);
fit2 <- survfit(formula = Surv(AGE,One) ~ SEX + Injury, data = FromDustToDust);fit1<-glm(PrematureDeath~PM25Aver, data = FromDustToDust, na.action = na.omit, family = binomial)
summary(fit1)## Call: survfit(formula = Surv(AGE, One) ~ SEX, data = FromDustToDust)
##
## SEX=1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 1 1416857 6718 9.95e-01 5.77e-05 9.95e-01 9.95e-01
## 2 1410139 4598 9.92e-01 7.48e-05 9.92e-01 9.92e-01
## 3 1405541 3008 9.90e-01 8.40e-05 9.90e-01 9.90e-01
## 4 1402533 2228 9.88e-01 9.03e-05 9.88e-01 9.88e-01
## 5 1400305 1858 9.87e-01 9.51e-05 9.87e-01 9.87e-01
## 6 1398447 1398 9.86e-01 9.86e-05 9.86e-01 9.86e-01
## 7 1397049 1143 9.85e-01 1.01e-04 9.85e-01 9.85e-01
## 8 1395906 1153 9.84e-01 1.04e-04 9.84e-01 9.85e-01
## 9 1394753 1040 9.84e-01 1.06e-04 9.83e-01 9.84e-01
## 10 1393713 1156 9.83e-01 1.09e-04 9.83e-01 9.83e-01
## 11 1392557 1016 9.82e-01 1.11e-04 9.82e-01 9.82e-01
## 12 1391541 800 9.82e-01 1.13e-04 9.81e-01 9.82e-01
## 13 1390741 822 9.81e-01 1.15e-04 9.81e-01 9.81e-01
## 14 1389919 991 9.80e-01 1.17e-04 9.80e-01 9.81e-01
## 15 1388928 1332 9.79e-01 1.19e-04 9.79e-01 9.80e-01
## 16 1387596 1834 9.78e-01 1.23e-04 9.78e-01 9.78e-01
## 17 1385762 2366 9.76e-01 1.28e-04 9.76e-01 9.77e-01
## 18 1383396 3227 9.74e-01 1.33e-04 9.74e-01 9.74e-01
## 19 1380169 3754 9.71e-01 1.40e-04 9.71e-01 9.72e-01
## 20 1376415 4259 9.68e-01 1.47e-04 9.68e-01 9.69e-01
## 21 1372156 4317 9.65e-01 1.54e-04 9.65e-01 9.66e-01
## 22 1367839 5024 9.62e-01 1.61e-04 9.62e-01 9.62e-01
## 23 1362815 4579 9.59e-01 1.67e-04 9.58e-01 9.59e-01
## 24 1358236 4871 9.55e-01 1.74e-04 9.55e-01 9.56e-01
## 25 1353365 4690 9.52e-01 1.80e-04 9.52e-01 9.52e-01
## 26 1348675 4655 9.49e-01 1.86e-04 9.48e-01 9.49e-01
## 27 1344020 4856 9.45e-01 1.91e-04 9.45e-01 9.46e-01
## 28 1339164 4542 9.42e-01 1.96e-04 9.42e-01 9.42e-01
## 29 1334622 4995 9.38e-01 2.02e-04 9.38e-01 9.39e-01
## 30 1329627 4872 9.35e-01 2.07e-04 9.35e-01 9.35e-01
## 31 1324755 4858 9.32e-01 2.12e-04 9.31e-01 9.32e-01
## 32 1319897 5101 9.28e-01 2.17e-04 9.28e-01 9.28e-01
## 33 1314796 5128 9.24e-01 2.22e-04 9.24e-01 9.25e-01
## 34 1309668 5322 9.21e-01 2.27e-04 9.20e-01 9.21e-01
## 35 1304346 4811 9.17e-01 2.32e-04 9.17e-01 9.18e-01
## 36 1299535 5328 9.13e-01 2.36e-04 9.13e-01 9.14e-01
## 37 1294207 5325 9.10e-01 2.41e-04 9.09e-01 9.10e-01
## 38 1288882 5684 9.06e-01 2.46e-04 9.05e-01 9.06e-01
## 39 1283198 6282 9.01e-01 2.51e-04 9.01e-01 9.02e-01
## 40 1276916 6261 8.97e-01 2.56e-04 8.96e-01 8.97e-01
## 41 1270655 7027 8.92e-01 2.61e-04 8.91e-01 8.92e-01
## 42 1263628 7608 8.86e-01 2.67e-04 8.86e-01 8.87e-01
## 43 1256020 8147 8.81e-01 2.72e-04 8.80e-01 8.81e-01
## 44 1247873 8665 8.75e-01 2.78e-04 8.74e-01 8.75e-01
## 45 1239208 9089 8.68e-01 2.84e-04 8.68e-01 8.69e-01
## 46 1230119 10378 8.61e-01 2.91e-04 8.60e-01 8.61e-01
## 47 1219741 11465 8.53e-01 2.98e-04 8.52e-01 8.53e-01
## 48 1208276 12559 8.44e-01 3.05e-04 8.43e-01 8.45e-01
## 49 1195717 13587 8.34e-01 3.12e-04 8.34e-01 8.35e-01
## 50 1182130 15121 8.24e-01 3.20e-04 8.23e-01 8.24e-01
## 51 1167009 17088 8.12e-01 3.29e-04 8.11e-01 8.12e-01
## 52 1149921 17963 7.99e-01 3.37e-04 7.98e-01 8.00e-01
## 53 1131958 19578 7.85e-01 3.45e-04 7.84e-01 7.86e-01
## 54 1112380 20308 7.71e-01 3.53e-04 7.70e-01 7.71e-01
## 55 1092072 21713 7.55e-01 3.61e-04 7.55e-01 7.56e-01
## 56 1070359 23391 7.39e-01 3.69e-04 7.38e-01 7.40e-01
## 57 1046968 24629 7.22e-01 3.77e-04 7.21e-01 7.22e-01
## 58 1022339 24679 7.04e-01 3.83e-04 7.03e-01 7.05e-01
## 59 997660 25858 6.86e-01 3.90e-04 6.85e-01 6.87e-01
## 60 971802 25539 6.68e-01 3.96e-04 6.67e-01 6.69e-01
## 61 946263 27025 6.49e-01 4.01e-04 6.48e-01 6.50e-01
## 62 919238 27918 6.29e-01 4.06e-04 6.28e-01 6.30e-01
## 63 891320 28167 6.09e-01 4.10e-04 6.08e-01 6.10e-01
## 64 863153 27650 5.90e-01 4.13e-04 5.89e-01 5.90e-01
## 65 835503 28105 5.70e-01 4.16e-04 5.69e-01 5.71e-01
## 66 807398 27934 5.50e-01 4.18e-04 5.49e-01 5.51e-01
## 67 779464 27786 5.31e-01 4.19e-04 5.30e-01 5.31e-01
## 68 751678 27654 5.11e-01 4.20e-04 5.10e-01 5.12e-01
## 69 724024 26307 4.92e-01 4.20e-04 4.92e-01 4.93e-01
## 70 697717 27614 4.73e-01 4.19e-04 4.72e-01 4.74e-01
## 71 670103 27727 4.53e-01 4.18e-04 4.53e-01 4.54e-01
## 72 642376 27627 4.34e-01 4.16e-04 4.33e-01 4.35e-01
## 73 614749 27526 4.14e-01 4.14e-04 4.14e-01 4.15e-01
## 74 587223 27341 3.95e-01 4.11e-04 3.94e-01 3.96e-01
## 75 559882 27416 3.76e-01 4.07e-04 3.75e-01 3.77e-01
## 76 532466 28831 3.55e-01 4.02e-04 3.55e-01 3.56e-01
## 77 503635 28829 3.35e-01 3.97e-04 3.34e-01 3.36e-01
## 78 474806 29604 3.14e-01 3.90e-04 3.13e-01 3.15e-01
## 79 445202 30436 2.93e-01 3.82e-04 2.92e-01 2.93e-01
## 80 414766 30782 2.71e-01 3.73e-04 2.70e-01 2.72e-01
## 81 383984 32027 2.48e-01 3.63e-04 2.48e-01 2.49e-01
## 82 351957 31964 2.26e-01 3.51e-04 2.25e-01 2.27e-01
## 83 319993 32112 2.03e-01 3.38e-04 2.03e-01 2.04e-01
## 84 287881 31509 1.81e-01 3.23e-04 1.80e-01 1.82e-01
## 85 256372 31217 1.59e-01 3.07e-04 1.58e-01 1.60e-01
## 86 225155 30143 1.38e-01 2.89e-04 1.37e-01 1.38e-01
## 87 195012 29089 1.17e-01 2.70e-04 1.17e-01 1.18e-01
## 88 165923 27748 9.75e-02 2.49e-04 9.70e-02 9.80e-02
## 89 138175 25737 7.94e-02 2.27e-04 7.89e-02 7.98e-02
## 90 112438 22930 6.32e-02 2.04e-04 6.28e-02 6.36e-02
## 91 89508 20258 4.89e-02 1.81e-04 4.85e-02 4.92e-02
## 92 69250 16869 3.70e-02 1.59e-04 3.67e-02 3.73e-02
## 93 52381 14111 2.70e-02 1.36e-04 2.67e-02 2.73e-02
## 94 38270 10659 1.95e-02 1.16e-04 1.93e-02 1.97e-02
## 95 27611 8243 1.37e-02 9.76e-05 1.35e-02 1.39e-02
## 96 19368 6202 9.29e-03 8.06e-05 9.14e-03 9.45e-03
## 97 13166 4517 6.10e-03 6.54e-05 5.98e-03 6.23e-03
## 98 8649 2756 4.16e-03 5.41e-05 4.05e-03 4.27e-03
## 99 5893 2201 2.61e-03 4.28e-05 2.52e-03 2.69e-03
## 100 3692 1505 1.54e-03 3.30e-05 1.48e-03 1.61e-03
## 101 2187 855 9.40e-04 2.57e-05 8.91e-04 9.92e-04
## 102 1332 620 5.03e-04 1.88e-05 4.67e-04 5.41e-04
## 103 712 318 2.78e-04 1.40e-05 2.52e-04 3.07e-04
## 104 394 154 1.69e-04 1.09e-05 1.49e-04 1.92e-04
## 105 240 109 9.25e-05 8.08e-06 7.79e-05 1.10e-04
## 106 131 48 5.86e-05 6.43e-06 4.72e-05 7.26e-05
## 107 83 37 3.25e-05 4.79e-06 2.43e-05 4.33e-05
## 108 46 12 2.40e-05 4.12e-06 1.71e-05 3.36e-05
## 109 34 4 2.12e-05 3.87e-06 1.48e-05 3.03e-05
## 110 30 6 1.69e-05 3.46e-06 1.14e-05 2.53e-05
## 111 24 12 8.47e-06 2.44e-06 4.81e-06 1.49e-05
## 112 12 9 2.12e-06 1.22e-06 6.83e-07 6.57e-06
## 116 3 3 0.00e+00 NaN NA NA
##
## SEX=2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 1 1337168 5375 9.96e-01 5.47e-05 9.96e-01 9.96e-01
## 2 1331793 3696 9.93e-01 7.10e-05 9.93e-01 9.93e-01
## 3 1328097 2629 9.91e-01 8.05e-05 9.91e-01 9.91e-01
## 4 1325468 1682 9.90e-01 8.61e-05 9.90e-01 9.90e-01
## 5 1323786 1609 9.89e-01 9.11e-05 9.89e-01 9.89e-01
## 6 1322177 1300 9.88e-01 9.49e-05 9.88e-01 9.88e-01
## 7 1320877 812 9.87e-01 9.72e-05 9.87e-01 9.87e-01
## 8 1320065 872 9.87e-01 9.96e-05 9.86e-01 9.87e-01
## 9 1319193 785 9.86e-01 1.02e-04 9.86e-01 9.86e-01
## 10 1318408 812 9.85e-01 1.04e-04 9.85e-01 9.86e-01
## 11 1317596 647 9.85e-01 1.06e-04 9.85e-01 9.85e-01
## 12 1316949 610 9.84e-01 1.07e-04 9.84e-01 9.85e-01
## 13 1316339 736 9.84e-01 1.09e-04 9.84e-01 9.84e-01
## 14 1315603 632 9.83e-01 1.10e-04 9.83e-01 9.84e-01
## 15 1314971 839 9.83e-01 1.13e-04 9.83e-01 9.83e-01
## 16 1314132 899 9.82e-01 1.15e-04 9.82e-01 9.82e-01
## 17 1313233 1015 9.81e-01 1.17e-04 9.81e-01 9.82e-01
## 18 1312218 1210 9.80e-01 1.20e-04 9.80e-01 9.81e-01
## 19 1311008 1306 9.79e-01 1.23e-04 9.79e-01 9.80e-01
## 20 1309702 1524 9.78e-01 1.26e-04 9.78e-01 9.79e-01
## 21 1308178 1575 9.77e-01 1.29e-04 9.77e-01 9.77e-01
## 22 1306603 1624 9.76e-01 1.33e-04 9.76e-01 9.76e-01
## 23 1304979 1791 9.75e-01 1.36e-04 9.74e-01 9.75e-01
## 24 1303188 1829 9.73e-01 1.40e-04 9.73e-01 9.73e-01
## 25 1301359 1792 9.72e-01 1.43e-04 9.72e-01 9.72e-01
## 26 1299567 1940 9.70e-01 1.46e-04 9.70e-01 9.71e-01
## 27 1297627 2083 9.69e-01 1.50e-04 9.69e-01 9.69e-01
## 28 1295544 1928 9.67e-01 1.54e-04 9.67e-01 9.68e-01
## 29 1293616 2054 9.66e-01 1.57e-04 9.66e-01 9.66e-01
## 30 1291562 2226 9.64e-01 1.61e-04 9.64e-01 9.65e-01
## 31 1289336 2481 9.62e-01 1.65e-04 9.62e-01 9.63e-01
## 32 1286855 2640 9.60e-01 1.69e-04 9.60e-01 9.61e-01
## 33 1284215 2758 9.58e-01 1.73e-04 9.58e-01 9.59e-01
## 34 1281457 2783 9.56e-01 1.77e-04 9.56e-01 9.57e-01
## 35 1278674 2924 9.54e-01 1.81e-04 9.54e-01 9.54e-01
## 36 1275750 2905 9.52e-01 1.85e-04 9.52e-01 9.52e-01
## 37 1272845 3277 9.49e-01 1.89e-04 9.49e-01 9.50e-01
## 38 1269568 3429 9.47e-01 1.94e-04 9.47e-01 9.47e-01
## 39 1266139 3906 9.44e-01 1.99e-04 9.44e-01 9.44e-01
## 40 1262233 3829 9.41e-01 2.04e-04 9.41e-01 9.41e-01
## 41 1258404 4314 9.38e-01 2.09e-04 9.37e-01 9.38e-01
## 42 1254090 4693 9.34e-01 2.14e-04 9.34e-01 9.35e-01
## 43 1249397 5380 9.30e-01 2.20e-04 9.30e-01 9.31e-01
## 44 1244017 5307 9.26e-01 2.26e-04 9.26e-01 9.27e-01
## 45 1238710 5863 9.22e-01 2.32e-04 9.22e-01 9.22e-01
## 46 1232847 6471 9.17e-01 2.38e-04 9.17e-01 9.18e-01
## 47 1226376 7631 9.11e-01 2.46e-04 9.11e-01 9.12e-01
## 48 1218745 8049 9.05e-01 2.53e-04 9.05e-01 9.06e-01
## 49 1210696 9164 8.99e-01 2.61e-04 8.98e-01 8.99e-01
## 50 1201532 10188 8.91e-01 2.70e-04 8.90e-01 8.91e-01
## 51 1191344 11171 8.83e-01 2.78e-04 8.82e-01 8.83e-01
## 52 1180173 11175 8.74e-01 2.87e-04 8.74e-01 8.75e-01
## 53 1168998 11915 8.65e-01 2.95e-04 8.65e-01 8.66e-01
## 54 1157083 12897 8.56e-01 3.04e-04 8.55e-01 8.56e-01
## 55 1144186 13445 8.46e-01 3.12e-04 8.45e-01 8.46e-01
## 56 1130741 14657 8.35e-01 3.21e-04 8.34e-01 8.35e-01
## 57 1116084 15170 8.23e-01 3.30e-04 8.23e-01 8.24e-01
## 58 1100914 15803 8.11e-01 3.38e-04 8.11e-01 8.12e-01
## 59 1085111 16269 7.99e-01 3.46e-04 7.99e-01 8.00e-01
## 60 1068842 15817 7.88e-01 3.54e-04 7.87e-01 7.88e-01
## 61 1053025 16727 7.75e-01 3.61e-04 7.74e-01 7.76e-01
## 62 1036298 17890 7.62e-01 3.68e-04 7.61e-01 7.62e-01
## 63 1018408 19673 7.47e-01 3.76e-04 7.46e-01 7.48e-01
## 64 998735 19668 7.32e-01 3.83e-04 7.31e-01 7.33e-01
## 65 979067 19774 7.17e-01 3.89e-04 7.17e-01 7.18e-01
## 66 959293 20106 7.02e-01 3.95e-04 7.02e-01 7.03e-01
## 67 939187 20767 6.87e-01 4.01e-04 6.86e-01 6.88e-01
## 68 918420 20410 6.72e-01 4.06e-04 6.71e-01 6.72e-01
## 69 898010 20707 6.56e-01 4.11e-04 6.55e-01 6.57e-01
## 70 877303 21603 6.40e-01 4.15e-04 6.39e-01 6.41e-01
## 71 855700 22010 6.23e-01 4.19e-04 6.23e-01 6.24e-01
## 72 833690 22394 6.07e-01 4.22e-04 6.06e-01 6.08e-01
## 73 811296 22526 5.90e-01 4.25e-04 5.89e-01 5.91e-01
## 74 788770 23582 5.72e-01 4.28e-04 5.71e-01 5.73e-01
## 75 765188 25066 5.53e-01 4.30e-04 5.53e-01 5.54e-01
## 76 740122 26522 5.34e-01 4.31e-04 5.33e-01 5.35e-01
## 77 713600 26504 5.14e-01 4.32e-04 5.13e-01 5.15e-01
## 78 687096 29157 4.92e-01 4.32e-04 4.91e-01 4.93e-01
## 79 657939 31289 4.69e-01 4.32e-04 4.68e-01 4.69e-01
## 80 626650 32186 4.45e-01 4.30e-04 4.44e-01 4.45e-01
## 81 594464 34261 4.19e-01 4.27e-04 4.18e-01 4.20e-01
## 82 560203 37000 3.91e-01 4.22e-04 3.90e-01 3.92e-01
## 83 523203 36768 3.64e-01 4.16e-04 3.63e-01 3.65e-01
## 84 486435 39297 3.34e-01 4.08e-04 3.34e-01 3.35e-01
## 85 447138 41565 3.03e-01 3.98e-04 3.03e-01 3.04e-01
## 86 405573 42186 2.72e-01 3.85e-04 2.71e-01 2.73e-01
## 87 363387 42589 2.40e-01 3.69e-04 2.39e-01 2.41e-01
## 88 320798 42435 2.08e-01 3.51e-04 2.07e-01 2.09e-01
## 89 278363 40791 1.78e-01 3.31e-04 1.77e-01 1.78e-01
## 90 237572 39295 1.48e-01 3.07e-04 1.48e-01 1.49e-01
## 91 198277 35761 1.22e-01 2.83e-04 1.21e-01 1.22e-01
## 92 162516 31736 9.78e-02 2.57e-04 9.73e-02 9.83e-02
## 93 130780 27719 7.71e-02 2.31e-04 7.66e-02 7.75e-02
## 94 103061 22766 6.00e-02 2.05e-04 5.96e-02 6.05e-02
## 95 80295 19362 4.56e-02 1.80e-04 4.52e-02 4.59e-02
## 96 60933 16184 3.35e-02 1.56e-04 3.32e-02 3.38e-02
## 97 44749 12712 2.40e-02 1.32e-04 2.37e-02 2.42e-02
## 98 32037 9467 1.69e-02 1.11e-04 1.67e-02 1.71e-02
## 99 22570 6847 1.18e-02 9.32e-05 1.16e-02 1.19e-02
## 100 15723 5525 7.63e-03 7.52e-05 7.48e-03 7.78e-03
## 101 10198 3628 4.91e-03 6.05e-05 4.80e-03 5.03e-03
## 102 6570 2294 3.20e-03 4.88e-05 3.10e-03 3.29e-03
## 103 4276 1676 1.94e-03 3.81e-05 1.87e-03 2.02e-03
## 104 2600 1038 1.17e-03 2.95e-05 1.11e-03 1.23e-03
## 105 1562 721 6.29e-04 2.17e-05 5.88e-04 6.73e-04
## 106 841 333 3.80e-04 1.69e-05 3.48e-04 4.14e-04
## 107 508 250 1.93e-04 1.20e-05 1.71e-04 2.18e-04
## 108 258 142 8.68e-05 8.05e-06 7.23e-05 1.04e-04
## 109 116 72 3.29e-05 4.96e-06 2.45e-05 4.42e-05
## 110 44 25 1.42e-05 3.26e-06 9.06e-06 2.23e-05
## 111 19 3 1.20e-05 2.99e-06 7.33e-06 1.95e-05
## 112 16 6 7.48e-06 2.36e-06 4.02e-06 1.39e-05
## 113 10 9 7.48e-07 7.48e-07 1.05e-07 5.31e-06
## 114 1 1 0.00e+00 NaN NA NAsummary(fit2)## Call: survfit(formula = Surv(AGE, One) ~ SEX + Injury, data = FromDustToDust)
##
## SEX=1, Injury=1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 1 171172 754 9.96e-01 1.60e-04 9.95e-01 9.96e-01
## 2 170418 543 9.92e-01 2.10e-04 9.92e-01 9.93e-01
## 3 169875 503 9.89e-01 2.47e-04 9.89e-01 9.90e-01
## 4 169372 496 9.87e-01 2.78e-04 9.86e-01 9.87e-01
## 5 168876 412 9.84e-01 3.02e-04 9.84e-01 9.85e-01
## 6 168464 357 9.82e-01 3.21e-04 9.81e-01 9.83e-01
## 7 168107 254 9.81e-01 3.33e-04 9.80e-01 9.81e-01
## 8 167853 262 9.79e-01 3.46e-04 9.78e-01 9.80e-01
## 9 167591 294 9.77e-01 3.60e-04 9.77e-01 9.78e-01
## 10 167297 232 9.76e-01 3.70e-04 9.75e-01 9.77e-01
## 11 167065 266 9.74e-01 3.81e-04 9.74e-01 9.75e-01
## 12 166799 180 9.73e-01 3.89e-04 9.73e-01 9.74e-01
## 13 166619 275 9.72e-01 4.00e-04 9.71e-01 9.73e-01
## 14 166344 420 9.69e-01 4.17e-04 9.69e-01 9.70e-01
## 15 165924 581 9.66e-01 4.38e-04 9.65e-01 9.67e-01
## 16 165343 1143 9.59e-01 4.78e-04 9.58e-01 9.60e-01
## 17 164200 1712 9.49e-01 5.30e-04 9.48e-01 9.50e-01
## 18 162488 2516 9.35e-01 5.98e-04 9.33e-01 9.36e-01
## 19 159972 2990 9.17e-01 6.66e-04 9.16e-01 9.18e-01
## 20 156982 3269 8.98e-01 7.32e-04 8.97e-01 8.99e-01
## 21 153713 3414 8.78e-01 7.91e-04 8.77e-01 8.80e-01
## 22 150299 3907 8.55e-01 8.50e-04 8.54e-01 8.57e-01
## 23 146392 3549 8.34e-01 8.98e-04 8.33e-01 8.36e-01
## 24 142843 3725 8.13e-01 9.43e-04 8.11e-01 8.15e-01
## 25 139118 3346 7.93e-01 9.79e-04 7.91e-01 7.95e-01
## 26 135772 3523 7.73e-01 1.01e-03 7.71e-01 7.75e-01
## 27 132249 3470 7.52e-01 1.04e-03 7.50e-01 7.54e-01
## 28 128779 3253 7.33e-01 1.07e-03 7.31e-01 7.35e-01
## 29 125526 3381 7.14e-01 1.09e-03 7.11e-01 7.16e-01
## 30 122145 3263 6.95e-01 1.11e-03 6.92e-01 6.97e-01
## 31 118882 3264 6.75e-01 1.13e-03 6.73e-01 6.78e-01
## 32 115618 3131 6.57e-01 1.15e-03 6.55e-01 6.59e-01
## 33 112487 3070 6.39e-01 1.16e-03 6.37e-01 6.42e-01
## 34 109417 2765 6.23e-01 1.17e-03 6.21e-01 6.25e-01
## 35 106652 2530 6.08e-01 1.18e-03 6.06e-01 6.11e-01
## 36 104122 2688 5.93e-01 1.19e-03 5.90e-01 5.95e-01
## 37 101434 2576 5.78e-01 1.19e-03 5.75e-01 5.80e-01
## 38 98858 2490 5.63e-01 1.20e-03 5.61e-01 5.65e-01
## 39 96368 2486 5.48e-01 1.20e-03 5.46e-01 5.51e-01
## 40 93882 2453 5.34e-01 1.21e-03 5.32e-01 5.37e-01
## 41 91429 2753 5.18e-01 1.21e-03 5.16e-01 5.20e-01
## 42 88676 2546 5.03e-01 1.21e-03 5.01e-01 5.06e-01
## 43 86130 2600 4.88e-01 1.21e-03 4.86e-01 4.90e-01
## 44 83530 2577 4.73e-01 1.21e-03 4.71e-01 4.75e-01
## 45 80953 2454 4.59e-01 1.20e-03 4.56e-01 4.61e-01
## 46 78499 2742 4.43e-01 1.20e-03 4.40e-01 4.45e-01
## 47 75757 2658 4.27e-01 1.20e-03 4.25e-01 4.29e-01
## 48 73099 2936 4.10e-01 1.19e-03 4.08e-01 4.12e-01
## 49 70163 2820 3.93e-01 1.18e-03 3.91e-01 3.96e-01
## 50 67343 2786 3.77e-01 1.17e-03 3.75e-01 3.79e-01
## 51 64557 2905 3.60e-01 1.16e-03 3.58e-01 3.62e-01
## 52 61652 2891 3.43e-01 1.15e-03 3.41e-01 3.46e-01
## 53 58761 3034 3.26e-01 1.13e-03 3.23e-01 3.28e-01
## 54 55727 2775 3.09e-01 1.12e-03 3.07e-01 3.12e-01
## 55 52952 2628 2.94e-01 1.10e-03 2.92e-01 2.96e-01
## 56 50324 2673 2.78e-01 1.08e-03 2.76e-01 2.81e-01
## 57 47651 2866 2.62e-01 1.06e-03 2.60e-01 2.64e-01
## 58 44785 2504 2.47e-01 1.04e-03 2.45e-01 2.49e-01
## 59 42281 2686 2.31e-01 1.02e-03 2.29e-01 2.33e-01
## 60 39595 2488 2.17e-01 9.96e-04 2.15e-01 2.19e-01
## 61 37107 2170 2.04e-01 9.74e-04 2.02e-01 2.06e-01
## 62 34937 2004 1.92e-01 9.53e-04 1.91e-01 1.94e-01
## 63 32933 1671 1.83e-01 9.34e-04 1.81e-01 1.84e-01
## 64 31262 1580 1.73e-01 9.15e-04 1.72e-01 1.75e-01
## 65 29682 1496 1.65e-01 8.96e-04 1.63e-01 1.66e-01
## 66 28186 1406 1.56e-01 8.78e-04 1.55e-01 1.58e-01
## 67 26780 1291 1.49e-01 8.60e-04 1.47e-01 1.51e-01
## 68 25489 1162 1.42e-01 8.44e-04 1.40e-01 1.44e-01
## 69 24327 1077 1.36e-01 8.28e-04 1.34e-01 1.37e-01
## 70 23250 1240 1.29e-01 8.09e-04 1.27e-01 1.30e-01
## 71 22010 1069 1.22e-01 7.92e-04 1.21e-01 1.24e-01
## 72 20941 1035 1.16e-01 7.75e-04 1.15e-01 1.18e-01
## 73 19906 853 1.11e-01 7.60e-04 1.10e-01 1.13e-01
## 74 19053 1055 1.05e-01 7.41e-04 1.04e-01 1.07e-01
## 75 17998 810 1.00e-01 7.26e-04 9.90e-02 1.02e-01
## 76 17188 916 9.51e-02 7.09e-04 9.37e-02 9.65e-02
## 77 16272 977 8.94e-02 6.89e-04 8.80e-02 9.07e-02
## 78 15295 969 8.37e-02 6.69e-04 8.24e-02 8.50e-02
## 79 14326 1012 7.78e-02 6.47e-04 7.65e-02 7.91e-02
## 80 13314 985 7.20e-02 6.25e-04 7.08e-02 7.33e-02
## 81 12329 1089 6.57e-02 5.99e-04 6.45e-02 6.68e-02
## 82 11240 997 5.98e-02 5.73e-04 5.87e-02 6.10e-02
## 83 10243 969 5.42e-02 5.47e-04 5.31e-02 5.53e-02
## 84 9274 939 4.87e-02 5.20e-04 4.77e-02 4.97e-02
## 85 8335 1083 4.24e-02 4.87e-04 4.14e-02 4.33e-02
## 86 7252 1001 3.65e-02 4.53e-04 3.56e-02 3.74e-02
## 87 6251 1033 3.05e-02 4.16e-04 2.97e-02 3.13e-02
## 88 5218 958 2.49e-02 3.77e-04 2.42e-02 2.56e-02
## 89 4260 818 2.01e-02 3.39e-04 1.95e-02 2.08e-02
## 90 3442 714 1.59e-02 3.03e-04 1.54e-02 1.65e-02
## 91 2728 652 1.21e-02 2.65e-04 1.16e-02 1.27e-02
## 92 2076 520 9.09e-03 2.29e-04 8.65e-03 9.55e-03
## 93 1556 438 6.53e-03 1.95e-04 6.16e-03 6.92e-03
## 94 1118 342 4.53e-03 1.62e-04 4.23e-03 4.86e-03
## 95 776 254 3.05e-03 1.33e-04 2.80e-03 3.32e-03
## 96 522 169 2.06e-03 1.10e-04 1.86e-03 2.29e-03
## 97 353 133 1.29e-03 8.66e-05 1.13e-03 1.47e-03
## 98 220 64 9.11e-04 7.29e-05 7.79e-04 1.07e-03
## 99 156 76 4.67e-04 5.22e-05 3.75e-04 5.82e-04
## 100 80 34 2.69e-04 3.96e-05 2.01e-04 3.59e-04
## 101 46 30 9.35e-05 2.34e-05 5.73e-05 1.53e-04
## 102 16 4 7.01e-05 2.02e-05 3.98e-05 1.23e-04
## 103 12 6 3.51e-05 1.43e-05 1.57e-05 7.80e-05
## 104 6 3 1.75e-05 1.01e-05 5.65e-06 5.43e-05
## 111 3 3 0.00e+00 NaN NA NA
##
## SEX=1, Injury=2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 1 1245685 5964 9.95e-01 6.18e-05 9.95e-01 9.95e-01
## 2 1239721 4055 9.92e-01 8.00e-05 9.92e-01 9.92e-01
## 3 1235666 2505 9.90e-01 8.94e-05 9.90e-01 9.90e-01
## 4 1233161 1732 9.89e-01 9.53e-05 9.88e-01 9.89e-01
## 5 1231429 1446 9.87e-01 1.00e-04 9.87e-01 9.88e-01
## 6 1229983 1041 9.87e-01 1.03e-04 9.86e-01 9.87e-01
## 7 1228942 889 9.86e-01 1.06e-04 9.86e-01 9.86e-01
## 8 1228053 891 9.85e-01 1.08e-04 9.85e-01 9.85e-01
## 9 1227162 746 9.85e-01 1.11e-04 9.84e-01 9.85e-01
## 10 1226416 924 9.84e-01 1.13e-04 9.84e-01 9.84e-01
## 11 1225492 750 9.83e-01 1.15e-04 9.83e-01 9.83e-01
## 12 1224742 620 9.83e-01 1.17e-04 9.82e-01 9.83e-01
## 13 1224122 547 9.82e-01 1.18e-04 9.82e-01 9.82e-01
## 14 1223575 571 9.82e-01 1.20e-04 9.82e-01 9.82e-01
## 15 1223004 751 9.81e-01 1.22e-04 9.81e-01 9.81e-01
## 16 1222253 691 9.81e-01 1.23e-04 9.80e-01 9.81e-01
## 17 1221562 654 9.80e-01 1.25e-04 9.80e-01 9.80e-01
## 18 1220908 711 9.80e-01 1.27e-04 9.79e-01 9.80e-01
## 19 1220197 764 9.79e-01 1.29e-04 9.79e-01 9.79e-01
## 20 1219433 990 9.78e-01 1.31e-04 9.78e-01 9.78e-01
## 21 1218443 903 9.77e-01 1.33e-04 9.77e-01 9.78e-01
## 22 1217540 1117 9.77e-01 1.36e-04 9.76e-01 9.77e-01
## 23 1216423 1030 9.76e-01 1.38e-04 9.75e-01 9.76e-01
## 24 1215393 1146 9.75e-01 1.41e-04 9.74e-01 9.75e-01
## 25 1214247 1344 9.74e-01 1.43e-04 9.73e-01 9.74e-01
## 26 1212903 1132 9.73e-01 1.46e-04 9.72e-01 9.73e-01
## 27 1211771 1386 9.72e-01 1.49e-04 9.71e-01 9.72e-01
## 28 1210385 1289 9.71e-01 1.51e-04 9.70e-01 9.71e-01
## 29 1209096 1614 9.69e-01 1.54e-04 9.69e-01 9.70e-01
## 30 1207482 1609 9.68e-01 1.58e-04 9.68e-01 9.68e-01
## 31 1205873 1594 9.67e-01 1.61e-04 9.66e-01 9.67e-01
## 32 1204279 1970 9.65e-01 1.64e-04 9.65e-01 9.66e-01
## 33 1202309 2058 9.64e-01 1.68e-04 9.63e-01 9.64e-01
## 34 1200251 2557 9.61e-01 1.72e-04 9.61e-01 9.62e-01
## 35 1197694 2281 9.60e-01 1.76e-04 9.59e-01 9.60e-01
## 36 1195413 2640 9.58e-01 1.81e-04 9.57e-01 9.58e-01
## 37 1192773 2749 9.55e-01 1.85e-04 9.55e-01 9.56e-01
## 38 1190024 3194 9.53e-01 1.90e-04 9.52e-01 9.53e-01
## 39 1186830 3796 9.50e-01 1.96e-04 9.49e-01 9.50e-01
## 40 1183034 3808 9.47e-01 2.01e-04 9.46e-01 9.47e-01
## 41 1179226 4274 9.43e-01 2.07e-04 9.43e-01 9.44e-01
## 42 1174952 5062 9.39e-01 2.14e-04 9.39e-01 9.40e-01
## 43 1169890 5547 9.35e-01 2.21e-04 9.34e-01 9.35e-01
## 44 1164343 6088 9.30e-01 2.29e-04 9.29e-01 9.30e-01
## 45 1158255 6635 9.24e-01 2.37e-04 9.24e-01 9.25e-01
## 46 1151620 7636 9.18e-01 2.45e-04 9.18e-01 9.19e-01
## 47 1143984 8807 9.11e-01 2.55e-04 9.11e-01 9.12e-01
## 48 1135177 9623 9.04e-01 2.64e-04 9.03e-01 9.04e-01
## 49 1125554 10767 8.95e-01 2.75e-04 8.94e-01 8.95e-01
## 50 1114787 12335 8.85e-01 2.86e-04 8.84e-01 8.86e-01
## 51 1102452 14183 8.74e-01 2.98e-04 8.73e-01 8.74e-01
## 52 1088269 15072 8.62e-01 3.09e-04 8.61e-01 8.62e-01
## 53 1073197 16544 8.48e-01 3.21e-04 8.48e-01 8.49e-01
## 54 1056653 17533 8.34e-01 3.33e-04 8.34e-01 8.35e-01
## 55 1039120 19085 8.19e-01 3.45e-04 8.18e-01 8.20e-01
## 56 1020035 20718 8.02e-01 3.57e-04 8.02e-01 8.03e-01
## 57 999317 21763 7.85e-01 3.68e-04 7.84e-01 7.85e-01
## 58 977554 22175 7.67e-01 3.79e-04 7.66e-01 7.68e-01
## 59 955379 23172 7.48e-01 3.89e-04 7.48e-01 7.49e-01
## 60 932207 23051 7.30e-01 3.98e-04 7.29e-01 7.31e-01
## 61 909156 24855 7.10e-01 4.07e-04 7.09e-01 7.11e-01
## 62 884301 25914 6.89e-01 4.15e-04 6.88e-01 6.90e-01
## 63 858387 26496 6.68e-01 4.22e-04 6.67e-01 6.69e-01
## 64 831891 26070 6.47e-01 4.28e-04 6.46e-01 6.48e-01
## 65 805821 26609 6.26e-01 4.34e-04 6.25e-01 6.26e-01
## 66 779212 26528 6.04e-01 4.38e-04 6.03e-01 6.05e-01
## 67 752684 26495 5.83e-01 4.42e-04 5.82e-01 5.84e-01
## 68 726189 26492 5.62e-01 4.45e-04 5.61e-01 5.63e-01
## 69 699697 25230 5.41e-01 4.46e-04 5.41e-01 5.42e-01
## 70 674467 26374 5.20e-01 4.48e-04 5.19e-01 5.21e-01
## 71 648093 26658 4.99e-01 4.48e-04 4.98e-01 5.00e-01
## 72 621435 26592 4.78e-01 4.48e-04 4.77e-01 4.78e-01
## 73 594843 26673 4.56e-01 4.46e-04 4.55e-01 4.57e-01
## 74 568170 26286 4.35e-01 4.44e-04 4.34e-01 4.36e-01
## 75 541884 26606 4.14e-01 4.41e-04 4.13e-01 4.15e-01
## 76 515278 27915 3.91e-01 4.37e-04 3.90e-01 3.92e-01
## 77 487363 27852 3.69e-01 4.32e-04 3.68e-01 3.70e-01
## 78 459511 28635 3.46e-01 4.26e-04 3.45e-01 3.47e-01
## 79 430876 29424 3.22e-01 4.19e-04 3.21e-01 3.23e-01
## 80 401452 29797 2.98e-01 4.10e-04 2.98e-01 2.99e-01
## 81 371655 30938 2.74e-01 3.99e-04 2.73e-01 2.74e-01
## 82 340717 30967 2.49e-01 3.87e-04 2.48e-01 2.49e-01
## 83 309750 31143 2.24e-01 3.73e-04 2.23e-01 2.24e-01
## 84 278607 30570 1.99e-01 3.58e-04 1.98e-01 2.00e-01
## 85 248037 30134 1.75e-01 3.40e-04 1.74e-01 1.76e-01
## 86 217903 29142 1.52e-01 3.21e-04 1.51e-01 1.52e-01
## 87 188761 28056 1.29e-01 3.00e-04 1.28e-01 1.30e-01
## 88 160705 26790 1.08e-01 2.78e-04 1.07e-01 1.08e-01
## 89 133915 24919 8.75e-02 2.53e-04 8.70e-02 8.80e-02
## 90 108996 22216 6.97e-02 2.28e-04 6.92e-02 7.01e-02
## 91 86780 19606 5.39e-02 2.02e-04 5.35e-02 5.43e-02
## 92 67174 16349 4.08e-02 1.77e-04 4.05e-02 4.11e-02
## 93 50825 13673 2.98e-02 1.52e-04 2.95e-02 3.01e-02
## 94 37152 10317 2.15e-02 1.30e-04 2.13e-02 2.18e-02
## 95 26835 7989 1.51e-02 1.09e-04 1.49e-02 1.53e-02
## 96 18846 6033 1.03e-02 9.04e-05 1.01e-02 1.05e-02
## 97 12813 4384 6.77e-03 7.35e-05 6.62e-03 6.91e-03
## 98 8429 2692 4.61e-03 6.07e-05 4.49e-03 4.73e-03
## 99 5737 2125 2.90e-03 4.82e-05 2.81e-03 3.00e-03
## 100 3612 1471 1.72e-03 3.71e-05 1.65e-03 1.79e-03
## 101 2141 825 1.06e-03 2.91e-05 1.00e-03 1.12e-03
## 102 1316 616 5.62e-04 2.12e-05 5.22e-04 6.05e-04
## 103 700 312 3.11e-04 1.58e-05 2.82e-04 3.44e-04
## 104 388 151 1.90e-04 1.24e-05 1.68e-04 2.16e-04
## 105 237 109 1.03e-04 9.08e-06 8.64e-05 1.22e-04
## 106 128 48 6.42e-05 7.18e-06 5.16e-05 8.00e-05
## 107 80 37 3.45e-05 5.26e-06 2.56e-05 4.65e-05
## 108 43 12 2.49e-05 4.47e-06 1.75e-05 3.54e-05
## 109 31 4 2.17e-05 4.17e-06 1.49e-05 3.16e-05
## 110 27 6 1.69e-05 3.68e-06 1.10e-05 2.59e-05
## 111 21 9 9.63e-06 2.78e-06 5.47e-06 1.70e-05
## 112 12 9 2.41e-06 1.39e-06 7.77e-07 7.47e-06
## 116 3 3 0.00e+00 NaN NA NA
##
## SEX=2, Injury=1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 1 74478 594 0.992024 3.26e-04 9.91e-01 0.992664
## 2 73884 464 0.985794 4.34e-04 9.85e-01 0.986645
## 3 73420 334 0.981310 4.96e-04 9.80e-01 0.982283
## 4 73086 279 0.977564 5.43e-04 9.77e-01 0.978628
## 5 72807 248 0.974234 5.81e-04 9.73e-01 0.975373
## 6 72559 229 0.971159 6.13e-04 9.70e-01 0.972362
## 7 72330 143 0.969239 6.33e-04 9.68e-01 0.970480
## 8 72187 136 0.967413 6.51e-04 9.66e-01 0.968689
## 9 72051 148 0.965426 6.69e-04 9.64e-01 0.966739
## 10 71903 105 0.964016 6.82e-04 9.63e-01 0.965355
## 11 71798 164 0.961814 7.02e-04 9.60e-01 0.963192
## 12 71634 130 0.960069 7.17e-04 9.59e-01 0.961476
## 13 71504 231 0.956967 7.44e-04 9.56e-01 0.958426
## 14 71273 231 0.953866 7.69e-04 9.52e-01 0.955373
## 15 71042 230 0.950777 7.93e-04 9.49e-01 0.952332
## 16 70812 416 0.945192 8.34e-04 9.44e-01 0.946828
## 17 70396 479 0.938760 8.79e-04 9.37e-01 0.940484
## 18 69917 680 0.929630 9.37e-04 9.28e-01 0.931469
## 19 69237 779 0.919171 9.99e-04 9.17e-01 0.921130
## 20 68458 775 0.908765 1.06e-03 9.07e-01 0.910835
## 21 67683 863 0.897178 1.11e-03 8.95e-01 0.899362
## 22 66820 861 0.885617 1.17e-03 8.83e-01 0.887906
## 23 65959 883 0.873761 1.22e-03 8.71e-01 0.876150
## 24 65076 1045 0.859730 1.27e-03 8.57e-01 0.862228
## 25 64031 985 0.846505 1.32e-03 8.44e-01 0.849098
## 26 63046 955 0.833682 1.36e-03 8.31e-01 0.836361
## 27 62091 874 0.821947 1.40e-03 8.19e-01 0.824699
## 28 61217 860 0.810400 1.44e-03 8.08e-01 0.813220
## 29 60357 888 0.798477 1.47e-03 7.96e-01 0.801364
## 30 59469 817 0.787508 1.50e-03 7.85e-01 0.790451
## 31 58652 900 0.775424 1.53e-03 7.72e-01 0.778426
## 32 57752 877 0.763648 1.56e-03 7.61e-01 0.766706
## 33 56875 924 0.751242 1.58e-03 7.48e-01 0.754353
## 34 55951 792 0.740608 1.61e-03 7.37e-01 0.743762
## 35 55159 1017 0.726953 1.63e-03 7.24e-01 0.730160
## 36 54142 865 0.715339 1.65e-03 7.12e-01 0.718587
## 37 53277 743 0.705363 1.67e-03 7.02e-01 0.708644
## 38 52534 887 0.693453 1.69e-03 6.90e-01 0.696772
## 39 51647 812 0.682551 1.71e-03 6.79e-01 0.685902
## 40 50835 775 0.672145 1.72e-03 6.69e-01 0.675525
## 41 50060 868 0.660490 1.74e-03 6.57e-01 0.663900
## 42 49192 878 0.648702 1.75e-03 6.45e-01 0.652139
## 43 48314 1028 0.634899 1.76e-03 6.31e-01 0.638366
## 44 47286 966 0.621929 1.78e-03 6.18e-01 0.625421
## 45 46320 956 0.609093 1.79e-03 6.06e-01 0.612607
## 46 45364 967 0.596109 1.80e-03 5.93e-01 0.599643
## 47 44397 1285 0.578856 1.81e-03 5.75e-01 0.582412
## 48 43112 1123 0.563777 1.82e-03 5.60e-01 0.567350
## 49 41989 1271 0.546712 1.82e-03 5.43e-01 0.550299
## 50 40718 1354 0.528532 1.83e-03 5.25e-01 0.532129
## 51 39364 1381 0.509990 1.83e-03 5.06e-01 0.513592
## 52 37983 1223 0.493569 1.83e-03 4.90e-01 0.497172
## 53 36760 1227 0.477094 1.83e-03 4.74e-01 0.480695
## 54 35533 1258 0.460203 1.83e-03 4.57e-01 0.463796
## 55 34275 1183 0.444319 1.82e-03 4.41e-01 0.447902
## 56 33092 999 0.430906 1.81e-03 4.27e-01 0.434477
## 57 32093 945 0.418217 1.81e-03 4.15e-01 0.421775
## 58 31148 991 0.404912 1.80e-03 4.01e-01 0.408452
## 59 30157 924 0.392505 1.79e-03 3.89e-01 0.396028
## 60 29233 723 0.382798 1.78e-03 3.79e-01 0.386304
## 61 28510 855 0.371318 1.77e-03 3.68e-01 0.374804
## 62 27655 814 0.360388 1.76e-03 3.57e-01 0.363853
## 63 26841 677 0.351298 1.75e-03 3.48e-01 0.354744
## 64 26164 710 0.341765 1.74e-03 3.38e-01 0.345189
## 65 25454 650 0.333038 1.73e-03 3.30e-01 0.336440
## 66 24804 616 0.324767 1.72e-03 3.21e-01 0.328148
## 67 24188 541 0.317503 1.71e-03 3.14e-01 0.320864
## 68 23647 577 0.309756 1.69e-03 3.06e-01 0.313095
## 69 23070 537 0.302546 1.68e-03 2.99e-01 0.305863
## 70 22533 484 0.296047 1.67e-03 2.93e-01 0.299344
## 71 22049 543 0.288756 1.66e-03 2.86e-01 0.292030
## 72 21506 723 0.279049 1.64e-03 2.76e-01 0.282289
## 73 20783 518 0.272094 1.63e-03 2.69e-01 0.275309
## 74 20265 490 0.265515 1.62e-03 2.62e-01 0.268705
## 75 19775 573 0.257821 1.60e-03 2.55e-01 0.260982
## 76 19202 575 0.250101 1.59e-03 2.47e-01 0.253230
## 77 18627 556 0.242635 1.57e-03 2.40e-01 0.245734
## 78 18071 749 0.232579 1.55e-03 2.30e-01 0.235633
## 79 17322 772 0.222213 1.52e-03 2.19e-01 0.225219
## 80 16550 850 0.210801 1.49e-03 2.08e-01 0.213750
## 81 15700 822 0.199764 1.47e-03 1.97e-01 0.202656
## 82 14878 1018 0.186095 1.43e-03 1.83e-01 0.188911
## 83 13860 830 0.174951 1.39e-03 1.72e-01 0.177701
## 84 13030 938 0.162357 1.35e-03 1.60e-01 0.165027
## 85 12092 1297 0.144942 1.29e-03 1.42e-01 0.147493
## 86 10795 1165 0.129300 1.23e-03 1.27e-01 0.131732
## 87 9630 1127 0.114168 1.17e-03 1.12e-01 0.116475
## 88 8503 1019 0.100486 1.10e-03 9.83e-02 0.102669
## 89 7484 1189 0.084522 1.02e-03 8.25e-02 0.086543
## 90 6295 1192 0.068517 9.26e-04 6.67e-02 0.070355
## 91 5103 852 0.057077 8.50e-04 5.54e-02 0.058768
## 92 4251 1052 0.042952 7.43e-04 4.15e-02 0.044433
## 93 3199 830 0.031808 6.43e-04 3.06e-02 0.033094
## 94 2369 496 0.025148 5.74e-04 2.40e-02 0.026298
## 95 1873 405 0.019711 5.09e-04 1.87e-02 0.020735
## 96 1468 388 0.014501 4.38e-04 1.37e-02 0.015385
## 97 1080 340 0.009936 3.63e-04 9.25e-03 0.010674
## 98 740 273 0.006270 2.89e-04 5.73e-03 0.006864
## 99 467 127 0.004565 2.47e-04 4.11e-03 0.005076
## 100 340 135 0.002752 1.92e-04 2.40e-03 0.003156
## 101 205 78 0.001705 1.51e-04 1.43e-03 0.002029
## 102 127 26 0.001356 1.35e-04 1.12e-03 0.001648
## 103 101 53 0.000644 9.30e-05 4.86e-04 0.000855
## 104 48 31 0.000228 5.54e-05 1.42e-04 0.000367
## 105 17 5 0.000161 4.65e-05 9.15e-05 0.000284
## 106 12 3 0.000121 4.03e-05 6.29e-05 0.000232
## 110 9 9 0.000000 NaN NA NA
##
## SEX=2, Injury=2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 1 1262690 4781 9.96e-01 5.47e-05 9.96e-01 9.96e-01
## 2 1257909 3232 9.94e-01 7.07e-05 9.94e-01 9.94e-01
## 3 1254677 2295 9.92e-01 8.01e-05 9.92e-01 9.92e-01
## 4 1252382 1403 9.91e-01 8.53e-05 9.91e-01 9.91e-01
## 5 1250979 1361 9.90e-01 9.01e-05 9.89e-01 9.90e-01
## 6 1249618 1071 9.89e-01 9.37e-05 9.89e-01 9.89e-01
## 7 1248547 669 9.88e-01 9.58e-05 9.88e-01 9.88e-01
## 8 1247878 736 9.88e-01 9.81e-05 9.87e-01 9.88e-01
## 9 1247142 637 9.87e-01 1.00e-04 9.87e-01 9.87e-01
## 10 1246505 707 9.87e-01 1.02e-04 9.86e-01 9.87e-01
## 11 1245798 483 9.86e-01 1.04e-04 9.86e-01 9.86e-01
## 12 1245315 480 9.86e-01 1.05e-04 9.86e-01 9.86e-01
## 13 1244835 505 9.85e-01 1.07e-04 9.85e-01 9.86e-01
## 14 1244330 401 9.85e-01 1.08e-04 9.85e-01 9.85e-01
## 15 1243929 609 9.85e-01 1.09e-04 9.84e-01 9.85e-01
## 16 1243320 483 9.84e-01 1.11e-04 9.84e-01 9.84e-01
## 17 1242837 536 9.84e-01 1.12e-04 9.84e-01 9.84e-01
## 18 1242301 530 9.83e-01 1.14e-04 9.83e-01 9.84e-01
## 19 1241771 527 9.83e-01 1.15e-04 9.83e-01 9.83e-01
## 20 1241244 749 9.82e-01 1.17e-04 9.82e-01 9.83e-01
## 21 1240495 712 9.82e-01 1.19e-04 9.82e-01 9.82e-01
## 22 1239783 763 9.81e-01 1.21e-04 9.81e-01 9.81e-01
## 23 1239020 908 9.81e-01 1.23e-04 9.80e-01 9.81e-01
## 24 1238112 784 9.80e-01 1.25e-04 9.80e-01 9.80e-01
## 25 1237328 807 9.79e-01 1.27e-04 9.79e-01 9.80e-01
## 26 1236521 985 9.78e-01 1.29e-04 9.78e-01 9.79e-01
## 27 1235536 1209 9.78e-01 1.32e-04 9.77e-01 9.78e-01
## 28 1234327 1068 9.77e-01 1.34e-04 9.76e-01 9.77e-01
## 29 1233259 1166 9.76e-01 1.37e-04 9.76e-01 9.76e-01
## 30 1232093 1409 9.75e-01 1.40e-04 9.74e-01 9.75e-01
## 31 1230684 1581 9.73e-01 1.43e-04 9.73e-01 9.74e-01
## 32 1229103 1763 9.72e-01 1.47e-04 9.72e-01 9.72e-01
## 33 1227340 1834 9.71e-01 1.50e-04 9.70e-01 9.71e-01
## 34 1225506 1991 9.69e-01 1.54e-04 9.69e-01 9.69e-01
## 35 1223515 1907 9.67e-01 1.58e-04 9.67e-01 9.68e-01
## 36 1221608 2040 9.66e-01 1.62e-04 9.66e-01 9.66e-01
## 37 1219568 2534 9.64e-01 1.66e-04 9.64e-01 9.64e-01
## 38 1217034 2542 9.62e-01 1.71e-04 9.61e-01 9.62e-01
## 39 1214492 3094 9.59e-01 1.76e-04 9.59e-01 9.60e-01
## 40 1211398 3054 9.57e-01 1.81e-04 9.57e-01 9.57e-01
## 41 1208344 3446 9.54e-01 1.86e-04 9.54e-01 9.55e-01
## 42 1204898 3815 9.51e-01 1.92e-04 9.51e-01 9.52e-01
## 43 1201083 4352 9.48e-01 1.98e-04 9.47e-01 9.48e-01
## 44 1196731 4341 9.44e-01 2.04e-04 9.44e-01 9.45e-01
## 45 1192390 4907 9.40e-01 2.11e-04 9.40e-01 9.41e-01
## 46 1187483 5504 9.36e-01 2.18e-04 9.36e-01 9.37e-01
## 47 1181979 6346 9.31e-01 2.25e-04 9.31e-01 9.31e-01
## 48 1175633 6926 9.26e-01 2.34e-04 9.25e-01 9.26e-01
## 49 1168707 7893 9.19e-01 2.42e-04 9.19e-01 9.20e-01
## 50 1160814 8834 9.12e-01 2.52e-04 9.12e-01 9.13e-01
## 51 1151980 9790 9.05e-01 2.61e-04 9.04e-01 9.05e-01
## 52 1142190 9952 8.97e-01 2.71e-04 8.96e-01 8.97e-01
## 53 1132238 10688 8.88e-01 2.80e-04 8.88e-01 8.89e-01
## 54 1121550 11639 8.79e-01 2.90e-04 8.78e-01 8.80e-01
## 55 1109911 12262 8.69e-01 3.00e-04 8.69e-01 8.70e-01
## 56 1097649 13658 8.58e-01 3.10e-04 8.58e-01 8.59e-01
## 57 1083991 14225 8.47e-01 3.20e-04 8.47e-01 8.48e-01
## 58 1069766 14812 8.35e-01 3.30e-04 8.35e-01 8.36e-01
## 59 1054954 15345 8.23e-01 3.39e-04 8.23e-01 8.24e-01
## 60 1039609 15094 8.11e-01 3.48e-04 8.11e-01 8.12e-01
## 61 1024515 15872 7.99e-01 3.57e-04 7.98e-01 8.00e-01
## 62 1008643 17076 7.85e-01 3.65e-04 7.85e-01 7.86e-01
## 63 991567 18996 7.70e-01 3.74e-04 7.70e-01 7.71e-01
## 64 972571 18958 7.55e-01 3.83e-04 7.54e-01 7.56e-01
## 65 953613 19124 7.40e-01 3.90e-04 7.39e-01 7.41e-01
## 66 934489 19490 7.25e-01 3.98e-04 7.24e-01 7.25e-01
## 67 914999 20226 7.09e-01 4.04e-04 7.08e-01 7.09e-01
## 68 894773 19833 6.93e-01 4.11e-04 6.92e-01 6.94e-01
## 69 874940 20170 6.77e-01 4.16e-04 6.76e-01 6.78e-01
## 70 854770 21119 6.60e-01 4.21e-04 6.59e-01 6.61e-01
## 71 833651 21467 6.43e-01 4.26e-04 6.42e-01 6.44e-01
## 72 812184 21671 6.26e-01 4.31e-04 6.25e-01 6.27e-01
## 73 790513 22008 6.09e-01 4.34e-04 6.08e-01 6.09e-01
## 74 768505 23092 5.90e-01 4.38e-04 5.89e-01 5.91e-01
## 75 745413 24493 5.71e-01 4.40e-04 5.70e-01 5.72e-01
## 76 720920 25947 5.50e-01 4.43e-04 5.50e-01 5.51e-01
## 77 694973 25948 5.30e-01 4.44e-04 5.29e-01 5.31e-01
## 78 669025 28408 5.07e-01 4.45e-04 5.06e-01 5.08e-01
## 79 640617 30517 4.83e-01 4.45e-04 4.82e-01 4.84e-01
## 80 610100 31336 4.58e-01 4.43e-04 4.57e-01 4.59e-01
## 81 578764 33439 4.32e-01 4.41e-04 4.31e-01 4.33e-01
## 82 545325 35982 4.03e-01 4.37e-04 4.03e-01 4.04e-01
## 83 509343 35938 3.75e-01 4.31e-04 3.74e-01 3.76e-01
## 84 473405 38359 3.45e-01 4.23e-04 3.44e-01 3.45e-01
## 85 435046 40268 3.13e-01 4.13e-04 3.12e-01 3.13e-01
## 86 394778 41021 2.80e-01 4.00e-04 2.79e-01 2.81e-01
## 87 353757 41462 2.47e-01 3.84e-04 2.47e-01 2.48e-01
## 88 312295 41416 2.15e-01 3.65e-04 2.14e-01 2.15e-01
## 89 270879 39602 1.83e-01 3.44e-04 1.82e-01 1.84e-01
## 90 231277 38103 1.53e-01 3.20e-04 1.52e-01 1.54e-01
## 91 193174 34909 1.25e-01 2.95e-04 1.25e-01 1.26e-01
## 92 158265 30684 1.01e-01 2.68e-04 1.01e-01 1.02e-01
## 93 127581 26889 7.97e-02 2.41e-04 7.93e-02 8.02e-02
## 94 100692 22270 6.21e-02 2.15e-04 6.17e-02 6.25e-02
## 95 78422 18957 4.71e-02 1.89e-04 4.67e-02 4.75e-02
## 96 59465 15796 3.46e-02 1.63e-04 3.43e-02 3.49e-02
## 97 43669 12372 2.48e-02 1.38e-04 2.45e-02 2.51e-02
## 98 31297 9194 1.75e-02 1.17e-04 1.73e-02 1.77e-02
## 99 22103 6720 1.22e-02 9.76e-05 1.20e-02 1.24e-02
## 100 15383 5390 7.91e-03 7.89e-05 7.76e-03 8.07e-03
## 101 9993 3550 5.10e-03 6.34e-05 4.98e-03 5.23e-03
## 102 6443 2268 3.31e-03 5.11e-05 3.21e-03 3.41e-03
## 103 4175 1623 2.02e-03 4.00e-05 1.94e-03 2.10e-03
## 104 2552 1007 1.22e-03 3.11e-05 1.16e-03 1.29e-03
## 105 1545 716 6.57e-04 2.28e-05 6.13e-04 7.03e-04
## 106 829 330 3.95e-04 1.77e-05 3.62e-04 4.31e-04
## 107 499 250 1.97e-04 1.25e-05 1.74e-04 2.23e-04
## 108 249 142 8.47e-05 8.19e-06 7.01e-05 1.02e-04
## 109 107 72 2.77e-05 4.69e-06 1.99e-05 3.86e-05
## 110 35 16 1.50e-05 3.45e-06 9.60e-06 2.36e-05
## 111 19 3 1.27e-05 3.17e-06 7.76e-06 2.07e-05
## 112 16 6 7.92e-06 2.50e-06 4.26e-06 1.47e-05
## 113 10 9 7.92e-07 7.92e-07 1.12e-07 5.62e-06
## 114 1 1 0.00e+00 NaN NA NAPlotting:
ggsurvplot(fit1)
ggsurvplot(fit2)
The Grouping Variable are SEX (1 = male, 2 = female) and Injury (death result of injury: 1 = no, 2 = yes)
My hypothesis is that females have a higher survival probability than males.
The hypothesis seemed to be supported.