BACKGROUND
genetic data collected in case-control or cohort studies of well-defined disease phenotypes can be used to identify genetic variants, typically single nucleotide polymorphisms , associated with disease status.
in a genetic association study, the data consist of genotypes and corresponding phenotypes from large numbers of individuals with the disease phenotype of interest  and disease-free controls.
a significant difference in the frequency of a genetic variant between the case and control groups is taken to be indicative of an association between that variant  and the phenotype of interest.
a framework also exists for analysing a continuous phenotypes.
once one or more associated variant have been identified, they can be included in predictive models for the estimation of disease risk in individuals for whom the  disease status is unknown.
technological developments, including commercially available chips for typing millions of genetic variants simultaneously, and next-generation sequencing technologies such as those used in the  <dig> genomes project  <cit> , are enabling the collection of large amounts of genetic data from large numbers of individuals which means that analysis of contemporary genetic data sets involves the study of high-dimensional data.
a number of statistical approaches have successfully been used to investigate the strength of association between genetic variants and a phenotype of interest.
these methods include testing the significance of the association of each variant with the phenotype independently using standard univariate tests .
however, there are disadvantages to relying on univariate methods for the analysis of genetic data.
the data from a genome-wide association study typically consists of thousands or millions of snps, and this large number of predictors introduces a multiple testing problem.
univariate methods therefore require stringent significance thresholds due to the large number of tests being undertaken to prevent a high false-positive rate  <cit> .
further, univariate methods fail to take into account the combined effect of multiple snps.
it is plausible that many genetic variants together contribute to the phenotype being studied  <cit> , and one advantage of using multivariate methods is that they allow for the study of the combined effect of multiple snps.
multivariate methods allow for the control of confounding variables, as in the lung cancer replication study in this paper, where gender, smoking status and age were included as unpenalised covariates in the model.
using multivariate methods, the regression model can be extended to include, for example, interaction or higher order terms, and in such a case a penalised regression approach would be appropriate  <cit> .
penalised regression methods have been applied to genetic data  <cit> .
among a number of regression approaches used for prediction in high-dimensional data, ridge regression has been shown to perform best in terms of prediction error  <cit> .
ridge regression has successfully been used to analyse genetic data where snps were in high ld  <cit> ; it is the test of significance used by malo, libiger & schork that we evaluate here.
to begin, we consider two regression models commonly used in the analysis of genetic data - the linear and the logistic regression models, as follows.
the standard linear regression model is given by
y=xβ+ε
where y is a  vector of dependent variables, yi, i =  and x is a  matrix of predictors.
β is a  vector of regression coefficients βj, j =  and ε is a  vector of normally distributed random errors, with εi ~iidn.
an example would be a model of the relationship between a continuous phenotype  measured in n individuals, and the genotype of these n individuals at m snps.
the ordinary least squares estimator for β is given by
β^=-1x′y
the significance of individual ols regression coefficients β^j in a multiple regression model can be estimated using a wald test.
the test statistic is
t0=β^jse
where se is an estimate of the standard error of the jth regression coefficient.
under the null hypothesis h0:β^j= <dig>  t <dig> follows a student t distribution with n - m degrees of freedom.
binary outcomes commonly arise in biomedical data where they may represent, for example, cases and controls.
in the logistic regression model, y is an n-dimensional vector of response variables taking values  <dig>  or  <dig> , and x the n × m matrix of explanatory variables, as before.
for the ith individual we denote
xi=
the ith response yi is a bernoulli variable with probability of success equal to pi.
the logistic regression model relates the probability pi that the ith observation is a case to the predictor variables as
p=pi=exiβ1+exiβ
where β is a vector of parameters to be estimated.
the significance of individual logistic regression coefficients, β^j, can be estimated using the test statistic
t0=β^jse
where se is an estimate of the standard error of the jth regression coefficient.
under the null hypothesis that βj =  <dig>  t <dig> asymptotically follows a standard normal distribution.
genetic data often comprises more predictor variables, m, than observations, n. in such a situation, unique maximum likelihood estimates of regression parameters do not exist.
further, collinearity in the predictors, due to linkage disequilibrium  in genetic data, which typically increases with the increasing density of available markers, results in unstable maximum likelihood estimates of regression coefficients.
an extensive literature exists on the application of modified regression techniques to the analysis of high-dimensional data.
penalised regression constrains the magnitude of the estimated regression coefficients, allowing their estimation when ordinary least squares  estimates cannot be obtained.
in a bayesian context, these techniques are equivalent to the specification of a particular prior distribution on the coefficients.
for example, lasso regression  <cit>  constrains the sum of the absolute value of the regression coefficients to be less than a constant.
this is equivalent to imposing a double exponential prior centred at zero on the coefficients.
lasso regression can estimate some coefficients to be exactly zero, permitting dimension reduction in the model.
hoggart, whittaker, de iorio & balding  <cit>  considered a modified lasso regression approach for the identification of causal snps in genome-wide or resequenced data, with the aim of identifying regions of association whilst considering all snps simultaneously.
when snps are in high ld, their method offers improvement over both single-snp analysis and lasso regression in terms of the power to detect causal variants, and a notable improvement over single-snp analysis in terms of false-positive rate.
ridge regression  <cit>  is a another penalised regression approach, in which a penalty is applied to the sum of the squared parameter estimates.
ridge regression has been used in a number of large-scale data analysis scenarios, including marker-assisted selection  <cit> , expression data analysis  <cit> , and genetic association studies when snps are in high ld  <cit> .
from a bayesian viewpoint, ridge regression can be considered as standard multiple regression with the coefficients estimates having a prior distribution that is normal with mean zero and known variance  <cit> .
in genetic epidemiology, it is desirable to estimate the strength of the association between a variant and a phenotype.
this is problematic when using ridge regression which, unlike other penalised regression approaches, does not reduce the number of parameters in the model, nor estimate the significance of each fitted coefficient.
a test of significance for coefficients estimated using ridge regression, based on an approximation of their distribution under the null hypothesis, was proposed by halawa & el bassiouni  <cit> .
the test was originally developed and evaluated for data with continuous outcomes, when different methods were used to compute the shrinkage parameter, λ, which controls the degree of shrinkage of the regression coefficients and hence their distribution under the null hypothesis.
malo, libiger & schork  <cit>  used the same test in an evaluation of the applicability of ridge regression as a means of accommodating ld in association studies.
they used the test in a comparison of the performance of ridge regression, multiple regression and single-snp analysis when snps are in varying degrees of ld.
they found that ridge regression identified different snps as associated with phenotype compared to single-snp analysis or multiple regression.
however, they did not consider the performance of the test itself.
in this paper, we evaluate the performance of a test of significance for ridge regression coefficients.
our test is based on the test proposed by halawa & el bassiouni  <cit> .
we extend the test, making it applicable in the m > n scenario that is common in contemporary genetic data sets.
we evaluate the performance of the test in simulation studies, using scenarios representative of realistic high-desnity genetic marker data, considering a range of data set dimensions and degrees of shrinkage.
ridge regression has also been applied in the logistic regression framework  <cit> .
we extend the test proposed by halawa & el bassiouni  <cit>  to the logistic ridge regression model, and again evaluate the test in a range of simulation scenarios at different values of λ.
in both linear and logistic ridge regression, we compare the approximate test of significance to a permutation test.
we view the permutation test as a benchmark as it produces an estimate of the null distribution of the parameter estimates.
however, the permutation test is computationally intensive and becomes more so when data are high-dimensional.
the test we propose makes it feasible to estimate significance with a much lower computational burden.
we introduce the p-value trace, a plot of the negative logarithm of the p-values of the ridge regression coefficients with increasing shrinkage parameter.
this plot enables the visualisation of the relative change in significance of each coefficient, and facilitates the identification of predictors most affected by increased penalisation in terms of significance.
we apply the approximate test of significance for logistic ridge regression coefficients to a lung cancer data set, demonstrating the utility of the test when correlation exists among the predictors.
this paper is organised as follows.
we first describe the approximate test of significance and the permutation test to which it is compared.
we then describe the simulation studies used in this paper.
in the results section we evaluate the performance of the proposed test in a range of simulation scenarios.
further, we apply the test to a lung cancer case control data set.
in the discussion we discuss the results and potential applications for the test.
methods
significance testing in linear ridge regression
ridge regression replaces the ols estimator β^ ) with the ridge regression estimator  <cit> :
β^λ=-1x′y
where collinearity exists in x, the ols estimates of β can be unstable, having large variance.
hoerl & kennard  <cit>  demonstrate that there exists a value of λ for which the ridge regression estimates β^λ have smaller mean square error  than the ols estimates.
where m > n, ols estimates of β cannot be obtained because the matrix x'x is singular.
the addition of the constant λ to the diagonal of the x'x matrix makes it invertible, so ridge regression estimates can be obtained.
halawa & el bassiouni  <cit>  use simulation studies to investigate applications of a 'non-exact' t-type test for the individual coefficients of a linear regression model fitted using ridge regression, based on the t-test above.
the test statistic is
tλ=β^jλse
where β^jλ is the estimate of the jth regression coefficient under the ridge regression model, and se is an estimate of the standard error.
estimates of the standard error of the jth element of β^λ are obtained as the square root of the jth element of the diagonal of the covariance matrix
var=σ2-1x′x- <dig>
in practise, σ <dig> is replaced by its estimate, given by the residual mean square of the ridge model:
σ^2=′ν
ν is the residual effective degrees of freedom.
halawa & el bassiouni  <cit>  use ν = n - m. however, when m > n this gives a negative estimate of the residual mean square.
instead, we use the definition of residual effective degrees of freedom given in hastie & tibshirani  <cit> , which makes use of the "hat matrix", h:
y^=xx′x+λi-1x′y
=hy
h is termed the 'hat matrix', because it 'puts the hat on' y, transforming it to Ŷ. degrees of freedom for error are defined as
ν=n−tr
in linear regression, the hat matrix reduces to h = x - <dig> x' and n - tr  reduces to n - m. the test statistic tλ is assumed to follow a student t distribution as in standard multiple linear regression.
however, the effective number of parameters of the penalised regression fit is smaller than m. hastie and tibshirani define tr  as the degrees of freedom taken up by the penalised regression fit  <cit> .tr  reduces to n - m in ordinary linear regression.
then, tλ is assumed to follow a student t distribution with n - tr  degrees of freedom.
in the case of large sample size, as is typically the case in genetic data, the distribution of the test statistic is asymptotically normal, as noted by halawa & el bassiouni  <cit> .
we compared the significance levels of the approximate test assuming both a normal and a student t-distribution of the test statistic and found that the results were substantially identical.
therefore we assume that under h <dig>  tλ~n <dig>  and use the normal distribution to test the significance of ridge regression coefficients.
the results from the corresponding tests assuming that under h <dig>  tλ ~ tn - tr are provided in an additional file.
see additional file  <dig>  tables s <dig> and s <dig>
significance testing in ridge logistic regression
ridge regression has been applied to the logistic regression model  <cit> .
cessie & van houwelingen  <cit>  show how ridge regression can be used to improve the parameter estimates in logistic regression when the number of predictors is relatively large or highly correlated.
they discuss different ways of choosing the shrinkage parameter to minimize prediction error.
vago & kennedy  <cit>  apply ridge logistic regression to a clinical data set.
in logistic ridge regression, the log-likelihood function is penalised with the penalty applied to the l <dig> norm of β  <cit> .
maximum likelihood estimates of β are obtained by maximising the logarithm of the likelihood function  <cit> , typically using the newton-raphson algorithm.
the approximate test statistic is
tλ=β^jλse
standard errors of the coefficient estimates are obtained as the square roots of the jth element of the diagonal of the covariance matrix.
this matrix is estimated from the final newton-raphson iteration:
varβ^λ=varx′wx+2λi-1x′wz=∂2ℓ∂β∂β′-1iβ∂2ℓ∂β∂β′-1=x′wx+2λi-1x′wxx′wx+2λi- <dig>
where i  is the observed information matrix, ℓ is the  log-likelihood, and w is the weight matrix:
w=diagpi ^pi ^=exiβ^λ1+exiβ^λ
z is an n ×  <dig> column vector with elements
zi=logitp^i+yi-p^ip^1-p^i
again we assume that under h <dig>  tλ~n <dig> and use the normal distribution to test the significance of ridge regression coefficients.
permutation test
to evaluate the proposed test, we compare its performance to that of a permutation test, which we view as a benchmark.
in the permutation test, to obtain a null distribution of the regression coefficients, the elements of the outcome vector y are randomly permuted.
the ridge regression model is fitted using the permuted observations, to obtain ridge regression coefficients.
by performing  <dig> such permutations, a null distribution of the regression coefficients is generated.
the permutation test p-value is calculated as the proportion of regression coefficients from the null distribution greater than or equal in absolute value to the absolute value of the coefficient fitted to the true  data.
choice of shrinkage parameter
hoerl & kennard  <cit>  present an existence theorem for ridge regression.
they demonstrate the existence of a value of the shrinkage parameter λ in equation  which will give estimates β^λ with smaller mean squared error than the ols estimates β^given in .
however, to date no analytical method to find the 'best' value of λ in terms of minimising mse has been determined.
a number of data-driven methods have been proposed.
these methods aim to determine a value of λ based on the data that will result in estimates of β^λwith improved mean squared error properties.
for example, hoerl, kennard & baldwin  <cit>  propose the following as an estimate of λ, with σ^ <dig> and β^taken from the ols estimates:
λhkb=mσ^2β^′β^
an alternative estimate was suggested by lawless & wang  <cit> , based on the adoption of the bayesian perspective mentioned above.
cross-validation based methods have also been proposed in the literature  <cit> .
hoerl & kennard  <cit>  introduce the ridge trace, a plot of the estimates β^λas λ increases from zero - see for example figure  <dig>  they propose choosing λ corresponding to the region on the ridge trace at which estimates of β^λ no longer change significantly as λ increases further.
following the ridge trace of hoerl & kennard  <cit> , we introduce a plot of p-values of the regression coefficients against λ as λ increases from zero .
we refer to this plot as a 'p-value trace'.
this p-value trace enables the visualisation of the change in p-values of the regression coefficients with increasing shrinkage.
simulation study
the proposed test was evaluated using simulated genotype data.
fregene software  <cit>  was used to simulate a population of haplotypes.
fregene simulates the forwards-in-time evolution of sequence-like genetic data.
the forward-in-time simulation allows demographic and selection scenarios to be implemented and recombination to be modelled.
haplotypes used in this study are taken from a simulation representing the neutral evolution of  <dig>  individuals over  <dig>  generations, with recombination parameters that realistically model recombination in humans.
this results in a population of  <dig>  haplotypes.
the simulated haplotypes, together with details of the simulation, are available to download from http://www.ebi.ac.uk/projects/bargen/ <cit> .
based on these haplotypes, genotypes and corresponding phenotypes were generated as follows.
a range of scenarios were considered, comprising n individuals at m snps.
<dig>  designate one snp as the causal snp, selected at random from those with population minor allele frequency in the range  <dig>  -  <dig> .
a subregion of m snps from the  <dig> mb region, containing this causal snp, is used for the genotype simulation.
the subregion is selected at random to be of length m and to contain the designated causal snp.
thus the correlation structure among the subregion depends on the correlation structure of the haplotype region at that point.
in the resultant sample, correlation with the causal snp ranges from low to perfect .
<dig>  sample two haplotypes  from the population of  <dig>  haplotypes.
sum the minor allele count at each snp to form a genotype.
<dig>  simulate the phenotype for this individual.
continuous phenotypes were generated as yi~n.
case-control phenotypes were generated following the liability model used by  <cit> .
the penetrance function, fk, is the probability of being a case, pr  given having k copies of the minor allele at the causal snp.
the genotype relative risk, r, is f <dig> /f <dig> and k is the population prevalence.
then, with the population frequency of the minor allele of the causal snp as p, under an additive genetic model, f <dig> = k/, f <dig> = rf <dig> and f <dig> = 2r f <dig> - f <dig>  a sample of n/ <dig> cases and n/ <dig> controls was generated by generating an individual genotype as described above, then assigning the individual to be a case with probability fk and a control otherwise.
this process is repeated until n/ <dig> cases and n/ <dig> controls are obtained.
<dig>  record the minor allele count  at the m snps for the ith individual, giving rise to an n × m matrix of minor allele counts.
ridge regression coefficients were fitted to data with continuous outcomes using lm.ridge from the package mass in r  <cit>  for both the simulated data and the permutation test.
estimates of regression coefficients under logistic ridge regression models were computed using the newton-raphson algorithm.
in the case of both continuous and binary outcomes, snps that were invariant in the sampled genotypes were removed from the data, and genotypes were standardised, prior to analysis.
the two tests were evaluated using the true positive rate  and false positive rate , averaged over all the replicates for each simulation scenario.
we define tpr to be the proportion of causal snps, as designated in the data simulation, significantly associated with phenotype at the nominal threshold α =  <dig> .
tpr is not reported for the null simulations, as there is no causal snp associated with phenotype in these data.
we define fpr to be the proportion of non-causal snps significantly associated with phenotype at the same significance threshold.
RESULTS
null simulation
genotypes and corresponding phenotypes were generated in two different sized data sets:  n =  <dig>  m =  <dig> and  n =  <dig>  m =  <dig>  in generating the null data, no snp was designated the causal snp.
continuous phenotypes were generated as yi~n; binary phenotypes were generated as yi ~ binom.
false positive rates are reported at the nominal significance threshold α =  <dig> .
four values of the ridge parameter λ were used.
results of the null simulations are shown in table  <dig>  we find that approximate test gives similar results to the permutation test in terms of false positive rate, especially in the case of continuous outcomes.
false positive rates at the nominal significance threshold α =  <dig>  in null datasets.
in each scenario, results are averaged over ten replicates
continuous phenotypes
genotypes and corresponding phenotypes were generated as described above, for a range of data sets dimensions: n =  <dig>   <dig>   <dig> and m =  <dig>   <dig>   <dig> and all snps in the  <dig> mb region .
phenotypes were generated as yi~n with k being the minor allele count at the causal snp.
a range of values of the shrinkage parameter λ were used: λ =  <dig> ,  <dig>   <dig>   <dig>  in figure  <dig>  the left column shows the null distributions, generated in a permutation test, used to estimate the significance of a ridge regression coefficient for a significant snp  and for a snp that is not associated with phenotype .
the coefficient fitted to the original data is indicated.
in the right column, the null distribution of the test statistic used in the approximate test is shown, with the test statistic of the fitted coefficient indicated.
ridge regression models were fitted using the shrinkage parameter λ =  <dig>  these results are examples taken from a single simulation, and above each plot the p-value according to the permutation test  or the approximate test  is shown.
performance comparison between a permutation test and the approximate test.
data are simulated genotype data with continuous phenotypes.
reported are proportion of true positive and false positive results at significance threshold α =  <dig> .
tpr = true positive rate, fpr = false positive rate.
results for each simulation scenario are averaged over ten replicates.
when ranking the snps in order of significance, the approximate test and the permutation test ranked the snps identically or nearly so .
continuous phenotypes with multiple causal snps
for complex diseases, multiple causal snps are likely to affect the phenotype.
we investigated the performance of the test when more than one snp is associated with phenotype.
we simulated data from two different scenarios: n =  <dig>  m =  <dig> and n =  <dig>  m =  <dig>  in each region of m simulated genotypes, ten snps with minor allele frequency  <dig>  -  <dig>  were designated causal and given effect size 1; the non-causal snps had effect size  <dig>  phenotypes were simulated as y = xβ + ε, ε~n, σ <dig> =  <dig> where β is the vector of effect sizes.
results are presented in table  <dig>  we see that the conclusions drawn about the test, of adequate power at the cost of a slightly higher false positive rate, are equally valid when multiple snps in the data are associated with phenotype.
performance comparison between a permutation test and the approximate test in simulated genotype data with continuous phenotypes and multiple causal snps.
reported are proportion of true positive and false positive results at significance threshold α =  <dig> .
n =  <dig> and in each scenario ten snps are with maf  <dig>  -  <dig>  are designated causal with effect size 1; the rest of the snps have  <dig> effect size.
tpr = true positive rate, fpr = false positive rate.
results for each simulation scenario are averaged over ten replicates.
computational performance comparison
using an example simulation, we compared the computational time required to compute the approximate and the permutation tests.
a data set with dimensions n =  <dig>  m =  <dig> and λ =  <dig> was used.
approximate test and permutation test p - values were computed and the time taken to arrive at the p-values was recorded.
calculations were done using r version  <dig> .
<dig>  <cit>  on an imac running mac os x version  <dig> .
<dig>  fitted with an  <dig>  ghz intel core i <dig> processor and  <dig> gb  <dig> mhz ddr <dig> ram.
computational times are compared in table  <dig>  we see that the permutation test takes approximately  <dig> times longer to compute than the approximate test.
comparison of computation performance of approximate test and permutation test.
n =  <dig>  m =  <dig> and λ =  <dig>  approximate test and permutation test p - values were computed and the time taken to arrive at the p-values was recorded.
calculations were done using r version  <dig> .
<dig>  <cit>  on an imac running mac os x version  <dig> .
<dig>  fitted with an  <dig>  ghz intel core i <dig> processor and  <dig> gb  <dig> mhz ddr <dig> ram.
binary phenotypes
genotypes and corresponding binary phenotypes were generated for nine different data set dimensions: n =  <dig>   <dig> and m =  <dig>   <dig>   <dig>   <dig>  and n =  <dig>  m = all snps in the  <dig> mb region .
the genotype relative risk, r, was specified as  <dig>  for the largest data dimensions , results are not shown for the permutation test, due to the computational time required to fit ridge regression models on data of this size.
similarly, results are not shown when n =  <dig> and m ≈  <dig>   <dig>
in figure  <dig>  the left column shows the null distributions, generated in a permutation test, used to estimate the significance of a ridge regression coefficient for a significant snp  and for a snp that is not associated with  in data with binary outcomes.
the coefficient fitted to the original data is indicated.
in the right column, the null distribution of the test statistic used in the approximate test is shown, with the test statistic of the fitted coefficient indicated.
ridge regression models were fitted using the shrinkage parameter λ =  <dig>  these results are examples taken from a single simulation, and above each plot the p-value according to the permutation test  or the approximate test  is shown.
performance comparison between a permutation test and the approximate test in simulated genotype data with case-control phenotypes.
reported are proportion of true positive and false positive results at significance threshold α =  <dig> .
tpr = true positive rate, fpr = false positive rate.
results for each simulation scenario are averaged over ten replicates.
comparison with univariate tests of significance
the performance of tests of significance of ridge regression coefficients, in terms of true and false positive rates, was compared to the performance of univariate tests.
comparisons were made in each of the simulation settings: the null simulation, and the simulations with continuous and with binary outcomes.
the results are shown in additional file  <dig>  tables s <dig>  s <dig> and s <dig>  as would be expected when using a penalised regression approach such as ridge regression, the performance of the corresponding significance test depends on the degree of shrinkage, which is controlled by the shrinkage parameter λ. the performance of the approximate test is comparable to that of a permutation test, with the advantage of a much-reduced computational burden.
further, ridge regression has the advantage over univariate tests of significance that it results in a much lower false positive rate.
these advantages of ridge regression compared to univariate methods are further illustrated in the study of lung cancer data which follows.
lung cancer data
genome-wide association studies have identified snps associated with lung cancer disease status.
snps have been identified at chromosomal locations 15q <dig>  <cit> , 5p <dig>  <cit>  and 6p <dig>  <cit> .
the associations at 15q <dig> and 5p <dig> have been replicated in white populations, but the association at 6p <dig> has not  <cit> .
not all studies successfully replicated the associations at 15q <dig> and 5p <dig>  <cit> .
here, we use ridge regression to re-evaluate a set of  <dig> snps for association with lung cancer disease status.
whilst these data are not as high-dimensional as those from a genome-wide study, they allow us to illustrate the features of using ridge regression for genetic data.
we show that ridge regression is a useful technique when data are correlated, and illustrate that multivariate methods have advantages over univariate tests of significance.
data consist of genotypes and non-genetic predictors from approximately  <dig> individuals in the european prospective investigation into cancer and nutrition .
missing genotypes were imputed using mean imputation.
gender, smoking status and age were included as unconstrained parameters in the model.
for the purpose of comparison, univariate  p-values were calculated.
univariate p-values were calculated by fitting a logistic regression model for each snp independently, with gender, smoking status and age included in the model.
a wald test was then used to estimate the significance of the coefficient.
in this example, m < n and  multivariate regression does give rise to unique parameter estimates.
multivariate p-values  are also reported here.
in a data set of this size, it is computationally feasible to fit a ridge regression model with a range of values of λ. the ridge trace  is a plot of parameter estimates against λ  <cit> .
the ridge trace typically suggests a range of values of λ rather than a single best value.
in figure  <dig>  λ =  <dig> seems to be the point at which the ridge estimates stabilise, and table  <dig> reports the corresponding approximate p-values.
of the  <dig> snps in the lung cancer data set, four have been previously reported to be associated with lung cancer disease status.
for these four snps, this table reports univariate, multivariate, approximate and permutation test p-values.
further, with increasing shrinkage and stabilisation of the estimates, the approximate p-values for these snps become significant .
rs <dig>  which is not in ld with other snps in the data, is significant in a multiple logistic regression model even when no penalty is included in the model .
a further snp, rs <dig>  in chromosome  <dig>  has a nominally significant p-value in the approximate test , but again this snp has not previously been shown to be associated with lung cancer risk.
the p-value trace using permutation test p-values  shows good agreement with the approximate p-value trace .
figure 2b is much more computationally expensive to produce than figure 2a, thus the use of the approximate test makes the plotting of a p-value trace for a range of values of λ more feasible.
the results in table  <dig> demonstrate the advantage of multivariate tests, and specifically of ridge regression, over univariate tests of significance.
snp rs <dig>  which has previously been shown to be associated with lung cancer disease status  <cit> , was not found to be significant using the univariate test but was found as such by the multivariate methods.
the two snps that are correlated, rs <dig> and rs <dig>  were not significant in multiple regression but were significant in ridge regression, demonstrating the advantage of using ridge regression when snps are correlated.
these results demonstrate that this approximate test of significance for coefficients fitted using logistic ridge regression reproduces previously ascertained associations, at reduced computational cost compared to a permutation test, even when snps are highly correlated.
CONCLUSIONS
we present and evaluate the performance of a test of significance for coefficients estimated using ridge regression.
we evaluate the test as applied to both linear and logistic ridge regression models.
our evaluation is by means of simulation studies across a range of scenarios representative of genetic data.
we evaluate the test by comparing its performance to that of a permutation test.
we evaluate the performance of the test when it is applied to a real data set.
the data set comprises lifestyle data and genotypes together with lung cancer case-control status.
using the proposed test, we successfully replicate previously found associations at much reduced computational cost compared to a permutation test.
this demonstrates the utility of the test for detecting significant variables when predictor variables are highly correlated, as were two significant snps in the lung cancer data.
stability selection  <cit>  is a method for variable selection that has received attention in the literature in recent years.
it combines subsampling of the data with a dimension reduction technique, with the aim of finding consistently significant variables.
penalised regression methods that perform variable selection, such as the lasso  <cit>  and the elastic net  <cit> , have been used in stability selection.
to use ridge regression with stability selection, a way of determining which variables to select in each subsample of the data is required.
a permutation test would be computationally expensive in this context, rapidly becoming infeasible if for large n, if the guidelines of  <dig> subsamples of size n/ <dig> given by meinshausen & bühlmann  <cit>  were followed.
the test of significance proposed here, being much less computationally expensive than a permutation test, makes the combination of ridge regression and stability selection a feasible possibility.
the choice of shrinkage parameter in ridge regression is discussed in the literature, but no consensus method provides an universally optimum choice.
the proposed test performs well over a range of values of λ. we introduce the p-value trace, a plot of the change in the negative logarithm of the p-value of the regression coefficients as the shrinkage parameter increases from zero.
this trace can be used in combination with the ridge trace of  <cit>  to visualise how the significance of the regression parameters, as well as their value, changes with increasing shrinkage.
such a plot could aid the identification of significant coefficients in the ridge model.
authors' contributions
ec ran the simulations, analysed the data, and wrote the paper.
pv was involved in the collection of the lung cancer data, supervised the analysis, and read and commented on the draft.
mdi came up with the idea, supervised the work, and commented on drafts of the paper.
all authors read and approved the final manuscript.
supplementary material
additional file 1
table s <dig> - performance comparison in null simulation using t-type test.
table s <dig> - performance comparison in simulated data with continuous outcomes using t-type test.
table s <dig> - performance comparison in null simulation including comparison to univariate tests of significance.
table s <dig> - performance comparison with continuous outcomes including comparison to univariate tests of significance.
table s <dig> - performance comparison with binary outcomes including comparison to univariate tests of significance.
click here for file
acknowledgements
we thank the european prospective investigation into cancer and nutrition : elio riboli, anne tjønneland, kim overvad, francoise clavel-chapelon, rudolf kaaks, heiner boeing, antonia trichopoulou, franco berrino, domenico palli, salvatore panico, rosario tumino, paolo vineis, bas buneo de mesquita, petra peeters, eiliv lund, carlos a. gonzález, aurelio barricarte, miren dorronsoro, maria josé sánchez, carmen navarro, j. ramón quirós, jonas manjer, goran hallmans, nicholas e. day, kay-tee khaw, nicholas wareham and timothy j. key.

 genetic data collected in case-control or cohort studies of well-defined disease phenotypes can be used to identify genetic variants, typically single nucleotide polymorphisms, associated with disease status. in a genetic association study, the data consist of genotypes and corresponding phenotypes from large numbers of individuals with the disease phenotype of interest and disease-free controls. a significant difference in the frequency of a genetic variant between the case and control groups is taken to be indicative of an association between that variant and the phenotype of interest. a framework also exists for analysing a continuous phenotypes. once one or more associated variant have been identified, they can be included in predictive models for the estimation of disease risk in individuals for whom the disease status is unknown. technological developments, including commercially available chips for typing millions of genetic variants simultaneously, and next-generation sequencing technologies such as those used in the genomes project, are enabling the collection of large amounts of genetic data from large numbers of individuals which means that analysis of contemporary genetic data sets involves the study of high-dimensional data. a number of statistical approaches have successfully been used to investigate the strength of association between genetic variants and a phenotype of interest. these methods include testing the significance of the association of each variant with the phenotype independently using standard univariate tests. however, there are disadvantages to relying on univariate methods for the analysis of genetic data. the data from a genome-wide association study typically consists of thousands or millions of snps, and this large number of predictors introduces a multiple testing problem. univariate methods therefore require stringent significance thresholds due to the large number of tests being undertaken to prevent a high false-positive rate. further, univariate methods fail to take into account the combined effect of multiple snps. it is plausible that many genetic variants together contribute to the phenotype being studied, and one advantage of using multivariate methods is that they allow for the study of the combined effect of multiple snps. multivariate methods allow for the control of confounding variables, as in the lung cancer replication study in this paper, where gender, smoking status and age were included as unpenalised covariates in the model. using multivariate methods, the regression model can be extended to include, for example, interaction or higher order terms, and in such a case a penalised regression approach would be appropriate. penalised regression methods have been applied to genetic data. among a number of regression approaches used for prediction in high-dimensional data, ridge regression has been shown to perform best in terms of prediction error. ridge regression has successfully been used to analyse genetic data where snps were in high ld; it is the test of significance used by malo, libiger & schork that we evaluate here. to begin, we consider two regression models commonly used in the analysis of genetic data - the linear and the logistic regression models, as follows. the standard linear regression model is given by y=xβ+ε where y is a vector of dependent variables, yi, i = and x is a matrix of predictors. β is a vector of regression coefficients βj, j = and ε is a vector of normally distributed random errors, with εi ~iidn. an example would be a model of the relationship between a continuous phenotype measured in n individuals, and the genotype of these n individuals at m snps. the ordinary least squares estimator for β is given by β^=-1x′y the significance of individual ols regression coefficients β^j in a multiple regression model can be estimated using a wald test. the test statistic is t0=β^jse where se is an estimate of the standard error of the jth regression coefficient. under the null hypothesis h0:β^j= t follows a student t distribution with n - m degrees of freedom. binary outcomes commonly arise in biomedical data where they may represent, for example, cases and controls. in the logistic regression model, y is an n-dimensional vector of response variables taking values or, and x the n × m matrix of explanatory variables, as before. for the ith individual we denote xi= the ith response yi is a bernoulli variable with probability of success equal to pi. the logistic regression model relates the probability pi that the ith observation is a case to the predictor variables as p=pi=exiβ1+exiβ where β is a vector of parameters to be estimated. the significance of individual logistic regression coefficients, β^j, can be estimated using the test statistic t0=β^jse where se is an estimate of the standard error of the jth regression coefficient. under the null hypothesis that βj = t asymptotically follows a standard normal distribution. genetic data often comprises more predictor variables, m, than observations, n. in such a situation, unique maximum likelihood estimates of regression parameters do not exist. further, collinearity in the predictors, due to linkage disequilibrium in genetic data, which typically increases with the increasing density of available markers, results in unstable maximum likelihood estimates of regression coefficients. an extensive literature exists on the application of modified regression techniques to the analysis of high-dimensional data. penalised regression constrains the magnitude of the estimated regression coefficients, allowing their estimation when ordinary least squares estimates cannot be obtained. in a bayesian context, these techniques are equivalent to the specification of a particular prior distribution on the coefficients. for example, lasso regression constrains the sum of the absolute value of the regression coefficients to be less than a constant. this is equivalent to imposing a double exponential prior centred at zero on the coefficients. lasso regression can estimate some coefficients to be exactly zero, permitting dimension reduction in the model. hoggart, whittaker, de iorio & balding considered a modified lasso regression approach for the identification of causal snps in genome-wide or resequenced data, with the aim of identifying regions of association whilst considering all snps simultaneously. when snps are in high ld, their method offers improvement over both single-snp analysis and lasso regression in terms of the power to detect causal variants, and a notable improvement over single-snp analysis in terms of false-positive rate. ridge regression is a another penalised regression approach, in which a penalty is applied to the sum of the squared parameter estimates. ridge regression has been used in a number of large-scale data analysis scenarios, including marker-assisted selection, expression data analysis, and genetic association studies when snps are in high ld. from a bayesian viewpoint, ridge regression can be considered as standard multiple regression with the coefficients estimates having a prior distribution that is normal with mean zero and known variance. in genetic epidemiology, it is desirable to estimate the strength of the association between a variant and a phenotype. this is problematic when using ridge regression which, unlike other penalised regression approaches, does not reduce the number of parameters in the model, nor estimate the significance of each fitted coefficient. a test of significance for coefficients estimated using ridge regression, based on an approximation of their distribution under the null hypothesis, was proposed by halawa & el bassiouni. the test was originally developed and evaluated for data with continuous outcomes, when different methods were used to compute the shrinkage parameter, λ, which controls the degree of shrinkage of the regression coefficients and hence their distribution under the null hypothesis. malo, libiger & schork used the same test in an evaluation of the applicability of ridge regression as a means of accommodating ld in association studies. they used the test in a comparison of the performance of ridge regression, multiple regression and single-snp analysis when snps are in varying degrees of ld. they found that ridge regression identified different snps as associated with phenotype compared to single-snp analysis or multiple regression. however, they did not consider the performance of the test itself. in this paper, we evaluate the performance of a test of significance for ridge regression coefficients. our test is based on the test proposed by halawa & el bassiouni. we extend the test, making it applicable in the m > n scenario that is common in contemporary genetic data sets. we evaluate the performance of the test in simulation studies, using scenarios representative of realistic high-desnity genetic marker data, considering a range of data set dimensions and degrees of shrinkage. ridge regression has also been applied in the logistic regression framework. we extend the test proposed by halawa & el bassiouni to the logistic ridge regression model, and again evaluate the test in a range of simulation scenarios at different values of λ. in both linear and logistic ridge regression, we compare the approximate test of significance to a permutation test. we view the permutation test as a benchmark as it produces an estimate of the null distribution of the parameter estimates. however, the permutation test is computationally intensive and becomes more so when data are high-dimensional. the test we propose makes it feasible to estimate significance with a much lower computational burden. we introduce the p-value trace, a plot of the negative logarithm of the p-values of the ridge regression coefficients with increasing shrinkage parameter. this plot enables the visualisation of the relative change in significance of each coefficient, and facilitates the identification of predictors most affected by increased penalisation in terms of significance. we apply the approximate test of significance for logistic ridge regression coefficients to a lung cancer data set, demonstrating the utility of the test when correlation exists among the predictors. this paper is organised as follows. we first describe the approximate test of significance and the permutation test to which it is compared. we then describe the simulation studies used in this paper. in the results section we evaluate the performance of the proposed test in a range of simulation scenarios. further, we apply the test to a lung cancer case control data set. in the discussion we discuss the results and potential applications for the test. methods significance testing in linear ridge regression ridge regression replaces the ols estimator β^ ) with the ridge regression estimator: β^λ=-1x′y where collinearity exists in x, the ols estimates of β can be unstable, having large variance. hoerl & kennard demonstrate that there exists a value of λ for which the ridge regression estimates β^λ have smaller mean square error than the ols estimates. where m > n, ols estimates of β cannot be obtained because the matrix x'x is singular. the addition of the constant λ to the diagonal of the x'x matrix makes it invertible, so ridge regression estimates can be obtained. halawa & el bassiouni use simulation studies to investigate applications of a'non-exact' t-type test for the individual coefficients of a linear regression model fitted using ridge regression, based on the t-test above. the test statistic is tλ=β^jλse where β^jλ is the estimate of the jth regression coefficient under the ridge regression model, and se is an estimate of the standard error. estimates of the standard error of the jth element of β^λ are obtained as the square root of the jth element of the diagonal of the covariance matrix var=σ2-1x′x- in practise, σ is replaced by its estimate, given by the residual mean square of the ridge model: σ^2=′ν ν is the residual effective degrees of freedom. halawa & el bassiouni use ν = n - m. however, when m > n this gives a negative estimate of the residual mean square. instead, we use the definition of residual effective degrees of freedom given in hastie & tibshirani, which makes use of the "hat matrix", h: y^=xx′x+λi-1x′y =hy h is termed the'hat matrix', because it'puts the hat on' y, transforming it to Ŷ. degrees of freedom for error are defined as ν=n−tr in linear regression, the hat matrix reduces to h = x - x' and n - tr reduces to n - m. the test statistic tλ is assumed to follow a student t distribution as in standard multiple linear regression. however, the effective number of parameters of the penalised regression fit is smaller than m. hastie and tibshirani define tr as the degrees of freedom taken up by the penalised regression fit.tr reduces to n - m in ordinary linear regression. then, tλ is assumed to follow a student t distribution with n - tr degrees of freedom. in the case of large sample size, as is typically the case in genetic data, the distribution of the test statistic is asymptotically normal, as noted by halawa & el bassiouni. we compared the significance levels of the approximate test assuming both a normal and a student t-distribution of the test statistic and found that the results were substantially identical. therefore we assume that under h tλ~n and use the normal distribution to test the significance of ridge regression coefficients. the results from the corresponding tests assuming that under h tλ ~ tn - tr are provided in an additional file. see additional file tables s and s significance testing in ridge logistic regression ridge regression has been applied to the logistic regression model. cessie & van houwelingen show how ridge regression can be used to improve the parameter estimates in logistic regression when the number of predictors is relatively large or highly correlated. they discuss different ways of choosing the shrinkage parameter to minimize prediction error. vago & kennedy apply ridge logistic regression to a clinical data set. in logistic ridge regression, the log-likelihood function is penalised with the penalty applied to the l norm of β. maximum likelihood estimates of β are obtained by maximising the logarithm of the likelihood function, typically using the newton-raphson algorithm. the approximate test statistic is tλ=β^jλse standard errors of the coefficient estimates are obtained as the square roots of the jth element of the diagonal of the covariance matrix. this matrix is estimated from the final newton-raphson iteration: varβ^λ=varx′wx+2λi-1x′wz=∂2ℓ∂β∂β′-1iβ∂2ℓ∂β∂β′-1=x′wx+2λi-1x′wxx′wx+2λi- where i is the observed information matrix, ℓ is the log-likelihood, and w is the weight matrix: w=diagpi ^pi ^=exiβ^λ1+exiβ^λ z is an n × column vector with elements zi=logitp^i+yi-p^ip^1-p^i again we assume that under h tλ~n and use the normal distribution to test the significance of ridge regression coefficients. permutation test to evaluate the proposed test, we compare its performance to that of a permutation test, which we view as a benchmark. in the permutation test, to obtain a null distribution of the regression coefficients, the elements of the outcome vector y are randomly permuted. the ridge regression model is fitted using the permuted observations, to obtain ridge regression coefficients. by performing such permutations, a null distribution of the regression coefficients is generated. the permutation test p-value is calculated as the proportion of regression coefficients from the null distribution greater than or equal in absolute value to the absolute value of the coefficient fitted to the true data. choice of shrinkage parameter hoerl & kennard present an existence theorem for ridge regression. they demonstrate the existence of a value of the shrinkage parameter λ in equation which will give estimates β^λ with smaller mean squared error than the ols estimates β^given in. however, to date no analytical method to find the'best' value of λ in terms of minimising mse has been determined. a number of data-driven methods have been proposed. these methods aim to determine a value of λ based on the data that will result in estimates of β^λwith improved mean squared error properties. for example, hoerl, kennard & baldwin propose the following as an estimate of λ, with σ^ and β^taken from the ols estimates: λhkb=mσ^2β^′β^ an alternative estimate was suggested by lawless & wang, based on the adoption of the bayesian perspective mentioned above. cross-validation based methods have also been proposed in the literature. hoerl & kennard introduce the ridge trace, a plot of the estimates β^λas λ increases from zero - see for example figure they propose choosing λ corresponding to the region on the ridge trace at which estimates of β^λ no longer change significantly as λ increases further. following the ridge trace of hoerl & kennard, we introduce a plot of p-values of the regression coefficients against λ as λ increases from zero. we refer to this plot as a'p-value trace'. this p-value trace enables the visualisation of the change in p-values of the regression coefficients with increasing shrinkage. simulation study the proposed test was evaluated using simulated genotype data. fregene software was used to simulate a population of haplotypes. fregene simulates the forwards-in-time evolution of sequence-like genetic data. the forward-in-time simulation allows demographic and selection scenarios to be implemented and recombination to be modelled. haplotypes used in this study are taken from a simulation representing the neutral evolution of individuals over generations, with recombination parameters that realistically model recombination in humans. this results in a population of haplotypes. the simulated haplotypes, together with details of the simulation, are available to download from http://www.ebi.ac.uk/projects/bargen/. based on these haplotypes, genotypes and corresponding phenotypes were generated as follows. a range of scenarios were considered, comprising n individuals at m snps. designate one snp as the causal snp, selected at random from those with population minor allele frequency in the range -. a subregion of m snps from the mb region, containing this causal snp, is used for the genotype simulation. the subregion is selected at random to be of length m and to contain the designated causal snp. thus the correlation structure among the subregion depends on the correlation structure of the haplotype region at that point. in the resultant sample, correlation with the causal snp ranges from low to perfect. sample two haplotypes from the population of haplotypes. sum the minor allele count at each snp to form a genotype. simulate the phenotype for this individual. continuous phenotypes were generated as yi~n. case-control phenotypes were generated following the liability model used by. the penetrance function, fk, is the probability of being a case, pr given having k copies of the minor allele at the causal snp. the genotype relative risk, r, is f /f and k is the population prevalence. then, with the population frequency of the minor allele of the causal snp as p, under an additive genetic model, f = k/, f = rf and f = 2r f - f a sample of n/ cases and n/ controls was generated by generating an individual genotype as described above, then assigning the individual to be a case with probability fk and a control otherwise. this process is repeated until n/ cases and n/ controls are obtained. record the minor allele count at the m snps for the ith individual, giving rise to an n × m matrix of minor allele counts. ridge regression coefficients were fitted to data with continuous outcomes using lm.ridge from the package mass in r for both the simulated data and the permutation test. estimates of regression coefficients under logistic ridge regression models were computed using the newton-raphson algorithm. in the case of both continuous and binary outcomes, snps that were invariant in the sampled genotypes were removed from the data, and genotypes were standardised, prior to analysis. the two tests were evaluated using the true positive rate and false positive rate, averaged over all the replicates for each simulation scenario. we define tpr to be the proportion of causal snps, as designated in the data simulation, significantly associated with phenotype at the nominal threshold α =. tpr is not reported for the null simulations, as there is no causal snp associated with phenotype in these data. we define fpr to be the proportion of non-causal snps significantly associated with phenotype at the same significance threshold.  null simulation genotypes and corresponding phenotypes were generated in two different sized data sets: n = m = and n = m = in generating the null data, no snp was designated the causal snp. continuous phenotypes were generated as yi~n; binary phenotypes were generated as yi ~ binom. false positive rates are reported at the nominal significance threshold α =. four values of the ridge parameter λ were used. results of the null simulations are shown in table we find that approximate test gives similar results to the permutation test in terms of false positive rate, especially in the case of continuous outcomes. false positive rates at the nominal significance threshold α = in null datasets. in each scenario, results are averaged over ten replicates continuous phenotypes genotypes and corresponding phenotypes were generated as described above, for a range of data sets dimensions: n = and m = and all snps in the mb region. phenotypes were generated as yi~n with k being the minor allele count at the causal snp. a range of values of the shrinkage parameter λ were used: λ =, in figure the left column shows the null distributions, generated in a permutation test, used to estimate the significance of a ridge regression coefficient for a significant snp and for a snp that is not associated with phenotype. the coefficient fitted to the original data is indicated. in the right column, the null distribution of the test statistic used in the approximate test is shown, with the test statistic of the fitted coefficient indicated. ridge regression models were fitted using the shrinkage parameter λ = these results are examples taken from a single simulation, and above each plot the p-value according to the permutation test or the approximate test is shown. performance comparison between a permutation test and the approximate test. data are simulated genotype data with continuous phenotypes. reported are proportion of true positive and false positive results at significance threshold α =. tpr = true positive rate, fpr = false positive rate. results for each simulation scenario are averaged over ten replicates. when ranking the snps in order of significance, the approximate test and the permutation test ranked the snps identically or nearly so. continuous phenotypes with multiple causal snps for complex diseases, multiple causal snps are likely to affect the phenotype. we investigated the performance of the test when more than one snp is associated with phenotype. we simulated data from two different scenarios: n = m = and n = m = in each region of m simulated genotypes, ten snps with minor allele frequency - were designated causal and given effect size 1; the non-causal snps had effect size phenotypes were simulated as y = xβ + ε, ε~n, σ = where β is the vector of effect sizes. results are presented in table we see that the conclusions drawn about the test, of adequate power at the cost of a slightly higher false positive rate, are equally valid when multiple snps in the data are associated with phenotype. performance comparison between a permutation test and the approximate test in simulated genotype data with continuous phenotypes and multiple causal snps. reported are proportion of true positive and false positive results at significance threshold α =. n = and in each scenario ten snps are with maf - are designated causal with effect size 1; the rest of the snps have effect size. tpr = true positive rate, fpr = false positive rate. results for each simulation scenario are averaged over ten replicates. computational performance comparison using an example simulation, we compared the computational time required to compute the approximate and the permutation tests. a data set with dimensions n = m = and λ = was used. approximate test and permutation test p - values were computed and the time taken to arrive at the p-values was recorded. calculations were done using r version. on an imac running mac os x version. fitted with an ghz intel core i processor and gb mhz ddr ram. computational times are compared in table we see that the permutation test takes approximately times longer to compute than the approximate test. comparison of computation performance of approximate test and permutation test. n = m = and λ = approximate test and permutation test p - values were computed and the time taken to arrive at the p-values was recorded. calculations were done using r version. on an imac running mac os x version. fitted with an ghz intel core i processor and gb mhz ddr ram. binary phenotypes genotypes and corresponding binary phenotypes were generated for nine different data set dimensions: n = and m = and n = m = all snps in the mb region. the genotype relative risk, r, was specified as for the largest data dimensions, results are not shown for the permutation test, due to the computational time required to fit ridge regression models on data of this size. similarly, results are not shown when n = and m ≈ in figure the left column shows the null distributions, generated in a permutation test, used to estimate the significance of a ridge regression coefficient for a significant snp and for a snp that is not associated with in data with binary outcomes. the coefficient fitted to the original data is indicated. in the right column, the null distribution of the test statistic used in the approximate test is shown, with the test statistic of the fitted coefficient indicated. ridge regression models were fitted using the shrinkage parameter λ = these results are examples taken from a single simulation, and above each plot the p-value according to the permutation test or the approximate test is shown. performance comparison between a permutation test and the approximate test in simulated genotype data with case-control phenotypes. reported are proportion of true positive and false positive results at significance threshold α =. tpr = true positive rate, fpr = false positive rate. results for each simulation scenario are averaged over ten replicates. comparison with univariate tests of significance the performance of tests of significance of ridge regression coefficients, in terms of true and false positive rates, was compared to the performance of univariate tests. comparisons were made in each of the simulation settings: the null simulation, and the simulations with continuous and with binary outcomes. the results are shown in additional file tables s s and s as would be expected when using a penalised regression approach such as ridge regression, the performance of the corresponding significance test depends on the degree of shrinkage, which is controlled by the shrinkage parameter λ. the performance of the approximate test is comparable to that of a permutation test, with the advantage of a much-reduced computational burden. further, ridge regression has the advantage over univariate tests of significance that it results in a much lower false positive rate. these advantages of ridge regression compared to univariate methods are further illustrated in the study of lung cancer data which follows. lung cancer data genome-wide association studies have identified snps associated with lung cancer disease status. snps have been identified at chromosomal locations 15q, 5p and 6p. the associations at 15q and 5p have been replicated in white populations, but the association at 6p has not. not all studies successfully replicated the associations at 15q and 5p. here, we use ridge regression to re-evaluate a set of snps for association with lung cancer disease status. whilst these data are not as high-dimensional as those from a genome-wide study, they allow us to illustrate the features of using ridge regression for genetic data. we show that ridge regression is a useful technique when data are correlated, and illustrate that multivariate methods have advantages over univariate tests of significance. data consist of genotypes and non-genetic predictors from approximately individuals in the european prospective investigation into cancer and nutrition. missing genotypes were imputed using mean imputation. gender, smoking status and age were included as unconstrained parameters in the model. for the purpose of comparison, univariate p-values were calculated. univariate p-values were calculated by fitting a logistic regression model for each snp independently, with gender, smoking status and age included in the model. a wald test was then used to estimate the significance of the coefficient. in this example, m < n and multivariate regression does give rise to unique parameter estimates. multivariate p-values are also reported here. in a data set of this size, it is computationally feasible to fit a ridge regression model with a range of values of λ. the ridge trace is a plot of parameter estimates against λ. the ridge trace typically suggests a range of values of λ rather than a single best value. in figure λ = seems to be the point at which the ridge estimates stabilise, and table reports the corresponding approximate p-values. of the snps in the lung cancer data set, four have been previously reported to be associated with lung cancer disease status. for these four snps, this table reports univariate, multivariate, approximate and permutation test p-values. further, with increasing shrinkage and stabilisation of the estimates, the approximate p-values for these snps become significant. rs which is not in ld with other snps in the data, is significant in a multiple logistic regression model even when no penalty is included in the model. a further snp, rs in chromosome has a nominally significant p-value in the approximate test, but again this snp has not previously been shown to be associated with lung cancer risk. the p-value trace using permutation test p-values shows good agreement with the approximate p-value trace. figure 2b is much more computationally expensive to produce than figure 2a, thus the use of the approximate test makes the plotting of a p-value trace for a range of values of λ more feasible. the results in table demonstrate the advantage of multivariate tests, and specifically of ridge regression, over univariate tests of significance. snp rs which has previously been shown to be associated with lung cancer disease status, was not found to be significant using the univariate test but was found as such by the multivariate methods. the two snps that are correlated, rs and rs were not significant in multiple regression but were significant in ridge regression, demonstrating the advantage of using ridge regression when snps are correlated. these results demonstrate that this approximate test of significance for coefficients fitted using logistic ridge regression reproduces previously ascertained associations, at reduced computational cost compared to a permutation test, even when snps are highly correlated.  we present and evaluate the performance of a test of significance for coefficients estimated using ridge regression. we evaluate the test as applied to both linear and logistic ridge regression models. our evaluation is by means of simulation studies across a range of scenarios representative of genetic data. we evaluate the test by comparing its performance to that of a permutation test. we evaluate the performance of the test when it is applied to a real data set. the data set comprises lifestyle data and genotypes together with lung cancer case-control status. using the proposed test, we successfully replicate previously found associations at much reduced computational cost compared to a permutation test. this demonstrates the utility of the test for detecting significant variables when predictor variables are highly correlated, as were two significant snps in the lung cancer data. stability selection is a method for variable selection that has received attention in the literature in recent years. it combines subsampling of the data with a dimension reduction technique, with the aim of finding consistently significant variables. penalised regression methods that perform variable selection, such as the lasso and the elastic net, have been used in stability selection. to use ridge regression with stability selection, a way of determining which variables to select in each subsample of the data is required. a permutation test would be computationally expensive in this context, rapidly becoming infeasible if for large n, if the guidelines of subsamples of size n/ given by meinshausen & bühlmann were followed. the test of significance proposed here, being much less computationally expensive than a permutation test, makes the combination of ridge regression and stability selection a feasible possibility. the choice of shrinkage parameter in ridge regression is discussed in the literature, but no consensus method provides an universally optimum choice. the proposed test performs well over a range of values of λ. we introduce the p-value trace, a plot of the change in the negative logarithm of the p-value of the regression coefficients as the shrinkage parameter increases from zero. this trace can be used in combination with the ridge trace of to visualise how the significance of the regression parameters, as well as their value, changes with increasing shrinkage. such a plot could aid the identification of significant coefficients in the ridge model. authors' contributions ec ran the simulations, analysed the data, and wrote the paper. pv was involved in the collection of the lung cancer data, supervised the analysis, and read and commented on the draft. mdi came up with the idea, supervised the work, and commented on drafts of the paper. all authors read and approved the final manuscript. supplementary material additional file 1 table s - performance comparison in null simulation using t-type test. table s - performance comparison in simulated data with continuous outcomes using t-type test. table s - performance comparison in null simulation including comparison to univariate tests of significance. table s - performance comparison with continuous outcomes including comparison to univariate tests of significance. table s - performance comparison with binary outcomes including comparison to univariate tests of significance. click here for file acknowledgements we thank the european prospective investigation into cancer and nutrition: elio riboli, anne tjønneland, kim overvad, francoise clavel-chapelon, rudolf kaaks, heiner boeing, antonia trichopoulou, franco berrino, domenico palli, salvatore panico, rosario tumino, paolo vineis, bas buneo de mesquita, petra peeters, eiliv lund, carlos a. gonzález, aurelio barricarte, miren dorronsoro, maria josé sánchez, carmen navarro, j. ramón quirós, jonas manjer, goran hallmans, nicholas e. day, kay-tee khaw, nicholas wareham and timothy j. key.
