structure 2.3- ÖÐÎÄ

ÊÇ׿ÏȵÄʵ¼Ê±ÈÀýÊÇÓÐÆ«µÄ¡£ÄǸöÀý×ÓÒ²ÊÇÎÒÃÇÓÃÕæÕýµÄÊý¾ÝµÄ¸üнüµÄ¾­ÀúµÄ´ú±í¡£

Õâ·¢ÉúÒòΪÔÚȱ·¦Èκγ̶ȷDzôºÍµÄ¸öÌåµÄÇé¿öÏ£¬¿ÉÄÜÓÐһЩ²»¿Éʶ±ðÐÔ£¨nonidentifiability£©£¬ÔÚÄÇÀïÓпÉÄܰѵÈλ»ùÒòƵÂÊÍÆµÃ¸ü¿ª£¬°Ñ»ìºÏ±ÈÀý¼·Ñ¹µ½Ò»Æð£¨·´¹ýÀ´Ò²ÊÇÕâÑù£©£¬»ñµÃ¼¸ºõÒ»ÑùµÄÄ£ÐÍÄâºÏ¡£µ±ÓÐÇ¿ÁҷǶԳƵĻìºÏʱ£¬Ê¹ÓÃPOPALPHAS = 1£¨Ã¿¸öȺÌåµ¥¶ÀµÄ?£©ÄܰïÖúÒ»µã£¬µ«ÊDz»½â¾ö»ù±¾µÄÎÊÌâ¡£Òò´Ë£¬ÔÚÕâЩÇéÐÎϵĻìºÏÐԵĹÀ¼ÆÓ¦¸Ã±»½÷É÷¶Ô´ý¡£

7 ´ÓÃüÁîÐÐÔËÐÐStructure

ÓÐÐí¶àÓÉÓû§ÉèÖõijÌÐò²ÎÊý¡£ÕâЩÔÚÁ½¸öÎļþÖУ¨mainparamsºÍextraparams£©£¬Ã¿µ±³ÌÐòÖ´ÐеÄʱºò£¬ÕâЩÎļþ±»¶ÁÈ¡¡£mainparamsÖ¸¶¨Êý¾ÝÎļþµÄÊäÈëÐÎʽºÍ×î»ù±¾µÄÔËÐвÎÊý¡£extraparamsÖ¸¶¨¶àÖÖ³ÌÐòÑ¡Ïî¡£ÄãÐèÒªÉèÖÃmainparamsÀïµÄËùÓÐÖµ£¬¶øÒÔextraparamsÀïµÄĬÈÏÖµ¿ªÊ¼ÔËÐлòÐíûÓÐÎÊÌâ¡£×¢ÒâĬÈÏÄ£Ðͼٶ¨»ìºÏ£¬²¢ÇÒ²»Ê¹ÓÃÓû§¶¨ÒåµÄPopData¡£

ÔÚÕâÁ½¸öÎļþÖÐÿ¸ö²ÎÊý¶¼ÒÔ´óд×Öĸ´òÓ¡£¬ÔÚ²ÎÊýÇ°ÃæÓШD#define¡¬£¨ËüÃÇÔÚÕû¸öÕâ·Ý×ÊÁÏÖÐҲȫ²¿ÒÔ´óд×Öĸ´òÓ¡¡££© ÖµÊǽô¸úÔÚ²ÎÊýµÄÃû×ÖÖ®ºóÁ¢¼´±»È·¶¨µÄ£¨ÀýÈç¨D#define NUMREPS 1000¡¬°ÑMCMCµÄÖØ¸´ÊýÁ¿ÉèÖÃΪ1000£©¡£

¼Ìÿ¸ö²ÎÊý¶¨ÒåÖ®ºó£¬ÓÐÒ»¸ö¼ò¶ÌµÄ×¢ÊÍ£¨ÓèD//¡¬±êÃ÷£©£¬ÃèÊöÕâ¸ö²ÎÊý¡£Õâ°üÀ¨ÆÚÍûµÄÖµµÄÀàÐÍ¡£ÕâЩ°üÀ¨£º ¨D(str)¡¬£¬ÓÃÓÚ×Ö·û´®£¨ÓÃÓÚÊäÈëºÍÊä³öÎļþµÄÃû×Ö£©£»¨D(int)¡¬£¬ÓÃÓÚÕûÊý£»¨D(d)¡¬£¬ÓÃÓÚË«¾«¶ÈÊý£¨¼´£¬Ò»¸öʵÊý£¬±ÈÈç3.14£»¨D(B)¡¬£¬ÓÃÓÚ²¼¶ûÊý£¨Boolean£©£¨¼´£¬²ÎÊýȡֵΪTRUE»òFALSE£¬Í¨¹ý·Ö±ð°ÑÕâ¸öÉèÖÃΪ1»ò0£©¡£

