# Patent

## Title
Differentially private processing and database storage

## Abstract

A hardware database privacy device is communicatively coupled to a private database system. The hardware database privacy device receives a request from a client device to perform a query of the private database system and identifies a level of differential privacy corresponding to the request. The identified level of differential privacy includes privacy parameters (ε,δ) indicating the degree of information released about the private database system. The hardware database privacy device identifies a set of operations to be performed on the set of data that corresponds to the requested query. After the set of data is accessed, the set of operations is modified based on the identified level of differential privacy such that a performance of the modified set of operations produces a result set that is (ε,δ)-differentially private.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
This application is a continuation of U.S. patent application Ser. No. 16/810,708, filed Mar. 5, 2020, continuation of U.S. patent application Ser. No. 16/238,439, filed Jan. 2, 2019, which is a continuation of U.S. patent application Ser. No. 15/793,907, filed Oct. 25, 2017, which is a continuation of U.S. patent application Ser. No. 15/203,797, filed Jul. 7, 2016, which claims the benefit of U.S. Patent Application No. 62/249,938, filed Nov. 2, 2015, all of which are incorporated by reference herein.


BACKGROUND
Field of Disclosure
The present invention generally relates to a database system, and more specifically to responding to a database query by executing a differentially private version of the query on the database.
Description of the Related Art
Personally identifiable information, such as health data, financial records, telecom data, and confidential business intelligence, such as proprietary data or data restricted by contractual obligations, is valuable for analysis and collaboration. Yet, only a fraction of such sensitive information is used by organizations or analysts for statistical or predictive analysis. Privacy regulations, security concerns, and technological challenges suppress the full value of data, especially personally identifiable information and confidential and proprietary records.
Methods that attempt to solve this problem, such as access controls, data masking, hashing, anonymization, aggregation, and tokenization, are invasive and resource intensive, compromise analytical utility, or do not ensure privacy of the records. For example, data masking may remove or distort data, compromising the statistical properties of the data. As another example, many of the above mentioned methods are not effective when information is stored in disparate data sources. Technology which enables organizations or analysts to execute advanced statistical and predictive analysis on sensitive information across disparate data sources without revealing record-level information is needed.
SUMMARY
A hardware database privacy device is communicatively coupled to a private database system. The hardware database privacy device receives a request from a client device to perform a query of the private database system and identifies a level of differential privacy corresponding to the request. The identified level of differential privacy includes a privacy parameter ε indicating the degree of information released about data in the private database system.
The differentially private hardware database privacy device identifies a set of data stored in the private database system and a set of operations to be performed on the set of data that corresponds to the requested query. After the set of data is accessed, the set of operations is modified based on the identified level of differential privacy such that a performance of the modified set of operations produces a result set that is differentially private. The modified set of operations is performed on the accessed set of data to produce the differentially private result set. The differentially private result set is provided to the client device for display on a hardware display of the client device.



BRIEF DESCRIPTION OF DRAWINGS
FIG. 1 illustrates a system for receiving a query for a private database, and for responding to the query by executing a differentially private version of the query on the private database.
FIG. 2 illustrates an example database structure, according to one embodiment.
FIG. 3 is a block diagram illustrating the privacy system of the system in FIG. 1, according to one embodiment.
FIG. 4 illustrates displaying results of a differentially private count query, according to one embodiment.
FIG. 5 illustrates an example binary decision tree for use in a differentially private random forest query, according to one embodiment.
FIG. 6 illustrates perturbing the counts for a differentially private histogram query, according to one embodiment.
FIG. 7A illustrates a recursive process for identifying threshold points of the classification output vector for a differentially private model testing query, according to one embodiment.
FIG. 7B illustrates an example confusion matrix generated during a differentially private model testing query.
FIG. 8 illustrates a system-level modification to the system of FIG. 1 that allows the client to access to a differentially private synthetic database, according to one embodiment.
FIG. 9 illustrates the application of a clustering query to entries of a differentially private synthetic database, according to one embodiment.
FIG. 10 illustrates a process for responding to a database query by executing a differentially private version of the query on the database, according to one embodiment.
FIG. 11 is a block diagram illustrating components of an example machine able to read instructions from a machine-readable medium and execute them in a processor (or controller).



DETAILED DESCRIPTION
The Figures (FIGS.) and the following description describe certain embodiments by way of illustration only. One skilled in the art will readily recognize from the following description that alternative embodiments of the structures and methods illustrated herein may be employed without departing from the principles described herein. Reference will now be made in detail to several embodiments, examples of which are illustrated in the accompanying figures. It is noted that wherever practicable similar or like reference numbers may be used in the figures and may indicate similar or like functionality.
System Overview
FIG. 1 is a system 100 for receiving a query 108 for a private database 106, and responding to the query 108 by executing a differentially private (DP) version of the query 114 on the private database 106. The system 100 includes a differentially private security system 102 that receives the analytical query 108 from a client 104 and applies a DP version of the query 114 on the database 106. Subsequently, the differentially private security system 102 returns the response of the DP query 114 to the client 104 as the DP response 112.
The database 106 is one or more private databases managed by one or more entities that can only be accessed by authorized or trusted users. For example, the database 106 may contain health data of patients, financial records, telecom data, and confidential business intelligence, such as proprietary data or data restricted by contractual obligations. The information stored in the database 106 is of interest to one or more clients 104, but clients 104 may not have the necessary authorization to access to information contained in the databases 106.
FIG. 2 illustrates an example database structure, according to one embodiment. For the remainder of the application, a database, including one or more of the private databases 106, may be referred to as a matrix with a number of rows and columns. Each row is an entry of the database and each column is a feature of the database. Thus, each row contains a data entry characterized by a series of feature values for the data entry. For example, as shown in FIG. 2, the example database 200 contains 8 entries and 11 features, and illustrates a list of patient profiles. Each patient is characterized by a series of feature values that contain information on the patient's height (Feature 1), country of residence (Feature 2), age (Feature 10), and whether the patient has contracted a disease (Feature 11).
The feature values may be numerical in nature, e.g., Features 1 and 10, or categorical in nature, e.g., Features 2 and 11. In the case of categorical feature values, each category may be denoted as an integer. For example, in Feature 11 of FIG. 2, “0” indicates that the patient has not contracted a disease, and “1” indicates that the patient has contracted a disease.
Returning to FIG. 1, the client 104 may be a human analyst or an organization that does not have direct access to the database 106, but is interested in applying an analytical query 108 to the database 106. For example, the client 104 may be a data analyst, data scientist, or a health analyst that is interested in the profiles of the patients but does not have direct access to the database 106. Each client 104 of the system 100 is associated with a privacy budget and specifies a set of privacy parameters each time the client 104 submits a query 108. The privacy budget is a numerical value representative of a number and/or type of remaining queries 108 available to the client 104 in terms of the privacy parameters specified for each query 108.
The query 108 submitted by the client 104 may be simple queries, such as count queries that request the number of entries in the databases 106 that satisfy a condition specified by the client 104, or complicated queries, such as predictive analytics queries that request a data analytics model trained on the databases 106. Upon submitting a query 108 to the differentially private security system 102, the client 104 receives a DP response 112 to a differentially private version of the submitted query 114.
The client 104 specifies a set of privacy parameters each time the client 104 submits query 108. The privacy parameters indicate an amount of decrease in the privacy budget of the client 104 in return for a response to the query 108. As described below in more detail with reference to the privacy system 160 in FIG. 3, the privacy parameters specified by the client 104 also indicate the amount of information released about the database 106 to the client 104.
The differentially private security system 102 receives an analytical query 108 from the client 104 and applies a differentially private version of the query 114 on the database 106, such that it releases a degree of information about the database 106 indicated by the privacy parameters specified by the client 104, but also protects a degree of privacy of the databases 106 specified by the entities managing the database 106. For example, the entities managing the database 106 may also set a maximum threshold on the degree of information released about the database 106 for a given query 108 that the client 104 may not exceed. Thus, the differentially private security system balances privacy protection of the database 106 while releasing useful information on the database 106 to the client 104. The differentially private security system 102 may have complete or partial access to the databases 106.
Upon receiving a query 108, the differentially private security system 102 applies DP query 114 to the database 106 and returns a DP response 112 to the client 104. The DP query 114 is a differentially private version of the query 108 that satisfies a definition of differential privacy described in more detail with reference to the privacy system 160 in FIG. 3. The DP query 114 may include perturbing the response or output of the query 108 with noise, or perturbing the process for generating the output of the query 108 with noise. The resulting output of the DP query 114 is returned to the client 104 as DP response 112. Ideally, the DP response 112 correlates to the original output of the query 108 on the databases 106 but maintains the degree of privacy specified by the entities managing the database 106.
Differentially Private Security System
The differentially private security system 102 includes a user interface 150, a library 152, an account management system 154, a query handling engine 156, a data integration module 158, and a privacy system 160. Some embodiments of the differentially private security system 102 have different or additional modules than the ones described here. Similarly, the functions can be distributed among the modules in a different manner than is described here. Certain modules and functions can be incorporated into other modules of the differentially private security system 102.
The user interface 150 can generate a graphical user interface on a dedicated hardware device of the differentially private security system 102 or the client 104 in which the client 104 can submit an analytical query 108 and the desired privacy parameters, and view DP response 112 in the form of numerical values or images. The client 104 may also inspect database 106 schemata, view an associated privacy budget, or cache the DP response 112 to view the response later. The user interface 150 submits properly formatted query commands to other modules of the differentially private security system 102.
The library 152 contains software components that can be included in external programs that allow the client 104 to submit the analytical query 108, receive the DP response 112, and other functions within a script or program. For example, the client 104 may use the software components of the library 152 to construct custom data analytic programs. Each of the software components in the library 152 submits properly formatted query commands to other modules of the differentially private security system 102.
The account management system 154 receives properly formatted query commands (herein “query commands” or “QC”), parses the received query commands, and updates the account of the client 104 according to the received query command. For example, the account management system 154 may check the query commands for syntactic correctness, or check whether a client 104 has access to a requested resource. As another example, the account management system 154 may check whether the privacy parameters specified by the client 104 for a given analytical query 108 can be accommodated, and if so, decrement the privacy budget of the client 104 by the amount specified in the query 108. Query commands verified by the account management system 154 are provided to the query handling engine 156. Examples of query commands accommodated by the differentially private security system 102 are listed below.












