## ARTICLE

## Differential Privacy for Sparse Classification Learning

Puyu Wang a , Hai Zhang a,b a School of Mathematics, Northwest University, Xi'an; b Faculty of Information Technology &amp; State Key Laboratory of Quality Research in Chinese Medicines, Macau University of Science and Technology, Macau

## ARTICLE HISTORY

Compiled August 5, 2019

## ABSTRACT

In this paper, we present a differential privacy version of convex and nonconvex sparse classification approach. Based on alternating direction method of multiplier (ADMM) algorithm, we transform the solving of sparse problem into the multistep iteration process. Then we add exponential noise to stable steps to achieve privacy protection. By the property of the post-processing holding of differential privacy, the proposed approach satisfies the glyph[epsilon1] -differential privacy even when the original problem is unstable. Furthermore, we present the theoretical privacy bound of the differential privacy classification algorithm. Specifically, the privacy bound of our algorithm is controlled by the algorithm iteration number, the privacy parameter, the parameter of loss function, ADMM pre-selected parameter, and the data size. Finally we apply our framework to logistic regression with L 1 regularizer and logistic regression with L 1 / 2 regularizer. Numerical studies demonstrate that our method is both effective and efficient which performs well in sensitive data analysis.

## KEYWORDS

Differential privacy; Classification; ADMM; Unstable; Sparse

## 1. Introduction

With the development of the internet of things, a large amount of data has been collected which contains the individuals sensitive information, such as medical purchase records, hospital electronic medical records, web-site search information, home address and contact information. Inappropriate use of sensitive information will cause a serious threat to personal privacy, which will bring the risk of personal privacy leakage. For example, an attacker can infer personal information or predict individual behavior by capturing personal sensitive behavior data easily. In order to avoid the leakage of privacy, scholars have proposed many different privacy-preserving approaches. Sweeney et al.[1, 2] proposed K-anonymous method, Machanavajjhala et al. [3] proposed l -diversity principle, Li et al.[4] proposed T-closeness and so on. These technologies have been studied extensively and applied to privacy protecting. However, the reliability of these privacy protection models is related to the background

This work was partially supported by National Natural Sciences China-Guangdong Joint Fund under grant numbers U1811461, National Natural Science Foundation of China under grant numbers 11571011. CONTACT Hai Zhang. Email: zhanghai@nwu.edu.cn

knowledge of the attacker. When the attacker has enough background knowledge, these methods may fail[5][6]. Dwork [7] proposed differential privacy in 2006. Unlike traditional privacy-preserving methods, differential privacy defines a strict attack model which is independent of background knowledge, and gives a quantitative representation of the degree of privacy leakage, which can protect individual privacy information effectively. In the seminal paper [7], Dwork gave the definition of differential privacy, which can be stated formally as follows.

Definition Given a random algorithm M Range , ( ) denote a collection of all possible · outputs of algorithm. For any two data sets X and X ′ that differing on at most one data point, and for any S ⊆ Range M ( ), if

$$\frac { P ( M ( X ) ) \in S ) } { P ( M ( X ^ { ^ { \prime } } ) \in S ) } \leq \exp ( \epsilon ),$$

then random algorithm M satisfies glyph[epsilon1] -differential privacy.

As a privacy-preserving technology with rigorous mathematical theory, differential privacy is suitable for privacy preserving in the era of big data, and has a broad application [8-11]. In the framework of differential privacy, the degree of privacy information leakage is controlled by the parameter glyph[epsilon1] . As the glyph[epsilon1] increases, the more privacy is revealed, which means that the probability of individual indentification increases. The advantages of differnetial privacy mainly include: (1) Differential privacy has good privacy protection ability even if attacker has background knowledge, because the attack model of differential privacy is independent of the background knowledge. (2) Differential privacy has the post-processing property, specifically, if data published satisfies a given amount glyph[epsilon1] -differential privacy, the privacy leakage will not increase with data analyzing or re-releasing. The above advantages make differential privacy have good scalability. Laplace mechanism[12] and exponential mechanism[13] are two common differential privacy mechanisms. Although differential privacy has a wide range of applications, most of applications focus on data publishing [14, 15] and query processing [16, 17].

In this paper, we focus on the regularized empirical risk minimization (ERM) problems for binary classification, which is defined as follows

$$\min _ { w } \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } O ( y _ { i } w ^ { T } x _ { i } ) + P _ { \lambda } ( w ),$$

where O ( ) · is loss function, x i ∈ R p , y i ∈ {-1 1 , } , w ∈ R p is parameter coefficient vector, P λ ( w ) is penalty function, and λ &gt; 0 is the tunning parameter which controls the complexity of the model. Specifically, we focus on the L q regularizer, which have the form P λ ( w ) = ‖ w ‖ q q , where 0 &lt; q ≤ 1. When q = 1, it is L 1 Lasso type[18], and when q = 1 2, it is / L 1 2 / regularizer[19]. Our goal is to train a sparse classifier w ∗ over the dataset while protect the individual privacy.

Differential privacy requires the output of the random algorithm being stable under small perturbations of the input data[20]. However, the most sparse algorithms are unstable. Xu[21] studied the relationship between sparsity and stability. Their results showed that a sparse algorithm can not be stable and vice versa. Thus, the unstable property of sparse algorithm makes it hard to design a differential privacy sparse

algorithm. Under some assumptions, there are some works on the convex sparse regularization approaches. For example, Kifer et al.[22] gave the first results on private sparse regression in high dimension. They designed a computationally efficient algorithm, implicitly based on subsampling stability, for support recovery using the LASSO estimator. Smith et al.[23] extended and improved on the results of [22], based on the sufficient conditions for the LASSO estimator to be robust to small changes in the data set. They proposed the privacy algorithms for sparse linear regression. [24] designed a nearly optimal differentially private version of Lasso. As compared to the previous work, they assumed that the input data has bounded l ∞ norm. In this paper, we focus on differential privacy sparse classification algorithms. Chaudhuri et al.[25] proposed the regularized classification method with differential privacy. Their results hold for loss functions and regularizers satisfying certain differentiability and convexity conditions, specifically, when regularizer is l 2 norm. Zhang and zhu[26] proposed differential privacy regularized classification methods for distributed stored data. Their mechanisms lead to algorithms that can provide privacy guarantees under mild conditions. Zhang et al.[27] extended and improved the results of [26], and based on modified ADMM proposed a perturbation method to improve privacy without compromising accuracy. Wang and Zhang[28] proposed distributed differential privacy logistic regression based on three step ADMM by noising the output of each iteration, which can protect the local privacy. However, the above methods both require the regularizer being stronly-convex and differentiable.

In this paper, we design a differential privacy version of convex and nonconvex sparse classification approach with mild conditions on the regularizers. To solve the optimization classification problem, we transform the solving process into three-step sub-problems that can be easily solved by iteration process. And only at the second step, the algorithm access the raw data. Interestingly, the algorithm in this step is stable which fits to the tramework of differential privacy. So we design a privacy preserving algorithm. And by the property of the post-processing holding of differential privacy, the proposed approach satisfies the glyph[epsilon1] -differential privacy even when the original problem is unstable. Then, we present the theoretical privacy bound of the classification algorithm. At last, we apply our framework to logistic regression with L 1 regularizer and logistic regression with L 1 2 / regularizer. Numerical studies demonstrate that our method is both effective and efficient which performs well in sensitive data analysis.

