Crash R with larger number of subjects than requested

[mattfidler]

library(nlmixr2)
# the model 
PK_1cmt <- function() {
  ini({
    lka    <- 0.45 
    lcl    <- 1 
    lvc    <- 3.45 
    eta.cl ~ 0.3
    propSd <- 0.5 
  })
  model({
    ka <- exp(lka)
    cl <- exp(lcl + eta.cl)
    vc <- exp(lvc)
    linCmt() ~ prop(propSd)
  })
}
# fit the model
fit <- nlmixr(PK_1cmt, theo_sd, est="focei")

# simulate the model (more IDs in param then in events – will crash R)
parms <- data.frame(ID=1:12,lka=0.45,lcl=1,lvc=3.45,eta.cl=0.3)
mdsim <- rxSolve(fit$ui,params=parms,events = subset(theo_sd,ID%in%1:6))

@RichardHooijmaijers