QC1. Count


‘SELECT COUNT (<column>) FROM <database.table> WHERE <where_clause> BUDGET


<eps> <delta>.


QC2. Median


‘SELECT MEDIAN (<column>) FROM <database.table> WHERE <where_clause> BUDGET


<eps> <delta>.


QC3. Mean


‘SELECT MEAN (<column>) FROM <database.table> WHERE <where_clause> BUDGET


<eps> <delta>.


QC4. Variance


‘SELECT VARIANCE (<column>) FROM <database.table> WHERE <where_clause>


BUDGET <eps> <delta>.


QC5. Inter-Quartile Range


‘SELECT IQR (<column>) FROM <database.table> WHERE <where_clause> BUDGET <eps>


<delta>.


QC6. Batch Gradient Descent


‘SELECT <GLM> (<columns_x>,<column_y>,<params>) FROM <database.table> WHERE


<where_clause> BUDGET <eps> <delta>.


QC7. Stochastic Gradient Descent


‘SELECT SGD <GLM> (<column>) FROM <database.table> WHERE <where_clause>


BUDGET <eps> <delta>.


QC8. Random Forest


‘SELECT RANDOMFOREST (<columns_x>,<columns_y>) FROM <database.table> WHERE


<where_clause> BUDGET <eps> <delta>.


QC9. Histogram


‘SELECT HISTOGRAM (<column>) FROM <database.table> WHERE <where_clause_i>


BUDGET <eps> <delta>.









The query handling engine 156 transforms the received query commands into appropriate function calls and database access commands by parsing the query command string. The function calls are specific to the query 108 requested by the client 104, and the access commands allow access to the required database 106. Different databases 106 require different access commands. The access commands are provided to the database integrator 158.
The database integrator 158 receives the access commands to one or more databases 106 and collects the required databases 106 and merges them into a single data object. The data object has a structure similar to that of a database structure described in reference to FIG. 2. The data object is provided to the privacy system 160.
The privacy system 160 receives the data object from the database integrator 158, appropriate function calls from the query handling engine 156 indicating the type of query 108 submitted by the client 104, privacy parameters specified for the query 108, and produces a DP response 112 to a differentially private version of the query 108 with respect to the databases 106. The privacy system 160 will be described in further detail in reference to FIG. 3 below.
Privacy System
FIG. 3 is a block diagram illustrating the privacy system 160 of the system 100 shown in FIG. 1, according to one embodiment. The privacy system 160 includes a count engine 302, a mean engine 304, a median engine 306, a variance engine 308, an IQR engine 310, a batch gradient engine 312, a stochastic gradient engine 314, a random forest engine 316, a histogram engine 318, a model testing engine 320, and a synthetic database engine 322. Some embodiments of the privacy system 160 have different or additional modules than the ones described here. Similarly, the functions can be distributed among the modules in a different manner than is described here. Certain modules and functions can be incorporated into other modules of the privacy system 160.
Definition of Differential Privacy
For a given query 108, the privacy system 160 receives a data object X, function calls indicating the type of query 108, privacy parameters specified by the client 104, and outputs a DP response 112 to a differentially private version of the query 108 with respect to X. Each data object X is a collection of row vectors xi=1, 2, . . . n, in which each row vector x; has a series of p elements xij=2, . . . , p.
A query M satisfies the definition of e-differential privacy if for all:




  
    ∀
    X
  
  ,
  
    
      X
      ′
    
    ∈
    D
  
  ,
  
    ∀
    
      S
      ⊆
      
        
          
            Range
            ⁡
            
              (
              M
              )
            
          
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            
              Pr
              ⁡
              
                [
                
                  
                    M
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                  ∈
                  S
                
                ]
              
            
            
              Pr
              ⁡
              
                [
                
                  
                    M
                    ⁡
                    
                      (
                      
                        X
                        ′
                      
                      )
                    
                  
                  ∈
                  S
                
                ]
              
            
          
        
        ≤
        
          e
          ɛ
        
      
    
  




where is the space of all possible data objects, X, X′ neighboring data objects, S is an output space of query M, and neighboring databases are defined as two data objects X, X′ that have at most one different entry from one another. That is, given two neighboring data objects X, X′ in which one has an individual's data entry, and the other does not, there is no output of query M that an adversary can use to distinguish between X, X′. That is, an output of such a query M that is differentially private reveals no information about the data object X. The privacy parameter ε controls the amount of information that the query M reveals about any individual data entry in X, and represents the degree of information released about the entries in X. For example, in the definition given above, a small value of ε indicates that the probability an output of query M will disclose information on a specific data entry is small, while a large value of ε indicates the opposite.

As another definition of differential privacy, a query M is (ε,δ)-differentially private if for neighboring data objects X, X′:




  
    ∀
    X
  
  ,
  
    
      X
      ′
    
    ∈
    D
  
  ,
  
    ∀
    
      S
      ⊆
      
        
          
            Range
            ⁡
            
              (
              M
              )
            
          
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            
              Pr
              ⁡
              
                [
                
                  
                    M
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                  ∈
                  S
                
                ]
              
            
            
              Pr
              ⁡
              
                [
                
                  
                    M
                    ⁡
                    
                      (
                      
                        X
                        ′
                      
                      )
                    
                  
                  ∈
                  S
                
                ]
              
            
          
        
        ≤
        
          
            e
            ɛ
          
          +
          
            δ
            .
          
        
      
    
  




The privacy parameter δ measures the improbability of the output of query M satisfying ε-differential privacy. As discussed in reference to FIG. 1, the client 104 may specify the desired values for the privacy parameters (ε,δ) for a query 108.
There are three important definitions for discussing the privacy system 160: global sensitivity, local sensitivity, and smooth sensitivity. Global sensitivity of a query M is defined as




  
    
      GS
      M
    
    ⁡
    
      (
      X
      )
    
  
  =
  
    
      max
      
        X
        ,
        
          
            
              X
              ′
            
            ⁢
            
              :
            
            ⁢
            
              d
              ⁡
              
                (
                
                  X
                  ,
                  
                    X
                    ′
                  
                
                )
              
            
          
          =
          1
        
      
    
    ⁢
    
      
      
        
          M
          ⁡
          
            (
            X
            )
          
        
        -
        
          M
          ⁡
          
            (
            
              X
              ′
            
            )
          
        
      
      
    
  




where X, X′ are any neighboring data objects, such that d(X, X′)=1. This states that the global sensitivity is the most the output of query M could change by computing M on X and X′.

The local sensitivity of a query M on the data object X is given by:




  
    
      LS
      M
    
    ⁡
    
      (
      X
      )
    
  
  =
  
    
      max
      
        
          
            X
            ′
          
          ⁢
          
            :
          
          ⁢
          
            d
            ⁡
            
              (
              
                X
                ,
                
                  X
                  ′
                
              
              )
            
          
        
        =
        1
      
    
    ⁢
    
      
      
        
          M
          ⁡
          
            (
            X
            )
          
        
        -
        
          M
          ⁡
          
            (
            
              X
              ′
            
            )
          
        
      
      
    
  




where the set {X′: d(X, X′)=1} denotes all data objects that have at most one entry that is different from X. That is, the local sensitivity LSM(X) is the sensitivity of the output of the query M on data objects X′ that have at most one different entry from X, measured by a norm function.

Related to the local sensitivity LSM(X), the smooth sensitivity given a parameter/p is given by:




  
    
      S
      M
    
    ⁡
    
      (
      
        X
        ;
        β
      
      )
    
  
  =
  
    
      max
      
        
          X
          ⁢
          
            
          
          ⁢
          ′
        
        ∈
        D
      
    
    ⁢
    
      
      
        
          
            LS
            M
          
          ⁡
          
            (
            X
            )
          
        
        ·
        
          e
          
            
              -
              β
            
            ·
            
              d
              ⁡
              
                (
                
                  X
                  ,
                  
                    X
                    ′
                  
                
                )
              
            
          
        
      
      
    
  