it is therefore desirable to obtain an estimate of the significance of each ridge regression coefficient. we develop and evaluate a test of significance for ridge regression coefficients. we introduce the p-value trace, a plot of the negative logarithm of the p-values of ridge regression coefficients with increasing shrinkage parameter, which enables the visualisation of the change in p-value of the regression coefficients with increasing penalisation. the proposed test is a useful alternative to a permutation test for the estimation of the significance of ridge regression coefficients, at a much-reduced computational cost. the p-value trace is an informative graphical tool for evaluating the results of a test of significance of ridge regression coefficients as the shrinkage parameter increases, and the proposed test makes its production computationally feasible. using simulation studies, we demonstrate that the performance of the test is comparable to that of a permutation test, with the advantage of a much-reduced computational cost. we apply the proposed method to a lung cancer case-control data set from epic, the european prospective investigation into cancer and nutrition. however, snps typed using these methods can be highly correlated due to linkage disequilibrium among them, and standard multiple regression techniques fail with these data sets due to their high dimensionality and correlation structure. technological developments have increased the feasibility of large scale genetic association studies. ridge regression is one such penalised regression technique which does not perform variable selection, instead estimating a regression coefficient for each predictor variable.
 technological developments have increased the feasibility of large scale genetic association studies. densely typed genetic markers are obtained using snp arrays, next-generation sequencing technologies and imputation. however, snps typed using these methods can be highly correlated due to linkage disequilibrium among them, and standard multiple regression techniques fail with these data sets due to their high dimensionality and correlation structure. there has been increasing interest in using penalised regression in the analysis of high dimensional data. ridge regression is one such penalised regression technique which does not perform variable selection, instead estimating a regression coefficient for each predictor variable. it is therefore desirable to obtain an estimate of the significance of each ridge regression coefficient.  we develop and evaluate a test of significance for ridge regression coefficients. using simulation studies, we demonstrate that the performance of the test is comparable to that of a permutation test, with the advantage of a much-reduced computational cost. we introduce the p-value trace, a plot of the negative logarithm of the p-values of ridge regression coefficients with increasing shrinkage parameter, which enables the visualisation of the change in p-value of the regression coefficients with increasing penalisation. we apply the proposed method to a lung cancer case-control data set from epic, the european prospective investigation into cancer and nutrition.  the proposed test is a useful alternative to a permutation test for the estimation of the significance of ridge regression coefficients, at a much-reduced computational cost. the p-value trace is an informative graphical tool for evaluating the results of a test of significance of ridge regression coefficients as the shrinkage parameter increases, and the proposed test makes its production computationally feasible. 