Patent Publication Number: US-11664098-B2

Title: Determining journalist risk of a dataset using population equivalence class distribution estimation

Description:
CROSS REFERENCE 
     This application is a continuation of U.S. Ser. No. 14/953,195 filed Nov. 27, 2015, entitled, “Determining Journalist Risk of a Dataset Using Population Equivalence Class Distribution Estimation”, which claims priority from U.S. Provisional Application No. 62/085,307 filed Nov. 27, 2014; both of which are incorporated herein in their entirety. 
    
    
     TECHNICAL FIELD 
     The present disclosure relates to databases and particularly to protecting privacy by estimating risk of identification of personal data stored in the databases. 
     BACKGROUND 
     Personal information is being continuously captured in a multitude of electronic databases. Details about health, financial status and buying habits are stored in databases managed by public and private sector organizations. These databases contain information about millions of people, which can provide valuable research, epidemiologic and business insight. For example, examining a drugstore chain&#39;s prescriptions can indicate where a flu outbreak is occurring. To extract or maximize the value contained in these databases, data custodians must often provide outside organizations access to their data. In order to protect the privacy of the people whose data is being analyzed, a data custodian will “de-identify” or “anonymize” information before releasing it to a third-party. An important type of de-identification ensures that data cannot be traced to the person about whom it pertains, this protects against ‘identity disclosure’. 
     When de-identifying records, many people assume that removing names and addresses (direct identifiers) is sufficient to protect the privacy of the persons whose data is being released. The problem of de-identification involves those personal details that are not obviously identifying. These personal details, known as quasi-identifiers, include the person&#39;s age, sex, postal code, profession, ethnic origin and income, financial transactions, medical procedures (to name a few). To be able to de-identify data the assessment of the risk of re-identification is required to be determined. Therefore there is a need for improved risk assessment of data sets. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       Further features and advantages of the present disclosure will become apparent from the following detailed description, taken in combination with the appended drawings, in which: 
         FIG.  1    shows a method for determining journalist risk for a cross-sectional dataset using population equivalence class distribution estimation; 
         FIG.  2    shows a method for determining journalist risk for a longitudinal dataset using population equivalence class distribution estimation; and 
         FIG.  3    shows system for determining journalist risk of a dataset using population equivalence class distribution estimation. 
     
    
    
     It will be noted that throughout the appended drawings, like features are identified by like reference numerals. 
     DETAILED DESCRIPTION 
     In accordance with an aspect of the present disclosure there is provided a computer implemented method of determining journalist risk associated with a dataset, the method comprising: retrieving the dataset containing a plurality of records containing personal data, the dataset representing a sample of individuals and data associated with a larger population; determining sample equivalence class (EC) distribution of the dataset; equating a population EC distribution to the determined sample EC distribution; calculating probability that an EC in the dataset of size x came from population of size y for all x and y; and calculating the journalist risk measurement using calculated probability; wherein the equivalence classes define a collection of all records in the dataset containing identical values for all quasi-identifiers in the data. 
     In accordance with another aspect of the present disclosure there is provided a system for determining journalist risk associated with a dataset, the system comprising: a memory containing the dataset containing a plurality of records containing personal data, the dataset representing a sample of individuals and data associated with a larger population, dataset having equivalence classes defining a collection of all records in the dataset containing identical values for all quasi-identifiers in the data; and a processor coupled to the memory, the processor executing instructions for: determining sample equivalence class (EC) distribution of the dataset; equating a population EC distribution to the determined sample EC distribution; calculating probability that an EC in the dataset of size x came from population of size y for all x and y; and calculating the journalist risk measurement using calculated probability. 
     In accordance with yet another aspect of the present disclosure there is provided a non-transitory computer readable memory containing instructions for determining journalist risk associated with a dataset, the instructions which when executed by a processor perform: retrieve the dataset containing a plurality of records containing personal data, the dataset representing a sample of individuals and data associated with a larger population; determine sample equivalence class (EC) distribution of the dataset; equate a population EC distribution to the determined sample EC distribution; calculate probability that an EC in the dataset of size x came from population of size y for all x and y; and calculate the journalist risk measurement using calculated probability; wherein the equivalence classes define a collection of all records in the dataset containing identical values for all quasi-identifiers in the data. 
     Embodiments are described below, by way of example only, with reference to  FIGS.  1 - 3   . 
     The methodology presented can augment the risk measurement processes within risk measurement and de-identification software. It handles the case when the dataset on which risk is being measured, is a sample of a larger dataset; these will be referred to as the sample and the population, respectively. It is further assumed that the dataset is a random sample of the population or that it is sampled on some unknowable attribute. For example, measuring journalist risk would not be appropriate in the case of a sample dataset of all births of twins from a birth database. A dataset is longitudinal if it tracks the same type of information on the same subjects at multiple points in time. For example, part of a longitudinal dataset could contain specific patients and their medical results over a period of years. Each patient may have varying times and number of visits. A cross-sectional dataset comprises data collected on many subjects at the same point of time, or without regard to differences in time. Analysis of cross-sectional datasets usually consists of comparing the differences among the subjects. 
     In the description the following terms are used:
         Equivalence Class: collection of all individuals in a dataset containing identical values for all quasi-identifiers.   Similarity Measure: used in the context of longitudinal risk measurement, it is the number of individuals who look like a given individual in a dataset.   Journalist Risk: Risk measured when the dataset is a sample of a larger dataset.   Prosecutor Risk: Risk measured only in reference to the dataset itself.
 