where d(X, X′) denotes the number of entries that differ between X and X′.

Notation for Random Variables

The notation in this section is used for the remainder of the application to denote the following random variables.
1) G(σ2), denotes a zero-centered Gaussian random variable with the probability density function




  
    f
    ⁡
    
      (
      
        x
        ❘
        
          σ
          2
        
      
      )
    
  
  =
  
    
      1
      
        σ
        ⁢
        
          
            2
            ⁢
            π
          
        
      
    
    ⁢
    
      
        e
        
          -
          
            
              x
              2
            
            
              2
              ⁢
              
                σ
                2
              
            
          
        
      
      .
    
  




2) L(b) denotes a zero-centered Laplacian random variable with the probability density function





  
    f
    ⁡
    
      (
      
        x
        ❘
        b
      
      )
    
  
  =
  
    
      1
      
        2
        ⁢
        b
      
    
    ⁢
    
      
        e
        
          -
          
            
              
              x
              
            
            b
          
        
      
      .
    
  




3) C(γ) denotes a zero-centered Cauchy random variable with the probability density function





  
    f
    ⁡
    
      (
      
        x
        |
        γ
      
      )
    
  
  =
  
    
      1
      
        π
        ⁢
        
          γ
          ⁡
          
            (
            
              1
              +
              
                
                  (
                  
                    x
                    γ
                  
                  )
                
                2
              
            
            )
          
        
      
    
    .
  




Further, a vector populated with random variables R as its elements is denoted by ν(R). A matrix populated with random variables R as its elements is denoted by M(R).
Count Engine 302
The count engine 302 produces a DP response 112 responsive to the differentially private security system 102 receiving a query 108 for counting the number of entries in a column of the data object X that satisfy a condition specified by the client 104, given privacy parameters (ε,δ). An example query command for accessing the count engine 302 is given in QC1 above. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for the number of patients that are above the age of 30.
The count engine 302 retrieves the count q from X. If privacy parameter δ is equal to zero, the count engine 302 returns




  
    y
    ≈
    
      q
      +
      
        L
        ⁡
        
          (
          
            
              c
              1
            
            ·
            
              1
              ε
            
          
          )
        
      
    
  
  ,




as the DP response 112 for display on the user interface 150, where c1 is a constant. An example value for c1 may be 1. If the privacy parameter δ is non-zero, the count engine 302 returns





  
    y
    ≈
    
      q
      +
      
        G
        ⁡
        
          (
          
            
              
                c
                1
              
              ·
              2
              ·
              log
            
            ⁢
            
              
                2
                δ
              
              ·
              
                1
                
                  ε
                  2
                
              
            
          
          )
        
      
    
  
  ,




as the DP response 112 for display on the user interface 150, where c1 is a constant. An example value for c1 may be 1.

The client 104 may request visualization of entries in the data object X for analysis of trends or patterns that depend on the features of the entries. In one embodiment, the privacy system 160 generates a differentially private visualization of the requested data entries from X. FIG. 4 illustrates displaying results of a differentially private count query to the user interface of the client, according to one embodiment.
The privacy system 160 first maps the requested entries from X for the selected features specified by the client 104. For example, as shown in the visualization 410 of FIG. 4, a series of requested entries are plotted depending on their values for Feature 1 and Feature 2. The privacy system 160 then generates disjoint regions on the plot and retrieves the counts of entries in each of the disjoint regions. In visualization 410, the privacy system 160 divides the plot into disjoint squares and retrieves the count of entries in each square.
For each disjoint region, the privacy system 160 submits a differentially private count query to the count engine 302, and randomly plots a number of entries determined by the DP response 112 of the count engine 302 for that region. The resulting DP visualization plot is returned to the client 104 for display to a user by the user interface 150. For example, square 440 in visualization 410 contains 3 entries, while the same square in DP visualization 420 contains 4 randomly plotted entries determined by the DP response 112 of the count engine 302.
Median Engine 304
The median engine 304 produces a DP response 112 responsive to the differentially private security system 102 receiving a query 108 for generating the median of entries in a column of the data object X that satisfy a condition specified by the client 104, given privacy parameters (ε,δ). An example query command for accessing the median engine 304 is given in QC2 above. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for the median age of all patients in X.
The median engine 304 aggregates the values of entries satisfying the condition specified by the client 104 into a list U, and retrieves the median q from U. If privacy parameter δ is equal to zero, the median engine 304 returns




  
    
      
        y
        ≈
        
          q
          +
          
            
              c
              1
            
            ·
            
              
                S
                M
              
              ⁡
              
                (
                
                  U
                  ;
                  
                    
                      c
                      2
                    
                    ·
                    ε
                  
                
                )
              
            
            ·
            
              
                C
                ⁡
                
                  (
                  1
                  )
                
              
              ε
            
          
        
      
    
    
      
        
      
    
  




as the DP response 112 for display on the user interface 150, in which c1, c2 are constant factors. Example values for c1, c2 may be 6 and 1⁄6, respectively. If δ is non-zero, the median engine 304 returns





  
    
      
        y
        ≈
        
          q
          +
          
            
              c
              1
            
            ·
            
              
                S
                M
              
              (
              
                U
                ;
                
                  
                    c
                    2
                  
                  ·
                  
                    ε
                    
                      
                        2
                        ·
                        log
                      
                      ⁢
                      
                        1
                        δ
                      
                    
                  
                
              
              )
            
            ·
            
              
                L
                ⁡
                
                  (
                  1
                  )
                
              
              ε
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




as the DP response 112 for display on the user interface 150. Example values for c1, c2 may be 2 and 1, respectively.

Mean Engine 306

The mean engine 306 produces a DP response 112 responsive the differentially private security system 102 receiving a query 108 for generating the mean of entries in a column of the data object X that satisfy a condition specified by the client 104, given privacy parameters (ε,δ). An example query command for accessing the mean engine 306 is given in QC3 above. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for generating the mean age of patients that are above the age of 30.
The mean engine 306 aggregates the values of entries satisfying the condition specified by the client 104 into a list U. Assuming there are n values in U, the mean engine 306 further divides U into m sub-lists Vj=1, 2, . . . , m each with n m values. The mean engine 306 aggregates each mean rj of sub-list Vj into a list R. The mean engine 306 requests a differentially private median query of the values in R to the median engine 304. The resulting output from the median engine 304 is returned as the DP response 112 for display on the user interface 150.
Variance Engine 308
The variance engine 308 produces a DP response 112 responsive to the differentially private security system 102 receiving a query 108 for generating the variance of entries in a column of the data object X that satisfy a condition specified by the client 104, given privacy parameters (ε,δ). An example query command for accessing the variance engine 308 is given in QC4 above. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for generating the variance of the age of all patients in X.
The variance engine 308 aggregates the values of entries satisfying the condition specified by the client 104 into a list U. Assuming there are n values in U, the variance engine 308 further divides U into m sub-lists Vj=1, 2, . . . , m each with n m values. The variance engine 308 aggregates each variance r; of sub-list Vj into a list R. The variance engine 308 requests a differentially private median query of the values in R to the median engine 304. The resulting output from the median engine 304 is returned as the DP response 112 for display on the user interface 150.
IQR Engine 310
The IQR engine 310 produces a DP response 112 responsive to the differentially private security system 102 receiving a query 108 for generating the interquartile range (IQR) of entries in a column of the data object X that satisfy a condition specified by the client 104, given privacy parameters (ε,δ). An example query command for accessing the IQR engine 310 is given in QC5 above. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for generating the IQR of the age of all patients in X.
In one embodiment, the IQR engine 310 aggregates the values of entries satisfying the condition specified by the client 104 into a list U. Assuming there are n values in U, the sample IQR of U is denoted as IQR(U), and a log transform of IQR(U) is denoted as:




  
    
      H
      n
    
    ⁡
    
      (
      U
      )
    
  
  =
  
    
      log
      
        1
        +
        
          1
          
            l
            ⁢
            o
            ⁢
            g
            ⁢
            n
          
        
      
    
    ⁢
    IQ
    ⁢
    
      
        R
        ⁡
        
          (
          U
          )
        
      
      .
    
  




The IQR engine 310 further maps the quantity Hn(U) to an integer k0 such that Hn(U)∈[k0, k0+1). The IQR engine 310 extracts a value A0(U) indicating the number of entries in U required to change in order for the new list U to satisfy Hn(Ũ) ∉[k0, k0+1).

The IQR engine 310 then generates a value R0(U) given by:




  
    
      R
      0
    
    ⁡
    
      (
      U
      )
    
  
  ≈
  
    
      
        A
        0
      
      ⁡
      
        (
        U
        )
      
    
    +
    
      L
      ⁡
      
        (
        
          
            c
            1
          
          ε
        
        )
      
    
  




in which c1 is a constant factor. If R0(U) is greater than a predetermined threshold, the IQR engine 310 returns





  
    y
    =
    
      IQ
      ⁢
      
        
          R
          ⁡
          
            (
            U
            )
          
        
        ·
        
          
            (
            
              1
              
                1
                +
                
                  log
                  ⁢
                  
                    
                  
                  ⁢
                  n
                
              
            
            )
          
          
            L
            ⁡
            
              (
              
                
                  c
                  1
                
                ε
              
              )
            
          
        
      
    
  
  ,