## 2. Differential privacy Sparse Classification Framework

In this section, we present the differential privacy sparse classification framework based on ADMM algorithm, then we show that under mild assumptions for loss function and penalty function, our method satisfies glyph[epsilon1] -differential privacy.

## 2.1. ADMM algorithm

ADMM was proposed in the early 1970s, and has since been studied extensively. In recent years, it has been used in many areas such as computer vision, signal processing and networking[29]. Now, we use ADMM to transform the solving of sparse problem into the multistep iteration process. First, consider adding the auxiliary variable Z ∈ R p , then model (1) can be rewritten as

$$\min _ { w, z } \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } O ( y _ { i } w ^ { T } x _ { i } ) + P _ { \lambda } ( Z ),$$

$$s. t. \, w = Z.$$

It's obvious that model (1) and model (2) are equivalent. Then we have the quadratically augmented Lagrangian function

$$L ( Z, w, V ) = \frac { 1 } { n } \sum _ { i = 1 } ^ { n } O ( y _ { i } w ^ { T } x _ { i } ) + P _ { \lambda } ( Z ) + \frac { c } { 2 } ( \| Z - w + \frac { V } { c } \| _ { 2 } ^ { 2 } - \| \frac { V } { c } \| _ { 2 } ^ { 2 } ),$$

where V is Lagrange multipliers, and c &gt; 0 is a preselected penalty coefficient. ADMM algorithm entails three steps per iteration

Step 1: Z updates,

$$Z ( k + 1 ) & = \underset { Z } { \arg \min } \ L ( Z, w ( k ), V ( k ) ) \\ & = \underset { Z } { \arg \min } \ P _ { \lambda } ( Z ) + \frac { c } { 2 } ( \| Z - w ( k ) + \frac { V ( k ) } { c } \| _ { 2 } ^ { 2 } ).$$

Step 2: w updates,

$$w ( k + 1 ) & = \underset { w } { \arg } \min \, L ( Z ( k + 1 ), w, V ( k ) ) \\ & = \underset { w } { \arg } \min \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } O ( y _ { i } w ^ { T } x _ { i } ) + \frac { c } { 2 } ( \| Z ( k + 1 ) - w + \frac { V ( k ) } { c } \| _ { 2 } ^ { 2 } ). \quad ( 4 )$$

$$\underset { w } { r g } \min \, \frac { 1 } { n } \sum _ { i = 1 } O ( y _ { i } w ^ { T } x _ { i } ) + \frac { c } { 2 } ( \| Z ( k + 1 ) - w + \frac { v \, \{ \kappa \} } { c } \| _ { 2 } ^ { 2 } ). \quad ( 4 )$$

Step 3: Lagrange multiplier V updates,

$$V ( k + 1 ) = V ( k ) + c ( Z ( k + 1 ) - w ( k + 1 ) ).$$

For detials, in iteration k +1, we update Z k ( +1) via minimum L Z,w k ,V ( ( ) ( k )) with respect to Z , update w k ( +1) via minimum L Z k ( ( +1) , w, V ( k )) with respect to w , and update Lagrange multiplier via V ( k +1) = V ( k ) + c Z k ( ( +1) -w k ( +1)), until the algorithm converges.

To establish the convergence of ADMM, we have the following assumptions on the data, loss function and regularizer.

Assumption 1 Suppose that ‖ x i ‖ 2 ≤ 1 for all i = 1 , ..., n.

Assumption 2 Loss function O ( ) · is strongly convex and Lipschitz differentiable, with | O ′ | ≤ c 1 , and 0 ≤ O ′′ ≤ c 2 , where c 1 and c 2 are constants.

Assumption 3 Regularizer has the form P λ ( Z ) = ∑ p i =1 p λ ( Z i ), where p λ ( ) · is restricted prox-regular, that is, if for any M &gt; 0 and bounded set T ⊆ dom f ( ( )), there ·

exist γ &gt; 0 such that ∀ z 2 ∈ T \ S M , d ∈ ∂f z ( 1 ) , ‖ d ‖ ≤ M,

$$f ( z _ { 1 } ) + \frac { \gamma } { 2 } \| z _ { 2 } - z _ { 1 } \| _ { 2 } ^ { 2 } \geq f ( z _ { 2 } ) + \langle d, z _ { 1 } - z _ { 2 } \rangle \,,$$

where S M := { Z i ∈ dom f ( ( )) : · ‖ d ‖ &gt; M for all d ∈ ∂f ( ) · } .

Remark 1 In fact, the fitting loss functions include logistic function, and other smooth classification functions, and the regularizers can be some sparsity-inducing functions such as L 1 regularizer, and l q quasi-norms for 0 &lt; q &lt; 1.

Now, we present the convergence of the above ADMM algorithm.

Theorem 1 If preselected parameter c satisfies suitable conditions, the sequence generated by the ADMM algorithm (3), (4),(5) converges to the stationary point of L Z,w,V ( ).

Remark 2 Theorem 1 can be considered as a special case of [30]. This paper analyzed the convergence of ADMM for minimizing possible nonconvex objective problem. Under Assumption 1-3, our sparse classification problem satisfies the convergence conditions A1-A5 in [30], so Theorem 1 holds.

Although ADMM algorithm dose not provide faster convergence compared to other gradient-descent algorithms, it has other advantages. ADMM transform solving of the original optmization problem into solving two simple sub-optimization problems. In Step 1, we only need to solve an univariate sparse regression which can be solved by the existing sparse algorithm. In Step 2, we solve a L 2 regularization problem, which is easy to be handled. More interestingly, the properties of ADMM are useful for obtaining a private version of the sparse calssification model. Notice that in Step 1 of ADMM, the algorithm will never call the training data set, thus it will not cause any leakage of individual privacy, although in this step, the algorithm is unstable due to the sparse regularizer. In step 2, the algorithm will call the training data, but now, the objective function is convex and stable which is suitable for the framework of differential privacy. Step 3 is also data-independent, which will not call the training data in each iteration.

Based on the discussion of this subsection, it's no longer difficult to study the differential privacy for sparse algorithm. In the next subsection, we propose the differential privacy sparse classification method by noising the data-related sub-optimization problem (4). By the post-processing property of differential privacy, the update of Z and V in each iteration will not cause the extra leakage of privacy.

## 2.2. Sparse Classification method with Differential Privacy

In this subsection, we propose differential privacy sparse classification algorithm, and give theoretical privacy bound of the algorithm.

To achieve differential privacy, we perturb the results of each iteration of ADMM by adding noise. Notice that the original data will only be accessed when updating variable w at step 2, so we perturb the objective function (4) by adding cb T w during each iteration, where b is a random vector, with the probability density proportional to exp( - ‖ ‖ γ b 2 ), and γ is the privacy parameter which controls the degree of privacy

leakage. To generate this noisy vector, we choose the norm from the gamma distribution with shape p and scale 1 /γ and the direction uniformly, where p is the number of the feature.

Specifically, in iteration k +1, we first update variable Z via (3), then update variable w by the following perturbing optimization

