structure 2.3 - ÖÐÎÄ ÏÂÔØ±¾ÎÄ

Êý(?, ?, ..., ?)µÄDirichlet¡£Ð¡µÄ?°µÊ¾´ó¶àÊý¸öÌå»ù±¾ÉÏÀ´×ÔÒ»¸öȺÌå»òÕßÁíÒ»¸öȺÌ壬¶ø?>1 °µÊ¾´ó¶àÊý¸öÌå±»»ìºÍ¡££©

POPALPHAS (²¼¶ûÖµ) Ϊÿ¸öȺÌåÍÆ¶ÏÒ»¸öµ¥¶ÀµÄ?¡£¶àÊýÇé¿öϲ»ÍƼö£¬µ«ÊǶÔÓÚ¾ßÓÐ

·Ç¶Ô³Æ»ìºÏµÄÇé¿ö¿ÉÄÜÊÇÓÐÓõġ£

ALPHA (Ë«¾«¶È) »ìºÏ³Ì¶ÈµÄDirichlet²ÎÊý£¨?£©£¨Èç¹ûINFERALPHA ==1£¬Õâ¸ö¾ÍÊdzõʼ

Öµ£©¡£

INFERLAMBDA (²¼¶ûÖµ) Ϊ?ÍÆ¶ÏÒ»¸öÊʵ±µÄÖµ¡£²»Îª´ó¶àÊý·ÖÎöÍÆ¼ö¡£ POPSPECIFICLAMBDA (²¼¶ûÖµ) Ϊÿ¸öȺÌåÍÆ¶ÏÒ»¸öµ¥¶ÀµÄ?¡£

LAMBDA (Ë«¾«¶È) ²ÎÊý»¯µÈλ»ùÒòƵÂʵÄÏÈÑ飬¶ÔÓÚ´ó¶àÊýÊý¾ÝĬÈÏÖµ1ºÃÏñ¹¤×÷µÃÏ൱

ºÃ¡£Èç¹û´ó¶àÊý±ê¼ÇÉÏµÄÆµÂÊÑÏÖØÏòµÍ/¸ßƵÂÊÆ«Ð±£¬¸üСµÄ?Öµ¿ÉÄܵ¼ÖÂDZÔÚ¸üºÃµÄ±íÏÖ¡£ËüºÃÏñ¶Ô¹À¼Æ?ûÓÐ¶ÔÆäËü³¬²ÎÊý£¨?ºÍF£©µÄ¹À¼ÆÄÇôºÃ¡£

ÏÈÑé

ÕâЩֵÓÃÀ´²ÎÊý»¯¼Ù¶¨µÄ¸ÅÂÊÄ£ÐÍ¡£¶àÊýÇé¿öÏÂȱʡÉèÖÃÓ¦¸ÃÏ൱Ã÷ÖÇ£¬Äã¿ÉÄܲ»ÏëÒªµ£ÐÄÕâЩ¡£

FPRIORMEAN£¬FPRIORSD (Ë«¾«¶È) ¼ûFREQSCORR¡£FkµÄÏÈÑ鱻ȡ¾ßÓÐÆ½¾ùÊý

FPRIORMEANºÍ±ê׼ƫ²îFPRIORSDµÄ¦Ã¡£ÎÒÃǵÄȱʡÉèÖý«Ðí¶àÈ¨ÖØ·ÅÓÚFµÄСֵÉÏ¡£ÎÒÃÇ·¢ÏÖÕâʹµÃËã·¨¶Ôϸ΢µÄ½á¹¹Ãô¸Ð£¬µ«ÊÇÔö¼ÓÁ˹ý¸ß¹À¼ÆKµÄ·çÏÕ£¨Falush et al., 2003a£©¡£

UNIFPRIORALPHA (²¼¶ûÖµ) £¬ALPHAMAX (Ë«¾«¶È) ¶Ô?¼ÙÉèÒ»¸ö¾ùÔÈ·Ö²¼µÄÏÈÑ飬ÔË

ÐÐÔÚ0ºÍALPHAMAXÖ®¼ä¡£Õâ¸öÄ£ÐͺÃÏñ¹¤×÷ºÜºÃ£»ÆäËüµÄÄ£ÐÍ£¨µ±UNIFPRIORALPHA = 0ʱ£©ÊÇÈ¡?×÷Ϊ¾ßÓЦ÷ֲ¼µÄÏÈÑ飬¾ßÓÐÆ½¾ùÊýALPHAPRIORA¡ÁALPHAPRIORB£¬·½²îALPHAPRIORA¡ÁALPHAPRIORB2¡£