as the DP response 112 for display on the user interface 150. If R0(U) is equal to or less than the predetermined threshold, the IQR engine 310 returns “No Answer” as the DP response 112 for display on the user interface 150.

In another embodiment, the IQR engine 310 aggregates the values of entries satisfying the condition specified by the client 104 into an ordered list U. The IQR engine 310 retrieves the first quartile and the third quartile from U, given by q and q′, respectively. If δ is zero, the IQR engine 310 returns:




  
    
      
        y
        ≈
        
          
            (
            
              q
              +
              
                
                  c
                  1
                
                ·
                
                  
                    S
                    M
                  
                  ⁡
                  
                    (
                    
                      U
                      ;
                      
                        
                          c
                          2
                        
                        ·
                        ε
                      
                    
                    )
                  
                
                ·
                
                  
                    C
                    ⁡
                    
                      (
                      1
                      )
                    
                  
                  
                    ε
                    2
                  
                
              
            
            )
          
          -
          
            (
            
              
                q
                ′
              
              +
              
                
                  c
                  1
                
                ·
                
                  
                    S
                    M
                  
                  ⁡
                  
                    (
                    
                      U
                      ;
                      
                        
                          c
                          2
                        
                        ·
                        ε
                      
                    
                    )
                  
                
                ·
                
                  
                    C
                    ⁡
                    
                      (
                      1
                      )
                    
                  
                  
                    ε
                    /
                    2
                  
                
              
            
            )
          
        
      
    
    
      
        
      
    
  




as the DP response 112 for display on the user interface 150, in which c1, c2 are constant factors.

If δ is non-zero, the IQR engine 310 returns:




  
    
      
        y
        ≈
        
          
            (
            
              q
              +
              
                
                  c
                  1
                
                ·
                
                  
                    S
                    M
                  
                  (
                  
                    U
                    ;
                    
                      
                        c
                        2
                      
                      ·
                      
                        ε
                        
                          
                            2
                            ·
                            log
                          
                          ⁢
                          
                            1
                            δ
                          
                        
                      
                    
                  
                  )
                
                ·
                
                  
                    L
                    ⁡
                    
                      (
                      1
                      )
                    
                  
                  
                    ε
                    /
                    2
                  
                
              
            
            )
          
          -
          
            (
            
              
                q
                ′
              
              +
              
                
                  c
                  1
                
                ·
                
                  
                    S
                    M
                  
                  (
                  
                    U
                    ;
                    
                      
                        c
                        2
                      
                      ·
                      
                        ε
                        
                          
                            2
                            ·
                            log
                          
                          ⁢
                          
                            1
                            δ
                          
                        
                      
                    
                  
                  )
                
                ·
                
                  
                    L
                    ⁡
                    
                      (
                      1
                      )
                    
                  
                  
                    ε
                    /
                    2
                  
                
              
            
            )
          
        
      
    
    
      
        
      
    
  




as the DP response 112 for display on the user interface 150, in which c1, c2 are constant factors.

Batch Gradient Engine 312

The batch gradient engine 312 produces a DP response 112 responsive to the differentially private security system 102 receiving a valid query 108 for generating a set of parameters θ for a general linear model that captures the correlation between a series of observable features and a dependent feature, given privacy parameters (ε,δ). The general linear model is trained on the selected columns of X. An example query command for accessing the batch gradient engine 312 is given in QC6 above.
Given a row vector x that contains a series of observable features and a label feature y, the correlation between the observable features and the label feature in a general linear model may be given as:

y=xθT,

where θ is a row vector containing parameters of the model. That is, the label feature is modeled as a weighted sum of the observable features, where each value in θ is the weight given to a corresponding observable feature.

For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for generating a set of parameters θ for a general linear model that captures the correlation between the height of the patients (observable feature) and the age of the patients (label feature). As another example, the features may be categorical in nature, and the requested parameters θ may be for a general linear model that captures the correlation between the height, age, residence of the patients (observable features) and whether the patient will or has contracted a disease (label feature).
Examples of general linear models supported by the batch gradient engine 312 are, but not limited to, linear regression, logistic regression, and support vector machine (SVM) classifiers.
The optimal values for the set of parameters θ is found by training the general linear model on training data (Xtrain, ytrain) consisting of selected columns of data object X. Specifically, Xtrain is a matrix database in which each column corresponds to a selected observable feature, and y is a column vector of the selected label feature values. Each entry in Xtrain has a one-to-one correspondence with an entry in y. The optimal θ is generally found by minimizing a loss function on (Xtrain, ytrain) over possible values of θ. Mathematically, the minimization is given by:




  θ
  =
  
    
      argmin
      θ
    
    ⁢
    
      
    
    ⁢
    l
    ⁢
    
      
    
    ⁢
    
      
        (
        
          
            X
            
              t
              ⁢
              r
              ⁢
              a
              ⁢
              i
              ⁢
              n
            
          
          ,
          
            
          
          ⁢
          
            
              y
              
                t
                ⁢
                r
                ⁢
                a
                ⁢
                i
                ⁢
                n
              
            
            ;
            θ
          
        
        )
      
      .
    
  




The batch gradient engine 312 returns a DP response 112 θDP of a differentially private batch gradient query by perturbing the loss function to be minimized. Specifically, the perturbed minimization is given by:




  
    
      θ
      
        D
        ⁢
        P
      
    
    =
    
      
        
          argmin
          θ
        
        ⁢
        
          
        
        ⁢
        
          l
          ⁡
          
            (
            
              
                X
                
                  t
                  ⁢
                  r
                  ⁢
                  ain
                
              
              ,
              
                
                  y
                  
                    t
                    ⁢
                    r
                    ⁢
                    a
                    ⁢
                    i
                    ⁢
                    n
                  
                
                ;
                θ
              
            
            )
          
        
      
      +
      
        
          θ
          T
        
        ⁢
        
          v
          (
          
            G
            (
            
              
                4
                ·
                
                  K
                  2
                
                ·
                
                  R
                  2
                  2
                
                ·
                
                  (
                  
                    
                      l
                      ⁢
                      
                        
                      
                      ⁢
                      og
                      ⁢
                      
                        1
                        δ
                      
                    
                    +
                    ε
                  
                  )
                
              
              
                ε
                2
              
            
            )
          
          )
        
      
    
  
  ,




in which K is the Lipschitz constant for loss function (⋅). If j is the index of the columns in Xtrain, xij denotes the value of entry i and column j in Xtrain, and it is publicly known that for each column j, aj≤xij≤bj, R2 may be given by:

R2=max(∥u∥2|aj≤uj≤bj)

where a is a vector having elements uj. The DP response 112 θDP may be provided for display on the user interface 150.

Stochastic Gradient Engine 314

Similarly to the batch gradient engine 312, the stochastic gradient engine 314 produces a DP response 112 responsive to the differentially private security system 102 receiving a valid query 108 for generating a set of parameters θ for a general linear model that captures the correlation between a series of observable features and a label feature, given privacy parameters (ε,δ). An example query command for accessing the stochastic gradient engine 314 is given in QC7 above.
Similar to the batch gradient engine 312, examples of general linear models supported by the stochastic gradient engine 314 are, but not limited to, linear regression, logistic regression, and support vector machine (SVM) classifiers.
The stochastic gradient engine 314 also minimizes a loss function on training data (Xtrain, ytrain) over possible values of 0 to find the optimal vales of parameter vector θ. However, the stochastic gradient engine 314 may minimize the loss function based on individual points or a subset of the training data, instead of the entire training data.
As discussed in reference to the batch gradient engine 312, a general minimization problem for finding the optimal values for θ over training data (Xtrain, ytrain) is given by:




  θ
  =
  
    
      argmin
      θ
    
    ⁢
    
      
    
    ⁢
    
      l
      ⁡
      
        (
        
          
            X
            
              t
              ⁢
              r
              ⁢
              a
              ⁢
              i
              ⁢
              n
            
          
          ,
          
            
              y
              
                t
                ⁢
                r
                ⁢
                a
                ⁢
                i
                ⁢
                n
              
            
            ;
            θ
          
        
        )
      
    
  




where (⋅) is a loss function. The minimization is solved by applying stochastic gradient descent on the loss function (⋅) with respect to θ. This involves the steps of identifying an initial set of values for θ, calculating the gradient of the loss function with respect to θ, and updating θ based on the calculated gradient. The steps are repeated until the algorithm reaches convergence, and an optimal set of values for θ that minimize the loss function are identified.

Specifically, given the estimate for the parameter θt at time t, stochastic gradient descent generates a new estimate θt+1 at the next time step t+1 by the following equation:

θt+1θt−ηt·n·∇θt(Xtrain,ytrain;θ),

in which ∇θt(Xtrain, ytrain; θ) is the gradient of the loss function with respect to θ, and ηt is the learning rate. The algorithm is repeated until the estimate for θ converges.

