In [`cFit`](https://github.com/cdriveraus/ctsem/blob/master/R/ctFit.R#L207) it looks like you are hardcoding ```R idcol='id' obsTpoints=max(unlist(lapply(unique(dat[,idcol]),function(x) length(dat[dat[,idcol]==x, idcol]) ))) ``` But the user specifies that in the model function so I think you can just make it ```R idcol <- ctmodelobj$id ```