³ÌÐò¶Ô²ÎÊýµÄ´ÎÐò²»Ãô¸Ð£¬Òò´ËÄã¿ÉÒÔÖØÅÅËüÃÇ»òÕßÔö¼Ó˵Ã÷£¬µÈµÈ¡£Ò»´Î¸ø¶¨µÄÔËÐÐËùʹÓõÄÈ«²¿²ÎÊýµÄÖµ±»´òÓ¡ÔÚÊä³öÎļþµÄĩβ¡£

7.1 ³ÌÐò²ÎÊý

ÔÚÕâÒ»½ÚÀïÎÒÃÇÁоٿÉÒÔ±»Óû§È·¶¨µÄËùÓвÎÊý¡£ÕâЩ±»¸ù¾ÝÓÃÔÚStructureµÄÃüÁîÐа汾µÄ²ÎÊýÎļþÅÅÐò¡£

7.2 ÎļþmainparamsÀïµÄ²ÎÊý

ÔÚÔËÐгÌÐò֮ǰ£¬Óû§ÐèҪȷ¶¨ËùÓÐÕâЩ²ÎÊý¡£ÕâЩ²ÎÊýµÄһЩ£¨LABEL, POPDATA,

POPFLAG, PHENOTYPE, EXTRACOLS£©±íÃ÷ÊÇ·ñÌØÊâÀàÐ͵ÄÊý¾Ý´æÔÚÓÚÊäÈëÎļþÖУ»ÕâЩÔÚµÚ2½ÚÖÐÃèÊö¡£

»ù±¾µÄ³ÌÐò²ÎÊý

MAXPOPS (ÕûÊý) Ϊ³ÌÐòµÄÒ»´ÎÌØ¶¨µÄÔËÐмÙÉèµÄȺÌåµÄÊýÄ¿¡£PritchardµÈ£¨2000a£©³ÆÖ®

ΪK¡£ÓÐʱ£¨È¡¾öÓÚÊý¾ÝµÄÐÔÖÊ£©ÓпÉÒÔ±»Ê¹ÓõÄKµÄ×ÔȻֵ£¬·ñÔòK¿ÉÒÔͨ¹ýÔÚKµÄ²»Í¬Öµ¼ì²éÄ£Ð͵ÄÅäºÏ±»¹À¼Æ£¨²ÎÔĵÚ5½Ú£©¡£

BURNIN (ÕûÊý) ÔÚ¿ªÊ¼Êý¾ÝÊÕ¼¯Ö®Ç°µÄburninʱÆÚµÄ³¤¶È¡££¨²ÎÔĵÚ3.3½Ú£© NUMREPS (ÕûÊý) ÔÚburninÖ®ºóµÄMCMCÖØ¸´µÄÊýÄ¿¡££¨²ÎÔĵÚ3.3½Ú£©

ÊäÈë/Êä³öÎļþ

INFILE (×Ö·û´®) ÊäÈëÊý¾ÝÎļþµÄÃû×Ö¡£×î´ó³¤¶È30¸ö×Ö·û£¨»òÕßÒ²Ðí¸üÉÙ£¬È¡¾öÓÚ²Ù×÷

ϵͳ£©¡£

OUTFILE (×Ö·û´®) ³ÌÐòÊä³öÎļþµÄÃû×Ö£¨ ºó׺¨D_1¡¬, ¨D_2¡¬, ...,¨D_m¡¬£¨ÓÃÓÚÖмäµÄ½á¹û£©ÒÔ

¼°¨D_f¡¬£¨×îÖÕ½á¹û£©±»Ìí¼Óµ½Õâ¸öÃû×Ö£©¡£ÓÐÕâЩÃû×ÖµÄÏÖÓеÄÎļþ½«±»¸²¸Ç¡£Ãû×ÖµÄ×î´óµÄ³¤¶ÈΪ30¸ö×Ö·û£¨»òÕßÒ²Ðí¸üÉÙ£¬È¡¾öÓÚ²Ù×÷ϵͳ£©¡£

Êý¾ÝÎļþ¸ñʽ

NUMINDS (ÕûÊý) Êý¾ÝÎļþÀïµÄ¸öÌåµÄÊýÄ¿¡£ NUMLOCI (ÕûÊý) Êý¾ÝÎļþÀïµÄλµãµÄÊýÄ¿¡£ PLOIDY (ÕûÊý) ÉúÎïÌåµÄ±¶ÐÔ¡£Ä¬ÈÏÊÇ2£¨¶þ±¶Ì壩¡£