$$w ( k + 1 ) & = \underset { w } { \arg \min } \, l _ { w } ^ { p r i v } ( k + 1 ) \\ & = \underset { w } { \arg \min } \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } O ( y _ { i } w ^ { T } x _ { i } ) + \frac { c } { 2 } ( \| Z ( k + 1 ) - w + \frac { V ( k ) } { c } \| _ { 2 } ^ { 2 } ) + c b ^ { T } w, \quad ( 6 )$$

$$i = 1$$

At last, update Lagrange multiplier V via (5). Differential privacy sparse classification algorithm(DPSC) is tabulated as Algorithm 1.

## Algorithm 1 Differential Privacy Sparse Classification Algorithm(DPSC)

Input: D = { x , y i i } n i =1 , parameter c , maximum number of iterations K , and privacy degree glyph[epsilon1] .

Initialize: Generate

Z

(0),

w

(0) randomly, let

V

(0) = 0.

Generate noise

b

∼

exp(

- ‖ ‖

γ

b

2

).

for

k

= 0 to

K

-

1

do

Update Z k ( +1) via (3).

Update w k

(

+1) via (6).

Update

V

(

k

+1) via (5).

end for

$$\text{Output: } w ^ { * } = w ( K ).$$

We now state our main result on the privacy property of DPSC. First, we give the theoretical privacy bound of each iteration of DPSC.

Lemma 2 Suppose that Assumption 1-3 hold, and the preselected coefficient c glyph[greaterorequalslant] 2 c 2 n , where c 2 is the upper bound of the second derivative of the loss function O , then the iteration k of DPSC satisfies the glyph[epsilon1] k -DP with

$$\epsilon _ { k } \, = \, \frac { 2 \gamma c _ { 1 } + 2. 8 c _ { 2 } } { n c },$$

where γ is a differential privacy parameter, and c 1 is the upper bound of the first derivative of the loss function.

Proof We prove that for any two data sets D and D ′ differing on at most one data point, given { Z r ( ) } k -1 r =1 , { w r ( ) } k -1 r =1 and { V ( r ) } k -1 r =1 , iteration k of DPSC satisfies

$$\frac { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ) } { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ^ { \prime } ) } \leq \exp ( \epsilon _ { k } ),$$

where S k ( ) is the set of possible outputs in iterations k .

We have

$$& \frac { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ) } { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ^ { \prime } ) } \\ & = \frac { P ( Z ( k ) = Z ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ) } { P ( Z ( k ) = Z ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ^ { \prime } ) } \cdot \frac { P ( w ( k ) = w ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, Z ^ { * } ( k ), D ) } { P ( w ( k ) = w ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, Z ^ { * } ( k ), D ^ { \prime } ) }.$$

Next, we analyze the two parts of (8) separately.

For the first part of (8), because the update of Z is independent of the data, this step will not leak privacy even though no noise is added, and we have

$$\frac { P ( Z ( k ) = Z ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ) } { P ( Z ( k ) = Z ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ^ { \prime } ) } = 1.$$

For the second part of (8), by the KKT condition, we have ∇ w k ( ) l priv w ( k ) = 0, which implies

$$b = - \frac { 1 } { c n } \sum _ { i = 1 } ^ { n } y _ { i } O ^ { ^ { \prime } } ( y _ { i } w ^ { T } ( k ) x _ { i } ) x _ { i } + ( Z ( k ) - w ( k ) + \frac { V ( k - 1 ) } { c } ). \quad \ ( 1 0 )$$

Given { Z r , w r , V ( ) ( ) ( r ) } k -1 r =1 and Z ∗ ( k ), B and W will be bijective, where B and W are random variables, of which the realizations are b and W k ( ). Now, let g k ( · , D ) : R p → R p denote the one-to-one mapping from B to W using dataset D , then

$$& \frac { P ( w ( k ) = w ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, Z ^ { * } ( k ), D ) } { P ( w ( k ) = w ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, Z ^ { * } ( k ), D ^ { \prime } ) } \\ & = \frac { P ( g _ { k } ^ { - 1 } ( w ( k ), D ) ) } { P ( g _ { k } ^ { - 1 } ( w ( k ), D ^ { \prime } ) ) } \cdot \frac { | \det ( J ( g _ { k } ^ { - 1 } ( w ( k ), D ) ) ) | } { | \det ( J ( g _ { k } ^ { - 1 } ( w ( k ), D ^ { \prime } ) ) ) | }.$$

