search

mclements / rstpm2

An R package for generalised survival models
28 stars 11 forks source link

fix invalid check #13

Closed boennecd closed 4 years ago

boennecd commented 4 years ago

Before the new test gave 

dat <- data.frame(
    y     = c(1, 2, 1, 2, 1, 1), 
    event = c(0, 1, 1, 0, 1, 1),
    grp   = c(1, 1, 2, 2, 3, 3))

fit     <- gsm(Surv(y, event) ~ 1, dat, df = 1, cluster = dat$grp)
fit_new <- gsm(Surv(y, event) ~ 1, dat, df = 1, cluster = dat$grp, 
               init = head(coef(fit), -1), logtheta = tail(coef(fit), 1))
#R> Error in gsm(Surv(y, event) ~ 1, dat, df = 1, cluster = dat$grp, init = head(coef(fit),  : 
#R>   length(init) + frailty == length(coef(lm.obj)) is not TRUE

but the stopifnot check was not correct.

mclements commented 4 years ago

Thanks for this.

I had already committed a change for the stopifnot(), so I removed that part.