The stochastic gradient engine 314 returns a DP response 112 θDP of a differentially private stochastic gradient query by perturbing the update of θ at one or more time steps of the stochastic gradient descent algorithm. Specifically, a perturbed update at time t to t+1 is given by:




  
    θ
    
      t
      +
      1
    
  
  =
  
    
      θ
      t
    
    -
    
      
        η
        t
      
      ·
      n
      ·
      
        
          ∇
          
            θ
            t
          
        
        ⁢
        
          l
          ⁡
          
            (
            
              
                X
                
                  t
                  ⁢
                  r
                  ⁢
                  ain
                
              
              ,
              
                
                  y
                  
                    t
                    ⁢
                    r
                    ⁢
                    a
                    ⁢
                    i
                    ⁢
                    n
                  
                
                ;
                θ
              
            
            )
          
        
      
    
    -
    
      
        η
        t
      
      ·
      
        v
        (
        
          G
          (
          
            
              
                c
                1
                2
              
              ·
              
                n
                2
              
              ·
              
                (
                
                  log
                  ⁢
                  
                    
                      n
                      δ
                    
                    ·
                    log
                  
                  ⁢
                  
                    1
                    δ
                  
                
                )
              
            
            
              ε
              4
            
          
          )
        
        )
      
    
  




    
    
        where ηt is the learning rate.
    
    


The stochastic gradient engine 314 may output the perturbed update at each time step as the DP response 112 for display on the user interface 150, or the converged parameter vector θDP as the DP response 112 for display on the user interface 150.
Random Forest Engine 316
The random forest engine 316 produces a DP response 112 responsive to the differentially private security system 102 receiving a valid query 108 for generating a trained random forest classifier that bins a series of feature values into one among multiple categories, given privacy parameters (ε,δ). The random forest classifier is trained on the selected columns of X. An example query command for accessing the random forest engine 316 is given in QC8 above. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for generating a trained random forest classifier that receives values for the height and age of a patient and determines whether the patient has contracted the disease or not.
The random forest classifier, is trained on training data (Xtrain, ytrain) to learn the correlation between the selected features of an entry and the category the entry belongs to. Specifically, Xtrain is a matrix database in which each column corresponds to a selected feature of interest to the client 104, and y is a column vector of already known labels indicating the category of a corresponding entry. Each entry in Xtrain has a one-to-one correspondence with a label entry in y. Upon being trained, the random forest classifier, or a classifier in general, receives a new data entry with selected feature values and generates an estimate of the category for the new entry.
The random forest classifier is an ensemble of individual binary decision tree classifiers, in which each binary decision tree generates an estimate for the category of an entry. Given a new data entry, the random forest classifier aggregates the category estimates from each binary decision tree and produces a final estimate for the category of the data entry.
FIG. 5 is an example diagram of a trained binary decision tree, according to one embodiment. Each decision tree includes a hierarchical structure with a plurality of T nodes ti=1, 2, . . . , T and a plurality of directed edges between a parent node and a child node. A parent node is a node that has outgoing edges to a child node, and a child node is a node that receives edges from a parent node. In the particular embodiment of a binary decision tree, each parent node has two child nodes. The nodes are one among a root node, in which the node has no incoming edges, an internal node, in which the node has one incoming edge with two outgoing edges, and a leaf node, in which the node has one incoming edge with no outgoing edges. For example, the example decision tree in FIG. 5 has seven nodes t1, t2, . . . , t7 and six edges. t1 is the root node, t2 and t3 are internal nodes, and t4-t7 are leaf nodes.
For each trained binary decision tree, each node except the root node corresponds to a partition of training data entries formed by a split s at a parent node. The split s at the parent node is based on a test condition of a feature of the training data (Xtrain, ytrain) that compares the feature value of an entry to a reference value, and verifies whether the feature value meets that condition or not. Returning to the example shown in FIG. 5, node t1 creates a split si into two child nodes t2 and t3 based on the test condition x5≥5, which checks if an entry contains a fifth feature value equal to or greater than 5. The training data (Xtrain, ytrain) is thus split at si into one partition that contains entries with x5≥5, and another partition that contains entries with x5<5. The former partition is directed to child node t1 and the latter partition is directed to child node t2. The partitioning process is repeated until the leaf nodes of the binary decision tree are determined.
At the end of the training process, each leaf node is associated with a category that has a dominant proportion in the corresponding partition at the leaf node. In FIG. 5, leaf node t4 is assigned label “1,” since the proportion of “1” labels in leaf node t4, denoted by p(1|t4), is greater than the proportion of “0” labels in leaf node t4, denoted by p(0|t4). Given a new data entry with an unknown category, the trained decision tree generates a label estimate by checking the appropriate feature values of the data entry at each node as it propagates through the tree to a destination leaf node. Upon arriving at the leaf node, the data entry is assigned the category label associated with the leaf node.
The random forest engine 316 returns a DP response 112 of a differentially private random forest query by perturbing the proportion of training data entries at leaf nodes of each trained binary decision tree. Specifically, the random forest engine 316 trains a random forest classifier T with an ensemble of Ntrees binary decision trees Bj=1, 2, . . . , Ntrees using training data (Xtrain, ytrain) from the data object X. Assuming a binary classification problem with two labels “0” and “1,” the random forest engine 316 perturbs the proportion of data entries associated with each category for each leaf node tL. The perturbed proportion pDP(⋅|tL) is given by:




  
    
      
        p
        DP
      
      ⁡
      
        (
        
          0
          |
          
            t
            L
          
        
        )
      
    
    ≈
    
      
        p
        ⁡
        
          (
          
            0
            |
            
              t
              L
            
          
          )
        
      
      +
      
        L
        ⁡
        
          (
          
            ε
            
              l
              ⁢
              o
              ⁢
              g
              ⁢
              
                N
                trees
              
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        p
        
          D
          ⁢
          P
        
      
      ⁡
      
        (
        
          1
          |
          
            t
            L
          
        
        )
      
    
    ≈
    
      
        p
        ⁡
        
          (
          
            1
            |
            
              t
              L
            
          
          )
        
      
      +
      
        
          L
          ⁡
          
            (
            
              ε
              
                l
                ⁢
                o
                ⁢
                g
                ⁢
                
                  N
                  trees
                
              
            
            )
          
        
        .
      
    
  




The random forest engine 316 returns the random forest classifier TDP containing an ensemble of perturbed binary decision trees BDPj=1, 2, . . . , Ntrees as the DP response 112. Moreover, the random forest engine 316 may display the perturbed proportion of data entries for leaf nodes of each binary decision tree BDPj=1, 2, . . . , Ntrees for display on the user interface 150.
Histogram Engine 318
The histogram engine 318 produces a DP response 112 responsive to the differentially private security system 102 receiving a query 108 for generating a histogram of a selected column in X, given privacy parameters (ε,δ). The histogram engine 318 creates one or more bins corresponding to one or more disjoint ranges of the selected feature values, and indicates the number or proportion of entries that belong to each bin. An example query command for accessing the histogram engine 318 is given in QC9 above. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for generating a histogram of the age of all patients for bins age 0-10, 11-20, 21-30, and so on.
The histogram engine 318 returns a DP response 112 of a differentially private histogram query by perturbing the counts for each bin.
In one embodiment, the histogram engine 318 generates the requested histogram from the selected column of X, and perturbs the counts of each bin by submitting a request to the count engine 302. FIG. 6 illustrates perturbing the counts for a differentially private histogram query, according to one embodiment. As shown in FIG. 6, the histogram engine 318 generates histogram 600 by counting the number of entries corresponding to each bin bi=1, 2, . . . , B. The histogram engine 318 then requests the count engine 302 to perturb the counts qi for each bin to generate a perturbed histogram 602. As shown in FIG. 6, the count 608 for bin b5 in the perturbed histogram 602 is a perturbed version of count 604 in the original histogram 600.
In another embodiment, the histogram engine 318 generates the requested histogram from the selected column of X, and perturbs the counts of each bin by decomposing the counts using a private wavelet decomposition algorithm. In such an embodiment, the histogram engine 318 aggregates the counts qi=1, 2, . . . , B for each bin bi=1, 2, . . . , B into a matrix (or vector) Q. The histogram engine 318 decomposes Q into a tree structure that is representative of a wavelet decomposition. Each leaf node of the tree corresponds to a count qi, and each parent node of the tree corresponds to one of multiple wavelet coefficients Cj=1, 2, . . . , m. The value of a wavelet coefficient cj is calculated based on the counts qt incorporated in the leaf nodes of the tree. This allows a count qi to be reconstructed as a function ƒi of the wavelet coefficients cj=1, 2, . . . , m. That is, for each count qi:

qi=ƒi(c0,C1, . . . ,Cm).

The histogram engine 318 generates a perturbed histogram by perturbing the wavelet coefficients, and reconstructing the counts using the perturbed wavelet coefficients. Specifically, the perturbed wavelet coefficients cDPi=1, 2, . . . , m are given by:



 
  
   c
   i
   DP
  
  =
  
   
    c
    i
   
   +
   
    
     L
     (
     
      
       2
       
        l
        +
        1
       
      
      
       ε
       ·
       m
      
     
     )
    
    .
   
  
 



The reconstructed counts from the perturbed wavelet coefficients is now given by:

qiDP=ƒi(C0DP,C1DP, . . . ,CmDP.

The histogram engine 318 outputs the perturbed histogram as the DP response 112 for display on the user interface 150.
In one embodiment, the histogram engine 318 may also be used to generate a differentially private visualization of data entries as described above in reference to the count engine 302 and FIG. 4. For example, the histogram module 318 may construct a multi-dimensional histogram corresponding to counts of the requested data entries in each region, perturb the histogram using mechanisms described above (e.g., private wavelet decomposition algorithm), and display the differentially private plot of the requested data entries on the user interface 150.
Model Testing Engine 320
The model testing engine 320 produces a DP response 112 responsive to the differentially private security system 102 receiving a query 108 for testing the performance of a classification model, given privacy parameters (ε,δ). The classification model is trained and tested on selected columns of X. As such, the model testing engine 320 may be appended to any other module that trains a classifier on X, such as the batch gradient engine 312, the stochastic gradient engine 314, or the random forest engine 316. For the example data object X shown in FIG. 2, the client 104 may submit a query 108 to return a DP response 112 for generating a performance evaluation of a support vector machine classifier that was trained using the stochastic gradient engine 314.
As discussed in reference to the random forest engine 316, classification models in general is trained on training data (Xtrain, ytrain) to learn the correlation between selected features of an entry and the category the entry belongs to. The training data (Xtrain, ytrain) may be extracted from a subset of entries contained in the data object X. Upon being trained, the classifier is able to receive a new data entry containing values for the selected features and generate an estimate of the category for the new entry.
Often times, the estimate of the category for an entry is determined by applying a cutoff threshold to a numerical, not categorical, output of a classifier. For example, in the random forest classifier described in reference to the random forest engine 316, the category associated with a leaf node tL is determined by the proportion of training data entries associated with each category, which is a numerical value. The random forest engine 316 may determine that a leaf node is associated with category “0” if the proportion of entries associated with label “0” is above a cutoff threshold of 0.5, 0.6, or 0.7. As another example, logistic regression classifiers output a numerical value in the range of [0, 1] given an entry of feature values. The entry may be classified into category “0” if the associated output is below a cutoff threshold of 0.5, 0.4, or 0.3. Regardless of the example, the cutoff threshold for determining the boundary between each category is a critical parameter depending on the context the classifier is applied to.
The model testing engine 320 receives a trained classifier and tests the performance of the trained classifier a series of cutoff thresholds, and generates a confusion matrix for each threshold indicating the performance of the classifier. The model testing engine 320 may test the performance of the classifier on testing data (Xtest, ytest). Similarly to training data, Xtest contains a set of entries with selected feature values, and ytest contains a vector of already known labels for each corresponding entry in Xtest. However, in contrast to training data, testing data (Xtest, ytest) comprises entries that are not present in training data (Xtrain, ytrain). That is, testing data comprises entries that the classifier has not “seen” yet.
The model testing engine 320 generates a series of cutoff thresholds based on the numerical values of p. FIG. 7A illustrates a recursive process for identifying threshold points of the classification output vector for the model testing engine 320, according to one embodiment. As shown in FIG. 7A, sample values of p are plotted on a range of 0 to 1. A series of k cutoff thresholds, or a series of intervals, are recursively identified by the median engine 304 such that the number of elements of p in each interval is approximately equal. Specifically, the median engine 304 recursively identifies the perturbed median for an interval and subsequently, its corresponding sub-intervals generated by dividing the interval by the identified perturbed median, until k thresholds are identified.
For example, as shown in FIG. 7, the median engine 304 identifies m1 as the perturbed median for the first interval [0, 1]. Subsequently, the median engine 304 identifies m2 as the perturbed median for the resulting sub-interval [0, m1], and m5 as the perturbed median for the sub-interval [m1, 1]. This process is repeated for sub-intervals [0, m2], [m2, m1], [m1, m5], [m5, 1] and for its sub-intervals until k thresholds, mi=1, 2, . . . , k are identified.
For each threshold mi, the model testing engine 320 generates corresponding category label estimates from p, and compares the estimates to the vector of known labels ytest. Given the comparisons, the model testing engine 320, constructs a confusion matrix that evaluates the performance of the classifier.
FIG. 7B illustrates an example confusion matrix 700 generated by the model testing engine 320, according to one embodiment. As shown in FIG. 7B, the confusion matrix 700 contains the number of testing data entries for 4 categories: i) entries that have an actual category of “1” and an estimate category of “1” (“True Positive” entries), ii) entries that have an actual category of “0” and an estimate category of “0” (“True Negative” entries), iii) entries that have an actual category of “0” and an estimate category of “1” (“False Positive” entries), and iv) entries that have an actual category of “1” and an estimate category of “0” (“False Negative” entries). For a given threshold, a data entry only contributes to one among the 4 categories.
For each threshold mi, the model testing engine 320 generates a perturbed confusion matrix by using the histogram engine 318. This is because each entry contributes to only one among the 4 disjoint categories, and thus, the entries in the confusion matrix 700 can be viewed as a histogram. The model testing engine 320 outputs each threshold mi, and the corresponding perturbed confusion matrix as the DP response 112 for display on the user interface 150.
Synthetic Database Engine 322
The synthetic database engine 322 produces a DP response 112 responsive to the differentially private security system 102 receiving a query 108 for transforming X into a synthetic database S, given privacy parameters (ε,δ). The resulting synthetic database S has a number of entries corresponding to that in X, but a fewer number of columns or features than X. Moreover, the spatial relationship between a pair of entries in X is retained in S. The transformation of X to S is (ε,δ)-differentially private with respect to a neighboring data object X′ with a 1-element difference from X.
The synthetic database engine 322 produces a DP response 112 of a differentially private synthetic database query by projecting the elements of X to S using a projection matrix. Assuming that data object X is a n×p matrix having n rows and p columns, the transformation by the synthetic database engine 322 is given by:



 
  S
  =
  
   X
   *
   
    J
    (
    
     G
     (
     
      
       
        4
        ·
        log
       
       ⁢
       δ
      
      
       
        n
        2
       
       ·
       
        ε
        2
       
      
     
     )
    
    )
   
  
 