MISSING (ÕûÊý) ΪÒÅ©»ùÒòÐÍÊý¾Ý¸øµÄ¼ÛÖµ¡£±ØÐëÊÇÒ»¸öÕûÊý£¬²¢ÇÒ²»¿ÉÔÚÊý¾Ý¼¯ÀïµÄ

ÆäËûµØ·½³öÏÖ¹ý¡£Ä¬ÈÏÊÇ-9¡£

ONEROWPERIND (²¼¶ûÖµ) ÿ¸ö¸öÌåµÄÊý¾Ý±»°²ÅÅÔÚÒ»µ¥¸öµÄÐÐÖС£ÀýÈ磬¶ÔÓÚ¶þ±¶Ìå

Êý¾ÝÀ´Ëµ£¬Õ⽫Òâζ×Åÿ¸öλµãµÄÁ½¸öµÈλ»ùÒòÔÚÏàͬµÄÐÐÖа´ÕÕÁ¬ÐøµÄ´ÎÐò£¬¶ø²»ÊDZ»°²ÅÅÔÚÏàͬµÄ×ÝÐÐÀïµÄÁ¬ÐøÁ½ÐÐÖС£¹ØÓÚÊäÈë¸ñʽµÄϸ½ÚÇë²ÎÔĵÚ2½Ú¡£ LABEL (²¼¶ûÖµ) ÊäÈëÎļþ°üº¬Ã¿¸ö¸öÌåµÄ±êÇ©£¨Ãû×Ö£©¡£1 = Êǵģ»0 = ²»¡£

POPDATA (²¼¶ûÖµ) ÊäÈëÎļþΪÿ¸ö¸öÌå°üº¬Ò»¸öÓû§¶¨ÒåµÄȺÌåÀ´Ô´¡£1 = Êǵģ»0 =

²»¡£

POPFLAG (²¼¶ûÖµ) ÊäÈëÎļþ°üº¬Ò»¸öָʾ±äÁ¿£¬Ëü±íÃ÷µ±USEPOPINFO==1£¨¼ûÏÂÃæ£©

ʱÊÇ·ñʹÓÃpopinfo¡£1 = Êǵģ»0 = ²»¡£

LOCDATA (²¼¶ûÖµ) ÊäÈëÎļþΪÿ¸ö¸öÌå°üº¬Ò»¸öÓû§¶¨ÒåµÄÈ¡ÑùλÖá£1 = Êǵģ»0 =

²»¡£ÓÃÓÚLOCPRIORÄ£ÐÍ¡£¿ÉÒÔÉèÖÃLOCISPOP = 1ÒÔ±ãÔÚLOCPRIORÄ£ÐÍÀï¸ÄΪʹÓÃPOPDATA¡£

PHENOTYPE (²¼¶ûÖµ) ÊäÈëÎļþ°üº¬Ò»¸ö±íÏÖÐÍÐÅÏ¢µÄÁС£1 = Êǵģ»0 = ²»¡£ EXTRACOLS (ÕûÊý) ÔÚ»ùÒòÐÍÊý¾Ý¿ªÊ¼Ö®Ç°£¬ÔÚ±íÏÖÐÍÖ®ºóµÄ¸½¼ÓµÄÊý¾ÝÁеÄÊýÄ¿¡£Õâ

Щ±»³ÌÐòºöÂÔ¡£0 = ûÓжîÍâµÄÁС£

MARKERNAMES (²¼¶ûÖµ) Êý¾ÝÎļþµÄ¶¥Ðаüº¬¶ÔÓ¦ÓÚʹÓõıê¼ÇµÄÒ»¸öLÃû×ÖµÄÁÐ±í¡£ RECESSIVEALLELES (²¼¶ûÖµ) Êý¾ÝÎļþµÄÏÂÒ»Ðаüº¬Ò»¸öLÕûÊýµÄÁÐ±í£¬±íÃ÷ÿ¸öλµã

ÉϵĵÈλ»ùÒòÊÇÒþÐԵġ£°ÑÕâÉèÖÃΪ1Òâζ×ÅʹÓõÄÊÇÏÔÐÔ±ê¼ÇÄ£ÐÍ¡£

MAPDISTANCES (²¼¶ûÖµ) Êý¾ÝÎļþµÄÏÂÒ»ÐУ¨»òÕßÈç¹ûMARKERNAMES==0£¬ÔòΪµÚ