Average Risk Measurement:
       

     The method by which average risk can be measured for both cross-sectional and longitudinal data can be equivalent to calculating the risk of re-identification of every patient in the database and then taking the average. 
     The re-identification risk of a record is 1/F where F is either the similarity measure of the record (for longitudinal data) or the size of the equivalence class the record belongs to (for cross-sectional data). Because it is possible to construct an alias equivalence class distribution from the similarity measure distribution, risk can measured once the appropriate equivalence class distribution is known. 
     When prosecutor risk is being measured and the dataset does not represent a sample, the equivalence class distribution of the dataset is used to measure the risk However, when the dataset represents a sample, F would be the size of the equivalence class to which the patient belongs in the population. Generally, it is not possible to know the equivalence class distribution of the population or the population equivalence class size from which the record was sampled. 
     Maximum Risk Measurement 
     Maximum journalist risk can be found by inverting the minimum similarity measure in the population for longitudinal data or inverting the minimum equivalence class size in the population for cross-sectional data. 
     Cross-Sectional Dataset 
     To begin, the estimation of the population equivalence class distribution and the subsequent risk measurement will be considered within the context of a cross-sectional dataset. Referring to  FIG.  1   , journalist risk for cross-sectional dataset using population equivalence class distribution estimation can be assessed by retrieving the data set from a store ( 102 ). The dataset may be accessed in whole or in segments either stored local or remotely by a computing device. 
     Estimating Population Equivalence Class Distribution 
     The sample equivalence class (EC) distribution is determined from dataset ( 104 ). In order measure journalist risk for cross-sectional data, the population equivalence class distribution must be estimated. Following a method developed by Zayatz, L. (Estimation of the percent of unique population elements on a microdata file using the sample. Washington: US Bureau of the Census; 1991.) the initial assumption that the population equivalence class distribution can be equated to the sample equivalence class distribution is made ( 106 ). 
     From there a combinatorics calculation gives the probability that an equivalence class of size x in the sample came from an equivalence class of size y in the population. First, the probability that an equivalence class of size y in the population will be sampled to an equivalence class of size x in the sample, is calculated: 
               Pr   ⁢           ⁢     (       size   ⁢           ⁢   x   ⁢           ⁢   in   ⁢           ⁢   sample     ❘     size   ⁢           ⁢   y   ⁢           ⁢   in   ⁢           ⁢   population       )       =         
where n is the sample size, N is the population size and
 
     
       
         
           
             = 
             
           
           
         
       
     
     The probability that an equivalence class is of size y in the population and is an equivalence class of size x in the population will be: 
             Pr   ⁢           ⁢     (       size   ⁢           ⁢   x   ⁢           ⁢   in   ⁢           ⁢   sample     ❘     size   ⁢           ⁢   y   ⁢           ⁢   in   ⁢           ⁢   population       )     ×     Pr   ⁡     (   y   )             