where J is a p×k projection matrix, with k<p. The resulting synthetic database matrix S is a n×k matrix containing equal number of entries or rows as data object matrix X, but containing a smaller number of features or columns than the original data object X.

As discussed above, the transformation using projection matrix J is (ε,δ)-differentially private. Moreover, the spatial relationship between a pair of entries in X is retained in S. That is, the distance between a pair of entries (xi,xj) in the p-dimensional feature space of X is approximately equal to the distance between a pair of entries (si,sj) in the k-dimensional feature space of S. The synthetic database engine 322 outputs S as the DP response 112 for display on the user interface 150.
FIG. 8 is a modification 800 of the system 100 in FIG. 1 that allows the client 104 access to synthetic database 806 generated by the synthetic database engine 322, according to one embodiment. As shown in FIG. 8, the modified system 800 may allow the client 104 to access the synthetic database 806 generated by the synthetic database engine 322 through the differentially private security system 102. Since the transformation from X to S is (ε,δ)-differentially private, the privacy of X is retained.
FIG. 9 illustrates applying a clustering query to entries of the synthetic database, according to one embodiment. The various modules of the privacy system 160 and other analytical modules may be applied to the synthetic database 806 instead of the data object X to extract patterns or trends related to the information in X. The results of the analysis on the synthetic database 806 may also be provided for display on the user interface 150. For example, the client 104 may submit a clustering query on the entries of the synthetic database 806 using the batch gradient engine 312 or the stochastic gradient engine 314.
In the example of FIG. 9, clusters 902 and 904 are results of a non-differentially private clustering query on a data object X having two columns or features. Each cluster 902 and 904 are associated with a group of entries. Since the clustering query is not differentially private on X, the results of the query are not shown to the client 104. Clusters 906 and 908 are results of a non-differentially private clustering query on the synthetic database S having one column or feature due to the transformation by the synthetic database engine 322. Since the transformation preserves the spatial relationship between a pair of entries, cluster 906 is largely associated with the same entries in cluster 902, and cluster 908 is largely associated with the same entries in cluster 904. Since the synthetic database S is (ε,δ)-differentially private, the results of the clustering query may be displayed to the client 104 using the user interface 150.
Validation Engine 324
The validation engine 324 produces a DP response 112 responsive to the differentially private security system 102 receiving a request for whether a query 108 satisfies the definition of (ε,δ)-differential privacy for privacy parameters (ε,δ). In one embodiment, the validation engine 324 may receive a function call from the client 104 that points to the query 108. The query 108 may be, for example, an analytical model or an algorithm that can be applied to a data object X.
The validation engine 324 certifies whether the received query 108 satisfies the definition of (ε,δ)-differential privacy by applying the query 108 to example pairs of neighboring data objects (Z, Z′). Specifically, the validation engine 324 generates pairs of neighboring data objects (Z, Z′), having at most 1 entry different from each other. The validation engine 324 applies the received query 108 to each example pair of neighboring data objects (Z, Z′) and determines whether an estimate of the quantity Pr[M(X)∈S]/Pr[M(X′)∈S] satisfies the definition of (ε,δ)-differential privacy a sampling of outputs from S of the query M and over the randomness of the query M.
In one embodiment, the validation engine 324 may output a binary value to the client 104 as the DP response 112 for display on the user interface 150 that indicates whether or not the query 108 is (ε,δ)-differentially private. In some embodiments, the validation engine 324, in response to a determination that the query 108 is not (ε,δ)-differentially private, can reject or deny the query.
FIG. 10 illustrates a process for responding to a database query by executing a differentially private version of the query on the database, according to one embodiment.
A request from a client device to perform a query is received 1010 and a level of differential privacy corresponding to the request is identified. A set of data in the private database system and a set of operations to be performed based on the received request is identified 1012. The set of identified data in the private database system is accessed 1014. The set of operations is modified 1016 based on the received level of differential privacy. The set of modified operations is performed 1018 on the set of data to produce a differentially private result set. The differentially private result set is provided 1020 to the client device for display on the client device.
FIG. 11 is a block diagram illustrating components of an example machine able to read instructions from a machine-readable medium and execute them in a processor (or controller). Specifically, FIG. 11 shows a diagrammatic representation of a machine in the example form of a computer system 1100. The computer system 1100 can be used to execute instructions 1124 (e.g., program code or software) for causing the machine to perform any one or more of the methodologies (or processes) described herein. In alternative embodiments, the machine operates as a standalone device or a connected (e.g., networked) device that connects to other machines. In a networked deployment, the machine may operate in the capacity of a server machine or a client machine in a server-client network environment, or as a peer machine in a peer-to-peer (or distributed) network environment.
The machine may be a server computer, a client computer, a personal computer (PC), a tablet PC, a set-top box (STB), a smartphone, an internet of things (IoT) appliance, a network router, switch or bridge, or any machine capable of executing instructions 1124 (sequential or otherwise) that specify actions to be taken by that machine. Further, while only a single machine is illustrated, the term “machine” shall also be taken to include any collection of machines that individually or jointly execute instructions 1124 to perform any one or more of the methodologies discussed herein.
The example computer system 1100 includes one or more processing units (generally processor 1102). The processor 1102 is, for example, a central processing unit (CPU), a graphics processing unit (GPU), a digital signal processor (DSP), a controller, a state machine, one or more application specific integrated circuits (ASICs), one or more radio-frequency integrated circuits (RFICs), or any combination of these. The computer system 1100 also includes a main memory 1104. The computer system may include a storage unit 1116. The processor 1102, memory 1104 and the storage unit 1116 communicate via a bus 1108.
In addition, the computer system 1106 can include a static memory 1106, a display driver 1110 (e.g., to drive a plasma display panel (PDP), a liquid crystal display (LCD), or a projector). The computer system 1100 may also include alphanumeric input device 1112 (e.g., a keyboard), a cursor control device 1114 (e.g., a mouse, a trackball, a joystick, a motion sensor, or other pointing instrument), a signal generation device 1118 (e.g., a speaker), and a network interface device 1120, which also are configured to communicate via the bus 1108.
The storage unit 1116 includes a machine-readable medium 1122 on which is stored instructions 1124 (e.g., software) embodying any one or more of the methodologies or functions described herein. The instructions 1124 may also reside, completely or at least partially, within the main memory 1104 or within the processor 1102 (e.g., within a processor's cache memory) during execution thereof by the computer system 1100, the main memory 1104 and the processor 1102 also constituting machine-readable media. The instructions 1124 may be transmitted or received over a network 1126 via the network interface device 1120.
While machine-readable medium 1122 is shown in an example embodiment to be a single medium, the term “machine-readable medium” should be taken to include a single medium or multiple media (e.g., a centralized or distributed database, or associated caches and servers) able to store the instructions 1124. The term “machine-readable medium” shall also be taken to include any medium that is capable of storing instructions 1124 for execution by the machine and that cause the machine to perform any one or more of the methodologies disclosed herein. The term “machine-readable medium” includes, but not be limited to, data repositories in the form of solid-state memories, optical media, and magnetic media.
Other Considerations
Some portions of the above description describe the embodiments in terms of algorithmic processes or operations. These algorithmic descriptions and representations are commonly used by those skilled in the data processing arts to convey the substance of their work effectively to others skilled in the art. These operations, while described functionally, computationally, or logically, are understood to be implemented by computer programs comprising instructions for execution by a processor or equivalent electrical circuits, microcode, or the like. Furthermore, it has also proven convenient at times, to refer to these arrangements of functional operations as modules, without loss of generality. The described operations and their associated modules may be embodied in software, firmware, hardware, or any combinations thereof.
As used herein any reference to “one embodiment” or “an embodiment” means that a particular element, feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment. The appearances of the phrase “in one embodiment” in various places in the specification are not necessarily all referring to the same embodiment.
Some embodiments may be described using the expression “coupled” and “connected” along with their derivatives. It should be understood that these terms are not intended as synonyms for each other. For example, some embodiments may be described using the term “connected” to indicate that two or more elements are in direct physical or electrical contact with each other. In another example, some embodiments may be described using the term “coupled” to indicate that two or more elements are in direct physical or electrical contact. The term “coupled,” however, may also mean that two or more elements are not in direct contact with each other, but yet still co-operate or interact with each other. The embodiments are not limited in this context.
As used herein, the terms “comprises,” “comprising,” “includes,” “including,” “has,” “having” or any other variation thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, article, or apparatus that comprises a list of elements is not necessarily limited to only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Further, unless expressly stated to the contrary, “or” refers to an inclusive or and not to an exclusive or. For example, a condition A or B is satisfied by any one of the following: A is true (or present) and B is false (or not present), A is false (or not present) and B is true (or present), and both A and B are true (or present).
In addition, use of the “a” or “an” are employed to describe elements and components of the embodiments herein. This is done merely for convenience and to give a general sense of the disclosure. This description should be read to include one or at least one and the singular also includes the plural unless it is obvious that it is meant otherwise.
Upon reading this disclosure, those of skill in the art will appreciate still additional alternative structural and functional designs for a system and a process for receiving a query for a private database, and responding to the query by executing a differentially private version of the query on the private database. Thus, while particular embodiments and applications have been illustrated and described, it is to be understood that the described subject matter is not limited to the precise construction and components disclosed herein and that various modifications, changes and variations which will be apparent to those skilled in the art may be made in the arrangement, operation and details of the method and apparatus disclosed herein.

## Claims

1. A database system configured to implement differential privacy, comprising:
a processor configured to execute computer program instructions; and
a non-transitory computer-readable storage medium storing computer program instructions executable by the processor to perform actions comprising:
receiving a database query requesting a differentially private response to the query;
determining a privacy parameter ε associated with the query, wherein ε describes a degree of information to release about a set of data stored by the database system that is responsive to the query;
identifying a privacy budget associated with the query, the privacy budget specified in terms of ε and representing a degree of information available to be released about data by the database system;
decrementing the privacy budget by a first ε spend determined responsive to the query; and
applying the query to the database system by performing a differentially private set of operations on the set of data stored by the database system that is responsive to the query to produce a differentially private result set that releases the degree of information about the set of data described by the privacy parameter ε, wherein the set of operations comprises a count operation on a subset of the set of data, and wherein the set of operations comprises perturbing results of the count operation on the subset by a factor defined by a Gaussian random variable G( ) to produce the differentially private result set.

2. The database system of claim 1, wherein the set of operations comprises 1) identifying a loss function for the set of data, the loss function comprising a function of a set of parameters θ describing correlations in the set of data, and 2) minimizing a perturbed loss function over the set of parameters θ, and wherein the set of operations comprises perturbing the loss function by a factor defined by:



 
  
   
    θ
    T
   
   ⁢
   
    v
    (
    
     G
     (
     
      
       4
       ·
       
        K
        2
       
       ·
       
        R
        2
        2
       
       ·
       
        (
        
         
          log
          ⁢
          
           1
           δ
          
         
         +
         ε
        
        )
       
      
      
       ε
       2
      
     
     )
    
    )
   
  
  ,
 



