Hosmer-Lemeshow适合度(GOF)测试

时间:2015-11-05 08:56:39

标签: r survival-analysis goodness-of-fit

我想在我的生存分析中应用HL GoF测试。

可以找到测试说明here

根据一些现成的数据,我们可以用95%CI

绘制cdf
library(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)

0 个答案:

没有答案