Lampiran: Summary (Survfit (Surv (Data$Survt, Data$Relapse 1) 1) )

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 2

LAMPIRAN

> summary(survfit(Surv(Data$Survt,Data$Relapse==1)~1))
Call: survfit(formula = Surv(Data$Survt, Data$Relapse == 1) ~ 1)

time n.risk n.event survival std.err lower 95% CI upper 95% CI


1 42 2 0.952 0.0329 0.8901 1.000
2 40 2 0.905 0.0453 0.8202 0.998
3 38 1 0.881 0.0500 0.7883 0.985
4 37 2 0.833 0.0575 0.7279 0.954
5 35 2 0.786 0.0633 0.6709 0.920
6 33 3 0.714 0.0697 0.5899 0.865
7 29 1 0.690 0.0715 0.5628 0.845
8 28 4 0.591 0.0764 0.4588 0.762
10 23 1 0.565 0.0773 0.4325 0.739
11 21 2 0.512 0.0788 0.3783 0.692
12 18 2 0.455 0.0796 0.3227 0.641
13 16 1 0.426 0.0795 0.2958 0.615
15 15 1 0.398 0.0791 0.2694 0.588
16 14 1 0.369 0.0784 0.2437 0.560
17 13 1 0.341 0.0774 0.2186 0.532
22 9 2 0.265 0.0765 0.1507 0.467
23 7 2 0.189 0.0710 0.0909 0.395

> plot(kmfit2, lty = c('solid', 'dashed'), col=c('blue','orange'),


xlab='survival time (weeks)', ylab='survival probabilities')
> survdiff(Surv(Survt,Relapse)~Rx, data=Data)
Call:
survdiff(formula = Surv(Survt, Relapse) ~ Rx, data = Data)

N Observed Expected (O-E)^2/E (O-E)^2/V


Rx=0 21 9 19.3 5.46 16.8
Rx=1 21 21 10.7 9.77 16.8

Chisq= 16.8 on 1 degrees of freedom, p= 4e-05


> modpar1=survreg(Surv(Data$Survt,Data$Relapse) ~ log.WBC +
Rx,data=Data,dist='exponential')
> modpar1
Call:
survreg(formula = Surv(Data$Survt, Data$Relapse) ~ log.WBC +
Rx, data = Data, dist = "exponential")

Coefficients:
(Intercept) log.WBC Rx
5.9578732 -0.8843946 -1.0930739

Scale fixed at 1

Loglik(model)= -100.8 Loglik(intercept only)= -116.8


Chisq= 31.99 on 2 degrees of freedom, p= 1.13e-07
n= 42

You might also like