This article has Open Peer Review reports available.
An investigation of gene-gene interactions in dose-response studies with Bayesian nonparametrics
- Andrew L Beam^{1}Email author,
- Alison A Motsinger-Reif^{2, 3} and
- Jon Doyle^{4}
https://doi.org/10.1186/s13040-015-0039-3
© Beam et al.; licensee BioMed Central. 2015
Received: 7 July 2014
Accepted: 18 January 2015
Published: 6 February 2015
Abstract
Background
Best practice for statistical methodology in cell-based dose-response studies has yet to be established. We examine the ability of MANOVA to detect trait-associated genetic loci in the presence of gene-gene interactions. We present a novel Bayesian nonparametric method designed to detect such interactions.
Results
MANOVA and the Bayesian nonparametric approach show good ability to detect trait-associated genetic variants under various possible genetic models. It is shown through several sets of analyses that this may be due to marginal effects being present, even if the underlying genetic model does not explicitly contain them.
Conclusions
Understanding how genetic interactions affect drug response continues to be a critical goal. MANOVA and the novel Bayesian framework present a trade-off between computational complexity and model flexibility.
Keywords
Background
Understanding the genetic factors underlying differential drug response continues to be a primary goal in pharmecogenomics and drug development. Association studies based on the use of in vitro cell lines, such as lymphoblastoid cell lines (LCLs), are becoming an attractive alternative to traditional, human based clinical trials [1-3]. Cell lines offer increased sample sizes relative to traditional studies, resulting in higher statistical power to detect genetic variants that potentially drive drug response. However, the challenges offered by these unique data has yet to be comprehensively studied and evaluated.
Recent work has shown that considering the full set dose-responses instead of summary statistics can greatly increase the power to detect trait-associated genetic variants [3-7]. This new approach is based on a statistical method known as multivariate analysis of variance or MANOVA, which is an extension to the well known analysis of variance (ANOVA) framework. Recent studies using a MANOVA-based framework have revealed genetic loci associated with differential responses in anti-cancer agents [3,4]. However, little study to date has been done to investigate the potential effect of gene-gene interactions, or epistasis, which is thought to be a critical piece in the genetic architecture of many complex phenotypes [8-10].
The role of epistasis in complex genetic disease continues to be debated in the literature, as well as the potential of epistasis to solve the so-call missing heritability problem [11-13]. Recent work has validated the role of epistasis in complex disease in humans [14] and drosophila [15]. To date, nearly all of the work on epistasis has focused on case-control studies, or single-valued QTLs (e.g. gene expression). No work has yet been done to understand how epistasis may affect the unique, multiresponse data that comes from dose-response studies.
To investigate how genetic interactions may affect the power of MANOVA we performed a simulation study across several plausible models of SNP-driven dose-response involving multiple loci. In addition to investigating the MANOVA approach, we develop a novel Bayesian nonparametric model that is capable of automatically accounting for genetic interactions in multi-response data. We have previously developed a Bayesian neural network testing framework [16] for case-control studies that is capable of identifying trait-associated loci in the presence of genetic interactions in a computationally tractable way. In this study, we extend this framework to model dose-response data that contain multiple continuous responses. Through comparisons of several possible genetic models, we hope to gain insight as to how deviation from an additive model (i.e. the presence of epistasis) may affect the statistical power of each method.
Methods
Multivariate analysis of variance (MANOVA)
where Max,Min are the upper and lower asymptotes, and w is the hill-slope. The IC_{50} parameter is the concentration at which the response is 50% of maximum and is usually given special importance because it a measure of chemical potency. The IC_{50} is then treated as quantitative trait and quantitative trait locus (QTL) mapping techniques attempt to identify loci that appear to be associated with this trait. In general, this approach results in a drastic loss of power relative to potential alternatives [7,17]. Recent work has shown that multivariate analysis of variance (MANOVA) has high power across a wide variety of possible dose-response relationships, making it an attractive option for genome-wide association mapping.
The least-squares estimator for B_{p×k}=<β_{1},…,β_{ p }> is given by (X^{ T }X)^{−1}X^{ T }Y. where X_{n×p}=<x_{ i },…,x_{ n }>^{ T } and Y_{n×k}=<y_{1},…,y_{ n }>. For the case when n<p or when X is not full rank, as is often the case for GWASs, smaller sub-models may be fit and analyzed. For example, if 1 million markers are genotyped for only 1,000 individuals, each marker may be fit separately along with any possible confounding covariates such as age, gender, or sub-population status.
This model has several advantages including interpretability and a well understood theoretical foundation that allows for null-hypothesis significance testing. However, the linearity assumption may be somewhat restrictive and if a trait is influenced or determined by interactions between markers, this model may only partially capture the true relationship. One approach to lessen this restriction would be to explicitly include all interactions as terms in the model, and build the linear model on this expanded set of covariates. However, this quickly becomes infeasible for even 2^{ n d } order interactions. In a study containing one million markers, there are approximately 5∗10^{9} possible 2^{ n d } order interactions. The computational issues associated with this approach make it infeasible, as do the corresponding multiple testing issues associated with testing billions or trillions of simultaneous hypotheses.
A different approach to account for interactions would be to use a class of models that relax the linearity assumption and impose very little structural constraints on the relationship between x_{ i } and y_{ ik }. Neural networks are one such approach and have a rich history of success in the machine learning and genetic epidemiology [18,19]. Neural nets come with several theoretical guarantees that make them appealing for modeling potentially nonlinear functions. Perhaps the most germane property is given by the universal function approximation theorem [20], which states that a sufficiently complex neural network is capable of modeling any smooth function on a compact set to an arbitrary degree of precision. This theorem ensures that a neural network is capable in principle of modeling a rich class of functions between input and output. Subsumed in the class of functions neural networks can represent is the linear model used by MANOVA. Thus, if the true model is indeed linear, the relationship learned by the neural network would automatically collapse to represent the linear mapping between input and output, while if the true relationship is more complex, it will include any nonlinearities without having to specify them a priori. In the next section we present and develop a Bayesian neural network framework for dose-response studies.
Bayesian neural networks
Bayesian neural networks represent an extension to the familiar neural network framework. Recent work in [16] has shown they can successfully model genetic interactions in case-control studies. Here we present an augmented design capable of modeling multi-response data, such as the kind observed in dose-response studies.
where N(·,·) and IG(·,·) represent Normal and Inverse-Gamma densities, respectively. This prior allows the network to automatically determine which inputs are the most relevant. If SNP j seems to be related to the response, then \({\sigma _{j}^{2}}\) will be concentrated around large values in the posterior. We have previously developed a Bayesian framework that allows for testing of variable importance relative to a null model [16]. This framework provides a posterior probability of each SNP’s involvement in the response and is capable of being applied to continuous, dose-response data in addition to the case-control scenario for which it was originally developed.
This framework comes with several advantages relative to the MANOVA approach. As discussed previously, the relaxation of the linear assumption allows for a broader class of possible genetic architectures to be detected. Additionally, the BNN approach provides a nonparameteric procedure for variable selection while MANOVA relies on normality of errors for its testing procedure to be valid. For instance, p-values from a MANOVA analysis may not be valid for scenarios in which the number of individuals for one genotype is rare (e.g. very few homozygotes for the minor allele). Additionally, normalization procedures must often be performed in order to ensure the error distributions conform as closely as possible to the assumed distribution. In contrast, due the lack of strict assumptions required by the Bayesian neural network framework, there is no reason in principle why a BNN would be affected by any of these issues.
Results and discussion
We investigated three plausible relationships between genotype and drug-response. For all simulations a Bayesian neural network with 10 hidden units was used. The HMC simulation was performed for 375 iterations, with the first 25 discarded as burn-in. For the HMC sampler a step-size value of 0.02, a momentum-persistence value of 0.75, and a temperature value of 1000 were used. The ARD hyper-parameter values were set to α_{0}=3 and β_{0}=1 while the output layer used α_{0}=0.1 and β_{0}=0.1, and we used an ARD cut-off value of 0.4. The ARD hyper-parameters, α_{0}, β_{0} control the shape and scale (respectively) of the inverse-gamma prior distribution in the ARD hierarchy. This prior distribution controls the degree to which weights in the hidden layer are constrained, relative to a standard neural network. In our experience, the network is relatively insensitive to these values but please see [16] for the explicit role these parameters play. Each response was normalized to have mean 0, unit variance. All BNN code was written in Python and is avilable at https://github.com/beamandrew/BNN.
We used a Bonferonni cut-off value of p<0.05 for the MANOVA procedure. MANOVA was conducted using the manova() function in R. Processing each file file took the BNN approximately 1 minute, while the MANOVA procedure took approximately 5 seconds for all simulation settings.
Additive model
Both BNNs and MANOVA display good power across a variety of parameter combinations, with MANOVA having a slight edge for a few parameter combinations. This is expected as MANOVA is testing the linear hypothesis directly, while the BNN is testing a much more general hypothesis.
Additive model with interactions
Again, both methods appear to have good power across a wide spectrum of parameter values. However, MANOVA does not perform as well as for larger values of MAF. The BNN also experiences a loss of power for the highest level of MAF. For most parameter combinations, the BNN outperformed MANOVA.
Purely interactive model
Perhaps surprisingly, MANOVA displays good power to detect the causal loci for several parameter combinations, despite the true model ostensibly lacking any marginal effects. Both methods struggled for smaller values of MAF, but this is expected due the relatively rare nature of positive responses for small levels of MAF.
Further analysis of simulated models
The results of the preceding section motivated a more in-depth analysis of the reasons behind MANOVA’s apparent ability to detect models that only contain interactions and the reduced ability by both methods to detect the model containing main effects and an interaction term. The apparent loss of power by both methods for higher values of MAF in Figure 4 may be an artifact of the specific simulation parameters used. Notice for example, that both MANOVA and BNNs have good power to detect one locus for all simulated scenarios, as evidenced by the dotted line. This suggests that the loss of power under this model is due to the decreased ability to detect the second locus (S_{2}), as examination of the raw simulation results indicates this locus was rarely identified as significant. To examine this hypothesis, we will assume a more general population genetics framework and then reexamine the results of the simulated models in the previous section as specific instances of this broader viewpoint.
A marginal effect is the expected value of y_{ ik } given the status at one locus, averaged over all possible values for all other loci. Since the MANOVA procedure only detects marginal effects, a deeper understanding of how they may manifest in these putative models of genetic influence may be valuable. What we intend to show is that the discrete nature of the minor allele coding can introduce artifacts such as main effects in models where none are explicitly present and how marginal effects are attenuated by an interaction term in models containing both. First, we will examine the marginal effect of having a specified genotype at one locus (e.g. S_{2}=1) independently of the status at the other locus, and determine generally how this marginal effect changes under different simulated models as a function of effect size and MAF.
Relationship beteween μ _{ AI } and μ _{ A } for each possible level of S _{ 2 }
S_{2}Status | μ _{ AI } | μ _{ A } | |μ_{ AI }| ? |μ_{ A }| |
---|---|---|---|
0 | \(\frac {\theta }{2}(p_{10} + 2p_{20})\) | \(\frac {\theta }{2}(p_{10} + 2p_{20})\) | = |
1 | \(-\frac {\theta }{2}(p_{01} - p_{11} - 2p_{21})\) | \(-\frac {\theta }{2}(p_{01} - p_{21})\) | < |
2 | \(-\frac {\theta }{2}(2p_{02} - p_{12} + 2p_{22})\) | \(-\frac {\theta }{2}(2p_{02} + p_{12})\) | < |
Note that all of the marginal effects for S_{2} are smaller than or equal to the additive models marginal effect size. This explains why there was a decrease in power to detect S_{2} by MANOVA in the simulations.
Relationship beteween μ _{ I } and μ _{ A } for each possible level of S _{ 2 }
S_{2}Status | μ _{ I } | μ _{ A } | |μ_{ I }| ?|μ_{ A }| |
---|---|---|---|
0 | 0 | \(\frac {\theta }{2}(p_{10} + 2p_{20})\) | < |
1 | \(\frac {\theta }{2}(p_{11} + 2 p_{21})\) | \(-\frac {\theta }{2}(p_{01} - p_{21})\) | ? |
2 | \(\frac {\theta }{2}(p_{12} + 4p_{22})\) | \(-\frac {\theta }{2}(2p_{02} + p_{12})\) | ? |
Here there is no clear, strict inequality, but assessments can be made for specific values of p_{01},p_{11},p_{21}. Note that in general |μ_{ I }|>0 in all instances so there will always be some amount of marginal effect present for this model. For many values of MAF, |μ_{ I }| will be a non-trivial amount relative to |μ_{ A }| resulting in high power for a model can only detect additive effects, such as MANOVA. Thus the reason MANOVA in Figure 5 has good power for high values of MAF is because there actually are marginal effects present, even though they were not explicitly included in the model.
The analysis in this section is meant as a small complement to the vast literature on quantitative trait loci (QTL) and the role of epistasis in settings other than the dose-response framework. Several comprehensive investigations have been made using model organisms such as Drosophila melanogaster [15,26] for which there is considerable evidence of the role of epistasis in quantitative traits. How the results of these simple, 2 loci models might translate to larger epistatic networks in other contexts involving many more loci is yet unclear. However, the simulated models in this study suggest that epistatic interactions in a dose-response framework may manifest as marginal effects for each loci involved in the interaction, at least for some configurations of the minor allele frequency.
Analysis of the anticancer agent etoposide
In this section we apply the the Bayesian neural network to real dose-response data originally presented in [3]. We focus on the anticancer agent Etoposide which is a topoisomerase inhibitor used in the treatment of wide variety of cancers. Etoposide was chosen for analysis due to the high heritability of cytotoxicity for this compound (approximately 40% heritable) [27] obvserved in previous studies, yet estimated marginal effects have only been able to account for a small fraction of this heritability.
The cytotoxicity of Eptoposide at six concentrations was assessed via cell viability counts in lymphoblastoid cell lines (LCLs) of 520 individuals of European descent. Genotype status for each individual was determined for either 314,621 or 620,901 SNPs, using HumanHap300 bead chip or HumanQuad610 bead chip platforms [3]. Status for 2.5 million total SNPs was imputed using the approach in [28]. Quality control on the 2.5 million SNPs was performed as in [3]. SNPs failing a test for Hardy-Weinberg equilibrium at the p=0.01 level were removed, which resulted in approximately 30,000 SNPs being discarded. QC and initial data preparation were performed using the software package plink [29]. In the previous section it was shown that even SNPs contributing to genetic risk through an interaction term only will still most likely result in observed marginal effect. Using this rationale we first screened for marginal SNP effects using MAGWAS [17]. Following [3], we included temperature, growth rate, the first 3 principal components, and the experimental batch as covariates in the model. After screening for SNPs with statistically significant marginal effects at the p<10^{−5} level, we were left with 41 SNPs for analysis by the Bayesian neural network. Genome-wide screening using MAGWAS took approximately an hour on an Intel i5, quad-core desktop CPU.
We fit a neural network model with 10 hidden units, a logistic activation function and six Gaussian output units (one for each concentration). Sampling was done for 25 burn-in iterations, followed by 2,000 iterations to be used for inference. The sampler settings were as follows: a of step-size 4∗10^{−2}, L=20 leap-frog steps per iteration, and an initial annealing temperature of T_{0}=1000. The ARD prior parameters were set to α_{0}=5,β_{0}=2. Please see [16] for more detail and discussion on these parameter settings. Sampling took approximately 10 minutes on a GeForce GTX 650 desktop GPU. The HMC analysis procedure was performed for 3 independent replications and trace plots for each variable were inspected to assess if parameter estimates had converged.
The top SNP from [3] (rs2076112) has a Bayesian posterior probability which ranks it as the 20th most important SNP of the 41 SNPs passing the MAGWAS filter. As was found in [3], several experimental condition variables were found to be very important by the BNN as well. The top 7 most important variables were all related to experimental and cellular growth conditions, underscoring the importance of accounting for these conditions as covariates when performing cell-based gene-mapping. The top SNP according the BNN was rs12650820 (Bayesian probability 0.20, MAGWAS p-value 7.878∗10^{−6}), a SNP located on chromosome 4 within the follistatin-like 5 (FSTL5, location 4q32.3) gene. Dysfunction of FSTL5 has been implicated as biomarker in medulloblastoma, a highly malignant primary brain cancer [30]. All of the remaining SNPs have marginal importance scores which are suggestive of involvement, but less definitive.
Conclusions
In this study we have examined how gene-gene interactions can affect the ability to detect trait-associated loci in cell-based, dose-response association studies. We have presented a novel nonparametric procedure in the form of a Bayesian neural network and compared its performance to the MANOVA framework. Using a simulation study of plausible genetic models and a population genetics based analysis, we have shown that MANOVA may be able to detect causal loci even in the presence of genetic interactions. Additionally, we have shown that the BNN is also very capable of detecting causal loci across a range of possible genetic models. Each approach comes with trade-offs - the MANOVA approach is conceptually more simple and computationally less expensive while the BNN approach is more flexible and built upon fewer assumptions. In light of the results of this study, we recommend a two-stage approach (as was used in the Etoposide analysis) as a viable analysis strategy. MAGWAS is able to screen millions of SNPs quickly for marginal effects, which will mostly likely be present even if the underlying genetic architecture is epistatic. Next, the Bayesian neural network can be employed on a smaller subset of SNPs to explore a richer model space.
It’s also worth discussing some of the modelling benefits offered by the Bayesian neural network approach. Due to flexibility afforded by the BNN framework, new constraints or data types can be easily accommodated. Note that to model a new type of data, one only must be able to write down the likelihood and incorporate this into the output layer of the network. Thus, when analysing a new type of data, the only change that must be made is the type of output layer used, which reflects the new data’s likelihood function. Everything else, such as input layer, hidden layers, and ARD-testing framework, remain the same. There is no assumed model and no distributional assumptions, making this proposed framework very widely applicable to a variety of association mapping scenarios. More work remains to be done, however, before this approach can be used in a more turn-key manner. The most immediate improvement that can be made is to the HMC-based sampler. As with any MCMC method, sub-optimal parameter values will result in a chain failing to converge. The HMC-MCMC sampler used here is no exception and requires several tuning parameters that can dramatically impact the performance of the approach as a whole. Most often, several pilot runs are performed to find good values for these parameters. Recently, several groups have improved upon the generic HMC sampling scheme we used [31-33] which alleviate some of these issues. Incorporating these recent advances in HMC sampling techniques will help to make the method more robust and user-friendly.
In this study we have investigated how genetic interactions may affect the MANOVA analysis procedure, provided a deeper understanding of how these interactions may affect dose-response data, and presented a novel Bayesian nonparametric analysis technique. Cell-based dose-response studies hold much promise for pharmacogenomics, and the work presented here will help future research efforts utilizing this type of data.
Software availability
Software implementing the approaches outlined in this paper is available at https://github.com/beamandrew/BNN.
Declarations
Acknowledgements
AMR received funding from 1R01CA161608 from the National Cancer Institute. ALB received funding form NIEHS Bioinformatics Training Grant 5T32ES007329-14 - Ruth L. Kirschstein National Research Service Award (NRSA) Institutional Research Training Grants (T32).
Authors’ Affiliations
References
- Welsh M, Mangravite L, Medina MW, Tantisira K, Zhang W, Huang RS, et al. Pharmacogenomic discovery using cell-based models. Pharmacol Rev. 2009; 61(4):413–29.PubMedPubMed CentralView ArticleGoogle Scholar
- Wheeler HE, Dolan ME. Lymphoblastoid cell lines in pharmacogenomic discovery and clinical translation. Pharmacogenomics. 2012; 13(1):55–70.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown CC, Havener TM, Medina MW, Jack JR, Krauss RM, McLeod HL, et al. Genome-wide association and pharmacological profiling of 29 anticancer agents using lymphoblastoid cell lines. Pharmacogenomics. 2014; 15(2):137–46.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown CC, Havener TM, Medina MW, Auman JT, Mangravite LM, Krauss RM, et al. A genome-wide association analysis of temozolomide response using lymphoblastoid cell lines reveals a clinically relevant association with mgmt. Pharmacogenetics Genomics. 2012; 22(11):796.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown C, Havener TM, Everitt L, McLeod H, Motsinger-Reif AA. A comparison of association methods for cytotoxicity mapping in pharmacogenomics. Front Genet. 2011; 2:86.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown CC, Havener TM, Medina MW, Krauss RM, McLeod HL, Motsinger-Reif AA, et al. Multivariate methods and software for association mapping in dose-response genome-wide association studies. BioData Min. 2012; 5(1):21.PubMedPubMed CentralView ArticleGoogle Scholar
- Beam A, Motsinger-Reif A. Beyond ic50s: Towards robust statistical methods for in vitro association studies. J Pharmacogenom Pharmacoproteomics. 2013; 2(120):2153–645.Google Scholar
- Moore JH. The ubiquitous nature of epistasis in determining susceptibility to common human diseases. Human Heredity. 2003; 56(1-3):73–82.PubMedView ArticleGoogle Scholar
- Moore JH. A global view of epistasis. Nat Genet. 2005; 37(1):13–4.PubMedView ArticleGoogle Scholar
- Carlborg Ö, Haley CS. Epistasis: too often neglected in complex trait studies?Nat Rev Genet. 2004; 5(8):618–25.PubMedView ArticleGoogle Scholar
- Zuk O, Hechter E, Sunyaev SR, Lander ES. The mystery of missing heritability: Genetic interactions create phantom heritability. Proc Nat Acad Sci USA. 2012; 109(4):1193–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Wood AR, Tuke MA, Nalls MA, Hernandez DG, Bandinelli S, Singleton AB, et al. Another explanation for apparent epistasis. Nature. 2014; 514(7520):3–5.View ArticleGoogle Scholar
- Bahcall O. Global epistasis. Nat Genet. 2014; 46(8):811.View ArticleGoogle Scholar
- Hemani G, Shakhbazov K, Westra H-J, Esko T, Henders AK, McRae AF, et al. Detection and replication of epistasis influencing transcription in humans. Nature. 2014.Google Scholar
- Huang W, Richards S, Carbone MA, Zhu D, Anholt RR, Ayroles JF, et al. Epistasis dominates the genetic architecture of drosophila quantitative traits. Proc Nat Acad Sci USA. 2012; 109(39):15553–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Beam AL, Motsinger-Reif A, Doyle J. Bayesian neural networks for detecting epistasis in genetic association studies. BMC Bioinf. 2014; 15(1):368.View ArticleGoogle Scholar
- Brown CC, Havener TM, Medina MW, Krauss RM, McLeod HL, Motsinger-Reif AA. Multivariate methods and software for association mapping in dose-response genome-wide association studies. BioData Min. 2012; 5(1):1–15.View ArticleGoogle Scholar
- Motsinger-Reif AA, Ritchie MD. Neural networks for genetic epidemiology: past, present, and future. BioData Min. 2008; 1(3):3.PubMedPubMed CentralView ArticleGoogle Scholar
- Motsinger-Reif AA, Dudek SM, Hahn LW, Ritchie MD. Comparison of approaches for machine-learning optimization of neural networks for detecting gene-gene interactions in genetic epidemiology. Genet Epidemiol. 2008; 32(4):325–40.PubMedView ArticleGoogle Scholar
- Hornik K, Stinchcombe M, White H. Multilayer feedforward networks are universal approximators. Neural Networks. 1989; 2(5):359–66.View ArticleGoogle Scholar
- Neal RM. Assessing relevance determination methods using delve. Nato Asi Ser F Comput Syst Sci. 1998; 168:97–132.Google Scholar
- Wipf DP, Nagarajan SS. A new view of automatic relevance determination. In: Advances in Neural Information Processing Systems. Curran Associates, Inc.: 2007. p. 1625–32.Google Scholar
- Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH, Teller E. Equation of state calculations by fast computing machines. J Chem Phys. 1953; 21:1087.View ArticleGoogle Scholar
- Hastings WK. Monte carlo sampling methods using markov chains and their applications. Biometrika. 1970; 57(1):97–109.View ArticleGoogle Scholar
- Beam AL, Ghosh SK, Doyle J. Fast hamiltonian monte carlo using gpu computing. ArXiv e-prints. 2014; 1402:4089. Provided by the SAO/NASA Astrophysics Data System.Google Scholar
- Mackay TF, Stone EA, Ayroles JF. The genetics of quantitative traits: challenges and prospects. Nat Rev Genet. 2009; 10(8):565–77.PubMedView ArticleGoogle Scholar
- Peters EJ, Motsinger-Reif A, Havener TM, Everitt L, Hardison NE, Watson VG, et al. Pharmacogenomic characterization of us fda-approved cytotoxic drugs. Pharmacogenomics. 2011; 12(10):1407–15.PubMedPubMed CentralView ArticleGoogle Scholar
- Li Y, Willer C, Sanna S, Abecasis G. Genotype imputation. Annu Rev Genomics Human Genet. 2009; 10:387.View ArticleGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. Plink: a tool set for whole-genome association and population-based linkage analyses. Am J Human Genet. 2007; 81(3):559–75.View ArticleGoogle Scholar
- Remke M, Hielscher T, Korshunov A, Northcott PA, Bender S, Kool M, et al. Fstl5 is a marker of poor prognosis in non-wnt/non-shh medulloblastoma. J Clin Oncol. 2011; 2011:3852 - 61.View ArticleGoogle Scholar
- Girolami M, Calderhead B. Riemann manifold langevin and hamiltonian monte carlo methods. J R Stat Soc: Ser B (Stat Methodology). 2011; 73(2):123–214.View ArticleGoogle Scholar
- Sohl-Dickstein J, Mudigonda M, DeWeese M. Hamiltonian monte carlo without detailed balance. In: Proceedings of the 31st International Conference on Machine Learning. 2014. The Journal of Machine Learning Research. p. 719–26.Google Scholar
- Shahbaba B, Lan S, Johnson WO, Neal RM. Split hamiltonian monte carlo. Stat Comput. 2014; 24(3):339–49.View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.