Issues surrounding data security and privacy of personal information are of growing concern to the public, governments and commercial sectors. Privacy concerns can erode confidence in the ability of organisations to store data securely, with a consequence that individuals may be reticent to contribute their personal information to scientific studies or to commercial organisations (Sanderson et al., 2015, Naveed et al., 2015, Gymrek et al., 2013).
In this paper we demonstrate that statistical regression methods can be applied directly to encrypted data without compromising privacy. Our work involves adapting existing methodology for regression analysis in such a way as to enable full computation within the mathematical and computational constraints of recently developed fully homomorphic encryption schemes(Gentry, 2010, Fan and Vercauteren, 2012). We show empirically that traditional state-of-the-art convergence acceleration techniques can under-perform when such constraints are taken into account.
Fully homomorphic encryption (FHE) differs from differential privacy (DP) in that it provides for exact computation with cryptographically strong privacy of all the data during the fitting process itself, at the expense of much greater restrictions on the operations and computational cost (Gentry, 2010, Aslett et al., 2015a). However, it is complementary to DP (as opposed to competing with it) in the sense that FHE maintains privacy during model fitting and prediction, while DP can ensure privacy post-processing of the data if the model itself is to be made public. For an overview of DP see Dwork and Roth (2014).
FHE allows for secure operations to be performed on data, and statistical analysis and machine learning is the major reason why people want to perform mathematical operations on data. Thus, there is a real opportunity for machine learning scientists to be involved in shaping the research agenda in FHE (United States EOP, 2014). The applications of encrypted statistics and machine learning include general purpose cloud computing when privacy concerns exist, and are especially important in e-health and clinical decision support (Basilakis et al., 2015, McLaren et al., 2016).
§2 is a brief accessible introduction to FHE and §3 recaps regression to fix notation and our method of representing data prior to encryption. A detailed examination of coordinate and gradient descent methods in an encrypted context follows in §4, including encrypted scaling for correctness, computational considerations, prediction, inference, regularisation and theoretical proofs for parameters in a popular FHE scheme. §5 discusses acceleration methods optimal for encrypted computation, with examples provided in §6 and discussion in §7.
2 Homomorphic encryption
A common technique for ensuring the privacy of data is to encrypt it, but typically once one wishes to fit a model it is necessary to first decrypt and risk exposing the data. However, recent advances in cryptography enable a very limited amount of computation to be performed directly on the encrypted content, rendering the correct result upon decryption.
2.1 Standard cryptography
A public key encryption scheme is one which has two algorithms, and , to perform encryption and decryption respectively, together with two keys: the public key, , can be widely distributed and used by anyone to encrypt a message; the secret key, , is required to decrypt any message encrypted using and so is kept private. The fundamental identity is:
The data to be encrypted, , is referred to as the message or plaintext and after encryption is referred to as the ciphertext. In conventional encryption algorithms manipulation of the ciphertext does not typically lead to meaningful modification of the message.
2.2 Fully Homomorphic encryption
An encryption scheme is said to be fully homomorphic if it also possesses two operations, and , which satisfy the following identities:
for all , which can be applied a theoretically arbitrary number of times. In other words, a homomorphic encryption scheme allows computation directly on ciphertexts which will correctly decrypt the result as if the corresponding operations had been applied to the original messages.
However there are many constraints in practical implementation, reviewed in Aslett et al. (2015a). For the purposes of this work they may be synopsised as:
Typically can only easily represent binary or integers.
Data size grows substantially after encryption.
The computational cost of and is orders of magnitude higher than standard and .
Operations such as comparisons () and division () are not possible.
Implementation of current schemes necessitate a highly computationally expensive ‘bootstrap’ operation (unrelated to statistical bootstrap) which must be applied frequently between operations to control the noise in the ciphertext.
Consequently, a crucial property of any FHE scheme is the multiplicative depth111Having only addition and multiplication operations mean all computations form polynomials. Simply put, the multiplicative depth corresponds to the degree of the maximal degree term of the polynomial minus one. e.g., has depth 1 and has depth . which is possible before a ‘bootstrap’ is required. Typically, one can select the parameters of the FHE scheme to support a pre-specified multiplicative depth, but this is a trade-off because parameters supporting greater depth between ‘bootstraps’ also result in larger ciphertexts and slower homomorphic operations (). Therefore, it is essential to consider the Maximum Multiplicative Depth (MMD) required to evaluate an algorithm encrypted, since this dramatically affects speed and memory usage. Indeed, cryptographic parameters are typically chosen to match the MMD of the target algorithm being run so as to avoid ‘bootstrap’ altogether (which is simply deemed too computationally costly).
The above constraints often mean that standard statistical methodology cannot be applied unmodified on encrypted content.
2.3 Privacy preserving statistics
There has been some work towards using FHE schemes for statistics and machine learning. Often this has involved identifying existing algorithms which can be run with minimal modifications (Wu and Haven, 2012, Graepel et al., 2013), or fitted on unencrypted data for later prediction on new encrypted data (Dowlin et al., 2016). However, some recent work has also begun on developing new methodology inspired by traditional techniques, specifically tailored to homomorphic computation so that the whole analysis — model fitting and prediction — can be computed encrypted (Aslett et al., 2015b).
In particular, the topic of linear regression in the context of FHE has not been covered systematically in the literature thus far.Hall et al. (2011) propose protocols for regression analysis which involve substantial communication and intermediate decryption between multiple parties: it takes two days to complete a problem of size observations and predictors. In this work, we want to develop methods capable of fitting and prediction without any intermediate communication or decryption. Wu and Haven (2012) were the first to tackle linear regression in the context of FHE by using Cramer’s rule for matrix inversion. Unfortunately, this approach calls for the computation of the determinant of , which quickly becomes intractable for even low dimensional problems (e.g., ). Specifically, the multiplicative depth is unbounded with growing dimension and, consequently, bootstrapping seems unavoidable. Graepel et al. (2013) mention that regression can be done by gradient descent, but do not implement the method or give further details.
The approach we present in this work enjoys several notable properties: (i) estimation and prediction can both be performed in the encrypted domain; (ii) bootstrapping can be avoided even for moderately large problems, (iii) scales linearly with the number of predictors; (iv) and admits the option of(ridge) regularisation.
3 The linear regression model
The standard linear regression model assumes
is a response vector of length; is a design matrix of size ; is a parameter vector of length ; and is a vector of length . Provided that
is invertible, the ordinary least squares (OLS) solution to
Regularisation techniques trade the unbiasedness of OLS for smaller parameter variance by adding a constraint of the form , for some (Hastie et al., 2009). Bounding the norm of the parameter vector imposes a penalty for complexity, resulting in shrinkage of the regression coefficients towards zero. We focus on regularisation (ridge, henceforth RLS; Hoerl and Kennard, 1970), where . Other options are available, although no one method seems to dominate the others (Zou and Hastie, 2005, Tibshirani, 1996, Fu, 1998). The standard solution to the regularised problem
revealing that is also key in converting ill-conditioned problems into well-conditioned ones (Hoerl and Kennard, 1970).
3.1 Data representation and encoding
Because FHE schemes can only naturally encrypt representations of integers, dealing with non-integer data requires a special encoding. We use the transformation , for and . Here, represents the desired level of accuracy or, more precisely, the number of decimal places to be retained. By construction, , so that smooth relative distances between elements encoded in this way are approximately preserved. The encoding can accommodate both discrete (integer-valued) data and an approximation to continuous (real-valued) data. Throughout, we assume that covariates are standardised and responses centred before integer encoding and encryption.
4 Least squares via descent methods
Conventional solutions to (3) and (5) can be found in closed form using standard matrix inversion techniques with runtime complexity for the matrix inversion and for the matrix product. This direct approach was the one advocated for encrypted regression by Wu and Haven (2012), though dimensionality is constrained by school-book matrix inversion to enable homomorphic computation.
4.1 Iterative methods
We analyse two variants of an iterative algorithm, one which updates all parameters in simultaneously at each iteration , using the vector (ELS-GD); and another which updates these parameters sequentially, using always the most current estimate for each parameter (ELS-CD).
The sequential update mode is related to the Gauss–Seidel method (coordinate descent), while the simultaneous update mode is related to the Jacobi method (gradient descent). See Varga (2000, chapter 3); Björck (1996, chapter 7).
As we will see below, there are two competing concerns here: optimisation efficiency and computational tractability within the cryptographic constraints, with these two types of updates having different properties. We demonstrate that the properties of these methods in the encrypted domain means some standard optimality results no longer apply.
4.1.1 Sequential updates via coordinate descent
Standard coordinate descent (see e.g. Wright, 2015) for linear regression has the update:
where contains the components updated on the previous iteration or this one as appropriate, there being many variants for the schedule of coordinates to update. However, this cannot be computed encrypted because the required data-dependent division is not feasible. In an alternative variant we can replace this with a generic step-size and at each iteration choose one variable, say , and update the corresponding regression parameter:
Two things are noteworthy: first, these equations require only evaluation of polynomials; and second, a rescaling is necessary since only integer polynomial functions can be computed homomorphically, and the parameter is usually not an integer. We show how to perform this type of rescaling in the context of §4.1.2.
As each update uses the most recent estimate of every parameter, the MMD grows by with each parameter update (due to the term ). This implies that for iterations over covariates the MMD is equal to . This renders the algorithm very expensive computationally, as it requires bootstrapping of ciphertexts for problems with even moderately large .
4.1.2 Simultaneous updates via gradient descent
In the case of gradient descent, for the objective function with , the update equations are:
As with coordinate descent, these updates can be computed homomorphically, although the update equations must again be rescaled. This is strictly necessary to accommodate the transformed data, in order to overcome the fact that we cannot divide (see §3.1). Letting for , the transformed equations for the simultaneous updates are:
where now all transformed variables are represented with tildes, e.g., and similarly for the other variables, except the coefficients as their scaling is iteration dependent (see supplementary materials, §1). The rescaling factors are independent of the data and known a priori, and so can be grouped during computation, e.g., in (10) can be encrypted as a single value.
Retrieval of the coefficients can be done by the secret key holder by computing . Note the important difference between coordinate and gradient descent: for iterations, in CD each coefficient is updated times while in GD each coefficient is updated times.
There is a crucial difference between CD and GD in the encrypted domain: GD reduces the multiplicative depth from to , which is independent of , enabling scalability to higher dimensional models without bootstrapping or having to select parameters which support greater MMD. As discussed in §2.2, bootstrapping in current FHE schemes is to be avoided wherever possible because it is very computationally expensive. This is an interesting and important result, because it means that in the specific setting of encrypted computation, Nesterov’s faster rates of convergence (Nesterov, 2012) compared to gradient descent in a randomised coordinate descent setting will not apply, as we will show.
For these computational reasons we focus primarily on gradient descent hereinafter. For convergence in a regression setting recall:
Lemma 1 (Convergence of Els-Gd).
Define and let . Then, for where is the spectral radius of .
The optimal choice of step size, in the sense that it minimises the spectral radius, is , implying an optimal spectral radius , where and
denote the largest and smallest eigenvalues of, respectively (see Ryaben’kii and Tsynkov, 2006, Theorem 6.3 for all proofs).
Lemma 2 (Oscillatory nature of Els-Gd).
The iterative process (8) can be written as an oscillating sum:
This lemma is proved in the supplementary materials (§3). We show in §5.2 that it is possible to improve the convergence rate by using acceleration methods that exploit the oscillatory nature of the GD algorithm to accelerate the convergence of the series .
Note that the form of the GD equation (10) implies a common scaling factor, , for all parameters. Performing encrypted prediction is then straightforward as it requires only the computation of the dot product , where denotes the predicted values in the space of the original data and denotes the corresponding transformed version. Upon decryption, rescaling can be done as before by the secret key holder. The procedure increases the MMD of the algorithm by 1.
The situation is more complex for coordinate descent since at the end of the final iteration each element of will have different scaling. Therefore the scaling must be unified before prediction, adding additional overhead.
However, the matrix inversion is intractable under homomorphic computation except for very small . An alternative is to estimate the standard errors by bootstrapping the data and using the variability in the parameter estimates obtained.
-regularised (ridge) least squares is easy to implement using the well known data augmentation procedure (Allen, 1974):
OLS estimates when using the augmented data, , are equivalent to RLS estimates when using the original data, ; that is,
Because the augmentation terms and are independent of the data and , the iterative methods already developed for least squares (§4.1) can be used with the augmented data. Also note that the maximal eigenvalue is easily updated, and so a new step size can be chosen without additional computation.
4.5 Theoretical parameter requirements for the Fan and Vercauteren scheme
We provide results to guide the choice of cryptographic parameters for the encryption scheme of Fan and Vercauteren (2012) — hereinafter FV. This is implemented in the HomomorphicEncryption R package (Aslett et al., 2015a) and used in the examples.
FV represents data as a polynomial of fixed degree with coefficients in a finite integer ring.222The maximal degree and maximal ring element are tunable parameters, but even small increases make ciphertexts bigger and homomorphic operations slower. For example, is represented by where are the binary decomposition of , such that . Addition and multiplication operations on the ciphertext result in polynomial addition and multiplication on this representation.
The transformed regression coefficients grow substantially during computation, so that we must ensure (i) the maximal degree of the FV message polynomial is large enough to decrypt ; and (ii) that the coefficient ring is large enough to accommodate the worst case growth in coefficients.
Lemma 3 (FV parameter requirements for GD).
If data is represented in binary decomposed polynomial form, then after running the ELS-GD algorithm the degree and coefficient value of the encrypted regression coefficients is bound by:
where and ;
This lemma is proved in the supplementary materials (§2) and provides lower bounds on the choice of parameters and in the FV scheme.
Recall the MMD for GD is . Theoretical bounds on security level (Lindner and Peikert, 2011) and multiplicative depth (Lepoint and Naehrig, 2014) — together with the polynomial bounds and algorithmic MMD proved here — then enable full selection of encryption parameters to guarantee security and computational correctness of the encrypted GD algorithm.
Although ELS-GD is guaranteed to converge to the OLS solution, the rate of convergence can be slow, for instance when predictors are highly correlated. Here we analyse some classic acceleration methods.
A preconditioning matrix is often used to accelerate convergence of iterative methods (Björck, 1996, chapter 7) by solving the preconditioned problem
with the same solution as the original problem, but having more favourable spectral properties. A simple preconditioning is diagonal scaling of the columns of . Let be a diagonal matrix with diagonal entries where . Since for all as a result of standardisation (see §3.1) the preconditioning matrix becomes . The preconditioned update equation is then
5.2 van Wijngaarden transformation
The van Wijngaarden (1965) transformation (VWT) is a variant of the Euler transform for accelerating the convergence of an alternating series. Given the partial sums of an alternating series, , we can compute averages of these partial sums,
(and averages of the averages) to form a matrix of averaged partial sums. For a finite number of terms (), averages of partial sums are often closer to the limiting value, , than any single partial sum of the original series, (), so that the alternating nature of the sequence is averaged out, damping its oscillatory behaviour and speeding-up convergence.
As shown in (11), the values computed with ELS-GD form an alternating series, making the VWT a candidate for accelerating the convergence to the true regression coefficients.
The implementation of (17) has a simple, closed form solution. For a stopping column (as van Wijngaarden suggests) we compute the average partial sum , and take this as our best approximation to the value of the series at convergence. Notably, this can be implemented homomorphically with the exception of the division, but since the factor is independent of the data, we can compute instead and incorporate the appropriate correction upon decryption; i.e., once ELS-GD is completed, the final VWT estimate is:
The computational cost of the VWT is minimal, involving approximately additions and multiplications, and increasing the MMD by only 1.
5.3 Nesterov’s accelerated gradient
Nesterov’s accelerated gradient (NAG) achieves a convergence rate of as opposed to the achieved by regular GD (Nesterov, 1983). The NAG algorithm can be written as follows:
These equations must also be rescaled for homomorphic computation (similarly to GD in §4.1.2):
where is the transformed vector of responses, according to §3.1, and similarly for the remaining variables, except which have iteration dependent scaling factors.
All scaling constants are independent of the data and known a priori, and so can be incorporated into the scaling by the secret key holder to obtain the final parameter estimates as: . Because of the extra acceleration step, ELS-NAG has a MMD equal to (see Table 1). This is particularly interesting, because although Nesterov’s method is state-of-the-art for unencrypted GD, the increase in MMD makes it costly for encrypted analysis.
|Preconditioned gradient descent||2K|
|van Wijngaarden transformation||2K+1|
|Nesterov’s accelerated gradient||3K|
In this section we empirically analyse the methods proposed for encrypted linear regression using simulated and real data (see supplementary materials, §4, for details). We use the implementation of the Fan and Vercauteren (2012) cryptosystem provided by Aslett, Esperança, and Holmes (2015a). The runtimes reported for encrypted analysis are on a 48-core server. All error norms are root mean squared deviations w.r.t. OLS, and we use throughout.
For simulations with independent data we generate , and . For simulations with correlated data we use Normal copulas and generate predictors whose pairwise correlations are all equal to .
Figure 4 (left) illustrates the computational properties of the coordinate and gradient descent methods for encrypted regression for a fixed MMD. Since the MMD supported in the encryption scheme is a prime determinant of the computational cost of homomorphic operations, this serves as a proxy for the error as a function of encrypted computational complexity. ELS-GD clearly outperforms ELS-CD for a fixed encrypted computational cost, as expected from §5. Furthermore, Figure 4, (right) shows the VWT provides additional acceleration in convergence relative to GD.
In general, higher correlation among predictors implies less favourable spectral properties for , which in turn makes convergence slower for both ELS-GD-VWT and ELS-NAG (Figure 4). A fair comparison must control for the fact that the two algorithms have different encrypted computational complexities. Using the MMD as a proxy again, ELS-GD-VWT typically outperforms ELS-NAG for fixed level of complexity (Figure 4). In very high correlation settings this relationship can be reversed, but only for large numbers of iterations, which it is desirable to avoid.
We stress that this choice is conditional on the encrypted computational framework considered here. It is particularly interesting that when working unencrypted, NAG is the state-of-the-art; but in the restricted framework of FHE, VWT empirically appears to be a better choice.
Convergence is affected by the number of predictors. In particular, the maximum number of iterations required to reduce the norm of the initial error vector by a factor (reciprocal of the average convergence rate; Varga, 2000, p.69) gives us an idea of the relationship between number of predictors and speed of convergence. For any level of correlation, this measure of complexity increases linearly with (see Figure 1 in the supplementary materials).
Finally, the computational costs of ELS-GD are given in Figure 5. Runtime grows quickly with the algorithm’s multiplicative depth, which increases with the number of iterations. However, for a fixed multiplicative depth the runtime is roughly linear in both and . Memory requirements grow in a similar fashion.
Mood stability data
The first application is to mood stability in bipolar patients (Bonsall et al., 2012). Of interest in this application is the characterisation of the stochastic process governing the resulting time series, pre and post treatment, which we model as an autoregressive process of order two (). Convergence is achieved within 2 iterations (; Figure 8). The algorithm runs encrypted in 12 seconds and requires under 15 MBs of memory, excluding overheads.
, effective degrees of freedom. 
Prostate cancer data
The second application is to prostate cancer (Stamey et al., 1989). The model here is a standard linear regression (). Although not all parameters have completely converged by iteration 4 with unregularised ELS-GD-VWT (; Figure 8), the predictions are close to those produced by RLS (Figure 8). The algorithm runs encrypted in 30 minutes and requires 3.5 GBs of memory ().
For runtimes and memory requirements in these applications see Figure 2 in the supplementary materials.
We demonstrated that in the restricted framework of FHE, traditional state-of-the-art methods can perform poorly. Statistical and computational methods tailored for homomorphic computation are therefore required, which may differ from the state-of-the-art in an unrestricted framework.
For optimal convergence speed, the step size can be provided by the data holder, who can use the inequality to approximate to arbitrary precision, since as .
Choosing the penalty is less straightforward. Traditional methods involve cross-validation which is impossible under strict FHE. Alternatively, it is possible to do rounds of communication and decryption between two parties to achieve this, in which case Differential Privacy can used as a way to guarantee security during the intermediate communication steps.
P.M. Esperança: LSI-DTC doctoral studentship (EPSRC grant EP/F500394/1). L.J.M. Aslett and C.C. Holmes: i-like project (EPSRC grant EP/K014463/1).
- Allen (1974) D. Allen. The relationship between variable selection and data agumentation and a method for prediction. Technometrics, 16(1):125–127, 1974.
- Aslett et al. (2015a) L. J. M. Aslett, P. M. Esperança, and C. C. Holmes. A review of homomorphic encryption and software tools for encrypted statistical machine learning. arXiv:1508.06574: arxiv.org/abs/1508.06574, 2015a.
- Aslett et al. (2015b) L. J. M. Aslett, P. M. Esperança, and C. C. Holmes. Encrypted statistical machine learning: new privacy preserving methods. arXiv:1508.06845: arxiv.org/abs/1508.06845, 2015b.
- Basilakis et al. (2015) J. Basilakis, B. Javadi, and A. Maeder. The potential for machine learning analysis over encrypted data in cloud-based clinical decision support — background and review. In Health Informatics and Knowledge Management (HIKM’15), volume 164 of Conferences in Research and Practice in Information Technology, pages 3–13. ACS, 2015.
- Björck (1996) A. Björck. Numerical methods for least squares problems. SIAM, 1996.
- Bonsall et al. (2012) M. Bonsall, S. Wallace-Hadrill, J. Geddes, G. Goodwin, and E. Holmes. Nonlinear time-series approaches in characterizing mood stability and mood instability in bipolar disorder. Proceedings of the Royal Society B (Biological Sciences), 279(1730):916–924, 2012.
Dowlin et al. (2016)
N. Dowlin, R. Gilad-Bachrach, K. Laine, K. Lauter, M. Naehrig, and
Cryptonets: Applying neural networks to encrypted data with high throughput and accuracy.Technical Report MSR-TR-2016-3.: www.microsoft.com/en-us/research/publication/cryptonets-applying-neural-networks-to-encrypted-data-with-high-throughput-and-accuracy/, 2016.
- Dwork and Roth (2014) C. Dwork and A. Roth. The Algorithmic Foundations of Differential Privacy. Foundations and Trends in Theoretical Computer Science. Now Publisher, 2014.
- Fan and Vercauteren (2012) J. Fan and F. Vercauteren. Somewhat practical fully homomorphic encryption. IACR Cryptology ePrint Archive, Report 2012/144: eprint.iacr.org/2012/144, 2012.
- Fu (1998) W. Fu. Penalized regressions: The bridge versus the lasso. Journal of Computational and Graphical Statistics, 7(3):397–416, 1998.
- Gentry (2010) C. Gentry. Computing arbitrary functions of encrypted data. Communications of the ACM, 53(3):97–105, 2010.
- Graepel et al. (2013) T. Graepel, K. Lauter, and M. Naehrig. ML Confidential: Machine learning on encrypted data. In Information Security and Cryptology (ICISC’12), volume 7839 of Lecture Notes in Computer Science, pages 1–21. Springer, 2013.
- Gymrek et al. (2013) M. Gymrek, A. L. McGuire, D. Golan, E. Halperin, and Y. Erlich. Identifying personal genomes by surname inference. Science, 339(6117):321–324, 2013.
- Hall et al. (2011) R. Hall, S. E. Fienberg, and Y. Nardi. Secure multiple linear regression based on homomorphic encryption. Journal of Official Statistics, 27(4):669–691, 2011.
- Hastie et al. (2009) T. Hastie, R. Tibshirani, and J. Friedman. The elements of statistical learning: data mining, inference and prediction. Springer, 2009.
- Hoerl and Kennard (1970) A. Hoerl and R. Kennard. Ridge regression: Biased estimation for nonorthogonal problems. TheJournal, Technometrics(12):1, 1970.
- Lepoint and Naehrig (2014) T. Lepoint and M. Naehrig. A comparison of the homomorphic encryption schemes FV and YASHE. In Progress in Cryptology (AFRICACRYPT’14), volume 8469 of Lecture Notes in Computer Science, pages 318–335. Springer, 2014.
- Lindner and Peikert (2011) R. Lindner and C. Peikert. Better key sizes (and attacks) for LWE-based encryption. In Topics in Cryptology (CT-RSA’11), volume 6558 of Lecture Notes in Computer Science, pages 319–339. Springer, 2011.
- McLaren et al. (2016) P. J. McLaren, J. L. Raisaro, M. Aouri, M. Rotger, E. Ayday, I. Bartha, M. B. Delgado, Y. Vallet, H. F. Günthard, M. Cavassini, H. Furrer, T. Doco-Lecompte, C. Marzolini, P. Schmid, C. Di Benedetto, L. A. Decosterd, J. Fellay, J. Hubaux, A. Telenti, and the Swiss HIV Cohort Study. Privacy-preserving genomic testing in the clinic: a model using HIV treatment. Genetics in Medicine, 2016. Advance online publication.
- Naveed et al. (2015) M. Naveed, E. Ayday, E. W. Clayton, J. Fellay, C. A. Gunter, J.-P. Hubaux, B. A. Malin, and X. Wang. Privacy in the genomic era. ACM Computing Surveys, 48(1):Article 6, 2015.
- Nesterov (1983) Y. Nesterov. A method for solving a convex programming problem with convergence rate . Soviet Mathematics Doklady, 27(2):372–376, 1983.
- Nesterov (2012) Y. Nesterov. Efficiency of coordinate descent methods on huge-scale optimization problems. SIAM Journal on Optimization, 22(2):341–362, 2012.
- Ryaben’kii and Tsynkov (2006) V. S. Ryaben’kii and S. V. Tsynkov. A Theoretical Introduction to Numerical Analysis. Chapman and Hall/CRC, 2006.
- Sanderson et al. (2015) S. C. Sanderson, M. D. Linderman, S. A. Suckiel, G. A. Diaz, R. E. Zinberg, K. Ferryman, M. Wasserstein, A. Kasarskis, and E. E. Schadt. Motivations, concerns and preferences of personal genome sequencing research participants: baseline findings from the HealthSeq project. European Journal of Human Genetics, 24:14–20, 2015.
- Stamey et al. (1989) T. Stamey, J. Kabalin, J. McNeal, I. Johnstone, F. Freiha, E. Redwine, and N. Yang. Prostate specific antigen in the diagnosis and treatment of adenocarcinoma of the prostate. II. radical prostatectomy treated patients. Journal of Urology, 141(5):1076–1083, 1989.
- Tibshirani (1996) R. Tibshirani. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society. Series B (Methodological), 58(1):267–288, 1996.
- United States Executive Office of the President (2014) United States Executive Office of the President. Big data: Seizing opportunities, preserving values, 2014.
- van Wijngaarden (1965) A. van Wijngaarden. In Cursus: Wetenschappelijk Rekenen B: Process Analyse, pages 51–60. Stichting Mathematisch Centrum, 1965.
- Varga (2000) R. S. Varga. Matrix iterative analysis, volume 27 of Springer Series in Computational Mathematics. Springer, 2nd edition, 2000.
- Wright (2015) S. Wright. Coordinate descent algorithms. Mathematical Programming, 151(1):3–34, 2015.
- Wu and Haven (2012) D. Wu and J. Haven. Using homomorphic encryption for large scale statistical analysis. Technical Report: cs.stanford.edu/people/dwu4/papers/FHE-SI_Report.pdf., 2012.
- Zou and Hastie (2005) H. Zou and T. Hastie. Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society. Series B (Methodological), 67(2):301–320, 2005.