 Descriptive analysis will be calculated using mean and standard deviation for continuous variables, percentage and frequency for categorical variables, and median and quartile range when the continuous variables are not normally distributed. Independent t-tests, Mann-Whitney U test, chi-square test or Fishers exact test will be used to compare any differences in socio-demographic data or outcome variables between groups. Significant results are indicated if p-value of all statistical analysis is less than 0.05 for two-tailed test. To determine the differences or changes between study group and control group between group effects, within group effects time, and the interaction effects group x time, two way repeated measures Analysis of variances ANOVA or co-variances ANCOVA will be applied when assumptions for parametric test are fulfilled. When the data cannot fulfill the assumptions of parametric test, Mann-Whitney U test will be applied to compare the differences between groups at each time interval while Friedman test will be used to determine the time effects within each group. 