library(survival) ## Read in the data: ## if( TRUE ){ DF = read.table('http://lib.stat.cmu.edu/datasets/csb/ch7.dat', header=FALSE) colnames(DF) = c('ID', 'Time', 'Cens', 'Trawl_Time', 'Max_Minus_Min', 'Fork_Length', 'Handling_Time', 'LN_Total_Catch') DF$ID = NULL } res = coxph(Surv(Time, Cens) ~ ., data=DF) print(res)