search

graemeleehickey / joineRML

R package for fitting joint models to time-to-event data and multivariate longitudinal data
GNU General Public License v3.0
30 stars 11 forks source link

Cannot fit univariate random-intercepts only JMs #27

Closed graemeleehickey closed 7 years ago

graemeleehickey commented 7 years ago

Matrix multiplication error for univariate JM with random-intercepts only

pbc2$log.bilirubin <- log(pbc2$serBilir)

fit.pbc <- mjoint(
  formLongFixed = list("log.bil" = log.bilirubin ~ year),
  formLongRandom = list("log.bil" = ~ 1 | id),
  formSurv = Surv(years, status2) ~ drug + age,
  data = pbc2,
  timeVar = "year",
  control = list(earlyPhase = 50, convCrit = "sas", rav = 0.01),
  verbose = TRUE)