Ò»ÐУ©°üº¬ÏàÁÚµÄλµãÖ®¼äµÄmapdistancesµÄÒ»¸öÁÐ±í¡£

¸ß¼¶Êý¾ÝÎļþÑ¡Ïî

PHASED (²¼¶ûÖµ) ¹©Á¬ËøÄ£ÐÍʹÓᣱíÃ÷Êý¾Ý´¦ÓÚÕýÈ·µÄÁ¬ËøÏà¡£Èç¹û£¨LINKAGE=1,

PHASED=0£©£¬ÔòPHASEINFO¿ÉÒÔ±»Ê¹Ó᪡ªÕâÊÇÊäÈëÎļþÀïµÄÒ»¸ö¶îÍâµÄÐУ¬¸ø³öÁ¬ËøÏàµÄ¸ÅÂÊ¡£µ±PHASEINFO = 0ʱÿ¸öÖµ±»ÉèÖÃΪ0.5£¬°µÊ¾Ã»ÓÐÁ¬ËøÏàµÄÐÅÏ¢¡£µ±Á¬ËøÄ£Ðͱ»ÓÃÓÚ¶à±¶Ìåʱ£¬ÒªÇóPHASED=1¡£

PHASEINFO (²¼¶ûÖµ) ÿ¸ö¸öÌåµÄ»ùÒòÐÍÊý¾ÝµÄÐкóÃæÓÐÒ»ÐйØÓÚµ¥±¶ÐÍÁ¬ËøÏàµÄÐÅÏ¢¡£

ÕâÖ»¹©Á¬ËøÄ£ÐÍʹÓá£ÓйØÏêϸÇéÐΣ¬²ÎÔĵÚ2½ÚºÍ3.1½Ú¡£

MARKOVPHASE (²¼¶ûÖµ) Á¬ËøÏàÐÅÏ¢¸úËæÒ»¸öMarkovÄ£ÐÍ¡£ÓûÁ˽âÏêϸÐÅÏ¢£¬²Î¼ûµÚ

2.2½ÚºÍ9.6½Ú¡£

NOTAMBIGUOUS (ÕûÊý) µ±RECESSIVEALLELES=1ʱ£¬¹©¶à±¶ÌåʹÓ᣶¨Òå´úÂ룬±í

Ã÷Ò»¸ö±ê¼ÇÉϵĻùÒòÐÍÊý¾ÝÊDz»º¬ºýµÄ¡£²»¿ÉÓëÊý¾ÝÀïµÄMISSING»òÕßÈκεÈλ»ùÒòÖµÏàÆ¥Åä¡£

7.3 ÎļþextraparamsÀïµÄ²ÎÊý

ÕâЩѡÏîÔÊÐíÓû§ÒÔ¸÷ÖÖ¸÷ÑùµÄ·½Ê½À´¾«Á¶Ä£ÐÍ£¬²¢ÇÒ×ö¸´ÔӵķÖÎö¡£Ä¬ÈÏÖµÊ×ÏÈ»òÐíÊǺõġ£¶ÔÓÚBooleanÑ¡Ï¼üÈë1ÓÃÓÚ¨DYes¡¬£¬»òÕߨDUse this option¡¬£»0ÓÃÓÚ¨DNo¡¬»òÕߨDDon¡¯t use this option¡¬¡£

³ÌÐòÑ¡Ïî

NOADMIX (²¼¶ûÖµ) ¼Ù¶¨Ä£ÐÍûÓлìºÏ£¨Pritchard et al., 2000a£©¡££¨Ã¿¸ö¸öÌå±»¼ÙÉèΪÍê

È«À´×ÔK¸öȺÌåÖ®Ò»£©¡£ÔÚÊä³öÖУ¬³ÌÐò²»ÊÇÈçͬÔÚ»ìºÏµÄÇé¿öÏ´òÓ¡QµÄƽ¾ùÖµ£¬¶øÊÇ´òӡÿ¸ö¸öÌåÀ´×Ôÿ¸öȺÌåµÄºóÑé¸ÅÂÊ¡£1 = ûÓлìºÏ£»0 = ÓлìºÏµÄÄ£ÐÍ¡£ LINKAGE (²¼¶ûÖµ) ʹÓÃÁ¬ËøÄ£ÐÍ¡£²ÎÔĵÚ3.1½Ú¡£RLOG10STARTÉèÖÃÿµ¥Î»¾àÀëµÄÖØ