to produce the set of parameters θ that minimize the perturbed loss function private set, where K, R2 are constants, wherein δ is another privacy parameter and ν( ) denotes a vector populated with random variables.

3. The database system of claim 1, wherein the set of operations comprises a median operation on elements associated with a second subset of the set of data, and wherein the set of operations comprises perturbing results of the median operation on the elements by a factor defined by:



 
  
   
    c
    1
   
   ·
   
    
     S
     M
    
    (
    
     U
     ;
        
     
      
       c
       2
      
      ·
      
       ε
       2
      
     
    
    )
   
   ·
   
    
     C
     ⁡
     (
     1
     )
    
    
     ε
     2
    
   
  
  ,
 



to produce the differentially private result set, where c1, c2 are constants, C(1) is a Cauchy random variable, ε2 is a second ε spend, and SM(U; c2·ε) is a function that computes a smooth sensitivity of the elements with parameter c2·ε2.

4. The database system of claim 1, wherein the set of operations comprises an interquartile range operation on a second subset of the set of data, and wherein the set of operations comprises perturbing a result of the interquartile range operation by a factor defined by:



 
  
   
    (
    
     1
     
      1
      +
      
       log
       ⁢
       n
      
     
    
    )
   
   
    L
    (
    
     
      c
      1
     
     
      ε
      2
     
    
    )
   
  
  ,
 