LOG10RMIN£¬LOG10RMAX£¬LOG10PROPSD£¬LOG10RSTART (Ë«¾«¶È) µ±Á¬ËøÄ£ÐÍ

±»Ê¹ÓÃʱ£¬Çл»ËÙÂÊr±»È¡ÎªÔÚ¶ÔÊý³ß¶ÈÉϾßÓÐÒ»¸ö¾ùÔȵÄÏÈÑ飬ÔÚLOG10RMINºÍLOG10RMAXÖ®¼ä¡£ÕâЩֵÐèÒªÓÉÓû§È·¶¨£¬ÒÔ±ã¾Í±»Ê¹ÓõÄͼ¾àµ¥Î»µÄ³ß¶È¶øÑÔÓÐÒâÒå¡£

ʹÓÃÏÈÑéµÄȺÌåÐÅÏ¢£¨USEPOPINFO£©

GENSBACK (ÕûÊý) Õâ¶ÔÓ¦ÓÚG£¨Pritchard et al., 2000a£© ¡£µ±¶Ô¸öÌåʹÓÃÏÈÑéµÄȺÌåÐÅÏ¢

£¨USEPOPINFO = 1£©Ê±£¬³ÌÐò²âÊÔÊÇ·ñÿ¸öÌåÔÚ×îºóG´úÓÐÒ»¸öÒÆÃñ׿ÏÈ£¬ÆäÖÐG = 0¶ÔÓ¦ÓÚ¸öÌå±¾Éí×÷Ϊһ¸öÒÆÃñ¡£ÎªÁËÓÐÊʵ±µÄ¹¦Ð§£¬GÓ¦¸Ã±»ÉèÖõÃÏ൱С£¨±ÈÈç2£©£¬³ý·ÇÊý¾ÝµÄÐÅÏ¢·Ç³£·á¸»¡£

MIGRPRIOR (Ë«¾«¶È) ±ØÐëÊÇÔÚ[0£¬1]ÖС£ÕâÊÇPritchardµÈ(2000a)ÖеÄv¡£Ã÷ÖǵÄÖµ¿ÉÄÜÔÚ

0.001¡ª0.1µÄ·¶Î§ÄÚ¡£

PFROMPOPFLAGONLY (²¼¶ûÖµ) Õâ¸öÑ¡ÏîÊÇ2.0°æÐÂÓеģ¬Ê¹½ö½öʹÓÃÔ¤ÏÈÖ¸¶¨µÄ¸öÌå

µÄ×Ó¼¯À´¸üеÈλ»ùÒòƵÂÊP³ÉΪ¿ÉÄÜ¡£ÎªÁËʹÓÃÕâ¸öÑ¡Ï°üÀ¨Ò»¸öPOPFLAGÁУ¬¶ÔÓ¦¸ÃÓÃÀ´¸üÐÂPµÄ¸öÌåÉèÖÃPOPFLAG=1£¬¶Ô²»Ó¦¸ÃÓÃÀ´¸üÐÂPµÄ¸öÌåµÄÉèÖÃPOPFLAG = 0¡£USEPOPINFO´ò¿ª»òÕß²»´ò¿ª£¬¶¼¿ÉÒÔʹÓÃÕâ¸öÑ¡Ïî¡£

Õâ¸öÑ¡ÏÊÇÓÐÓõģ¬ÀýÈ磬Èç¹ûÄãÓÐÀ´×ÔÒÑ֪ȺÌåµÄÒ»×é×¼²Î¿¼¸öÌ壬ȻºóÄãÏëÒª¹À¼ÆÒ»Ð©Î´ÖªµÄ¸öÌåµÄ׿ÏÈ¡£Ê¹ÓÃÕâ¸öÑ¡Ïÿ¸öδ֪¸öÌåµÄq¹À¼ÆÖµÖ»ÒÀÀµÓڲο¼¸öÌå×飬¶ø²»ÒÀÀµÓÚÑù±¾ÖÐµÄÆäËüδ֪¸öÌå¡£Õâ¸öÌØµãÓÐʱÊǺϺõÐèÒªµÄ¡£

ʹÓõصãÐÅÏ¢µÄLOCPRIORÄ£ÐÍ

LOCISPOP (²¼¶ûÖµ) µ±LOCPRIORÄ£ÐÍ´ò¿ªÊ±£¬Õâ¸öÑ¡ÏîÖ¸µ¼³ÌÐòʹÓÃÊäÈëÎļþÖеÄ

PopDataÁÐ×÷ΪµØµãÊý¾Ý¡£µ±LOCISPOP = 0ʱ£¬³ÌÐòÐèÒªÒ»¸öLocDataÁÐÀ´Ê¹ÓÃLOCPRIOR¡£

