как я могу получить график соотношения рисков, используя rstpm2 и ggplot2

Question

• 1000 1007 *

haz <- predict(fit, **type="haz"**, newdata= data.frame(group=0:1), grid=TRUE,var="group", se.fit=TRUE,keep=TRUE)


haz <- cbind(haz,attr(haz,"newdata"))

haz <- transform(haz,group=factor(group,labels=c("Drug","Surgery")))

ggplot(haz, aes(duration,group=group)) + geom_line(aes(y=Estimate,colour=group)) + 
  geom_ribbon(aes(ymin=lower,ymax=upper), alpha=0.2) + xlab("Time since surgery (years)") + ylab("Hazard")

Мне удалось получить хороший сюжет.

Однако, когда я изменил «haz» на «hr», это не сработало. даже оценки не соответствовали. Во всех случаях имеют одинаковые оценки

, например:

Оценка нижнего верхнего

1 0,7664429 0,5020507 1,170071

2 0,7664429 0,5020507 1,170071

3 0,7664429 0,5020507 1,170071

4 0,7664429 0,5020507 1,170071

5 0,7664429 0,5020507 1,170071

6 0,7664429 0,5020507 1,170071

* 1030700 0,5020507 1,170071 *1029*

7

8 0,7664429 0,5020507 1,170071

9 0,7664429 0,5020507 1,170071

10 0,7664429 0,5020507 1,170071

11 0,7664429 0,5020507 1,170071

10407 1,170071 *1039* *1040* 1040

13 0,7664429 0,5020507 1,170071

14 0,7664429 0,5020507 1,170071

Что мне делать ???

помогите пожалуйста!