Fix bug from Richard Jackson

Question

Fix bug from Richard Jackson

Closed this issue 8 years ago · 2 comments

graemeleehickey commented 8 years ago

Answer 1 · 2017-01-04T17:11:12.000Z

See email from RJ

Answer 2 · 2017-01-09T11:55:21.000Z

Example:

data(heart.valve)
heart.surv <- UniqueVariables(heart.valve, 
                              var.col = c("fuyrs", "status"),
                              id.col = "num")
heart.surv <- heart.surv[, 3:1] # reverse order
heart.long <- heart.valve[, c("num", "time", "log.lvmi")]
heart.cov <- UniqueVariables(heart.valve, 
                             c("age", "hs", "sex"), 
                             id.col = "num")
heart.valve.jd <- jointdata(longitudinal = heart.long, 
                            baseline = heart.cov, 
                            survival = heart.surv, 
                            id.col = "num", 
                            time.col = "time")
fit <- joint(data = heart.valve.jd, 
             long.formula = log.lvmi ~ 1 + time + hs, 
             surv.formula = Surv(fuyrs, status) ~ hs, 
             model = "intslope")
jointSE(fit, nboot = 10) # fails