where g -1 k ( w k , D ( ) ) is the mapping from w k ( ) to b using data D and J g ( -1 k ( w k , D ( ) ) is the Jacobian matrix of it.

For the first part in (11), notice that Assumption 2 holds, and for any data point, ‖ x i ‖ 2 ≤ 1 , y i ∈ {-1 1 , } . We choose b ∼ exp( - ‖ ‖ γ b 2 ), and there is only one data difference between D and D ′ , say ( x , y 1 1 ) and ( x , y ′ 1 ′ 1 ) respectively, so from (10) we have

$$\| b ^ { ^ { \prime } } - b \| _ { 2 } = \frac { 1 } { c n } \| y _ { 1 } O ^ { ^ { \prime } } ( y _ { 1 } w ^ { T } ( k ) x _ { 1 } ) - y _ { 1 } ^ { ^ { \prime } } O ^ { ^ { \prime } } ( y _ { 1 } ^ { ^ { \prime } } w ^ { T } ( k ) x _ { 1 } ^ { ^ { \prime } } ) \| _ { 2 } \leq \frac { 2 c _ { 1 } } { c n }. \quad \quad ( 1 2 )$$

Thus, for fixed γ , we have

$$\frac { P ( g _ { k } ^ { - 1 } ( w ( k ), D ) ) } { P ( g _ { k } ^ { - 1 } ( w ( k ), D ^ { \prime } ) ) } & = \exp ( \gamma \| b ^ { ^ { \prime } } \| _ { 2 } - \gamma \| b \| _ { 2 } ) \\ & \leq \exp ( \gamma \| b ^ { ^ { \prime } } - b \| _ { 2 } ) \leq \exp ( \frac { 2 \gamma c _ { 1 } } { c n } ).$$

Consider the second part in (11), the Jacobian matrix is

$$J ( g _ { k } ^ { - 1 } ( w ( k ), D ) = - \frac { 1 } { c n } \sum _ { i = 1 } ^ { n } ( O ^ { ^ { \prime \prime } } ( y _ { i } w ^ { T } ( k ) x _ { i } ) x _ { i } x _ { i } ^ { T } ) - I _ { p \times p },$$

let G k ( ) = 1 cn ( O ′′ ( y w ′ 1 T ( k x ) ′ 1 ) x x ′ 1 ′ T 1 -O ′′ ( y w 1 T ( k x ) 1 ) x x 1 T 1 ), and H k ( ) = -J g ( -1 k ( w k , D ( ) ), then

$$\begin{array} { c } \Phi ^ { 2 } _ { \Phi } & \\ & \frac { | \det ( J ( g _ { k } ^ { - 1 } ( w ( k ), D ) ) ) | } { | \det ( J ( g _ { k } ^ { - 1 } ( w ( k ), D ^ { \prime } ) ) ) | } = \frac { | \det ( H ( k ) ) | } { | \det ( H ( k ) + G ( k ) ) | } \\ & = \frac { 1 } { | \det ( I + H ( k ) ^ { - 1 } \cdot G ( k ) ) | } = \frac { 1 } { \prod _ { i = 1 } ^ { r } ( 1 + \lambda _ { i } ( H ( k ) ^ { - 1 } \cdot G ( k ) ) ) }, \\ & \Phi ^ { 2 } _ { \Phi } & \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ ` } \\ & \Phi ^ { 2 } _ { \Phi } & \Phi ^ { 2 } _ { \Phi } & \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phphi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Philon } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ * \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ < \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 }_ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ 2 } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phis ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } _ { \Phi } \Phi ^ { 2 } }, \end{array}$$

where λ i ( H k ( ) -1 · G k ( ))) denotes the i -th largest eigenvalue of H k ( ) -1 G k ( ). Since G k ( ) has rank at most 2, it implies that H k ( ) -1 G k ( ) also has rank at most 2. Notice that 0 &lt; O ′′ ≤ c 2 , thus λ i ( H k ( )) ≥ 1 &gt; 0 , and -c 2 cn ≤ λ i ( G k ( )) ≤ c 2 cn , which implies that

$$- \frac { c _ { 2 } } { c n } \leq \lambda _ { j } ( H ( k ) ^ { - 1 } G ( k ) ) \leq \frac { c _ { 2 } } { c n }.$$

By the assumption c &gt; 2 c 2 n , we have

$$- \frac { 1 } { 2 } \leq \lambda _ { j } ( H ( k ) ^ { - 1 } G ( k ) ) \leq \frac { 1 } { 2 },$$

since λ min ( H k ( ) -1 G k ( )) ≥ -1 , then

$$\frac { 1 } { | 1 + \lambda _ { m a x } ( H ( k ) ^ { - 1 } G ( k ) ) | ^ { 2 } } \leq \frac { 1 } { | \det ( I + H ( k ) ^ { - 1 } G ( k ) ) | } \leq \frac { 1 } { | 1 + \lambda _ { m i n } ( H ( k ) ^ { - 1 } G ( k ) ) | ^ { 2 } }.$$

thus

$$\frac { | \det ( J ( g _ { k } ^ { - 1 } ( w ( k ), D ) ) ) | } { | \det ( J ( g _ { k } ^ { - 1 } ( w ( k ), D ^ { \prime } ) ) ) | } \leq \frac { 1 } { | 1 - \frac { c _ { 2 } } { c n } | ^ { 2 } } \leq \exp ( - 2 \ln ( 1 - \frac { c _ { 2 } } { c n } ) ) \leq \exp ( \frac { 2. 8 c _ { 2 } } { n c } ). \quad ( 1 5 )$$

where in the last inequality, we use the fact that for any x ∈ [ 0 , 1 2 ] , -ln(1 -x ) &lt; 1 4 . x. Combining (13) and (15),

$$\frac { P ( w ( k ) = w ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, Z ^ { * } ( k ), D ) } { P ( w ( k ) = w ^ { * } ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, Z ^ { * } ( k ), D ^ { \prime } ) } \leq \exp ( \frac { 2 \gamma c _ { 1 } + 2. 8 c _ { 2 } } { c n } ). \quad ( 1 6 )$$

Combining (9) and (16) , the lemma has been proved.

Theorem 3 Suppose that Assumption 1-3 hold, and the preselected coefficient c ≥ 2 c 2 n , then DPSC satisfies the glyph[epsilon1] -differential privacy with glyph[epsilon1] = K · ( 2 γc 1 +2 8 . c 1 cn ), where K is the number of iteration.

Proof We prove that for any two datasets D and D ′ that differing on at most one data point, DPSC satisfies

$$\frac { P ( \{ S ( k ) \} _ { k = 1 } ^ { K } \in S | D ) } { P ( \{ S ( k ) \} _ { k = 1 } ^ { K } \in S | D ^ { ^ { \prime } } ) } \leq \exp ( \epsilon ).$$

where S is the set of possible outputs during K iterations.

We have

$$P$$

$$& \quad \text{ we have} \\ & \frac { P ( \{ Z ( k ), w ( k ) \} _ { k = 1 } ^ { K } \in S \ | \ D ) } { P ( \{ Z ( k ), w ( k ) \} _ { k = 1 } ^ { K } \in S \ | \ D ^ { \prime } ) } \\ & = \frac { P ( \{ Z ( 0 ), w ( 0 ) \} \in S ( 0 ) \ | \ D ) } { P ( \{ Z ( 0 ), w ( 0 ) \} \in S ( 0 ) \ | \ D ^ { \prime } ) } \cdot \prod _ { k = 1 } ^ { K } \frac { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ) } { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ^ { \prime } ) } \\ & = \prod _ { k = 1 } ^ { K } \frac { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ) } { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) \ | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ^ { \prime } ) }, \\ & \quad \text{and since was an array of $Z(\mathbb{M}$ and $\mathbb{M}$ and $\mathbb{M}$ and $\mathbb{M}$ and $\mathbb{M}$ but that at anolitrix which would $1$} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{$$

and since we generate Z (0) and w (0) randomly, the last equality achieved. Lemma 2 shows that for each k , k ≥ 1,

$$\frac { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ) } { P ( \{ Z ( k ), w ( k ) \} \in S ( k ) | \ \{ Z ( r ), w ( r ) \} _ { r = 1 } ^ { k - 1 }, D ^ { \prime } ) } \leq \exp ( \frac { 2 \gamma + 2. 8 c _ { 1 } } { c n } ),$$

Thus we have

$$\frac { P ( \{ S ( k ) \} _ { k = 1 } ^ { K } \in S | D ) } { P ( \{ S ( k ) \} _ { k = 1 } ^ { K } \in S | D ^ { \prime } ) } \leq \exp ( K \cdot ( \frac { 2 \gamma + 2. 8 c _ { 1 } } { c n } ) ).$$

Hence the theorem is proved.

From Theorem 3, we can see that the privacy bound of DPSC is controled by the iteration number K , the privacy parameter γ , the upper bound of the first and second derivative of loss function c 1 and c 2 , ADMM pre-selected parameter c , and the data size n . With increasing of γ , the privacy bound increase, which means that the algorithm leaks more privacy. Meanwhile, the noisy vector b decrease, that is, we perturb the objective function with a smaller scale. Thus, the utility of algorithm increase. In addition, as the size of training data increase, the privacy bound decrease, i.e. the possibility of individual privacy leakage decrease. It shows that the large traning data make individual more difficult to be distinguished.

Remark 3 Based on the idea of perturbing the optimization problem[25], we consider noising the objective function of Z and w in each itertaion of ADMM to achieve differential privacy. Notice that the update of Z is independent of the data, that is, the original data will not be accessed when updating variable Z . Thus we only perturb the objective function of w , then propose DPSC algorithm.

## 3. Applications

In this section, we apply the differential privacy sparse classification model to two typical models - logistic regression with L 1 regularizer and logistic regression with L 1 2 / regularizer. We propose privacy classification algorithms, and present the privacy bound of them. The applications show that our method is suitable for dealing with sparse classification problems.

## 3.1. Differential Privacy Logistic Regression with L 1 Regularizer

The logistic regression with L 1 regularizer [31] is a useful method to solve over-fitting problem and has well generalization [32, 33]. It selects variables according to the amount of penalization on the l 1 -norm of the coefficients. The logistic regression with L 1 regularizer has the following form

$$\min _ { w } \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } \log ( 1 + \exp ( - y _ { i } w ^ { T } x _ { i } ) ) + \lambda \| w \| _ { 1 },$$

where x i ∈ R p is input variable, y i ∈ {-1 1 , } is output variable. First, we add the auxiliary variable Z , and rewrite the model as

$$\min _ { w, Z } \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } \log ( 1 + \exp ( - y _ { i } w ^ { T } x _ { i } ) ) + \lambda \| Z \| _ { 1 }$$

$$s. t. \ w & = Z.$$

which has the quadratically augmented Lagrangian function

$$L ( Z, w, V ) = \frac { 1 } { n } \sum _ { i = 1 } ^ { n } \log ( 1 + \exp ( - y _ { i } w ^ { T } x _ { i } ) ) + \lambda \| Z \| _ { 1 }$$

$$+ \frac { c } { 2 } ( \| Z - w + \frac { V } { c } \| _ { 2 } ^ { 2 } - \| \frac { V } { c } \| _ { 2 } ^ { 2 } ),$$

where V is Lagrange multipliers. From Algorithm 1, Differential Privacy Logistic Regression with L 1 Regularizer algorithm (DPLL) has the following steps

Step 1: Z updates,

$$Z ( k + 1 ) = \underset { Z } { \arg \min } \ \lambda \| Z \| _ { 1 } + \frac { c } { 2 } ( \| Z - w ( k ) + \frac { V ( k ) } { c } \| _ { 2 } ^ { 2 } ),$$

Step 2: w updates,

$$w ( k + 1 ) = & \arg \min _ { w } \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } \log ( 1 + \exp ( - y _ { i } w ^ { T } x _ { i } ) ) \\ & + \frac { c } { 2 } ( \| Z ( k + 1 ) - w + \frac { V ( k ) } { c } \| _ { 2 } ^ { 2 } ) + c b ^ { T } w,$$

Step 3: Lagrange multiplier V updates,

$$V ( k + 1 ) = V ( k ) + c ( Z ( k + 1 ) - w ( k + 1 ) ).$$

Consider the updating of Z , let Q k ( ) = w k ( ) -V ( k ) c , for i = 1 , · · · , p , the coordinate Z i of vector Z is updated by

$$Z _ { i } ( k + 1 ) = \underset { Z _ { i } } { \arg \min } \ \lambda | Z _ { i } | + \frac { c } { 2 } ( Q _ { i } ( k ) - Z _ { i } ) ^ { 2 },$$

where Q k i ( ) is the i th coordinate of vector Q k ( ). It's easy to obtain the solutions of (23), which are

$$Z _ { i } ( k + 1 ) = \begin{cases} \begin{array} { c c c } Q _ { i } ( k ) - \frac { \lambda } { c } & \quad i f \ Q _ { i } ( k ) \geq \frac { \lambda } { c } \\ 0 & \quad i f \ - \frac { \lambda } { c } < Q _ { i } ( k ) \leq \frac { \lambda } { c } \\ Q _ { i } ( k ) + \frac { \lambda } { c } & \quad i f \ Q _ { i } ( k ) < \frac { \lambda } { c }. \end{array} \end{cases}$$

To update variable w , we consider using gradient descent algorithm. Let m = 1 2 , , · · · , M denote the inner iteration index for the gradient descent algorithm used to solve (21). For the minimization at step k + 1 of the (outer) consensus iteration, the sequence of iterates w 0 ( k +1) = w k ( ), at step m

$$w ^ { m } ( k + 1 ) = w ^ { ( m - 1 ) } ( k + 1 ) - \alpha \frac { \partial l _ { w } ^ { p r i v } ( w ^ { ( m - 1 ) } ( k + 1 ) ) } { \partial w },$$

where ∂l priv w ( w ( m -1) ( k +1)) ∂w = 1 n ∑ n i =1 1+exp( -y x i i y w i T x i ) -c Z k ( ( +1) -w ( m -1) ( k )+ V ( k ) c )+ c · b, and α is learning rate.

Now, we present the Differential Privacy Logistic regression with L 1 regularizer algorithm (DPLL). The overall operation of the algorithms can be described as follows. During iteration k + 1, we update each coordinate i of Z via (24), then update w cyclically via (25), at last, update Lagrange multiplier V via (22).

The theoretical privacy bound of DPLL is presented in next subsection.

## 3.2. Differential Privacy Logistic Regression with L 1 2 / Regularizer

In this subsection, we focus on logistic regression with nonconvex penalty. As a typical nonconvex regularization, L 1 2 / has many desirable properties and has been studied extensively. [34] proposed the logistic regression with L 1 2 / regularizer which enhances the variable selection capability and alleviates the over-fitting problem of the traditional

Algorithm 2 Differential Privacy Logistic Regression with L 1 Regularizer algorithm(DPLL)

Input: D = { x , y i i } n i =1 , parameter c , maximum number of iterations K , maximum number of inner iteration M , learning rate α , and privacy degree glyph[epsilon1] . Initialize: Generate Z (0) and w (0) randomly, and let V (0) = (0 , ..., 0). Generate noise b ∼ exp( - ‖ ‖ γ b 2 ).

$$\text{for } k = 0 \text{ to } K - 1 \text{ do }$$

$$\text{for } i = 1 \text{ to } p \text{ do }$$

Update

Z

i

(

k

+1) via (24).

end for for

m

= 1 to

Update

M

w

m

do

(

k

end for

Let w k

(

Update

+1) =

+1) via (25).

w

M

(

k

+1).

V

(

k

+1) via (22).

end for

$$\text{Output: } w ^ { * } = w ( K ).$$

model. It has the following form

$$\min _ { w } \, \frac { 1 } { n } \sum _ { i = 1 } ^ { n } \log ( 1 + \exp ( - y _ { i } w ^ { T } x _ { i } ) ) + \lambda \| w \| _ { 1 / 2 } ^ { 1 / 2 },$$

where ‖ w ‖ 1 2 / 1 2 / = ∑ p i =1 | w i | 1 2 / is quasi-norm. Same as before, we add auxiliary variable Z , and rewrite the augmented Lagrangian function as

$$L ( Z, w, V ) = \frac { 1 } { n } \sum _ { i = 1 } ^ { n } \log ( 1 + \exp ( - y _ { i } w ^ { T } x _ { i } ) ) + \lambda \| Z \| _ { 1 / 2 } ^ { 1 / 2 }$$

$$+ \frac { c } { 2 } ( \| Z - w + \frac { V } { c } \| _ { 2 } ^ { 2 } - \| \frac { V } { c } \| _ { 2 } ^ { 2 } ).$$

Based on Algorithm 1, we propose Differential Privacy Logistic regression with L 1 2 / regularizer (DPLH), which is tabulated as Algorithm 3. Here, the update of w and V in each iteration are same as DPLL.

In iteration k +1 of DPLH, to solve (27), we resort to reweighted algorithm, which has the following steps

Step 1: Set the initial value Z 0 and the maximum iteration step T , and let t = 0.

Step 2: Solve

$$Z ^ { t + 1 } = \arg \min _ { Z } \, \frac { c } { 2 } ( \| Z - w ( k ) + \frac { V ( k ) } { c } \| _ { 2 } ^ { 2 } ) + \lambda \sum _ { i = 1 } ^ { p } \frac { | Z _ { i } | } { | Z _ { i } ^ { t } | ^ { 1 / 2 } }$$

```
                                                                                                                                                                                                       
                                                                                                                                                                                                       

                                                                                                                                                                                                        }
                                                                                                                                                                                                      
	                                                                                                                                                                                                       	                                                                                                                                                                                                      

                                                                                                                                                                                                       }

                                                                                                                                                                                                      
---------------------------------------------------------------------------------------------------------------------                                                                                                                                                                                                        \
                                                                                                                                                                                                       \

                                                                                                                                                                                                       }                                                                                                                                                                                                        |                                                                                                                                                                                                        ]
                                                                                                                                                                                                       |
                                                                                                                                                                                                      


                                                                                                                                                                                                       ]

                                                                                                                                                                                                       |
