[R] Survival Regression with multiple events per subject

[Fabian Hefner]

Dear R users!

I want to process a maximum likelihood estimation for a parametric
regression survival time model with multiple events per subject.

the STATA command for this survival regression is:

use survreg
stset failure(exercise), id(optionid)
local regressors itm posret negret
streg `regressors', distribution(weibull)

explanation:
stset declares data to be survival-time data;
exercise is the indicator variable, which shows if the subject is dead or
alive;
optionid is the multiple-record ID variable which means every subject has a
unique id and one subject can have multiple events. (see example below)

streg computes a maximum likelihood estimation for parametric regression
survival time models with multiple record data

now I search an equivalent command in R
I found the survival package but I have no solution for the use with
multiple records per subject.

library(survival)
survRegData <- survreg(formula=Surv(time,exercise)~itm+posret+negret,
data=Data, dist="weibull")
summary(survRegData)

My Question is: how can I modify the above command for the use of multiple
events per subject when the optionid is used for indicating the subject?

The dataset look like:

data  | exercisedate | itm    | posret  | negret | optionid| exercise | time
1996  | 1996         | 1.4518 | 0.05487 |-0.4485 |    1    |   0      |   1
1997  | 1997         | 2.4535 | 0.00385 |-0.2525 |    1    |   1      |   5
1998  | 1998         | 1.2523 | 0.04486 |-0.1482 |    2    |   1      |   8
1999  | 1999         | 3.4257 | 0.15287 |-0.8615 |    3    |   0      |   4
2000  | 2000         | 1.1457 | 0.07487 |-0.5485 |    3    |   1      |   5
2001  | 2001         | 2.4418 | 0.09553 |-0.3772 |    3    |   0      |   2


Thank you,

Fabian Hefner