In vaccine research, proportional hazards model including only first event have been widely used for estimating vaccine efficacy because it is easy to interpret and convenient. However, this method causes not only loss of information but also biased result when heterogeneity of study subject in exposure and susceptibility exists. Furthermore, it is hard to ignore the possibility that each event is correlated with each other in the repeated events. Therefore, we compare various statistical models to estimate vaccine efficacy under various situations with heterogeneity and event dependency.