Nothing
library(testthat)
library(data.table)
context("Prediction error")
# {{{ "Brier score censored data order"
test_that("Brier score censored data order",{
library(riskRegression)
library(survival)
library(prodlim)
data(Melanoma)
setDT(Melanoma)
fit <- coxph(Surv(time,status!=0)~invasion+epicel+logthick,data=Melanoma,x=TRUE)
## many ties in Melanoma
setkey(Melanoma,age)
a <- Score(list(fit),data=Melanoma,Surv(time,status!=0)~invasion+epicel+logthick,cens.model="marginal",metric="Brier")
A <- Score(list(fit),data=Melanoma,Surv(time,status!=0)~invasion+epicel+logthick,cens.model="cox",metric="Brier")
setkey(Melanoma,logthick)
b <- Score(list(fit),data=Melanoma,Surv(time,status!=0)~invasion+epicel+logthick,cens.model="marginal",metric="Brier")
B <- Score(list(fit),data=Melanoma,Surv(time,status!=0)~invasion+epicel+logthick,cens.model="cox",metric="Brier")
a$call <- b$call <- A$call <- B$call <- NULL
## expect_error(expect_equal(ignore_attr=TRUE,a,b,tolerance = .002))
expect_equal(ignore_attr=TRUE,a,b,tolerance = .02)
expect_equal(ignore_attr=TRUE,A,B,tolerance=.02)
})
# }}}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.