#data(data5)
#modM1=lm(M1 ~ Exposure, data = data5)
#modM2=lm(M2 ~ Exposure, data = data5)
#modM3=lm(M3 ~ Exposure, data = data5)
#lmodel.m=list(modM1,modM2, modM3)
#model.y=aalen(Surv(surv_time, event) ~ const(Exposure) + const(M1) + const(M2)
#+ const(M3), data = data5, robust=T)
#multi.media=multimediate(lmodel.m,correlated=TRUE,model.y,treat='Exposure',
#treat.value=1,control.value=0,J=1000,conf.level=0.95,data=data5)
#summary(multi.media, opt='avg')
Run the code above in your browser using DataLab