to produce the differentially private result set, where L( ) denotes a Laplacian random variable, c1 and n are constants, and ε2 is a second ε spend.

5. The database system of claim 1, wherein the set of operations comprises 1) identifying a loss function for a second subset of the set of data, the loss function a function of a set of parameters θ describing correlations in the second subset, and 2) for each time step, generating a perturbed estimate for the set of parameters θ that minimize the loss function, and wherein the set of operations comprises for each time step, perturbing an estimate for the set of parameters θ by:



 
  
   
    η
    t
   
   ·
   
    v
    (
    
     G
     (
     
      
       
        c
        1
        2
       
       ·
       
        n
        2
       
       ·
       
        (
        
         log
         ⁢
         
          
           n
           δ
          
          ·
          log
         
         ⁢
         
          1
          δ
         
        
        )
       
      
      
       
        ε
        4
       
       2
      
     
     )
    
    )
   
  
  ,
 



to produce the perturbed estimate for the set of parameters θ as the differentially private result set, where ηt, n, and c1 are constants, and ε2 is a second ε spend, wherein ν( ) denotes a vector populated with random variables and δ is another privacy parameter.

6. The database system of claim 1, wherein each entry in a second subset of the set of data is labeled with a category chosen from a set of two or more categories, and wherein the set of operations comprises:
receiving a trained classifier and generating an output vector by applying the classifier to entries of the second subset, each element of the output vector corresponding to a numerical output of the classifier for a corresponding entry in the second subset;
identifying a threshold value and assigning categories for each of the elements of the output vector based on a perturbed threshold value; and
recording counts related to a performance of the classifier, the counts generated by comparing the assigned categories of the elements of the output vector to the corresponding label in the second subset, and

the set of operations comprises:
perturbing the threshold value based on a second ε spend ε2 to generate the perturbed threshold value; and
perturbing the counts relating to the performance of the classifier based once to produce the perturbed counts as the differentially private result set.

7. The database system of claim 1, wherein the differentially private set of operations is performed responsive to determining that the first ε spend can be accommodated by the privacy budget.

8. The database system of claim 1, wherein the privacy budget is associated with a user of the database system.

9. A method for implementing differential privacy, comprising:
receiving a database query requesting a differentially private response to the query;
determining a privacy parameter ε associated with the query, wherein describes a degree of information to release about a set of data stored by a database system that is responsive to the query;
identifying a privacy budget associated with the query, the privacy budget specified in terms of ε and representing a degree of information available to be released about data by the database system;
decrementing the privacy budget by a first ε spend determined responsive to the query; and
applying the query to the database system by performing a differentially private set of operations on the set of data stored by the database system that is responsive to the query to produce a differentially private result set that releases the degree of information about the set of data described by the privacy parameter ε, wherein the set of operations comprises a count operation on a subset of the set of data, and wherein the set of operations comprises perturbing results of the count operation on the subset by a factor defined by a Gaussian random variable G( ) to produce the differentially private result set.

10. The method of claim 9, wherein the set of operations comprises 1) identifying a loss function for the set of data, the loss function comprising a function of a set of parameters θ describing correlations in the set of data, and 2) minimizing a perturbed loss function over the set of parameters θ, and wherein the set of operations comprises perturbing the loss function by a factor defined by:



 
  
   
    θ
    T
   
   ⁢
   
    v
    (
    
     G
     (
     
      
       4
       ·
       
        K
        2
       
       ·
       
        R
        2
        2
       
       ·
       
        (
        
         
          log
          ⁢
          
           1
           δ
          
         
         +
         ε
        
        )
       
      
      
       ε
       2
      
     
     )
    
    )
   
  
  ,
 



to produce the set of parameters θ that minimize the perturbed loss function as the differentially private set, where K, R2 are constants, wherein δ is another privacy parameter and ν( ) denotes a vector populated with random variables.

11. The method of claim 9, wherein the set of operations comprises a median operation on elements associated with a second subset the set of data, and wherein the set of operations comprises perturbing results of the median operation on the elements by a factor defined by:



 
  
   
    c
    1
   
   ·
   
    
     S
     M
    
    (
    
     U
     ;
        
     
      
       c
       2
      
      ·
      
       ε
       2
      
     
    
    )
   
   ·
   
    
     C
     ⁡
     (
     1
     )
    
    
     ε
     2
    
   
  
  ,
 



to produce the differentially private result set, where c1, c2 are constants, C(1) is a Cauchy random variable, ε2 is a second ε spend, and SM(U; c2·ε) is a function that computes a smooth sensitivity of the elements with parameter c2·ε2.

12. The method of claim 9, wherein the set of operations comprises an interquartile range operation on a second subset of the set of data, and wherein the set of operations comprises perturbing a result of the interquartile range operation by a factor defined by:



 
  
   
    (
    
     1
     
      1
      +
      
       log
       ⁢
       n
      
     
    
    )
   
   
    L
    (
    
     
      c
      1
     
     
      ε
      2
     
    
    )
   
  
  ,
 



to produce the differentially private result set, where L( ) denotes a Laplacian random variable, c1 and n are constants and ε2 is a second ε spend.

13. The method of claim 9, wherein the differentially private set of operations is performed responsive to determining that the first ε spend can be accommodated by the privacy budget.

14. The method of claim 9, wherein the privacy budget is associated with a user of the database system.

15. A non-transitory computer-readable storage medium storing computer program instructions executable by a processor to perform actions for implementing differential privacy comprising:
receiving a database query requesting a differentially private response to the query;
determining a privacy parameter ε associated with the query, wherein describes a degree of information to release about a set of data stored by a database system that is responsive to the query;
identifying a privacy budget associated with the query, the privacy budget specified in terms of ε and representing a degree of information available to be released about data by the database system;
decrementing the privacy budget by a first ε spend determined responsive to the query; and
applying the query to the database system by performing a differentially private set of operations on the set of data stored by the database system that is responsive to the query to produce a differentially private result set that releases the degree of information about the set of data described by the privacy parameter ε, wherein the set of operations comprises a count operation on a subset of the set of data, and wherein the set of operations comprises perturbing results of the count operation on the subset by a factor defined by a Gaussian random variable G( ) to produce the differentially private result set.

16. The non-transitory computer-readable storage medium of claim 15, wherein the set of operations comprises 1) identifying a loss function for the set of data, the loss function comprising a function of a set of parameters θ describing correlations in the set of data, and 2) minimizing a perturbed loss function over the set of parameters θ, and wherein the set of operations comprises perturbing the loss function by a factor defined by:



 
  
   
    θ
    T
   
   ⁢
   
    v
    (
    
     G
     (
     
      
       4
       ·
       
        K
        2
       
       ·
       
        R
        2
        2
       
       ·
       
        (
        
         
          log
          ⁢
          
           1
           δ
          
         
         +
         ε
        
        )
       
      
      
       ε
       2
      
     
     )
    
    )
   
  
  ,
 



to produce the set of parameters θ that minimize the perturbed loss function, where K, R2 are constants, wherein δ is another privacy parameter and ν( ) denotes a vector populated with random variables.

17. The non-transitory computer-readable storage medium of claim 15, wherein the set of operations comprises a median operation on elements associated with a second subset of the set of data, and wherein the set of operations comprises perturbing results of the median operation on the elements by a factor defined by:



 
  
   
    c
    1
   
   ·
   
    
     S
     M
    
    (
    
     U
     ;
        
     
      
       c
       2
      
      ·
      
       ε
       2
      
     
    
    )
   
   ·
   
    
     C
     ⁡
     (
     1
     )
    
    
     ε
     2
    
   
  
  ,
 



to produce the differentially private result set, where c1, c2 are constants, C(1) is a Cauchy random variable, ε2 is a second ε spend, and SM(U; c2·ε) is a function that computes a smooth sensitivity of the elements with parameter c2·ε2.

18. The non-transitory computer-readable storage medium of claim 15, wherein the set of operations comprises an interquartile range operation on a second subset of the set of data, and wherein the set of operations comprises perturbing a result of the interquartile range operation by a factor defined by:



 
  
   
    (
    
     1
     
      1
      +
      
       log
       ⁢
       n
      
     
    
    )
   
   
    L
    (
    
     
      c
      1
     
     
      ε
      2
     
    
    )
   
  
  ,
 



to produce the differentially private result set, where L( ) denotes a Laplacian random variable, c1 and n are constants and ε2 is a second ε spend.

19. The non-transitory computer-readable storage medium of claim 15, wherein the differentially private set of operations is performed responsive to determining that the first ε spend can be accommodated by the privacy budget.

20. The non-transitory computer-readable storage medium of claim 15, wherein the privacy budget is associated with a user of the database system.