data(tyge_2011)
# compute log relative risks
TE <- log((tyge_2011$event_e / tyge_2011$n_e) /
(tyge_2011$event_c / tyge_2011$n_c))
seTE <- sqrt(
1 / tyge_2011$event_e - 1 / tyge_2011$n_e +
1 / tyge_2011$event_c - 1 / tyge_2011$n_c
)
head(data.frame(TE, seTE))
Run the code above in your browser using DataLab