defD <- defData(
varname = "effect", formula = 0, variance = 1,
dist = "normal"
)
defE <- defDataAdd(
varname = "died", formula = "-2.5 + 0.3*period + effect",
dist = "binary", link = "logit"
)
d <- genData(1000, defD)
d <- addPeriods(d, 10)
dx <- genNthEvent(d, defEvent = defE, nEvents = 3)
Run the code above in your browser using DataLab