LOCPRIORINIT (Ë«¾«¶È) LOCPRIOR²ÎÊýrµÄ³õʼֵ£¬È·¶¨ÈºÌåµÄÐÅÏ¢Óжà·á¸»

£¨citepHubiszEtAl09£©¡£ÎÒÃÇ·¢ÏÖLOCPRIORINIT = 1ÓÐÖúÓÚÈ¡µÃºÃµÄÊÕÁ²¡£

MAXLOCPRIOR (Ë«¾«¶È) rµÄ·¶Î§ÊÇÀ´×Ô£¨0£¬MAXLOCPRIOR£©¡£ÎÒÃǽ¨Òé

MAXLOCPRIOR = 20¡£

Êä³öÑ¡Ïî

PRINTNET (²¼¶ûÖµ) ´òÓ¡¾ÛÀàÖ®¼äµÄ¡°¾»ºËÜÕËá¾àÀ롱¡£ÈºÌåAºÍBÖ®¼äµÄÕâ¸ö¾àÀëDAB

±»¼ÆËãΪ

DA,B1L?Jl(l)(l)?(HA?HB)?A,jp?B,j??£¬ £¨5£© ???1??pLl?1?j?12?(l)?xÆäÖÐp,jÊÇȺÌåxÖÐλµãlÉϵĵÈλ»ùÒòjµÄµÈλ»ùÒòƵÂʵĺóÑ鯽¾ù¹À¼Æ£¬LÊÇλµãµÄÊýÄ¿£¬

JlÊÇλµãlÉϵĵÈλ»ùÒòµÄÊýÄ¿£¬ÆäÖÐ

1L?Jl(l)2??x,j?£¬ £¨6£© Hx???1??pLl?1?j?1?¾»ºËÜÕËá¾àÀëÊÇ·Ö±ðÀ´×ÔȺÌåAºÍBµÄÒ»¶ÔµÈλ»ùÒò²»Í¬µÄƽ¾ù¸ÅÂÊ£¬Ð¡ÓÚȺÌåÄÚÆ½¾ùÔÓºÏÐÔ¡£»òÐí¸üÖ±¾õµØ£¬Õâ¿ÉÒÔ±»ÈÏΪÊÇÀ´×Ô²»Í¬ÈºÌåµÄµÈλ»ùÒòÖ®¼ä³É¶Ô²îÒìµÄƽ¾ùÊýÁ¿£¬³¬³öÿ¸öȺÌåÄÚ´æÔڵıäÒìµÄÊýÁ¿¡£¸Ã¾àÀëÓкÏÊʵÄÌØÐÔ£¬¼´ÏàËÆµÄȺÌåÖ®¼äµÄ¾àÀë½Ó½ü0£¬²¢ÇÒÓÈÆä£¬DAA = 0¡£×¢Òâµ½¾àÀëÊǶԳƵģ¬ÒÔÖÁÓÚDAB = DBA¡£Õâ¸ö¾àÀëÊÊÓÚ»­ÈºÌåµÄÊ÷À´°ïÖúÖ±¹Û»¯¾ÛÀàÖ®¼äµÄ²îÒìµÄˮƽ£¨Falush et al., 2003b£©¡£ PRINTKLD (²¼¶ûÖµ) [²»ÔÞ³É]Õâ¸öÑ¡Ïî²»ÔÙ¿ÉÓᣠPRINTLAMBDA (²¼¶ûÖµ) ´òÓ¡?µÄµ±Ç°µ½ÆÁÄ»¡£

PRINTQSUM (²¼¶ûÖµ) ´òÓ¡µ±Ç°µÄQ¹À¼ÆµÄÕªÒªµ½ÆÁÄ»£»Õâ´òÓ¡PopDataµÄÿ¸öÖµµÄƽ¾ùÊý¡£

SITEBYSITE (²¼¶ûÖµ) £¨Á¬ËøÄ£ÐÍ£©ÎªÊý¾ÝÀïµÄÿÖÖ»ùÒòÐÍ´òÓ¡·ÖÅɸÅÂʵÄÍêÕûµÄÕªÒª¡£Õâ±»´òÓ¡µ½Ò»¸öµ¥¶ÀµÄÎļþ£¬¾ßÓкó׺¨Dss¡¬¡£Õâ¸öÎļþ¿ÉÄÜÊÇ´óµÄ£¡

