when a crime is committed and a biological trace such as blood stain or semen is found at the crime scene or the body of the victim , forensic scientists can link the case to the arrested suspect through the matching of the trace evidence to the suspect , using dna profiling . for cases in which no suspect can be identified based on non - dna evidences such as fingerprints or witness reports , the police force may search a database of dna profiles from previously convicted criminals or unsolved crime cases and open an investigation on individuals with perfectly matched profiles . since 1995 when the first national dna database came to operation , dna database search has become an important crime - solving tool for suspect identification , benefit from the amassing of large offender dna databases in many countries . for example , as of june 2011 , the us federal dna database codis has assisted in over 141,300 investigations in the usa by more than 147,200 hits produced . the evaluation of the evidentiary value of perfect matches from dna database search has been thoroughly discussed by balding and donnelly , stockmarr , and meester and sjerps , among many others . if no offender profile in the database perfectly matches the crime trace , an additional search can be performed , hoping that an individual in the database is a close relative of the perpetrator and can be identified through the search . the use of familial search on dna database starts to become popular in recent years . as of may 2011 , about 40 serious crimes has been solved with the aid of about 200 familial searches in the uk , showing that the familial database search has potential to be an effective forensic tool that can increase the number of suspects identified based on dna evidences . since many and even all individuals in the database may be qualified as the relative of the perpetrator , a scoring scheme is needed to rank the offender profiles , so that the police force can focus their investigation upon the top - listed candidates . a widely adopted familial search score is the likelihood ratio ( lr ) of the following hypotheses : 
 hp :  a relative of the individual is a contributor to the crime trace.hd :  a relative of the individual is not a contributor to the crime trace.the lr for various two - person relationships can be easily calculated using the formulae provided by evett and weir . the performance of using lr scores for familial database searching has been evaluated in various articles [ 510 ] . based on simulated and real dna databases , the familial searches were demonstrated to be able to identify the first - order relative such as parent / child or full sibling of the perpetrator in the 100 top - listed candidates in over 70% of the cases , provided that the database really contains the profile of the relative . hp :  a relative of the individual is a contributor to the crime trace . hd :  a relative of the individual is not a contributor to the crime trace . in some typical crime cases such as murder cases or rape cases , the biological trace is often observed as a dna mixture contributed by the victim and the perpetrator . over the years , the evaluation of the dna mixtures in probable cause cases has been studied extensively [ 1114 ] , and various formulae have been developed to handle situations when relatives are involved [ 1517 ] . in two previous articles [ 18 , 19 ] , we have reported our progress on the attempt to apply database search on cases with dna mixtures as part of the evidences and derived several formulae for evaluating the evidentiary values of  cold hits .  to go a further step , chung et  al . has extended the familial database search method so as to handle mixture cases , and the performance is shown to be as good as in single source cases . based on the results of a familial search , the number of individuals to be further investigated may depend on the regional policy . for example , according to the california familial search policy , after the list of top - ranked candidates is produced , lineage dna testing will be conducted on up to 168 candidates by using y - str typing . in addition to familial search policy , the scale of the investigation will also vary based on practical considerations including the limitation of manpower , the resources of the police offices , and the scale of the crime severity . the main cost factor is in the need to investigate the background of each candidate so that irrelevant candidates can be eliminated as being the possible relative of the perpetrator . investigating on thousands of candidates for every case is impractical as the law enforcement resources will be overwhelmed by the heavy case load . it would be more effective if a large number of candidates are investigated only for serious and high - profile criminal cases . it is therefore necessary to determine the scale of the investigation so that the effectiveness of the search can be guaranteed under reasonable cost . ge et  al .   had given some suggestions on the thresholds to the lr as well as the identity - by - state scores , with the aim of balancing the false - positive and false - negative rates . these thresholds , however , are specific to a particular database as they are determined based on a simulation study using caucasian population data on the 13 codis str loci . the primary aim of this work is to establish a general strategy on deciding the number of top - listed candidates to be investigated after the familial search , according to the statistical criteria on the true and false hit rates required by the police force . this article is organized as follows : first , we describe the use of lr scores for producing a ranked list of candidates from familial search on dna mixtures , on the basis of our previous work in chung et  al . . we then establish an estimate of the true hit rate of the investigation on a specific number of top - listed candidates , from which a novel strategy is developed for determining an appropriate scale of the investigation . the performance of the familial search on dna mixtures and the proposed strategy are demonstrated through an example using swedish data . given a crime trace observed as a dna mixture m contributed by the perpetrator and the victim , the familial search score for a particular individual j in the database d is defined as a lr of the following hypotheses : 
 hj :  the victim and a relative of individual j are contributors.hd :  the victim and one unknown person are contributors.note that the prosecution hypothesis hj here states that a relative of individual j , rather than the individual , is the contributor and therefore is not the same as the hypothesis formulated in chung et  al .   and chung and fung . in general , a lr can be calculated for each of the possible genetic relationships for every member of the database . however , the familial search is practically used only for parent / child and sibling relationships because it would become less effective for other degrees of relatedness below sibling that share less genetic similarity . hd :  the victim and one unknown person are contributors . at a particular locus l , denote ml as the set of alleles present in the mixture and vl and xjl as the genotype of the victim and individual j , respectively , for j    d and l  =  1, ... ,l . under linkage equilibrium assumption , presented the following formula for calculating the lr of hj versus hd : 
 1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : lr } { \rm lr}_j=\prod\limits_{l=1}^{l}{\frac{p(m_l|v_l , x_{j\,l},h_j)}{p(m_l|v_l , h_{\rm r } ) } } $ $ \end{document}where hr is the hypothesis that the victim and a random person contribute to the mixture . the computation of p(ml|vl , xjl , hj ) and p(ml|vl , hr ) is based on the q - function presented in hu and fung   and fung and hu : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{lll } & & { \kern-6pt } p(m_l|v_l , x_{j\,l},h_j)\\ & & { \kern6pt}=k_0q(2,u_l)+k_1\big(i_m(t_1)q(1,u_l\setminus\{t_1\ } ) \\ & & { \kern138pt}+i_m(t_2)q(1,u_l\setminus\{t_2\})\big ) \\ & & { \kern18pt}+k_2i_m(t_1)i_m(t_2)q(0,u_l\setminus\{t_1,t_2\ } ) \end{array}$$\end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ p(m_l|v_l , h_{\rm r})=q(2,u_l ) $ $ \end{document}where t1t2 is the genotype present in xjl , ul  =  ml    vl is the set of alleles present in ml but absent in vl , and im(t ) is the indicator function defined by im(t )  =  1 if t    m and 0 otherwise . the quantities ( k0,2k1,k2 ) are the kinship coefficients for the relationship considered in hj . in particular , ( k0,2k1,k2 ) take the values of ( 0.25,0.5,0.25 ) for full siblings and ( 0,1,0 ) for parent / child relationship . table  1 shows the computational formulae of q ( . , . ) for a two - person mixture at a particular autosomal locus l with k alleles a1,a2, ... ,ak and corresponding allele frequencies p1,p2, ... ,pk\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(\sum_{i=1}^{k}p_i=1)$\end{document } , under the assumption of hardy  weinberg equilibrium . table  1the calculating formulae of q(j , b ) for different combinations of mixture ml and arbitrary set of alleles b at a particular autosomal locus l with alleles a1,a2, ... ,ak and corresponding allele frequencies p1,p2, ... ,pk , under the assumption of hardy weinberg equilibriummlbq(0,b)q(1,b)q(2,b)ai1pi\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_i^2$\end{document}ai , aj1pi  +  pj\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(p_i+p_j)^2$\end{document}ai0pi\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_i^2 + 2p_ip_j$\end{document}aj0pj\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_j^2 + 2p_ip_j$\end{document}ai , aj , ak1pi  +  pj  +  pk\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(p_i+p_j+p_k)^2$\end{document}ai0pipi(pi  +  2pj  +  2pk)aj0pjpj(pj  +  2pi +  2pk)ak0pkpk(pk  +  2pi  +  2pj)ai , aj002pipjai , ak002pipkaj , ak002pjpkai , aj , ak , at1pi  +  pj  +  pk  +  pt\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(p_i+p_j+p_k+p_t)^2$\end{document}ai0pipi(pi  +  2pj  +  2pk  +  2pt)aj0pjpj(pj  +  2pi  +  2pk  +  2pt)ak0pkpk(pk  +  2pi  +  2pj  +  2pt)al0ptpt(pt  +  2pi  +  2pj  +  2pk)ai , aj002pipjai , ak002pipkai , at002piptaj , ak002pjpkaj , at002pjptak , at002pkptthe indices i , j , k , and t are pairwise distinct the calculating formulae of q(j , b ) for different combinations of mixture ml and arbitrary set of alleles b at a particular autosomal locus l with alleles a1,a2, ... ,ak and corresponding allele frequencies p1,p2, ... ,pk , under the assumption of hardy  weinberg equilibrium the indices i , j , k , and t are pairwise distinct the individuals in the database are ranked in descending order by the lr scores , and the individuals with highest scores will be preliminarily identified as the suspects that need further investigations by the police force . in case when rare alleles are present in the unexplained profiles ml    vl , the random match probability p(ml|vl , hr ) will become extremely small , and the unknown contributor s relative will be assigned a large lr score so that the perpetrator will be more likely to be identified . therefore , the performance of the familial search depends much on the observed mixture m  =  { ml , l  =  1, ... ,l } and the victim profile v  =  { vl , l  =  1, ... ,l}. to predict the performance of the search , cowen and thomson   suggested fitting the logistic regression model 
 2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : logistic_model } { \rm logit}(p(\theta_k=1))=\alpha-\beta \log_{10}{p(m|v , h_{\rm r } ) } $ $ \end{document}where k is an binary response variable taking the value of 1 if the relative of the unknown contributor is located within the top k profiles and 0 otherwise . the model links the outcome of a particular search that limits the investigations to the k top - listed individuals to the random match probability p(m|v , hr ) and therefore can be used to predict the performance of the search given the observed dna profiles ( m , v ) . in particular , if the estimated probability p(k  =  1 ) is too small even for large k , the familial search may fail to provide effective assistant for identifying the suspect . this information can practically guide the police force to make decisions on whether it is worthwhile to open an investigation just based on the familial search results . the use of logistic regression for evaluating the general performance of a familial database search will be demonstrated through an example . the logistic regression approach , however , uses only the random match probabilities but not the detail information provided in the dna evidence for a particular case . for practical crime cases , it may be also necessary to determine the most appropriate scale of the crime investigation by fully utilizing the available dna evidences . in the next section , we develop a novel strategy on deciding the number of individuals to be investigated after the familial search , based on an estimate of the hit rate of the crime investigation on a specific number of top - listed candidates . denote hd as the hypothesis that the sibling of the contributor is in the database , i.e. , hd  =    i    dhi . under the assumption that the relative of someone in the database contributes to the mixture , the posterior probability that the contributor is the relative of individual j , given the mixture profile m , victim profile v and database profiles xd  =  { xjl , j    d , l  =  1, ... ,l } can be expressed as 
 3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{lll } \label{eqn : posterior } & & { \kern-6pt } p(h_j|m , v,\mathbf{x_{\rm d}},h_{\rm d})\\ & & { \kern6pt } = \frac{p(m|v , x_j , h_j)p(h_j|h_{\rm d})}{\sum_{i\in d}p(m|v , x_i , h_i)p(h_i|h_{\rm d})}. \end{array } $ $ \end{document}whereexj  =  { xjl , l  =  1, ... ,l } for j    d. the proof of eq .  3 is given in the  appendix .  under most common scenarios , 3 becomes 
 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : posterior2 } p(h_j|m , v,\mathbf{x_{\rm d}},h_{\rm d})=\frac{p(m|v , x_j , h_j)}{\sum_{i\in d}p(m|v , x_i , h_i)}. $ $ \end{document}alternatively , other priors can also be used if the police authority has acquired more information about the family status details of each individual in the database . without loss of generality , suppose that the database is sorted in descending order according to the posterior probabilities such that 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{rll } p(h_1|m , v,\mathbf{x}_{\rm d},h_{\rm d})&\ge & p(h_2|m , v,\mathbf{x}_{\rm d},h_{\rm d})\\ & \ge & \cdots \ge p(h_n|m , v,\mathbf{x}_{\rm d},h_{\rm d } ) \end{array}$$\end{document}where n is the size of the database . for any integer 1    k    n , the function defined by 
 5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : hit_rate } q(k)=\sum\limits_{i=1}^{k}p(h_i|m , v,\mathbf{x_{\rm d}},h_{\rm d } ) $ $ \end{document}evaluates the probability that the relative of the contributor can be identified by investigating the top - k individuals with respect to their posterior probabilities , given that the relative is in the database . in other words , q(k ) represents the hit rate of the crime investigation on k top - listed candidates , thereby providing a means for developing a practical crime investigation strategy under different criteria . for instance , a simple strategy is to investigate the top - k individuals in the database sorted with respect to their lrs or posterior probabilities , with k determined by 
 6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : kp0 } k(p_0)= \min\{k>0:q(k)\ge p_0\ } $ $ \end{document}where p0 is the hit rate required by the police force . given a crime trace observed as a dna mixture m contributed by the perpetrator and the victim , the familial search score for a particular individual j in the database d is defined as a lr of the following hypotheses : 
 hj :  the victim and a relative of individual j are contributors.hd :  the victim and one unknown person are contributors.note that the prosecution hypothesis hj here states that a relative of individual j , rather than the individual , is the contributor and therefore is not the same as the hypothesis formulated in chung et  al .   and chung and fung . in general , a lr can be calculated for each of the possible genetic relationships for every member of the database . however , the familial search is practically used only for parent / child and sibling relationships because it would become less effective for other degrees of relatedness below sibling that share less genetic similarity . hd :  the victim and one unknown person are contributors . at a particular locus l , denote ml as the set of alleles present in the mixture and vl and xjl as the genotype of the victim and individual j , respectively , for j    d and l  =  1, ... ,l . under linkage equilibrium assumption , presented the following formula for calculating the lr of hj versus hd : 
 1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : lr } { \rm lr}_j=\prod\limits_{l=1}^{l}{\frac{p(m_l|v_l , x_{j\,l},h_j)}{p(m_l|v_l , h_{\rm r } ) } } $ $ \end{document}where hr is the hypothesis that the victim and a random person contribute to the mixture . the computation of p(ml|vl , xjl , hj ) and p(ml|vl , hr ) is based on the q - function presented in hu and fung   and fung and hu : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{lll } & & { \kern-6pt } p(m_l|v_l , x_{j\,l},h_j)\\ & & { \kern6pt}=k_0q(2,u_l)+k_1\big(i_m(t_1)q(1,u_l\setminus\{t_1\ } ) \\ & & { \kern138pt}+i_m(t_2)q(1,u_l\setminus\{t_2\})\big ) \\ & & { \kern18pt}+k_2i_m(t_1)i_m(t_2)q(0,u_l\setminus\{t_1,t_2\ } ) \end{array}$$\end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ p(m_l|v_l , h_{\rm r})=q(2,u_l ) $ $ \end{document}where t1t2 is the genotype present in xjl , ul  =  ml    vl is the set of alleles present in ml but absent in vl , and im(t ) is the indicator function defined by im(t )  =  1 if t    m and 0 otherwise . the quantities ( k0,2k1,k2 ) are the kinship coefficients for the relationship considered in hj . in particular , ( k0,2k1,k2 ) take the values of ( 0.25,0.5,0.25 ) for full siblings and ( 0,1,0 ) for parent / child relationship . table  1 shows the computational formulae of q ( . , . ) for a two - person mixture at a particular autosomal locus l with k alleles a1,a2, ... ,ak and corresponding allele frequencies p1,p2, ... ,pk\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(\sum_{i=1}^{k}p_i=1)$\end{document } , under the assumption of hardy  weinberg equilibrium . table  1the calculating formulae of q(j , b ) for different combinations of mixture ml and arbitrary set of alleles b at a particular autosomal locus l with alleles a1,a2, ... ,ak and corresponding allele frequencies p1,p2, ... ,pk , under the assumption of hardy weinberg equilibriummlbq(0,b)q(1,b)q(2,b)ai1pi\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_i^2$\end{document}ai , aj1pi  +  pj\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(p_i+p_j)^2$\end{document}ai0pi\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_i^2 + 2p_ip_j$\end{document}aj0pj\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_j^2 + 2p_ip_j$\end{document}ai , aj , ak1pi  +  pj  +  pk\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(p_i+p_j+p_k)^2$\end{document}ai0pipi(pi  +  2pj  +  2pk)aj0pjpj(pj  +  2pi +  2pk)ak0pkpk(pk  +  2pi  +  2pj)ai , aj002pipjai , ak002pipkaj , ak002pjpkai , aj , ak , at1pi  +  pj  +  pk  +  pt\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(p_i+p_j+p_k+p_t)^2$\end{document}ai0pipi(pi  +  2pj  +  2pk  +  2pt)aj0pjpj(pj  +  2pi  +  2pk  +  2pt)ak0pkpk(pk  +  2pi  +  2pj  +  2pt)al0ptpt(pt  +  2pi  +  2pj  +  2pk)ai , aj002pipjai , ak002pipkai , at002piptaj , ak002pjpkaj , at002pjptak , at002pkptthe indices i , j , k , and t are pairwise distinct the calculating formulae of q(j , b ) for different combinations of mixture ml and arbitrary set of alleles b at a particular autosomal locus l with alleles a1,a2, ... ,ak and corresponding allele frequencies p1,p2, ... ,pk , under the assumption of hardy  weinberg equilibrium the indices i , j , k , and t are pairwise distinct the individuals in the database are ranked in descending order by the lr scores , and the individuals with highest scores will be preliminarily identified as the suspects that need further investigations by the police force . in case when rare alleles are present in the unexplained profiles ml    vl , the random match probability p(ml|vl , hr ) will become extremely small , and the unknown contributor s relative will be assigned a large lr score so that the perpetrator will be more likely to be identified . therefore , the performance of the familial search depends much on the observed mixture m  =  { ml , l  =  1, ... ,l } and the victim profile v  =  { vl , l  =  1, ... ,l}. to predict the performance of the search , cowen and thomson   suggested fitting the logistic regression model 
 2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : logistic_model } { \rm logit}(p(\theta_k=1))=\alpha-\beta \log_{10}{p(m|v , h_{\rm r } ) } $ $ \end{document}where k is an binary response variable taking the value of 1 if the relative of the unknown contributor is located within the top k profiles and 0 otherwise . the model links the outcome of a particular search that limits the investigations to the k top - listed individuals to the random match probability p(m|v , hr ) and therefore can be used to predict the performance of the search given the observed dna profiles ( m , v ) . in particular , if the estimated probability p(k  =  1 ) is too small even for large k , the familial search may fail to provide effective assistant for identifying the suspect . this information can practically guide the police force to make decisions on whether it is worthwhile to open an investigation just based on the familial search results . the use of logistic regression for evaluating the general performance of a familial database search will be demonstrated through an example . the logistic regression approach , however , uses only the random match probabilities but not the detail information provided in the dna evidence for a particular case . for practical crime cases , it may be also necessary to determine the most appropriate scale of the crime investigation by fully utilizing the available dna evidences . in the next section , we develop a novel strategy on deciding the number of individuals to be investigated after the familial search , based on an estimate of the hit rate of the crime investigation on a specific number of top - listed candidates . denote hd as the hypothesis that the sibling of the contributor is in the database , i.e. , hd  =    i    dhi . under the assumption that the relative of someone in the database contributes to the mixture , the posterior probability that the contributor is the relative of individual j , given the mixture profile m , victim profile v and database profiles xd  =  { xjl , j    d , l  =  1, ... ,l } can be expressed as 
 3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{lll } \label{eqn : posterior } & & { \kern-6pt } p(h_j|m , v,\mathbf{x_{\rm d}},h_{\rm d})\\ & & { \kern6pt } = \frac{p(m|v , x_j , h_j)p(h_j|h_{\rm d})}{\sum_{i\in d}p(m|v , x_i , h_i)p(h_i|h_{\rm d})}. \end{array } $ $ \end{document}whereexj  =  { xjl , l  =  1, ... ,l } for j    d. the proof of eq .  3 is given in the  appendix .  under most common scenarios , it is sensible to assume a uniform prior for hi so that eq .  3 becomes 
 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : posterior2 } p(h_j|m , v,\mathbf{x_{\rm d}},h_{\rm d})=\frac{p(m|v , x_j , h_j)}{\sum_{i\in d}p(m|v , x_i , h_i)}. $ $ \end{document}alternatively , other priors can also be used if the police authority has acquired more information about the family status details of each individual in the database . without loss of generality , suppose that the database is sorted in descending order according to the posterior probabilities such that 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{rll } p(h_1|m , v,\mathbf{x}_{\rm d},h_{\rm d})&\ge & p(h_2|m , v,\mathbf{x}_{\rm d},h_{\rm d})\\ & \ge & \cdots \ge p(h_n|m , v,\mathbf{x}_{\rm d},h_{\rm d } ) \end{array}$$\end{document}where n is the size of the database . for any integer 1    k    n , the function defined by 
 5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : hit_rate } q(k)=\sum\limits_{i=1}^{k}p(h_i|m , v,\mathbf{x_{\rm d}},h_{\rm d } ) $ $ \end{document}evaluates the probability that the relative of the contributor can be identified by investigating the top - k individuals with respect to their posterior probabilities , given that the relative is in the database . in other words , q(k ) represents the hit rate of the crime investigation on k top - listed candidates , thereby providing a means for developing a practical crime investigation strategy under different criteria . for instance , a simple strategy is to investigate the top - k individuals in the database sorted with respect to their lrs or posterior probabilities , with k determined by 
 6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \label{eqn : kp0 } k(p_0)= \min\{k>0:q(k)\ge p_0\ } $ $ \end{document}where p0 is the hit rate required by the police force . we illustrate here the performance of the familial search and the use of eqs .  5 and 6 . an offender database consisted of 50,000 unrelated dna profiles is generated according to the swedish allele frequencies given in montelius et al . , at the 10 sgm plus str loci which are commonly used in european national dna databases . the first 1,000 profiles in the database are used to represent the profiles of the relatives of the perpetrators in 1,000 different crime cases . for each of the first 1,000 profiles , a related profile ( full sibling or parent / child ) is generated to represent the unobserved profile of the unknown perpetrator in that particular case . the observed dna mixtures are produced by mixing these relative profiles with independently generated victim profiles . for each case , a familial search is performed , and a list of 50,000 lr scores is calculated by using eq .  1 . the rank of the lr score of the true relative of the unknown contributor indicates at least how many top - listed individuals must be investigated in order to successfully identify the relative of the perpetrator . in general , such ranks recorded from the 1,000 cases can provide reasonable assessment of the effectiveness of the familial search . for example , there are 704 cases in which the true sibling of the unknown contributor is within the top 100 individuals in the database ranked by the lr scores , showing a chance of 70.4% that the sibling of the unknown contributor can be successfully identified if the 100 top - listed individuals are investigated , provided that the unknown contributor is a sibling of one of the individuals in the database . the chance of successful identification of the parent / child of the unknown contributor from the top 100 profiles is 88.0% . 1 shows the empirical probabilities of identifying the relative of the contributor using lr scores . table  2empirical probabilities of identifying the parent / child and full sibling of the contributor in the top k profiles ranked by likelihood ratio scoresk% of identification in top k profilesparent 1empirical probability of identifying the parent / child and full sibling of the contributor in the topk profiles ranked by likelihood ratio scores in a one - unknown mixture case empirical probabilities of identifying the parent / child and full sibling of the contributor in the top k profiles ranked by likelihood ratio scores empirical probability of identifying the parent / child and full sibling of the contributor in the topk profiles ranked by likelihood ratio scores in a one - unknown mixture case clearly as shown in fig .  1 , the empirical probabilities indicate a more effective search for the parent / child of the unknown contributor than the sibling . it is interesting to note that the empirical probabilities of successful identification are almost linearly related with logarithm of k , the number of top - listed individuals to be investigated . this suggests that enlarging the scale of the investigation will have little improvement to the performance of familial search if the original scale is already substantially large . in roughly 46% of the cases , the sibling of the unknown contributor can be found within the top 20 profiles , while in 56% of the cases , investigating the top 10 profiles is already enough for identifying the parent / child of the unknown contributor . , the performance of the search here is slightly inferior as we use a larger database and less loci for profiling in this work . however , it is comparable to the result of curran and buckleton   in a single - source case which also use a system of 10 str loci but a smaller database of about 24,000 profiles . therefore , in cases with dna mixture as part of the available clue , the familial dna database search can be still applied for identifying the suspect , with the false hit rate as low as in single - source cases . following cowen and thomson , logistic regression models in the form of eq . 2 are fitted to the data to estimate the dependence of the outcome with respect to the observed profiles ( m , v ) . table  3 listed the coefficients of the fitted models , with p values all smaller than 0.001 . the performance of the search is strongly associated with the value of the random match probability . as shown in fig .  2 , the smaller the random match probability , the more likely the relative of the unknown contributor can be identified . in particular , if p(m|v , hr ) is less than 10 , which is about the median of the random match probabilities from the simulated profiles , there will be a 41% chance to identify the sibling and a 63% chance to identify the parent / child within the top 10 profiles of the database . the performance would be worse if p(m|v , hr ) is as low as 10 . the estimated probabilities of successfully identifying the parent / child and sibling are 24% and 17% , respectively . in such cases , an investigation on just 10 top - listed individuals from the familial search is definitely not enough . table  3estimated parameters of the logistic regression models for the probability of identifying the relative of the contributor in the top k profilesrelationshipkparameter estimate ( se)parent / child5  6.0635  ( 0.5733)0.5114  ( 0.0492)10  6.1751  ( 0.5970)0.5593  ( 0.0522)20  6.1630  ( 0.6355)0.5969 ( 0.0565)full sibling10  5.1906  ( 0.5513)0.4014  ( 0.0468)50  4.7499  ( 0.5610)0.4477  ( 0.0492)100  3.8619  ( 0.5926)0.4173  ( 0.0528)fig . 2predicted probability of identifying the a parent / child and b full sibling of the unknown contributor in the top k profiles estimated parameters of the logistic regression models for the probability of identifying the relative of the contributor in the top k profiles predicted probability of identifying the a parent / child and b full sibling of the unknown contributor in the top k profiles to determine an appropriate scale of the investigation , we can apply the proposed strategy . for illustrative purpose , we consider a search that aims at identifying the parent / child or the full sibling of the contributor from the database , though the basic principle can be applied to any kind of relationship . figure  3 shows the average of the estimated hit rates evaluated by using the q(k ) function in eq .  5 and the empirical hit rates obtained as in fig .  1 . as can be seen , the average of the estimated hit rates are close to the empirical hit rates , suggesting reliable estimates of the hit rates by using q(k ) . the number of top - listed individuals that needs investigation can be determined by k(p0 ) defined in eq .  6 for a particular required hit rate p0 . the empirical distributions of k(p0 ) for various values of p0 the distributions are all skewed toward the right as the averages are all greater than the medians . for the search that aims at identifying the parent / child of the contributor , large proportions of k(p0 ) are found to be less than 100 with respect to p0    0.8 . when the required hit rate p0 is as high as 0.9 , there are still about 75% of the cases in which the values of k(p0 ) are less than 200 , i.e. , further investigations on less than 200 top - listed candidates will be suggested . therefore , the investigation can be controlled under a reasonable scale even though a high hit rate is usually required . 3probabilities of identifying the a parent / child and b full sibling of the unknown contributor in the top k profiles estimated using empirical hit rates ( solid line ) and the average of the q(k ) function ( dashed line)table  4summary statistics of k(p0 ) from simulation resultsp0empirical probabilityaverage of k(p0)percentiles of k(p0 ) distribution10th25thmedian75th90thparent / child0.50.55317.1141123440.60.62525.8261838690.70.72445.131031621100.80.81675.2519531041790.90.925137114197192314full sibling0.50.51241.42112962970.60.61287.0527661301970.70.74419018701522834130.80.851458572083946699120.90.9561,3933137871,3131,9602,516the empirical hit rates are the probabilities of identifying the parent / child and full sibling of the contributor within the top - k(p0 ) individuals probabilities of identifying the a parent / child and b full sibling of the unknown contributor in the top k profiles estimated using empirical hit rates ( solid line ) and the average of the q(k ) function ( dashed line ) summary statistics of k(p0 ) from simulation results the empirical hit rates are the probabilities of identifying the parent / child and full sibling of the contributor within the top - k(p0 ) individuals on the other hand , the values of k(p0 ) determined in searching for the full sibling of the contributor are relatively larger . if p0 =  0.9 is required , the values of k(p0 ) are greater than 1,000 in the majority of the cases , indicating that it may not be feasible to achieve such a high hit rate with limited law enforcement resources . the result is expected as the performance in searching for the full sibling should be inferior to the cases in searching for the parent / child . nevertheless , lowering the requirement on the hit rate to p0    0.7 can still reduce the values of k(p0 ) to less than 150 in about 50% of the cases , which are more feasible scales . based on this information , the decision to open an investigation can be taken by the police force after thorough consideration on the crime severity and the law enforcement resources . it is remarkable to note that the empirical hit rates based on investigating the top - k(p0 ) individuals are all greater than the required hit rate p0 , further justifying our approach in providing information for the police force to determine an appropriate scale of crime investigation based on the results from a familial database search . this article illustrates how the traditional familial database search methods that were used only for single - source samples can be extended so as to handle mixture cases . it is demonstrated that the familial search applied to two - person mixture cases can perform as good as in single - source cases by using the lr scoring scheme . the results presented in previous section are based on the most common scenario when the mixed stain is originated from the victim whose dna profile is available . for two - unknown mixture cases , the set of possible profiles for the unknown contributors will be less restrictive . as a result , the discriminatory power of the likelihood ratio score for database search in two - unknown mixture cases will become relatively lower , comparing to one - unknown mixture cases and single - source cases . to demonstrate this , consider the following hypotheses for a two - unknown mixture case : 
 hj :  a relative of individual j and an unrelated unknown person are contributors.hd :  two unrelated unknown persons are contributors.and the corresponding lr can be calculated as 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \rm lr}_j=\prod\limits_{l=1}^{l}{\frac{p(m_l|x_{j\,l},h_j)}{p(m_l|h_{\rm r } ) } } $ $ \end{document}where 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{lll } p(m_l|x_{j\,l},h_j)&=&k_0q(4,m_l)\\&&+k_1\big(i_m(t_1)q(3,m_l\setminus\{t_1\})\\ & & \qquad\;+\ ; i_m(t_2)q(3,m_l\setminus\{t_2\})\big ) \\ & & + k_2i_m(t_1)i_m(t_2)q(2,m_l\setminus\{t_1,t_2\ } ) , \end{array}$$\end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ p(m_l|h_{\rm r})=q(4,m_l ) . $ $ \end{document}figure  4 shows the empirical search results on a two - unknown case using the same generated database . as expected , the performance is substantially inferior to the case in which the victim profile is available . there is only a 36% chance to identify the sibling and a 54% chance to identify the parent / child within the top 100 profiles of the database . 4empirical probability of identifying the parent / child and full sibling of the contributor in the topk profiles ranked by likelihood ratio scores in a two - unknown mixture case hj :  a relative of individual j and an unrelated unknown person are contributors . hd :  two unrelated unknown persons are contributors . empirical probability of identifying the parent / child and full sibling of the contributor in the topk profiles ranked by likelihood ratio scores in a two - unknown mixture case in addition to the logistic regression approach that models the predicted hit rate as a function of the random match probability , we have derived a formula that can accurately estimate the hit rate for a specific case . using this formula , a simple strategy is developed to determine the least number of individuals who should be included in the crime investigation , according to the desired hit rate required by the police force . it should be noted that the formula for estimating the hit rate , as well as the proposed strategy , is also applicable in single - source cases . as shown by the numerical example , the proposed strategy can provide useful information aiding the police force for deciding on the scale of investigation , thereby achieving desirable crime - solving rate with reasonable cost . as commented by chakraborty and ge , it is necessary to point out that the familial dna database search is an auxiliary tool for crime solving when there is no clue on the source of the crime trace , rather than playing a decisive role to bring the suspects identified through database search to trials . the evidentiary value of the cold hits from familial search should serve as supporting reference to aid the jury to make their decision , provided that a court case is raised after the investigation on the cold - hit suspects . to clearly present our idea for the application of familial search to mixture cases , several key assumptions have been made in the works presented in this article . these include the linkage equilibrium that assumes independence of alleles across all loci and the hardy  weinberg equilibrium that assumes the independence between the two alleles of a genotype at a particular locus . the former is usually guaranteed by the proper choice of str loci while the latter may fail to apply in case when the contributors of the mixture come from a small population or from different ethnic groups . although it is not difficult to modify the formulae presented in earlier sections and table  1 to handle allele dependence , the impact of the deviation from hardy  weinberg equilibrium to the performance of the familial search is yet to be studied . besides , the presented method is based on the common scenario that there is only one perpetrator committing the crime and a two - person mixture is included as part of the dna evidence . more general approach that considers multiple perpetrators as in group rape cases can make familial search a more useful crime - solving tool . therefore , another possible direction for future work is to develop general methodologies that can handle the complication arises due to the presence of multiple perpetrators . under the assumption that the sibling of the contributor is someone in the database , the posterior probability of hj given the dna evidences m , v , xd can be evaluated by using the bayes rule : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{rll } & & { \kern-6pt } p(h_j|m , v,\mathbf{x_{\rm d}},h_{\rm d})\\ & & { \kern6pt}=\frac{p(m|v,\mathbf{x_{\rm d}},h_j , h_{\rm d})p(h_j|h_{\rm d})}{\sum_{i=1}^n p(m|v,\mathbf{x_{\rm d}},h_i , h_{\rm d})p(h_i|h_{\rm d } ) } \\ & & { \kern6pt}=\frac{p(m|v,\mathbf{x_{\rm d}},h_j)p(h_j|h_{\rm d})}{\sum_{i\in d } p(m|v,\mathbf{x_{\rm d}},h_i)p(h_i|h_{\rm d } ) } \end{array}$$\end{document}where the equality results from the fact that hi    hd  =  hi for i    d and p(hi|hd )  =  0 for i    d. when hi holds true , the mixture m will depend only on v and xi and is unrelated to the profiles of other individuals in the database . substituting p(m|v , xd , hi )  =  p(m|v , xi , hi ) into the above equation leads to eq .