Hello,
I've got problems with censboot in the boot library based on Davison&Hinkley
(I'm running S-Plus Prof 4.5 rel 2 in Windows 98 and Bootstrap S-Plus Functions
(Version 1.1; May 1998)).
When the examples in boot help are run I get the following
> aml.fun <- function(data) {
+ surv <- survfit(Surv(data$time,data$cens)~data$group)
+ out <- NULL
+ st <- 1
+ for (s in 1:length(surv$strata)) {
+ inds <- st:(st+surv$strata[s]-1)
+ md <- min(surv$time[inds[1-surv$surv[inds]>=0.5]])
+ st <- st+surv$strata[s]
+ out <- c(out,md)
+ }
+ }
> aml.case <- censboot(aml,aml.fun,R=499,strata=aml$group)
Error in t[r, ] <- statistic(bootda..: Invalid righ
t-hand side to replacement
Dumped
and also
> aml.s1 <- survfit(Surv(aml$time,aml$cens)~aml$group)
>
> aml.s2 <- survfit(Surv(aml$time-0.001*aml$cens,1-aml$cens)~1)
> aml.cond <- censboot(aml,aml.fun,R=499,strata=aml$group,
+ F.surv=aml.s1,G.surv=aml.s2,sim="cond")
Error in t[r, ] <- statistic(bootda..: Invalid righ
t-hand side to replacement
Dumped
and as well
> aml.fun1 <- function(data,str) {
+ surv <- survfit(Surv(data[,1],data[,2])~str)
+ out <- NULL
+ st <- 1
+ for (s in 1:length(surv$strata)) {
+ inds <- st:(st+surv$strata[s]-1)
+ md <- min(surv$time[inds[1-surv$surv[inds]>=0.5]])
+ st <- st+surv$strata[s]
+ out <- c(out,md)
+ }
+ }
> aml.wei <- censboot(cbind(aml$time,aml$cens),aml.fun1,R=499,
+ strata=aml$group, F.surv=aml.s1,sim="weird")
Error in rbind(temp, as.matrix(data[..: Number of co
lumns of matrices and lengths of names vectors must
match (see arg 2)
Dumped
The melanoma boot works OK but i have further problems with the envelope
> mel.env <- envelope(mel.mod$t)$point
Error in envelope(mel.mod$t): Object "n" not found
Dumped
Has anyone some help to offer??
With thanks in advance
Fredrik Lundgren
Cardiovascular centre
Linköping
Sweden
-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news
|