×éÂÊrµÄ³õʼֵ¡£RLOG10MINºÍRLOG10MAXÈ·¶¨ÓÃÓÚlog10rµÄ×îСºÍ×î´óÔÊÐíÖµ¡£RLOG10PROPSDÉèÖÃÿ´Î¸üÐÂÖÐlog10rµÄ½¨ÒéµÄ¸Ä±äµÄ´óС¡£Ç°¶Ë¶ÔÕâЩ×ö³öһЩ²Â²â£¬µ«ÊÇÓû§ÐèҪעÒ⣬ȷÐÅÄÇЩֵ¶ÔÓÚÌØ¶¨µÄÓ¦ÓÃÊǺÏÀíµÄ¡£

USEPOPINFO (²¼¶ûÖµ) ʹÓÃÏÈÑéµÄȺÌåÐÅÏ¢À´·ÖÅɸöÌ嵽Ⱥ£¨cluster£©¡£Ò²¼ûMIGRPRIOR

ºÍGENSBACK¡£±ØÐëÓÐPOPDATA = 1¡£

LOCPRIOR (²¼¶ûÖµ) ʹÓõصãÐÅÏ¢À´Ìá¸ß½á¹¹ÐÅÏ¢ÈõµÄÊý¾ÝµÄÐÔÄÜ¡£

FREQSCORR (Ë«¾«¶È) ʹÓèDF model¡¬£¬ÆäÖеÈλ»ùÒòƵÂÊ¿çԽȺÌåÏà¹Ø (Falush et al.,

2003a)¡£¸ü¾ßÌåµØËµ£¬²»ÊǼٶ¨Ò»¸öÏÈÑ飬ÔÚÆäÖÐÿ¸öȺÌåÖеĵÈλ»ùÒòƵÂÊÊÇÀ´×ÔÒ»¸ö¾ùÔȵÄDirichlet·Ö²¼µÄ¶ÀÁ¢³éÑù£¬¶øÊÇ´ÓÒ»¸ö·Ö²¼¿ªÊ¼£¬Õâ¸ö·Ö²¼ÒÔÑù±¾µÄƽ¾ùµÈλ»ùÒòƵÂÊΪÖÐÐÄ¡£Õâ¸öÄ£ÐͶÔÓÚÇ×Ôµ¹ØÏµ·Ç³£½üµÄȺÌå¸üÕæÊµ£¨ÆäÖÐÎÒÃÇÔ¤ÆÚµÈλ»ùÒò ƵÂÊ¿çԽȺÌåÏàËÆ£©£¬²¢ÇÒ¿ÉÒÔ²úÉú¸üºÃµÄ¾ÛÀࣨµÚ3.2½Ú£©¡£FkµÄÏÈÑéÊÇÓÃFPRIORMEAN ºÍFPRIORSDÈ·¶¨µÄ¡£µ±FREQSCORR±»¿ªÆôʱ£¬¿ÉÄÜÓÐÒ»ÖÖ¹ý¸ß¹À¼ÆKµÄÇ÷ÊÆ¡£ ONEFST (²¼¶ûÖµ) ¶ÔÈ«²¿ÈºÌå¼Ù¶¨ÏàͬµÄFkµÄÖµ£¨ÀàËÆÓÚWrightµÄ´«Í³µÄFST£©¡£ÕⲻΪ´ó

¶àÊýÊý¾ÝÍÆ¼ö£¬ÒòΪÔÚʵ¼ùÖÐÄã»òÐíÆÚÍûÿ¸öȺÌåÖеIJ»Í¬µÄ·ÖÆç£¨divergence£©Ë®Æ½¡£µ±K = 2ʱ£¬ÓÐʱ¿ÉÄÜÄÑÒÔ·Ö±ð¹À¼ÆFSTµÄ¶þ¸öÖµ£¨µ«ÊǼûHarterµÈ£¨2004£©£©¡£µ±ÄãÊÔͼ¹À¼ÆKʱ£¬ÄãÓ¦¸Ã°ÑÏàͬµÄÄ£ÐÍÓÃÓÚÈ«²¿K£¨ÎÒÃǽ¨ÒéONEFST = 0£©¡£

INFERALPHA (²¼¶ûÖµ) ¸ù¾ÝÊý¾ÝÍÆ¶ÏÄ£ÐͲÎÊý?µÄÖµ£»·ñÔò?±»¹Ì¶¨ÔÚÖµALPHAÉÏ£¬Ëü

ÊÇÓÉÓû§Ñ¡ÔñµÄ¡£ÔÚNOADMIXÄ£ÐÍÏÂÕâ¸öÑ¡Ïî±»ºöÂÔ¡££¨×æÏÈÏòÁ¿QµÄÏÈÑéÊǾßÓвÎ

ÁªÏµ¿Í·þ£º779662525#qq.com(#Ìæ»»Îª@)