where Pr(y)=num. of equivalence classes of size y in pop./total number of equivalence classes in population. Because of the assumption that the equivalence class distributions of the sample and population are the same, it may also be written that:
 
     
       
         
           
             
               P 
               ⁢ 
               
                 r 
                 ⁡ 
                 
                   ( 
                   y 
                   ) 
                 
               
             
             = 
             
               
                 num 
                 . 
                 
                     
                 
                 ⁢ 
                 of 
               
               ⁢ 
               
                   
               
               ⁢ 
               equivalence 
               ⁢ 
               
                   
               
               ⁢ 
               classes 
               ⁢ 
               
                   
               
               ⁢ 
               of 
               ⁢ 
               
                   
               
               ⁢ 
               size 
               ⁢ 
               
                   
               
               ⁢ 
               y 
               ⁢ 
               
                   
               
               ⁢ 
               in 
               ⁢ 
               
                   
               
               ⁢ 
               
                 sample 
                 / 
                 total 
               
               ⁢ 
               
                   
               
               ⁢ 
               number 
               ⁢ 
               
                 
                     
                 
                 ⁢ 
                 
                     
                 
               
               ⁢ 
               of 
               ⁢ 
               
                   
               
               ⁢ 
               equivalence 
               ⁢ 
               
                   
               
               ⁢ 
               classes 
               ⁢ 
               
                   
               
               ⁢ 
               in 
               ⁢ 
               
                   
               
               ⁢ 
               sample 
             
           
         
       
     
     Finally Bayes&#39; theorem can be applied in conjunction with the law of total probability to calculate the probability that an equivalence class of size x in the sample came from an equivalence class of size y in the population: 
               Pr   ⁢           ⁢     (       size   ⁢           ⁢   y   ⁢           ⁢   in   ⁢           ⁢   population     ❘     size   ⁢           ⁢   x   ⁢           ⁢   in   ⁢           ⁢   sample       )       =         Pr   ⁡     (   y   )       ×     Pr   (       size   ⁢           ⁢   x   ⁢           ⁢   in   ⁢           ⁢   sample     ❘     size   ⁢           ⁢   y   ⁢           ⁢   in   ⁢           ⁢   population       )         ⁢     Pr   ⁡     (   C   )       ×     Pr   (       size   ⁢           ⁢   x   ⁢           ⁢   in   ⁢           ⁢   sample     ❘     size   ⁢           ⁢   C   ⁢           ⁢   in   ⁢           ⁢   population       )               
where max size is the maximum equivalence class size in the sample and the population. The probability sample EC of size x came from population EC of size y for all x and y is then calculated ( 108 ). These probabilities will be used to calculate journalist risk ( 112 ) as described below.
 
Calculating Journalist Risk
 
     Without any consideration of a population, the average prosecutor risk may be calculated by dividing the proportion of records belong to equivalence classes of a given size, by the equivalence class size and summing the result: 
             Risk   =         
where Prop(i)=number of records in equivalence classes of size i/n
 
     This equation is modified to calculate journalist risk when the dataset is recognized as a sample: 
     
       
         
           
             Risk 
             = 
             
               ⁢ 
               
                 Pr 
                 ⁡ 
                 
                   ( 
                   
                     
                       size 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       • 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       in 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         pop 
                         . 
                       
                     
                     ❘ 
                     
                       size 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       • 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       in 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         samp 
                         . 
                       
                     
                   
                   ) 
                 
               
               × 
               
                 
                   •••• 
                   ⁡ 
                   
                     ( 
                     • 
                     ) 
                   
                 
                 / 
                 • 
               
             
           
         
       
     
     This journalist risk equation assumes the re-identification attack involves choosing someone from the sample and trying to re-identify them in the population. An alternate attack would see an adversary trying to match someone from the population in the sample. The risk in that case would be: 
             Risk   =     ⁢     Pr   ⁡     (       size   ⁢           ⁢   •   ⁢           ⁢   in   ⁢           ⁢     pop   .       ❘     size   ⁢           ⁢   •   ⁢           ⁢   in   ⁢           ⁢     samp   .         )       ×   ••   ⁢           ⁢   _   ⁢           ⁢     •••   ⁡     (   •   )       ×   •           
