R version 3.1.0 (2014-04-10) -- "Spring Dance"
Copyright (C) 2014 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin13.1.0 (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(survey)
Attaching package: 'survey'
The following object is masked from 'package:graphics':
dotchart
> library(survival)
Loading required package: splines
>
> pbc2<-rbind(pbc,pbc)
> pbc2$id<-rep(1:418,2)
>
> dpbc1<-svydesign(id=~1, data=pbc)
Warning message:
In svydesign.default(id = ~1, data = pbc) :
No weights or probabilities supplied, assuming equal probability
> dpbc2<-svydesign(id=~id, data=pbc2)
Warning message:
In svydesign.default(id = ~id, data = pbc2) :
No weights or probabilities supplied, assuming equal probability
>
> s1<-svykm(Surv(time,status>0)~1, subset(dpbc1, bili>6), se=TRUE)
> s2<-svykm(Surv(time,status>0)~1, subset(dpbc2, bili>6), se=TRUE)
>
> (c1<-confint(s1,(1:5)*365))
0.025 0.975
365 0.6446215 0.8594153
730 0.5410938 0.7766848
1095 0.2683127 0.5103356
1460 0.1444731 0.3722001
1825 0.1009672 0.3204713
> (c2<-confint(s2,(1:5)*365))
0.025 0.975
365 0.6446215 0.8594153
730 0.5410938 0.7766848
1095 0.2683127 0.5103356
1460 0.1444731 0.3722001
1825 0.1009672 0.3204713
> all.equal(c1, c2)
[1] TRUE
>
> m1<-svycoxph(Surv(time,status>0)~log(bili), design=dpbc1)
> m2<-svycoxph(Surv(time,status>0)~log(bili), design=dpbc2)
>
> d<-data.frame(bili=c(5,10))
> p1<-predict(m1, se=TRUE, newdata=d,type="curve")
> p2<-predict(m2, se=TRUE, newdata=d,type="curve")
>
> (pc1<-confint(p1[[1]],(1:5)*365))
0.025 0.975
365 0.8410027 0.9266263
730 0.7371114 0.8548312
1095 0.5517779 0.7018583
1460 0.4335073 0.5992819
1825 0.3260899 0.5046241
> (pc2<-confint(p2[[1]],(1:5)*365))
0.025 0.975
365 0.8409490 0.9267054
730 0.7370152 0.8549432
1095 0.5515848 0.7019513
1460 0.4332252 0.5992968
1825 0.3257172 0.5045795
> all.equal(pc1, pc2)
[1] "Mean relative difference: 0.0002070722"
>
> (q1<-quantile(p1[[2]]))
0.75 0.5 0.25
489 930 1492
> (q2<-quantile(p2[[2]]))
0.75 0.5 0.25
489 930 1492
> all.equal(q1,q2)
[1] TRUE
>
> proc.time()
user system elapsed
3.410 0.099 3.519
Generated by dwww version 1.15 on Wed Jun 26 01:28:41 CEST 2024.