从Survreg解释Weibull参数

作者: 狗头军师
发布时间: 2024-12-18 03:55:13 (4天前)

2 条回复

1#
回复此人
v-star*위위 | 2020-08-24 15-35

手册页?survreg（在“示例”部分中）对此进行了说明。 ``` library(survival) y <- rweibull(1000, shape=2, scale=5) r <- survreg(Surv(y)~1, dist="weibull") a <- 1/r$scale # Approximately 2 b <- exp( coef(r) ) # Approximately 5 y2 <- b * ( -ln( 1-runif(1000) ) ) ^(1/a) y3 <- rweibull(1000, shape=a, scale=5) # Check graphically that the distributions are the same plot(sort(y), sort(y2)) abline(0,1) ```

编辑

登录后才能参与评论