PRINTQHAT (²¼¶ûÖµ) µ±Õâ¸ö±»´ò¿ªÊ±£¬QµÄµã¹À¼Æ²»½ö±»´òÓ¡µ½Ö÷ÒªµÄ½á¹ûÎļþ£¬¶øÇÒ±»´òÓ¡µ½Ò»¸öµ¥¶ÀµÄºó׺Ϊ¨Dq¡¬µÄÎļþ¡£ÎªÁËÔËÐÐæ¢ÃóÌÐòSTRAT£¬ÐèÒªÕâ¸öÎļþ¡£ UPDATEFREQ (ÕûÊý) ¶ÔÆÁÄ»´òÓ¡¸üÐÂµÄÆµÂÊ¡£Èç¹ûÕâ= 0Ôò×Ô¶¯ÉèÖᣠPRINTLIKES (²¼¶ûÖµ) ÔÚÿ¸öµü´ú¹ý³ÌÖаÑËÆÈ»º¯ÊýµÄµ±Ç°Öµ´òÓ¡µ½ÆÁÄ»¡£

INTERMEDSAVE (ÕûÊý) Èç¹ûÄã¿ÊÍûÔÚÔËÐнáÊøÖ®Ç°¿´¼û³õ²½µÄ½á¹û£¬Äã¿ÉÒÔÈóÌÐòÔÚMCMCÔËÐÐÆÚ¼äÿ¸ôÒ»¶Îʱ¼ä½«½á¹û´òÓ¡µ½Îļþ¡£×ܹ²´òÓ¡INTERMEDSAVE¸öÕâÑùµÄÎļþ£¬ÔÚ¼ÌBURNINµÄÍê³ÉÖ®ºóÏàµÈµÄ¼ä¸ô¡£Í¨¹ýÉèÖÃΪ0¹ØµôÕâ¸öÑ¡ÏʹÓÃOUTFILEÃû×Ö´´ÔìµÄÕâЩÎļþµÄÃû×Ö¡£

ECHODATA (²¼¶ûÖµ) °ÑÊý¾Ý¼¯µÄ¸ÅÒª´òÓ¡µ½ÆÁÄ»ºÍÊä³öÎļþ¡££¨´òÓ¡ÊäÈëÎļþ×î¶¥ÐкÍ×îµ×ÐеĿªÍ·ºÍ½á⣬ÒÔ±ãÔÊÐíÓû§¼ì²éËüÒѾ­±»ÕýÈ·µØ¶ÁÈ¡¡££©

ANCESTDIST (²¼¶ûÖµ) ÊÕ¼¯¹ØÓÚÿ¸öÌåµÄQµÄ·Ö²¼µÄÐÅÏ¢£¬Ò²Ö»¹À¼ÆÆ½¾ùÖµ¡£µ±Õâ±»´ò¿ªÊ±£¬Êä³öÎļþ°üÀ¨Ã¿¸öq(i)µÄ¸ÅÂÊÇø¼äµÄ×ó¶ËºÍÓÒ¶Ë¡££¨¸ÅÂÊÇø¼äÊÇÖÃÐÅÇø¼äµÄ±´Ò¶Ë¹ÀàËÆ

Îï¡££©´òÓ¡µÄÖµÏÔʾ¸ÅÂÊÇø¼äµÄÖмä100p%£¬ÆäÖÐpÊÇ0.0µ½1.0µÄ·¶Î§ÄÚµÄÒ»¸öÊý×Ö£¬ÓÃANCESTPINTÉèÖá£QµÄ·Ö²¼Í¨¹ý¼Ç¼ÔÚ0ºÍ1Ö®¼äµÄÐí¶àÏä×ÓµÄÿһ¸öµã»÷£¨hit£©µÄÊýÄ¿±»¹À¼Æ£¬ÒÔ±ãÐγÉijÖÖÖ±·½Í¼¡£ÕâЩÏä×ӵĿí¶ÈÒ»Ñù´óС£¬ÊÇʹÓÃNUMBOXESÈ·¶¨µÄ¡£

ÔÓÏî

COMPUTEPROB (²¼¶ûÖµ) ÔÚÿ´Î¸üÐÂʱ´òÓ¡Êý¾ÝµÄ¶ÔÊýËÆÈ»£¬²¢¹À¼Æ¸ø¶¨KºÍÄ£Ð͵ÄÊý

¾ÝµÄ¸ÅÂÊ£¨²ÎÔĵÚ5½Ú£©¡£Õâ±»ÓÃÓÚ¹À¼ÆK£¬²¢ÇÒ¶ÔÓÚburnin£¨ÀÏ»¯£¿£©ÊÇ·ñ×ã¹»³¤Ò²ÊÇÒ»¸öÓÐÓõÄÕï¶Ï¡£¹ØµôÕâ¸öÑ¡ÏîµÄÖ÷ÒªÔ­Òò½«ÊǼÓËÙ³ÌÐò£¨~ 10¨C15%£©¡£

