unsa.R

In camsabathe/pseudoICD_bad: Pseudo-observations for interval-censored data from an illness-death model

library(SmoothHazard)

data("Paq1000")

fit.w= idm(formula02=Hist(time=t,event=death)~1,
           formula01=Hist(time=list(l,r),event=dementia)~1,
           formula12=Hist(time=t,event=death)~1,data=Paq1000)

fit.splines <-  idm(formula02=Hist(time=t,event=death)~1,
                    formula01=Hist(time=list(l,r),event=dementia)~1,
                    formula12=Hist(time=t,event=death)~1,
                    method="Splines", kappa=c(200000,200000,20000),
                    data=Paq1000, CV=FALSE)


fit.loo.w=idmloo(object=fit.w, data=Paq1000, parallel = FALSE, cpus = NULL, type = NULL)


fit.loo.w.par=idmloo(object=fit.w, data=Paq1000, parallel = TRUE, cpus = 3, type = NULL)


fit.loo.sp=idmloo(object=fit.splines, data=Paq1000, parallel = FALSE, cpus = NULL, type = NULL)


fit.loo.w.sp=idmloo(object=fit.splines, data=Paq1000, parallel = TRUE, cpus = 3, type = NULL)


test=pseudo(fit.loo.w, t=10,lifeExpect=FALSE, parallel = FALSE, cpus=NULL, type=NULL)

test.par=pseudo(fit.loo.w, t=10,lifeExpect=FALSE, parallel = TRUE, cpus=NULL, type=NULL)


remove(les_fit)

object=fit.w
data=Paq1000
parallel=TRUE
cpus=3
type=NULL


formula01=formula(fit.w$terms$Formula01)

formula02=formula(fit.w$terms$Formula02)
formula12=formula(fit.w$terms$Formula12)

data=Paq1000


m <- match.call()
m01 <- m02 <- m12 <- m[match(c("","data","subset","na.action"),names(m),nomatch=0)]
m01$formula <-
m02$formula <-  formula(fit.w$terms$Formula02)
m12$formula <-  formula(fit.w$terms$Formula12)
m01[[1]] <- m02[[1]] <- m12[[1]] <- as.name("model.frame")
m01 <- eval(m01,parent.frame())
m02 <- eval(m02,parent.frame())
m12 <- eval(m12,parent.frame())
responseTrans <- stats::model.response(m01)
responseAbs <- stats::model.response(m02)
truncated <- nchar(attr(responseAbs,"entry.type"))>1


isdkjfi=formula(fit.w$terms$Formula01)
attributes(isdkjfi)


fit.weib <- idm(formula02=Hist(time=t,event=death,entry=e)~certif,
                formula01=Hist(time=list(l,r),event=dementia)~certif,
                data=Paq1000, conf.int=FALSE)

fit.loo.w=idmloo(object=fit.weib, data=Paq1000, parallel = FALSE, cpus = NULL, type = NULL)

formula02=Hist(time=t,event=death,entry=e)

fit.weib <- idm(formula02=Hist(time=t,event=death,entry=e)~1,
                formula01=Hist(time=list(l,r),event=dementia)~1,
                data=Paq1000, conf.int=FALSE)

fit.loo.w=idmloo(object=fit.weib, data=Paq1000, parallel = FALSE, cpus = NULL, type = NULL)


test1=check_trunc(object=fit.w)