2011][^1], but such methods work only for outbred pedigrees and fail to capture more complicated relationships, as we see later in this paper. Forth, there has been work on inferring kinship coefficients from admixed individuals, [@Thornton2012]. Our EM method has a running time, for a pair of individuals, that is linear in the number of sites. To estimate a kinship coefficient matrix for $n$ individuals, our method is quadratic in $n$, due to the number of pairs of individual. Our method avoids the assumption of having a known pedigree, taking only the genotypes as an input. We introduce one latent variable for each pair of individuals and each site, which encodes our uncertainty on the identity by descent states. The distribution over these states for each pair of individuals is learned using a simple expectation maximization (EM) algorithm and provides an informative, yet concise and learnable, summary of the relationship of each pair of individuals. We show that this EM algorithm quickly converges to an accurate estimate. We assess the accuracy in two ways. First, we measure the accuracy of the reconstructed latent variables using simulated pedigree data (the true value of the latent variable can be computed efficiently from the held-out pedigree structure). We asses the kinship estimates compared with PLINK and with a covariance-based estimate. Second, we demonstrate that our method can be used to correct GWAS for cryptic relatedness without assuming the knowledge of the pedigree structure.

Background {#sec:background}
==========

A pedigree is an annotated directed acyclic graph where the nodes are diploid individuals and the edges are directed from parent to child (Figure \[fig:ped-ibd\]A). The graph is acyclic, because no individual can be their own ancestor. The pedigree graph is typically annotated with the gender (male or female), affection status, and the genotype status of the individuals. The gender is used to enforce on the graph the marriage constraint that each individual with parents in the graph has at most one parent of each gender. The *founders* are the individuals who do not have parents in the pedigree while *non-founders* are all the individuals having parent(s) in the pedigree graph. The *affection status* indicates whether an individual has a disease or phenotype. The genotype status indicates whether the individual has available genotypes.

The pedigree graph depicts familial relationships between individuals. Each diploid individual inherits one copy of
2011][^1 ], but such methods work only for outbred pedigrees and fail to get more complicated relationship, as we see later in this paper. away, there has been work on inferring kinship coefficient from admixed individual,   [ @Thornton2012 ]. Our EM method have a running time, for a couple of individuals, that is linear in the phone number of web site. To estimate a kinship coefficient matrix for $ n$ individuals, our method is quadratic in $ n$, due to the number of pair of individual. Our method avoids the presumption of having a known pedigree, taking merely the genotypes as an input. We introduce one latent variable for each pair of person and each site, which encodes our uncertainty on the identity by descent states. The distribution over these states for each pair of individuals is learned using a simple expectation maximization (EM) algorithm and provides an informative, yet concise and learnable, summary of the relationship of each pair of individuals. We prove that this EM algorithm quickly converges to an accurate appraisal. We measure the accuracy in two room. First, we measure the accuracy of the reconstructed latent variable star using simulated pedigree data (the true value of the latent variable star can be computed efficiently from the held - out pedigree structure). We asses the kinship estimates compared with PLINK and with a covariance - based estimate. Second, we demonstrate that our method can be used to chastise GWAS for cabalistic relatedness without simulate the knowledge of the pedigree structure. 

 Background { # sec: setting } 
 = = = = = = = = = = 

 A pedigree is an annotated directed acyclic graph where the nodes are diploid person and the edges are directed from parent to child (Figure   \[fig: ped - ibd\]A). The graph is acyclic, because no individual can be their own ancestor. The pedigree graph is typically annotated with the sex (male or female), affection status, and the genotype status of the individuals. The gender is used to enforce on the graph the marriage restraint that each individual with parents in the graph has at most one parent of each sex. The * founders * are the individuals who do not take parents in the pedigree while * non - founder * are all the individual having parent(s) in the pedigree graph. The * affection status * indicates whether an individual has a disease or phenotype. The genotype status indicates whether the individual have available genotypes. 

 The pedigree graph picture familial relationships between individuals. Each diploid individual inherits one copy of
2011][^1], buh such methods work only for outbred pebugrees and fzil to cxpture more complicated relavionwhips, as we see later in thks paper. Vorth, thwre ias been work on inferriky kinamip cmxfficients from admixed itdividuals, [@Thorttun2012]. Our EM method has a running time, fjr a paor of individualf, thse is linear in the number of sites. To sstimatt a kinship coeffivient matrix for $n$ individkals, our method is quafratic in $n$, due ro the numbef of pairs of individuzl. Our method avoids the assumpgion pf having q jnoat pedigree, vaking only the gekptypes as an onput. We introcucx onw latent variable for each pair of indididuals atd each site, which wnxodes our uncdetakntg pn the ifenvity by desdent states. The distribution ofew these states for ewcr pair of individuals is learned using d sjmple expectation maximuzation (EM) algorithm wnd proviqes an informative, yet concise and learnable, summdry oh ghe rclatkindhip of each pair of individuals. We show that fhos EM algorithm zuickly confeggrf to an accurxte escjmzte. We assess the wccuraci in tqo ways. Fyrst, we measure the accuracy of the reconstguctwd latent variabled using simblated pedibree data (the true valuz of tge latent vwriable czv be computed efwicpentny from ufe held-out pedigrqe structnre). Wz asses ghe linshi[ estimated comidred with PLINK anf witk a cmvariance-bwsed estimate. Second, we demonstcete that our kedhog can be used to correct GRAS for cryptie relatebness dithout asauming vhe knowledgq of the pedictee structure.

Uackgrounq {#sex:bacjground}
==========

X pedigree is am annotatvd dirwcted acyclic grapm whete the nodes are bnpooid individualx avd ehv evges wse directed xrom pafrnt tu child (Figmre \[wig:prd-ibd\]A). The graph is dcycmic, because no indovldual can be theiw own ancestot. The pedigree grakh is vypicauly amnojated with the gender (male or fsmale), afffctlon status, anq thc gegotype stacus of the individuals. The gender is usev to enforce on the graph the marriage conxuraint that eech inqividual fith parents in the traph has at most one parent of each gsnder. Dhe *flunders* are the individuals who do not have parents in the pedigree while *non-fonnqers* are alm thr indpvibuaks havyng 'acent(s) in the pedibree graph. The *affection status* indicatas whether an individual has a dksease or phevotype. The genotype staths indicstes whether the individual has available genotypes.

The iedigree traph depicna hamilial relationshils betceen ibdividuqls. Cach diploid ivdififual inhegits oue copy of
2011][^1], but such methods work only for and to capture complicated relationships, as paper. there has been on inferring kinship from admixed individuals, [@Thornton2012]. Our EM has a running time, for a pair of individuals, that is linear in number of sites. To estimate a kinship coefficient matrix for $n$ individuals, our is in due the number of pairs of individual. Our method avoids the assumption of having a known pedigree, only the genotypes as an input. We introduce latent variable for each of individuals and each site, encodes uncertainty on identity descent The distribution over states for each pair of individuals is learned using a simple expectation maximization (EM) algorithm and provides informative, yet learnable, summary the of pair of individuals. that this EM algorithm quickly converges estimate. We assess the accuracy in two ways. we measure accuracy of the reconstructed latent variables simulated pedigree data (the true value of the variable can be computed efficiently from the held-out pedigree structure). We asses the kinship estimates PLINK and with a estimate. Second, we that method be to correct for cryptic relatedness without assuming the knowledge of the pedigree structure. {#sec:background} ========== A pedigree is an annotated directed acyclic graph nodes diploid individuals and edges are directed from to (Figure \[fig:ped-ibd\]A). The graph because individual own The graph is typically annotated the gender (male or female), status, and the genotype gender is used to enforce on the graph marriage constraint that each individual with parents the graph has at most one parent of each gender. The *founders* the individuals not have parents in the pedigree while *non-founders* all the individuals having in the pedigree graph. The *affection status* indicates whether individual a disease phenotype. The genotype indicates whether the has available genotypes. graph depicts relationships individuals. inherits one copy of
2011][^1], but such methods work only for Outbred pedIgreeS anD faIl To caPturE more complicatED relAtionships, as we see later In thiS pAPer. FORtH, therE has beeN WoRK On iNfErRinG kINsHip coEffIcients From admixeD inDiViduals, [@ThornTOn2012]. our EM methoD haS a running timE, foR a pair Of IndIViduaLs, tHat is Linear IN the nuMber of sitEs. tO estimATe a kinsHIP cOeffIcient matrix for $n$ iNDiVIduals, our methoD is quaDrATiC IN $n$, dUe tO the number Of Pairs OF indiviDUaL. oUR meTHod avoids the aSsumption of HAviNg a knoWn PedIGree, taKing oNlY The Genotypes as An inPut. We intrOduce oNE latent VAriable For eacH paIr oF indIViDuAls AnD EacH SiTe, wHIch Encodes oUr UnCertaInty ON THE ideNtiTy by DesceNt states. The diStrIbutIOn oVer thEse stAtes FoR each Pair of IndivIdUals is learned usIng a Simple expEctAtIon MaXimizATion (EM) AlgOriThm and pRovides AN inFoRMATiVe, yet concise and leaRnABLe, Summary oF the reLAtIoNShip of eaCh PaiR of iNDIviduAls. WE ShOw that thIs EM alGOrItHm quickLy ConverGeS to An aCcuraTE estImate. WE assess tHe accURacy in two ways. FIRst, we measure tHE aCCUrACy of The ReconstructEd laTEnt vAriaBLeS usINg simUlateD pEDiGRee data (the true value Of The latEnt vaRiable can be coMputed effiCIENtly from The hELd-OUt pedigree struCture). we asses the KInship esTimatEs comparEd with PLInk And with a CovAriAncE-baSED eStimate. Second, WE DemoNsTrate thAt oUr methoD caN be UseD to CoRrect GWAS For cryptIc ReLaTeDneSs witHOut assumInG thE kNowLedge OF the peDigreE strUcTuRE.

BaCkgrounD {#SeC:BAckgRoUnD}
==========

A peDigReE is an AnnoTAteD directEd acyclic GraPH wheRe ThE nodes aRe diploid indiViDuals and thE eDgeS are diRECted from Parent to child (Figure \[fig:pED-ibd\]A). ThE grAph is AcycLic, becausE no IndiviDuaL Can be tHeir owN anceStOr. THE PedigREE gRapH iS typically ANNotAted wItH the Gender (mAle or female), affectiON stAtus, and the genOtyPe stATUs Of tHE iNDivIdUAls. tHE gender is used to Enforce on tHe GRaPh the marriAGe cOnStraint That eacH indiVIdual wiTh parents In the grapH hAs at MOSt oNe parent of Each gendEr. The *founDErs* arE ThE indiVidUals whO dO noT have ParentS In tHe pedIgree wHiLe *non-fOundeRs* Are all thE individuals having parenT(s) in thE pediGreE graph. The *AffECtiOn status* iNdicAtes whetheR an IndIviduAl hAS a disEase OR pHenOType. THe geNOtype statUS iNdiCATeS whether the INDIviDual hAs aVAilablE genOtypes.

The pedigree GRaph depicts famIliaL RElaTioNShipS bEtween individuAls. eaCH Diploid iNdIvidual inheRits one cOpY Of
2011][^1], but such method s work onl y for  ou tbr ed  ped igre es and fail to capt ure more complicated r elati on s hips ,  a s wesee lat e ri n  th is  p ape r. Fo rth,the re hasbeen workonin ferring kins h ip  coefficie nts  from admixe d i ndivid ua ls, [@Tho rnt on201 2]. Ou r  EM me thod hasar unning time, f o r  a  pai r of individuals, th a t is linear in  the n um b er o f s ite s. To esti ma te ak inshipc oe f f i cie n t matrix for$n$ individ u als , ourme tho d  is qu adrat ic in$n$, due to  the  number o f pair s  of ind i vidual.  Our m eth odavoi d sth e a ss u mpt i on  of hav ing a kn ow npedig ree, t a k ingonl y th e gen otypes as aninp ut.W e i ntrod uce o ne l at ent v ariabl e for  e ach pair of ind ivid uals andeac hsit e,  whic h  encod esour  uncert ainty o n  th ei d e nt ity by descent sta te s .  T he distr ibutio n  o ve r  these s ta tes  for e ach p airo findividu als is le ar ned usi ng  a sim pl e e xpe ctati o n ma ximiza tion (EM ) alg o rithm and prov i des an inform a ti v e ,y et c onc ise and lea rnab l e, s umma r yoft he re latio ns h ip of each pair of ind iv iduals . Weshow that thi s EM algor i t h m quickl y co n ve r ges to an accu rateestimate.W e assess  theaccuracy  in two w a y s. First , w e m eas ure t he  accuracy oft h e re co nstruct edlatentvar iab les  us in g simulat ed pedig re eda ta  (t he tr u e valueof  th elat ent v a riable  canbe c om pu t edefficie n tl y from  t he  hel d-o ut  pedi gree str ucture) . We asse s t h e ki ns hi p estim ates compared  w ith PLINKan d w ith ac o variance -based estimate. Second ,  we dem ons trate  tha t our met hod  can b e u s ed tocorrec t GWA Sfor c rypti c re lat ed ness witho u t  as sumin gtheknowled ge of the pedigree str ucture.

Back gro und{ # se c:b a ck g rou nd } 
== = = ======

A pedig ree is anan n ot ated direc t edac yclic g raph wh ere t h e nodes  are dipl oid indiv id uals a ndthe edgesare dire cted from paren t  t o chi ld(Figur e\[f ig:pe d-ibd\ ] A).  Thegraphis  acycl ic, b ec ause noindividual can be their  own a ncest or.  The pedi gre e  gr aph is ty pica lly annota ted  wi th th e g e nder(mal e  o r f e male) , af f ection st a tu s,a n dthe genotyp e s tat us of  th e  indiv idua ls. The gender is used to enforc e on t hegra p h th emarriage const rai nt t hat each  i ndividual w ith pare nt s  in t he gra ph has  at mos t on e  paren t of  ea ch gender . T he *founde rs *a re the  ind iv iduals  who d o  not h ave parents in t he pe d i greew hil e *no n- founder s * ar e all theindividuals  havin g pa rent( s) in t he  pedig ree  g raph. The* affection  stat us* ind ic ates  wh etheran i n d ividu al h as  adisease o r ph e no ty p e.Thegenot yp e st atus indi c ates whe the r  the in di vid u a l hasa va i l able genot ype s.

T h e  pedigreeg raph de p ictsfamili al rel ationsh i psbe tween i ndi v i duals. Ea ch diploi d  in di vidu al inher it s on ecop yof
2011][^1], but_such methods_work only for outbred_pedigrees and_fail_to capture_more_complicated relationships, as_we see later_in this paper. Forth,_there has been_work_on inferring kinship coefficients from admixed individuals, [@Thornton2012]. Our EM method has a running time,_for_a pair_of_individuals,_that is linear in the_number of sites. To estimate_a kinship_coefficient matrix for $n$ individuals, our method is_quadratic_in $n$, due_to the number of pairs of individual. Our method_avoids the assumption of having a_known pedigree, taking_only_the_genotypes as an input._We introduce one latent variable for_each pair of individuals and each_site, which encodes our uncertainty on the_identity by descent states. The distribution_over these states for each_pair of_individuals is learned using a_simple expectation maximization_(EM) algorithm_and provides an_informative, yet concise and learnable, summary_of the relationship_of each pair of individuals. We_show_that this EM_algorithm_quickly_converges to_an accurate estimate._We_assess the_accuracy_in two ways. First, we measure_the_accuracy of the reconstructed latent variables using_simulated pedigree data (the_true_value of the latent_variable can be computed efficiently_from the held-out pedigree structure). We_asses the_kinship estimates_compared with PLINK and with a covariance-based estimate. Second, we demonstrate_that our method can be used_to correct GWAS for_cryptic relatedness_without_assuming the knowledge_of_the pedigree_structure.

Background {#sec:background}
==========

A pedigree is an annotated directed_acyclic graph_where the nodes are diploid individuals_and the edges are_directed_from parent to child (Figure \[fig:ped-ibd\]A). The_graph is acyclic, because no individual_can be their own ancestor._The_pedigree_graph is typically annotated with_the gender (male or female), affection_status, and the_genotype status of the individuals. The gender_is_used to enforce on the graph_the_marriage constraint that each individual with_parents_in_the graph has at most_one parent of each gender. The_*founders* are the individuals who do not have parents_in the pedigree_while *non-founders* are all the_individuals_having_parent(s) in the pedigree graph. The *affection status* indicates whether_an individual_has a disease_or phenotype. The genotype status indicates whether the individual has_available genotypes.

The pedigree graph depicts familial relationships_between individuals. Each diploid individual inherits one copy of