ADMBURNIN (ÕûÊý) £¨µ±RECOMBINE=1ʱʹÓ㩵±Ê¹ÓÃÁ¬ËøÄ£ÐÍʱ£¬ÔÚ´ó¶àÊýÇéÐÎÏÂ

Ç¿ÁÒÍÆ¼öÀûÓûìºÏÎïÄ£Ð͵Ķ̵Äburnin£¨±È·½Ëµ500´Îµü´ú£©¡£Ã»ÓÐÕâÑùµÄburnin£¬Á¬ËøÄ£Ð;­³£Éú²ú¶ÀÌØµÄ½á¹û¡£ÉèÖÃADMBURNIN

ALPHAPROPSD (Ë«¾«¶È) ?µÄMetropolis-Hastings¸üв½ÖèÉæ¼°´ÓÒ»¸ö¾ßÓÐÆ½¾ùÊý#ºÍ±ê

×¼²îALPHAPROPSD>0µÄÕý̬·Ö²¼ÖÐÌôѡһ¸öÖµ??¡£ALPHAPROPSDµÄÖµ²»Ó°ÏìMarkovÁ´µÄ½¥½üµÄÐÐΪ£¬µ«ÊÇ¿ÉÄܶÔÊÕÁ²µÄËÙÂʾßÓÐÏÔÖøµÄÓ°Ïì¡£Èç¹ûÓÐÐí¶à¹ØÓÚ?µÄÐÅÏ¢£¬ÔòСµÄALPHAPROPSDÖµ¶ÔÓÚ»ñµÃºÏÀíµÄ½ÓÊÜËÙÂʸü¿ÉÈ¡¡£Èç¹ûûÓÐºÜ¶à¹ØÓÚ#µÄÐÅÏ¢£¬Ôò¸ü´óµÄÖµ²úÉú¸üºÃµÄ»ìºÏ¡£

STARTATPOPINFO (²¼¶ûÖµ) ʹÓøø¶¨µÄȺÌå×÷ΪȺÌåÆðÔ´µÄÆðʼÌõ¼þ¡££¨ÐèÒª

POPDATA==1£©¡£Õâ¸öÑ¡ÏîÌṩMarkovÁ´ÕýÔÚÕýÈ·µØÊÕÁ²µÄÒ»¸ö¼ì²é£¬Èç¹ûÄãÆÚÍûÍÆ¶ÏµÄ½á¹¹ÓëÊäÈë±êÇ©ÏàÅ䣬¶øËüûÓС£Õâ¸öÑ¡Ïî¼Ù¶¨ÊäÈëÎļþÀïµÄPopDataÔÚ1ºÍkÖ®¼ä£¬ÆäÖÐk ? MAXPOPS¡£PopData²»ÔÚÕâ¸ö·¶Î§ÄڵĸöÌåËæ»úÆðʼ¡£

RANDOMIZE (²¼¶ûÖµ) Ϊÿ´ÎÔËÐÐʹÓÃÒ»¸ö²»Í¬µÄËæ»úÊýÖÖ×Ó£¬È¡×ÔϵͳʱÖÓ¡££¨Ò²¼û

SEED£©

SEED£¨ÕûÊý£© Èç¹ûRANDOMIZE==0£¬ÔòÄ£ÄâÖÖ×Ó±»Ô¤ÖÃΪSEED¡£ÕâÔÊÐíÔËÐб»×¼È·µØ

ÖØ¸´¡£Èç¹ûRANDOMIZE==1£¬ÔòÔÚSEEDÀïÖ¸¶¨µÄÈκÎÖµ±»ºöÂÔ¡£×¢Ò⼴ʹµ±RANDOMIZE==1ʱ£¬×¢ÒâÄǸö==1£¬³ÌÐòÊä³öÈÔÈ»±íÃ÷¿ªÊ¼µÄÖÖ×ÓÖµ£¬ÕâÑùµÄ»°£¬Èç¹ûÐèÒª£¬¿ÉÒÔÖØ¸´Ìض¨µÄÔËÐС£

METROFREQ (ÕûÊý) ÔÚ»ìºÏÎïÄ£ÐÍÏÂʹÓÃMetropolis-Hastings²½ÖèÀ´¸üÐÂQµÄƵÂÊ¡£µ±Õâ

¸öÑ¡ÏʹÓÃʱ£¬ÎªÃ¿¸öq(i)Ñ¡ÔñÒ»¸öеÄÌáÒéÖµq(i)?¡£Õâ¸öÌáÒéÖµÊÇ´ÓÏÈÑ飨¼´q(i)? ~ D(?,