omuhaz pacote estima ofunção de risco decensurado à direita dados usando métodos de suavização do kernel. Minha pergunta é: existe alguma maneira de obter intervalos de confiança para a função de risco quemuhaz calcula?

options(scipen=999)
library(muhaz)
data(ovarian, package="survival")
attach(ovarian)
fit1 <- muhaz(futime, fustat)
plot(fit1, lwd=3, ylim=c(0,0.002))

No exemplo acima, omuhaz.object fit tem algumas entradasfit1$msemin, fit1$var.min, fit1$haz.est no entanto, seu comprimento é metade dafit1$haz.est.

Alguma idéia se é possível extrair intervalos de confiança para a função de risco?

EDIT: Eu tentei o seguinte com base no que @ user20650 sugeriu

options(scipen=999)
library(muhaz)
data(ovarian, package="survival")
fit1 <- muhaz(ovarian$futime, ovarian$fustat,min.time=0, max.time=744)


h.df<-data.frame(est=fit1$est.grid, h.orig=fit1$haz.est)

for (i in 1:10000){
d.s.onarian<-ovarian[sample(1:nrow(ovarian), nrow(ovarian), replace = T),]
d.s.muhaz<-muhaz(d.s.onarian$futime, d.s.onarian$fustat, min.time=0, max.time=744 )
h.df<-cbind(h.df, d.s.muhaz$haz.est)
}


h.df$upper.ci<-apply(h.df[,c(-1,-2)], 1,  FUN=function(x) quantile(x, probs = 0.975))
h.df$lower.ci<-apply(h.df[,c(-1,-2)], 1,  FUN=function(x) quantile(x, probs = 0.025))
plot(h.df$est, h.df$h.orig, type="l", ylim=c(0,0.003), lwd=3)
lines(h.df$est, h.df$upper.ci,  lty=3, lwd=3)
lines(h.df$est, h.df$lower.ci,  lty=3, lwd=3)

Definir o tempo máximo parece funcionar, toda amostra de bootstrap possui os mesmos pontos de grade de estimativa. No entanto, o IC obtido faz pouco sentido. Normalmente, eu esperaria que os intervalos fossem estreitos em t = 0 e aumentassem com o tempo (menos informações, mais incerteza), mas os intervalos obtidos parecem ser mais ou menos constantes com o tempo.