```

with an existing L 1 algorithm, and let t = t +1 .

Step 3: If t &lt; T , go to Step 2, otherwise, output Z k ( +1) = Z t .

The initial value Z 0 normally can be taken as Z 0 = (1 1 , , ..., 1). And with such a setting, the first iteration ( t = 0) in Step 2 is exactly corresponding to solving a lasso problem. When t = 1, with easy linear transformation, step 2 can be transformed into a lasso problem. It is possible that when t &gt; 1, some Z t i = 0. To guarantee the feasibly, we replace | Z t i | 1 2 / with | Z t i + µ | 1 2 / in Step 2, where µ is any fixed positive real number.

The overall operation of DPLH can be described as follows. During iteration k +1, we update variable Z via reweighted algorithm, then update w cyclically via (25), at last, update Lagrange multiplier V via (22).

Now, we present the theoretical privacy bound of the DPLL and DPLH.

Theorem 4 Suppose that preselected penalty coefficient c ≥ 1 2 n and γ ≤ cn -7 20 , then DPLL and DPLH satisfy the glyph[epsilon1] -differential privacy with glyph[epsilon1] = K · glyph[epsilon1] 1 , where glyph[epsilon1] 1 = 8 γ +2 8 . 4 cn , and K is the iteration number.

Proof Obviously, the logistic function is strongly convex, Lipschitz differentiable, with 0 ≤ O ′′ ≤ 1 4 , which satisfies Assumption 2, here c 1 = 1 4 . And L 1 regularizer and L 1 2 / regularizer satisfy Assumption 3. Thus, from the proof of Theorem 3, theorem is established.

## 4. Experiments

We use the simulation and real data to illustrate the effectiveness of our algorithms DPLL and DPLH. In all datasets, we preprocessed the data such that the input vectors had maximum norm 1.

## 4.1. Simulation

We study the classification model

$$p ( y _ { i } = 1 | x _ { i } ) = \frac { 1 } { 1 + \exp ( - w ^ { T } x _ { i } ) }, \ i = 1, \dots, n.$$

where x i ∈ R 100 is an input vector, x i ∼ N (0 , Σ) , covariance matrix Σ = 0 5 . | i - | j , 1 ≤ i, j ≤ n, coefficient vector w ∈ R 100 , where w (1 : 8) = (10 9 8 7 6 5 4 0 5) , and , , , , , , , . ′ w (9 : 100) = (0 , 0 , ..., 0). glyph[epsilon1] is standardized normal random error, and sample number n . If p y ( i = 1 | x i ) ≥ 0 5 . , let y i = 1, else, y i = -1 .

## 4.1.1. Privacy-Accuracy Tradeoff

In this case, we study the tradeoff between the privacy requirement on the classifier, and its classification accuracy, when the classifier is trained on data of a fixed size. The privacy bound of DPLL and DPLH is glyph[epsilon1] = K · 8 γ +2 8 . 4 cn . Thus, the privacy requirement is measured by the privacy parameter γ , iteration number K , pre-selected parameter c and training data size n . In this case, we set K = 100, c = 2 5, and . n = 10000. Then we run DPLL and DPLH with γ increases. To measure accuary of two algorithms, we record the classification error rate (CE) and mean square error (MSE). We also record the average number of correctly identified zero coeffcients (C.A.N of zero) and the average number of incorrectly identfied zero coe cients (IC.A.N of zero) over the datasets to measure the ability of variable selection of two algorithms.

We have simulated 50 datasets consisting of 11,000 observations and each data set was divided into two parts: a training set with 10,000 observations and a test set with 1000 observations. To demonstrate the effectiveness of DPLL and DPLH, we applied the logistic regression with L 1 algorithm (denote as LLA), the logistic regression with L 1 2 / algorithm (denote as LHA), DPLL and DPLH to the 50 datasets. For each data set, we trained classifiers for 5 fixed values of λ and tested the error of these classifiers. Specifically, for non-private algorithms, we chose the tunning parameter λ by 5-fold CV. And for private algorithms, we set the same value of λ with non-private algorithms for the same dataset. We take differential privacy degree glyph[epsilon1] as 0 1 . , 0 5 . , 1 1 5 2 2 5 3 and , . , , . , 4, and let gradient descent step α as 0.5. Due to the randomness of the noise, we performed 50 independent runs of two privacy algorithms for each parameter setting, and recorded the average results of 50 trials. Figure.1 and Figure.2 show the variable selection results of non-private algorithms and private algorithms. Figure.3 and Figure.4 show the classification results of four algorithms.

From Figure 1 and Figure 2, we can see that C.A.N of zero of LLA and LHA are larger than C.A.N of zero of DPLL and DPLH. And with the increase of privacy degree glyph[epsilon1] , C.A.N of zero of DPLL and DPLH both increase. The stricter privacy-preserving (smaller glyph[epsilon1] ), the lower variable selection accuracy (smaller C.A.N of zero). These results show that privacy and variable selection accuracy can not be simultaneously satisfied when training data size n is fixed and not large enough. What's more, the variable selection results of logistic regression with L 1 2 / regularizer are more sparse than logistic regression with L 1 regularizer, even though IC.A.N of zero of logistic regression with L 1 2 / regularizer is slightly larger than logistic regression with L 1 regularizer.

From Figure 3 and Figure 4, we can see that prediction error (MSE and CE) of DPLL and DPLH are larger than prediction error (MSE and CE) of LLA and LHA. And with the increase of privacy degree glyph[epsilon1] , prediction error (MSE and CE) of DPLL

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-9919c397-0b0a-4d5b-9145-d3270d472514.jpg)

Privacy degree Epsilon, n=10000, p=100

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-45faa7d1-83b4-4264-af1e-250845deba47.jpg)

Figure 1. With the increase of privacy degree glyph[epsilon1] , results of C.A.N of LLA, LHA,DPLL and DPLH.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-550b86e3-82a1-48be-8532-5a5ace91fcf5.jpg)

Figure 2. With the increase of privacy degree glyph[epsilon1] , results of IC.A.N of LLA, LHA,DPLL and DPLH.

Figure 3. With the increase of privacy degree glyph[epsilon1] , results of Classification Error Rate of LLA, LHA,DPLL and DPLH.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-db84f1cb-54ae-418a-9907-c1ae10c5d971.jpg)

Figure 4. With the increase of privacy degree glyph[epsilon1] , results of Mean Square Error of LLA, LHA,DPLL and DPLH.

and DPLH both decrease. The stricter privacy-preserving (smaller glyph[epsilon1] ), the lower prediction accuracy (larger MSE and CE ). These results show that privacy and prediction accuracy can not be simultaneously satisfied, when training data size n is fied and not large enough.

## 4.1.2. Accuracy-Training data size Tradeoff

In this case, we study the tradeoff between training data size and prediction accuracy. Specifically, we examine how prediction accuracy varies as the size of the training set increases when the privacy degree is fixed.

In this case, we set K = 150, c = 2 5, and . glyph[epsilon1] = 0 5 and 1 5 separately. We have . . simulated 6 datasets consisting of 5000 , 10000 , ..., 30000 training observations and have simulated a test set with 1000 observations. Then we run LLA, LHA, DPLL and DPLH with training data size n increases. To measure accuary of two algorithms, we record the classification error rate (CE) and mean square error (MSE). For each data set, we also trained classifiers for 5 fixed values of λ , and let gradient descent step α as 0.5. Due to the randomness of the noise, we performed 50 independent runs of two privacy algorithms for each parameter setting, and recorded the average results of 50 trials. Figure.5 and Figure.6 show the prediction results (MSE and CE) of four algorithms.

From Figure 5 and Figure 6, we can see that prediction error (MSE and CE) of DPLL and DPLH are larger than prediction error (MSE and CE) of LLA and LHA. And with the increase of the size of training data, prediction error (MSE and CE) of DPLL and DPLH both decrease for two glyph[epsilon1] . The privacy bound glyph[epsilon1] = K · 8 γ +2 8 . 4 cn of

Figure 5. With the increase of training data size n , results of Classification Error Rate of LLA, LHA,DPLL and DPLH.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-436b5c16-436b-4201-b4e4-89d6fbeae25d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-ef7434d0-d31e-4f83-80ad-30c386cf4f37.jpg)

Figure 6. With the increase of training data size n , results of Mean Square Error of LLA, LHA,DPLL and DPLH.

DPLL and DPLH shows that for fixed glyph[epsilon1] and other parameters, with the increase of n , γ which controls the norm of random variable b increases. Thus the prediction error decreases. And our experiment results are consistent with this conclusion. Interestingly, when the size of training data n is large enough, the results of the differential privacy algorithms are similar as the corresponding non-private algorithms. Big data offers more effectively protection to individual privacy.

## 4.2. Real data

We analyze the KDDCup99 data set from the UCI Machine Learning Repository [35]. The task is to build a network classifier, which can predict whether a network connection is a denial-of-service attack or not. This data set contains about 5,000,000 instances, and each instance includes 41 attributes. For this data the average fraction of negative labels is 0.20. To preprocess the data, first we converted each classiffication attribute to a binary vector, and converted lables { good connections, bad connections } to { 1 , - } 1 . Then we normalized each column to ensure that the maximum value is 1, and finally, we normalized each row to ensure that the norm of any example is at most 1. After preprocessing, each instance was represented by a 118-dimensional vector, of norm at most 1.

In this case, we study the privacy-accuracy tradeoff and the accuracy-training data size tradeoff. We applied LLA, LHA, DPLL and DPLH to the subset of KDDCup99 data set with different parameter settings. For non-private and private algorithms, we chose the tunning parameter λ by 5-fold CV. When applied classification algorithms, let gradient descent step α as 0.2, iteration number K as 150, and pre-selected parameter c as 1.

To study the privacy-accuracy tradeoff, we chose 60000 examples from the original data set, of which the number of positive label is 38420. We applied four algorithms to this data set, set glyph[epsilon1] as 0 1 . , 0 5 . , 1 1 5 2 2 5 , . , , . , and 3, and performed 50 independent runs of two private algorithms for each parameter setting. We recorded the average results of 50 trials. Figure.7 and Figure.8 show the prediction results (MSE and CE) of four algorithms.

To study the accuracy-training data size tradeoff, we chose 6 sub-datasets of size n = 10000 to n = 60000 from the original data, and applied LLA, LHA, DPLL and DPLH to these datasets. We also performed 50 independent runs of two private algorithms for each parameter setting, and recorded the average results of 50 trials. Figure.9 and Figure.10 show the prediction results.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-e5d1b87c-ebb4-49b7-8850-0b0a56223780.jpg)

Figure 7. With the increase of privacy degree glyph[epsilon1] , CE results of KDDCup99.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-bcff7c62-49fb-4b25-8101-c64149b38382.jpg)

Figure 8. With the increase of privacy degree glyph[epsilon1] , MSE results of KDDCup99.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-bed0ae5b-0560-4c02-9d40-3acf1d6b5ffd.jpg)

Figure 9. With the increase of training data size n , CE results of KDDCup99.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00780v1-figure-7ff76eb7-6b95-46fd-9551-2759da02e681.jpg)

Figure 10. With the increase of training data size n , MSE results of KDDCup99.

From Figure 7-8, we can see that with the increase of the privacy budget, the prediction error of private algorithms decrease. From Figure 9-10, we can see that for private algorithms, the prediction error of private algorithms decrease. The results are consistent with theoratical results in Theorem 4. This case shows that our methods are suitable for dealing with real data.

## 5. Conclusions

In this paper, we have proposed a differential privacy version of convex and nonconvex sparse classification approach without stable condition on regularizer. Based on ADMM algorithm, we transform the solving of sparse problem into a multistep iteration process. In each iteration of our algorithm, we first deal with an univariate sparse regression which can be solved by the existing sparse algorithm, then we solve a datarelated L 2 regularization problem which is easy to be handled. To achieve privacy protecting, we perturb the data-related step in each iteration. And by the property of the post-processing holding of differential privacy, the proposed approach satisfies the glyph[epsilon1] -differential privacy even when the original problem is unstable. The theoretical privacy bound of the classification algorithm is glyph[epsilon1] = K · ( 2 γc 1 +2 8 . c 2 cn ). The privacy bound of our algorithm is controlled by the algorithm iteration number K , the privacy parameter γ , the parameter of loss function c 1 and c 2 , ADMM pre-selected parameter c , and the data size n . With the increase of privacy parameter γ , the privacy bound glyph[epsilon1] increases, which means that the algorithm leaks more privacy. Meanwhile, we perturb the objective functions with a smaller scale. Thus, the utility of algorithm increases. In addition, as the size of training data n increases, the privacy bound glyph[epsilon1] decreases, i.e. the

possibility of individual privacy leakage decreases. It shows that the large traning data make it difficult to distinguish individuals. At last, we apply our framework to logistic regression with L 1 regularizer and logistic regression with L 1 2 / regularizer. Numerical studies demonstrate that our method is both effective and efficient which performs well in sensitive data analysis. The differential privacy sparse classification framework proposed in this work can be easily generalized to other regularization methods, such as SCAD, MCP and so on. Furthermore, ADMM algorithm is an effective method to solve distributed problems. Our framework can also be generalized to distributed sparse optimization problems. All these problems are under our current research.

## References

- [1] L. Sweeney, K-anonymity: a model for protecting privacy[J], International Journal of Uncertainty, Fuzziness and Knowledge-Based Systems. 10(05)(2002) 557-570.
- [2] L. Sweeney, Achieving k-anonymity privacy protection using generalization and suppression[J], International Journal of Uncertainty, Fuzziness and Knowledge-Based Systems. 10(05)(2008) 571-588.
- [3] A. Machanavajjhala, J. Gehrke, D. Kifer, et al., L-diversity: privacy beyond kanonymity[C], Proceedings of the 22nd International Conference on Data Engineering. IEEE, 2006, pp.24.
- [4] N. Li, T. Li, S. Venkatasubramanian, T-closeness: privacy beyond k-anonymity and ldiversity[C], Proceedings of the 23rd International Conference on Data Engineering. IEEE, 2007, pp.106-115.
- [5] R. GantaP. KasiviswanathanA. Smith, Composition attacks and auxiliary information in data privacy[C], Proceedings of the 14th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. ACM, 2008. pp. 265-273.
- [6] W. Wong, A. Fu, K. Wang, et al., Can the utility of anonymized data be used for privacy breaches[J], ACM Transactions on Knowledge Discovery from Data. 2011, 5(16)(2011)1-24.
- [7] C. Dwork, Differential privacy[J], Lecture Notes Computer Science . Sci. 26(2)(2006) 1-12.
- [8] C. Dwork, A. Roth, The algorithmic foundations of differential privacy[M], Now Publishers Inc. 2014, 1-300.
- [9] P. Liu, Y. Xu and Q. Jiang et al., Local differential privacy for social network publishing[J], Neurocomputing, https://doi.org/10.1016/j.neucom.2018.11.104.
- [10] P.Sui, X.Li, A privacy-preserving approach for multimodal transaction data integrated analysis[J], Neurocomputing, 253 (2017) 56-64.
- [11] J. Brinkrolf, C. Gopfert, B. Hammer, Differential privacy for learning vector quantization[J], Neurocomputing, 342 (2019) 125-136.
- [12] C. Dwork, F. Mcsherry, K. Nissim, Calibrating noise to sensitivity in private data analysis[J]. Conference on Theory of Cryptography, Springer, Verlag, 2006, pp. 265-284.
- [13] F Mcsherry, K Talwar. Mechanism Design via Differential Privacy[C]. In: Proceedings of the 48th Annual Symposium.
- [14] J. Bolot, N. Fawaz, S. Muthukrishnan, et al., Private decayed predicate sums on streams[C], Proceedings of the 16th International Conference on Database Theory. ACM, 2013, pp. 284-295.
- [15] V. Karwa, S. Raskhodnikova, G. Yaroslavtsev, Private analysis of graph structure[J]. Proceedings of the VLDB Endowment, 4(11)(2011)1146-1157.
- [16] G. Cormode, C. Procopiuc, D. Srivastava, et al., Differentially private spatial decompositions[C]. Proceedings of the 28th International Conference on Data Engineering. IEEE, 2012, pp. 41(4)20-31.
- [17] C. Li, M. Hay, G. Miklau, et al., A data- and workload-aware algorithm for range queries under differential privacy[J], Proceedings of the VLDB Endowment ,7(5)341-352.
- [18] R. Tibshirani, Regression shrinkage and selection via the LASSO[J], Journal of the Royal

Statistical Society, 58(1996) 267-288.

- [19] Z Xu, X Chang, X Xu, et al., L 1 / 2 regularization: A thresholding represention theory and a fast sover[J], IEEE Transaction on Neural Network and Learning Systems, 23(7)(2012) 1013-1027.
- [20] C Dwork, J Lei, Differential privacy and robust statistics[C], in: Annual ACM Symposium on Theory of Computing, 2009,pp. 371-380.
- [21] H. Xu, C. Caramanis, S. Mannor, Sparse algorithms are not stable: a no-free-lunch theorem[J], IEEE Transactions on Pattern Analysis and Machine Intelligence, 34(1)(2012)187193.
- [22] D. Kifer, A. Smith, A. Thakurta, Private convex empirical risk minimization and highdimensional regression[J], Journal of Machine Learning Research, 25(2012)1-40.
- [23] A. Smith, A. Thakurta, Differentially private feature selection via stability arguments, and the robustness of the Lasso[C], Annual Conference on Learning Theory. PMLR, 2013, vol.30, pp. 819-850.
- [24] K. Talwar, A. Thakurta, L. Zhang, Nearly-optimal private LASSO[C], Proceedings of the 28th International Conference on Neural Information Processing Systems. ACM,2015, vol. 2, pp. 3025-3033.
- [25] K. Chaudhuri, C. Monteleoni, A. Sarwate, Differentially private empirical risk minimization.[J], Journal of Machine Learning Research, 12(2011)1069-1109.
- [26] T. Zhang, Q. Zhu, Dynamic differential privacy for ADMM-based distributed classification learning[J], IEEE Transactions on Information Forensics and Security, 12(1)(2017)172-187.
- [27] X. Zhang, M. Khalili, M. Liu, Improving the privacy and accuracy of ADMM-based distributed algorithms[C], Proceedings of the 35th International Conference on Machine Learning. PMLR,2018, pp. 5796-5805.
- [28] P. Wang, H. Zhang, Distributed logistic regression with differential privacy (in Chinese)[J], Sci Sin Inform, doi10.1360/N112018-00214.
- [29] S. Boyd, N. Parikh, E. Chu: Distributed optimization and statistical learning via the alternating direction method of multipliers[M], Now Foundations and Trends, 2011,1-128.
- [30] Y. Wang, W. Yin, J. Zeng: Global convergence of ADMM in nonconvex nonsmooth optimization[J], Journal of Scientific Computing,78(1)(2019)29-63.
- [31] M. Park, T. Hastie: L 1 -regularization path algorithm for generalized linear models[J], Journal of Royal Statistical Society: Series B, 69(4)(2007)659-677.
- [32] A. Genkin, D. Lewis, D .Madigan, Large-scale bayesian logistic regression for text categorization[J], Technometrics, 49(2007)291-304.
- [33] K. Koh, S. Kim, S. Boyd, An interior-point method for large-scale L 1 regularized logistic regression[J], Journal of Machine Learning Research, 8(8)(2007)1519-1555.
- [34] Q. Zhao, D. Meng, Z. Xu, L 1 / 2 regularized logistic regression[J],Pattern Recognition and Artificial Intelligence, 25(5)(2012)721-728.
- [35] S. Hettich and S.D. Bay. The UCI KDD Archive. University of California, Irvine, Department of Information and Computer Science, 1999. URL http://kdd.ics.uci.edu.
Puyu Wang Hai Zhang
Differential Privacy for Sparse Classification Learning
In this paper, we present a differential privacy version of convex and
nonconvex sparse classification approach. Based on alternating direction method
of multiplier (ADMM) algorithm, we transform the solving of sparse problem into
the multistep iteration process. Then we add exponential noise to stable steps
to achieve privacy protection. By the property of the post-processing holding
of differential privacy, the proposed approach satisfies the
$\epsilon-$differential privacy even when the original problem is unstable.
Furthermore, we present the theoretical privacy bound of the differential
privacy classification algorithm. Specifically, the privacy bound of our
algorithm is controlled by the algorithm iteration number, the privacy
parameter, the parameter of loss function, ADMM pre-selected parameter, and the
data size. Finally we apply our framework to logistic regression with $L_1$
regularizer and logistic regression with $L_{1/2}$ regularizer. Numerical
studies demonstrate that our method is both effective and efficient which
performs well in sensitive data analysis.