where J is the number of equivalence classes in the sample and EC_num(i) is the number of equivalence classes of size i in the sample.
 
     Overall, journalist risk will be the maximum of these two risk calculations. 
     Longitudinal Dataset 
     Calculating journalist risk for longitudinal data will incorporate the methodology involved in calculating journalist risk for cross-sectional data as described in connection with  FIG.  1   . With reference to  FIG.  2   , when dealing with a longitudinal dataset, instead of grouping records into equivalence classes, a similarity measure is calculated for each individual in the dataset. The similarity measure is a description of how many other individuals that individual looks like. The dataset may be retrieved or accessed in whole or in segments either stored local or remotely by a computing device ( 202 ). 
     Similarity measures for the dataset are determined ( 204 ). The reason that equivalence classes are not formed is because the measure for defining what it means for a selected individual to be similar to another individual is non-commutative: patient A might be similar to patient B, but patient B might not be similar to patient A. As an example, consider two individuals with the same values for all demographic quasi-identifiers. Patient A has disease X, while patient B has disease X and Y. An adversary trying to re-identify Patient A will look for patients with disease X and will find two such patients: patient B is similar to patient A. An adversary trying to re-identify Patient B will look for patients with disease X and Y and will only find one such patient: in this case patient A is not similar to patient B. 
     An equivalence class distribution is required for the Zayatz method and so the similar measure distribution must be transformed. A consequence of the non-commutative nature of the similarity measure is that the corresponding equivalence class distribution might contain non-integer values. 
     For example it might be determined that there are 3 individuals with similarity measures of 4, in which case the number of equivalence classes of size 4 would be ¾=0.75. 
     There are no difficulties in using non-integer equivalence counts with the population equivalence class estimation methodology as it only uses the corresponding probability distribution created by dividing the equivalence class size counts by the total number of equivalence classes. 
     Before the similarity measure distribution is converted to an equivalence class distribution, the similarity measure for each individual in the data set must be rounded to an integer value ( 206 ). The reason it may not be an integer value, is because when determining how many individuals a particular individual looks like, multiple combinations of the patients&#39; records will be considered individually yielding multiple integers which must then be averaged. 
     A rounding step is necessary for two reasons. Firstly, the similarity measure directly translates into an equivalence class size in the process described above. A combinatoric approach is used where integer sizes are implicitly assumed. Secondly, the method considers individuals who have the same number of individuals who are similar to them as equivalence classes. So if there are individuals with similarity counts of 2.90, 3, 3.1, it is desirable to model them as one equivalence class of size 3. This would not be the case unless rounding is performed. The number of each similarity measure count is divided by the measure to obtain sample EC distribution ( 208 ). 
     Once the similarity measure distribution of a longitudinal data set is converted to an equivalence class distribution ( 210 ), the same risk measurement steps are followed that have been previously detailed for a cross-sectional data set in reference to  FIG.  1   . 
     Measuring Longitudinal Risk by Taking Sample 
     There may be datasets that are so large that it becomes impractical to determine a similarity count for every individual in a longitudinal dataset. The similarity counts for a certain percentage of individuals in the sample dataset can be determined and subsequently the similarity counts for the entire dataset could be estimated by dividing those similarity counts by the percentage of individuals used to construct the counts. 
     For example if 10% of individuals in a dataset had their similarity measures determined and aggregated into counts by comparing against all individuals in the dataset, the similarity counts for the dataset would be determined by dividing those counts by 0.10. These counts would represent an estimate for the similarity measure counts in the sample dataset. The probability sample EC of size x came from population EC of size y for all x and y is then calculated ( 212 ). These numbers would then be used to begin the process for measuring journalist risk ( 214 ) as discussed above. The calculated risk for the dataset can then be presented ( 216 ), and as described above. 
       FIG.  3    provides a system for journalist risk assessment using population equivalence class distribution estimation as used in connection with the above described method. A computer or server  310  providing at least a processor  312 , memory  314  and input/output interface  316 , implements the code for executing the de-identification process. A source dataset  302  is stored on non-transitory computer readable storage memory which may reside locally or remotely from processing unit  312 . The dataset is processed by the computer  310  to provide risk assessment which can be used for the optimal de-identification. Generalization strategies and levels of suppression can also be provided through template files, user selection or input through interaction with the computer  310 , either directly through input devices such a keyboard/mouse and display or remotely through a connected computing network  326 . External storage  322 , or computer readable memory such as compact disc, digital versatile disc or other removable memory devices  324  may be used to provide the instructions for execution of the risk assessment and de-identification methods or provide input for generalization or suppression parameters via I/O unit  316 . Execution of the method on processor  312  retrieves  306  and provides an assessment of risk or provide the resulting parameters which can be utilized in performing de-identification of the dataset to meet a desired risk threshold. The de-identification process may use optimization such as optimal lattice anonymization for determine a level of de-identification which meets desired risk threshold. 
     Each element in the embodiments of the present disclosure may be implemented as hardware, software/program, or any combination thereof. Software codes, either in its entirety or a part thereof, may be stored in a non-transitory computer readable medium or memory (e.g., as a RAM, ROM, for example a non-volatile memory such as flash memory, CD ROM, DVD ROM, Blu-Ray™, a semiconductor ROM, USB, or a magnetic recording medium, for example a hard disk). The program may be in the form of source code, object code, a code intermediate source and object code such as partially compiled form, or in any other form. 
     It would be appreciated by one of ordinary skill in the art that the system and components shown in  FIGS.  1 - 3    may include components not shown in the drawings. For simplicity and clarity of the illustration, elements in the figures are not necessarily to scale, are only schematic and are non-limiting of the elements structures. It will be apparent to persons skilled in the art that a number of variations and modifications can be made without departing from the scope of the invention as defined in the claims.