我想在我的生存分析中应用HL GoF测试。
可以找到测试说明here:
根据一些现成的数据,我们可以用95%CI
绘制cdflibrary(survival)
library(ResourceSelection)
data(ovarian)
s <- Surv(ovarian$futime, ovarian$fustat)
sWei <- survreg(s ~ age,dist='weibull',data=ovarian)
probs <- seq(0.001,0.999,0.001)
pred <- predict(sWei,type="uquantile",p=probs,se.fit=TRUE)
val <- cbind(pred$fit[1,],pred$fit[1,]-1.96*pred$se.fit[1,],pred$fit[1,]+1.96*pred$se.fit[1,])
val <- exp(val)
# cumulative density function
plot(val[,1], probs)
lines(val[,2], probs, col='red')
lines(val[,3], probs, col='red')
但是我们想要运行HL GoF测试来评估Weibull分布是否不正确。
因此:
hoslem(x,y)
x = vector of observations
y = fitted probabilities
我试试,但不正确。有人可以引导我找到正确的拟合概率,它们需要与观察的长度相同
hl <- hoslem.test(ovarian$fustat, val[,1], g=10)