Based on our results the effect of the 5-HTTLPR might be better captured in specific subgroups of the population than in large but heterogeneous population samples. Our Bayesian Odds Ratio calculations demonstrated that in these subgroups the genetic ORs are increased. Such stratified analyses should result in greater power in future GxE interaction studies (see Table A in S1 File). In addition, our study demonstrated that the Bayesian network-based methodology, which was developed to analyze relevance of predictors with respect to a set of phenotypic, clinical and environmental descriptors, is a powerful approach to identify highly relevant genetic risk factors even when traditional (frequentist) analysis methods are not able to provide significant results.
