library(survival)
# coxdata=read.table("https://folk.ntnu.no/bo/STK4080/cox-hand.txt",header=T)
# fit.c=coxph(Surv(Time,Status==1)~x, data=coxdata)
Time = c(5,10,40,80,120,400,600)
x = c(12,10,3,5,3,4,1)
Status = c(0,1,0,0,1,1,0)
fit.c=coxph(Surv(Time,Status==1)~x)
summary(fit.c)
# SCHOENFELD RESIDUAL PLOT
plot(cox.zph(fit.c), df=2)
#
fit.p=survreg(Surv(Time,Status==1)~x, dist="weibull")
#
# BATTERY DATA
TimeB = c(164,164,218,230,263,467,538,639,669,917,1148,1678,1678,1678,1678,76,82,210,315,385,412,491,504,522,646,678,775,884,1131,1446,1824,1827,2248,2385,3077)
xB = c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
StatusB = c(1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,1,1,1)
fit.coxB=coxph(Surv(TimeB,StatusB==1)~xB)
summary(fit.coxB)
fit.parB=survreg(Surv(TimeB,StatusB==1)~xB, dist="weibull")
summary(fit.parB)
# SCHOENFELD RESIDUAL PLOT
plot(cox.zph(fit.coxB))