060310391&ȃʼngŘȄŶȀ&0560565&ȃégŘȄ&&&
ǦDžNJǴǡǺǝǎǖ&ă6S&
ŏŝWU��&
2014/5/22&&&8:40-10:10&@1LŰ-ă4Łěk�
&
Ÿ�
Ŗ�eŏŝWU&
• Ŗ�DZȁNjȁƶŖ�eƩÂģ��Ƥ�ęǀúƖƝWU&• öôƖƣƊƼŖ�eƩ�ęƶƄƛƩPƻƨƉƼŖ�DZȁNjȁǀùļƤƐƼ&ȃ ƌƫƄïl ��ƩŖ�eǀŊ*ƩŖ�DZȁNjȁƤŕ�ƘƼȄ&
• Ŗ�eŏŝWUƨŚęƔƽƝŖ�DZȁNjȁƩaYƥƄŜó�Ńȃ ƌƫƄHŜƶĤ�ȄƩaYƥƩŠ#ǀ´ƺƎƨƘƼƓƥƤƄŜó�ŃŖ�e�ȃQTLȄƩ�ęƶ@Áǀ§hƤƐƼ&ȃă7SƄă8SƤgĜƖƲƘȄ�
Ź�
¿³ƩǍȁǾȁǢ�
• ŏŝWUƥãéWU&• �ªƌƥĐªƌ&• �ũ5{Ƅaũ5{&• q�Ƅºq×&• Đªƌ"&• n®NJǝǗ&• �5×&• WUŠ®&• bZIJƄs�IJ&• ĐNƚºŔBRŬƄuSǘȁǻǖDZǿRŬ&• Ě÷ĞƄĦǐǽǣȁºŔB&
ź�
ŏŝWU&vs.&ãéWU&
• ãéWU&– DNA^[n®ƧƦƄãéóƧņŧƨ[ƢƑWU&
• ŏŝWU&– Ŗ�eşƩŏŝƩ�ƔȃĐªƌ"ȄǀņŧƥƖƝWU&
• �ĵ2þŭƩņŧƪ�ƊƨÓ ƖƧƊ&– Ŗ�eşƩĐªƌƩū�ǀ�«ƘƼƨƪƄŏŝWUƩ�\Əśĭ&
Ż�
DNAƪƁÂģ�ƂƨǧǝǎƔƽƣƊƼ�
DžǤƩ12Đȃ24¿ȄƩÂģ��
ǨǡƩ]NƄMĐƩÂģ�ƩƋƟƄ1¿ǀÑİƎƺ&ƵƋ1¿ǀâİƎƺƷƙƻIƒƼ�
ǨǡƩȉȊĐȃ46¿ȄƩÂģ��
h:p://wkp.fresheye.com/wikipedia/Âģ��DNAƩ3ė�
Kurata&and&Omura&(1982)&Jpn&J&Breed&32:253&�
ż�
Âģ�Ʃ�ªƌƥŖ�eƩĐªƌ�
ѲƎƺƩÂģ��
â²ƎƺƩÂģ��
�ªƌ�
�6���D ��=0);���E�5"!�54/1$;���!2��%'�+;�
A� B� C�
���7 ��=�-1�&:�5�*�,<;#!
+63("�$5�)5��.;!���D�%8"B3CE7�5�5�%:<9.$! Ž�
DžǤă8Âģ�&
AF��Ŋƪ�ªƌƏ�:ƔƽƣƊƼƝƴ&ãéņŧȃËňȄƉƝƻƩĐªƌ"ȃĖňȄƏ�Ɗ&
ž�
0 2 4 6 8
0.00
0.05
0.10
0.15
0.20
0.25
0.30
k
P(X=k)
òèƏ3ƄńèƏ7.ƠƝĆƎƺƄ&8SƩ�,¡4ǀħƊƄòèƄńèǀƛƽƜƽ4SƙơƬƑùç�
€
P(X = 4) =84"
# $ %
& ' × 0.34 × 0.74
�ĢƨƄĆƩƩòèƩ=NƏpȃńèƩ=NƏ1-pȄƩƥƐƄ&nSƩ�,¡4ƤƄkSòèǀƬƑùç&
€
nk"
# $ %
& ' =
n!k!(n − k)!
ƓƓƤƄ�
€
P(X = k) =nk"
# $ %
& ' pk (1− p)n−k
�ũ5{�
�ũ#®�
p=0.3,&n=8�
ſ�
Ȏ�ĆƩƩòèƩ=NƏpƩƥƐƨƄ8SƩ�,¡4ƤƄòèƄńèǀƛƽƜƽ4Sƙơ�Ƒùç�
€
P(X = 4 | p) =84"
# $ %
& ' p4 (1− p)4
�Ƅ¾÷ƩpǀÔƴƝƊƥƘƼƅ&ƓƩùçǀºbƨƘƼpǀIJƥƖƣƪƦƋƎȎ&
ºq×�
q�ȋƉƼ�ľȃǵǠǻȄƩƵƥƤƄıÝƔƽƝ łƏêƗƼùç&ȃƨÓ Ȅ�
q�ǀºbƨƘƼǵǠǻǧǹǴȁǚǀÔƴƼ²×�
€
L(p | X = 4)∝ P(X = 4 | p) =84#
$ % &
' ( p4 (1− p)4 →&ºbB�
ƀ�
q�ƩºbB�
• q�ȃL(p)ȄƩºbBƪƄq�ƩĠán®ȃln&L(p)ȄƩºbBƘƼƓƥƥOƗŶȃln&L(p)ƏEŀ_?Š®ƩƝƴȄ�
€
lnL(p) = 4 ln p + 4 ln(1− p) + const
0.0 0.2 0.4 0.6 0.8 1.0
-18
-14
-10
-6
p
4 ln
(p) +
4 ln
(1-p
)
p=0.5ƩƥƐºb�
d lnL(p)dp
= 4 (1− 2p)p(1− p)
= 0
0.0 0.2 0.4 0.6 0.8 1.0
-4-3
-2-1
0
p
log(p)
Ÿŷ�
òè3Ƅńè5ƄŨè2Ə.ƠƝĆƎƺƄ&9SƩ�,¡4ƤƄòèƄńèƄŨèǀƛƽƜƽ3SƙơƬƑùç�
€
P(X1 = 3,X2 = 3,X3 = 3) =9!3!3!3!
× 0.33 × 0.53 × 0.23
�ĢƨƄĆƩƩiíôȃi=1,2,…,kȄƩģƩèƩ=NƏpiƩƥƐƄnSƩ�,¡4ƤƄiíôƩģƩèǀxiSƬƑùç&
€
P(X1 = x1,...,Xk = xk ) =n!
x1!x2!p1x1 pk
xk =n!
x1!x2!pixi
i
k
∏
ƔƐƱƦƩ5{Ə�ũ5{ƞƠƝƩƨnƖƄƓƽƪaũ5{�
1� 1� 1� 2� 2�
2� 2� 2� 3� 3�
ŸŸ�
aũ5{Ʃn®q��
€
lnL(p1, p2) = 3ln p1 + 3ln p2 + 3ln(1− p1 − p2) + const
p1=1/3,&p2=1/3,&p3=1/3ƩƥƐºb�
€
∂ lnL(p1, p2)∂p1
= 3 1p1−
11− p1 − p2
$
% &
'
( ) =
(1− 2p1 − p2)p1(1− p1 − p2)
= 0
€
∂ lnL(p1, p2)∂p2
= 3 1p2−
11− p1 − p2
$
% &
'
( ) =
(1− p1 − 2p2)p2(1− p1 − p2)
= 0
òèƄńèƄŨèƏƛƽƜƽp1,p2,p3ƤÙƕƠƣƊƼĆƎƺƄ9SƩ�,¡4ǀħƊƄòèƄńèƄŨèǀƛƽƜƽ3Sƙơ�ƑƥƊƋĒÁƏ�ƺƽƝƅƓƩƥƐƩp1,p2,p3Ʃºq§hŜǀÔƴƹƅ�
Ȏ�
ŸŹ�
&Ü®5ĪµƨƍƒƼ&Âģ�Ʃ�ªƌƥŖ�eƩĐªƌ�
Ÿź�
F2ƨƍƒƼ5ŧ�
Ŧþă1�ȃF1Ȅ&
İ1ȃP1Ȅ& İ2ȃP2Ȅ&×&�Ś&
Ś)�&
Ü®5Ī&
İ1Y&
A&
B&
A&
B&
A&
B&
a&
b&
a&
b&
a&
b&
A&
B&
0.5(1-r)&
a&
b&
0.5(1-r)&r:&Đªƌ"&
İ2Y&
ĠÐ&
F2ťTƨƍƒƼ5ŧȃÌǖǹDžǢȄ&
ĐªƌY& ĐªƌY&
a&
B&
A&
b&
0.5r& 0.5r&
ŸŻ�
F2ƨƍƒƼ5ŧȃĕƐȄ�AB&0.5(1-r)&
Ab&0.5r&
aB&0.5r&
ab&0.5(1-r)&
AB&0.5(1-r)&
AABB&(1-r)2&&
AABb&r(1-r)&
AaBB&r(1-r)&
AaBb&(1-r)2&
Ab&0.5r&
AABb&r(1-r)&
AAbb&r2&
AaBb&r2&
Aabb&r(1-r)&
aB&0.5r&
AaBB&r(1-r)&
AaBb&r2&
aaBB&r2&
aaBb&r(1-r)&
ab&0.5(1-r)&
AaBb&(1-r)2&
Aabb&r(1-r)&
aaBb&r(1-r)&
aabb&(1-r)2&
1/4ƪõƊƣƉƼ&
Ÿż�
AA&&
Aa& aa&
BB& (1-r)2& 2r(1-r)&&
r2&&
Bb& 2r(1-r)& 2{r2+(1-r)2}& 2r(1-r)&&
bb& r2& 2r(1-r)&&
(1-r)2&
½�ū�ƨƪ1/4ǀƎƒƼ�ĭƏƉƼ&
!!!!!!!!!!!!!!!!!30!
!!!!!!!!!!!!!!!!!!7!
!!!!!!!!!!!!!!!!2!
!!!!!!!!!!!!!!!!!!!5!
!!!!!!!!!!!!!!!!57!
!!!!!!!!!!!!!!!!6!
!!!!!!!!!!!!!!!!!!!1!
!!!!!!!!!!!!!!!!!!8!
!!!!!!!!!!!!!!!29!
q�Š®ƪaũ5{Ʃq�Š®&ȃ9þŭƩģèǀƄ145SƩ¡4ǀħƠƝ]NƥOƗȄ&
€
lnL(r) = 30ln(1− r)2 + 7ln2r(1− r) ++ 29ln(1− r)2 + const
€
=144 ln(1− r) + 32ln r + 57ln(1− 2r + 2r2) + const
ŸŽ�
€
lnL(r) =144 ln(1− r) + 32ln r + 57ln(1− 2r + 2r2) + const
0.0 0.1 0.2 0.3 0.4 0.5
-170
-150
-130
-110
r
ln L
(r)
- const
d lnL(r)dr
= −1441− r
+32r+57(4r − 2)1− 2r + 2r2
0.0 0.1 0.2 0.3 0.4 0.5
-2000
200
600
1000
r
LOD&(ȋƃÓŶŶŶŶŶŶŶŶŶŶŶŶŶŶŶŶȃNJǝǗoddsƥƹƮȄƩ}ën®ǀƥƠƝƵƩƅ&&&&&&&&&&ŏŝƖƣƊƧƊƥƐƩq�L(0.5)ƨÓưƄ§hƔƽƝĐªƌ"Ʃ&ƃƃƃƵƥƤƩq�Ə�&ƉƼƎǀĩƖƝ(&
€
L( ˆ r ) /L(0.5)€
ˆ r = 0.117
LOD&
€
LOD = log10L( ˆ r )
L(0.5)"
# $
%
& ' = log10 L(ˆ r ) − log10 L(0.5) = 26.8
Ÿž�
�5×&(bisec^on&method)&
€
∂ lnL(r)∂r
= −1441− r
+32r
+57(4r − 2)1− 2r + 2r2
0.0 0.1 0.2 0.3 0.4 0.5
-200
0200
600
1000
r
0ƨƧƼrǀ¥Ƙ&
rƩ§hƨƪƄ�ƨƵNewton-Raphson×ƧƦƏëƊƺƽƼ&
f(x)=0ƥƧƼxǀÔƴƼƨƪƄ&!1. f(a)ƥf(b)ƥƤĂLƏîƧƼƹƋƧ�Ţaƥ�ŢbǀÕƴƼ&
2. aƥbƩşßmǀÔƴƼ&3. f(m)ƩĂLƏf(a)ƥOƗƧƺaǀmƤęƐªƌƼƅf(b)ƥOƗƧƺbǀmƤęƐªƌƼƅ&
4. 2-3ǀĘƻŋƘƓƥƤƄmǀf(x)=0ƥƧƼßƨŊƢƒƣƊƑ&a& b&
1Sô&
2Sô&b&
b&
b&
3Sô&
4Sô&
a&
a&
a&
m&
m&
m&
Ÿſ�
Đªƌ"ǀWUņŧƨ`ªƘƼ�
• Đªƌ"ƨƪ?×�ƏàƊ&– ACCşƩĐªƌ"ƪƄABCşƩĐªƌ"ƥBCCşƩĐªƌ"ƩQƹƻƵoƔƑƧƼ(aś�ªƌƩƝƴ)&
• ƛƓƤÌƩƹƋƧ?×�ƩƉƼr�ǀhěƘƼ&– OƗÂģ�ƨ��ƘƼ2Ŗ�e�şƨƍƊƣƄ�ªƌƏ�XƖƣÂģ5�ƉƝƻxSêƗƼƥƐƄƛƩ2�şƩņŧǀxƤĩƖƄWUņŧƥƹƮƅE�ƪǵǻnjǿȃMȄƵƖƑƪǘǿǜǵǻnjǿ(cM)&
– HaldaneŠ®ƶKosambiŠ®ƧƦƄÇƆƧ`ªŠ®ƏƉƼ&
HaldaneŠ®&
KosambiŠ®& ȃ�ªƌ~Ûǀĝ�ƖƝǵǠǻȄ&
ȃæāƨ�ªƌƏŅƓƼǵǠǻȄ&
€
x = −12ln(1− 2r)
€
x =14ln 1+ 2r1− 2r#
$ %
&
' (
Ÿƀ�
Đªƌ"ƥWUņŧƩŠ#&
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.0
0.1
0.2
0.3
0.4
0.5
x (M)
r
KosambiŠ®&
HaldaneŠ®&r&=&x&
.%�" )/��������0$16#����� 2*(& !�7��������+�5��38���.%�'�4.%�'������,���
Źŷ�
ŏŝWUƩÆć�Ū&1. MŖ�e�Ʃ5ŧÓƩŔN�ÃhǀħƋ&ȃ5ŧÎƳƩÃ4ƄjűƏƋƲƑħƠƣƊƧƊDZȁNjȁƵÃ4ƤƐƼȄ&
2. Ŗ�e�şƩĐªƌ"ǀĴąƘƼ&3. LOD(ƶWUņŧǀĴąƘƼ&4. LOD(ƶWUņŧƩƉƼ(ǀš(ƥƖƣƄ/Ŗ�e�ǀŏŝĚ(linkage&groups)ƨǏǻȁǬ5ƒƘƼ&ȃD5ƧDZȁNjȁ®ƏƉƼ]NƄŏŝĚƪÂģ�ƨ�ġƘƼȄ&
5. ŏŝĚÒƨŖ�e�Ʃƭǀ§hƘƼ&6. ȃ]NƨƹƠƣƪȄWUņŧǀĴąƖƧƍƘ&7. ŏŝWUǀ¨Ƒ&8. ǏǹǪDŽNjǻǕLJǥǚDžǬǀ¨ƊƝƻƄǫȁǡǖǡǹǝǬÃhǀħƠƝƻƖƣƄWUƩùƎƺƖƔǀùļƘƼ&
ŹŸ�
Ě5ƒƩǃǻǑǺǗdzƩ &
• š(×&1. Đªƌ"ƶWUņŧƄLOD(ƨnƖƣš(ǀķhƘƼ&
2. Đªƌ"ƶWUņŧƩ]Nƪš(ƹƻoƔƊ(ǀúƖƝDZȁNjȁnǀƄLOD(Ʃ]NƪbƐƊ(ǀúƖƝDZȁNjȁnǀOƗŏŝĚƨtƘƼƥ8°ƖƣO�ǏǻȁǬƨ5ŭƘƼ&
3. 1ƥ2ǀƄ/DZȁNjȁƏ5ŭƔƽƼƲƤĘƻŋƘ&
ŹŹ�
Ŗ�e�ƩŪ��ƒ�
• wUƩƹƋƧWUǀÆćƘƼƓƥƏºďôÈ&
• ƖƎƖƄwUƩƹƋƧŖ�e�ƩƭȃL1-L2-L3-L4-L5-L6ȄǀÕhƘƼ<ƪƄŖ�e�şƩņŧƞƒƖƎ5ƎƺƧƊ&
• ƛƓƤƄņŧ�\ƎƺƄ1Ì,ƩÂģ��ƨƮŖ�e�ƩŚęǀÕhƖƝƊ&
• ƦƋƘƽƫƹƊƎȎ&
L1&
L2&
L3&
L4&
L5&
L6&
14.5&
0.0&
5.2&
8.3&
17.1&
20.3&
Źź�
Ū��ƒƩįÞ&• ähƩ[Þƨ�ƠƣºŔƧŪ�ǀƳơƒƼ&
– ºbq�ȃMaximum&LikelihoodȄ&
– ºoĐªƌ"QȃMinimumƃsumƃofƃadjacentƃrecombina^onƃfrac^onȋSARFȄ&
∑−
=+++++ −−+=
1
11,1,1,1,1, )]1log()1(log[)(
l
iiiiiiiiiiinRL θθθθ
ƓƓƤȅni,i+1ƪŤ¦ƘƼŖ�e�şƩǒǿǬǻ®&θi,i+1ƪŤ¦ƘƼŖ�e�şƩĐªƌ".&&
∑−
=+=
1
11,
l
iiiSARF θ
ŹŻ�
ƌƫƄ&Đªƌ"ƩQǀºoƨƖƣƳƹƋȂ&
L1&
L2&
L3&
L4&
ƓƩƭƏÍƖƊƥƐƨ…&
L1&
L2&
L3&
L4&
ƌƫƄ&ƭǀL1-L3-L2-L4ƥ&§hƖƝƥƘƼƥ…&
L1&
L2&
L3&
L4&
ĒÁóƨƄĐªƌ"ƩQƏbƐƑƧƼ& Źż�
ºŔƧŪ�Ʃ¥Ď&! [ÞƪEčƞƏƄƉƻƋƼŪ�®ƪğbȂƃ&
" 10Ŗ�e�&…&1,804,400&" 100&&&&&&&&&&&& &…&4.7&x&10158ƃƃƃƃƃƃ→ƃ@çóƧ¥ĎƏ�ĭ&
" 1000&&&&&&&&& &…&2.0&x&102567ƃ&
ƌƫƄMAPMAKER&ƤƪƄ&
c&
t&
k&
q&
p&
5ơý�ƩŖ�e�ƤŲĐƳǀ�ƻ…&&ȃ5Ŗ�e�Ʃư¹ƌƪ5!/2&=&60ōƻȄ&
ÏƻƩŖ�e�ǀŪÌƄºŔƧ]�ƨŚęƖƣƊƑ&ȃq�ǀ[ÞƨƘƼȄ&
a&
c&
t&
k&
q&
p&
r&v&
RŬßȂ&&ºďóƨ�ƺƽƼŪ�ƏƄº7ƨŕƫƽƼŖ�e�ƶƄŚęƔƽƼŖ�eƩŪƨ!fƘƼƅ&ƓƩƹƋƧ²×ƤƪƄ�ƙƖƵƹƊIJƏ�ƺƽƼƥƪŢƺƧƊ&& ŹŽ�
bZóºŔIJƥs�IJ&(global&op^mal&solu^on,&local&op^mal&solu^on)&
s�IJƨ&ñƻơƴƣƖƲƋKĞ�Ƶ&
pƧƑƧƊ&
]�ƝƻóƧ¥ĎƤƪ…&
[Þȃ ƌƫƄq�Ȅ&
IJĀşȃ ƌƫƄŖ�e�ƩŪ�Ȅ&
IJĀşƏ�ƑƄaƑƩs�IJƏƉƼƹƋƧ]NƄ&s�IJƨ£ƿƽƙƄbZóƧºŔIJǀ@çóƨÔƴƼ²×&
ƲƝƪ&ºŔIJƤƉƼ$ĹƪàƊƏƄjë�RŬƩàƊIJȃŊ�ºŔIJȄǀ@çóƨÔƴƼ²×&
Ə�ĭ&
bZóºŔIJ&
Źž�
uSǘȁǻǖDZǿRŬ&Traveling&salesman&problem&(TSP)&
• ĐNƚºŔBRŬȃcombinatorial&op^miza^on&problemȄƩƬƥơ&
ĈV13,509řzȃ�J500���ȄǀuƼŪŇ&
h:p://www.crpc.rice.edu/newsArchive/tsp.html&
a®ƩnłƨơƊƣƄnłşƩņŧħ6Ə�ƌƺƽƝƥƐƨƄ/nłǀºøņŧƤuƼŪŇǀÔƴƼRŬ&ȃMnłƪ1�ƞƒĶRƔƽƼȄ&
Źſ�
ĦǐǽǣȁºŔB&(ant&colony&op^miza^on:&ACO)�
ĦƩ¤ůħAƨƍƒƼǻȁǻŶ�• ŶǪLJǽǵǿƏƧƊ]Nȃ7½å�ȄƪȅǹǿǛdzƨ¥Ď�• ŶůǀðĮƘƼƥ¬ð�ƩǪLJǽǵǿǀ5ÖƖƧƏƺvƨ|ŗ�• ŶǪLJǽǵǿƏ�ƑÏìƖƣƊƼĄ�ǀ+-óƨ¥Ď�
�%ƆƩĦƏƓƩƹƋƧEčƧǻȁǻƨ[ƢƑƓƥƤȅĦƪƇĚƽƥƖƣƈů]ƲƤƩºøǻȁǡǀ=ƻ4ƘƓƥƏƤƐƼ�
ĦǐǽǣȁºŔBƃȀȀȀƃƓƩƹƋƧĦƩħAǻȁǻǀǐǿǩǶȁǚ�Ƥ�óƨ9ëƖƄĐNƚºŔBRŬǀIJƑ� Źƀ�
ACOƩºŔBƩǬǽǘǖ�1. ĦǀřzȃŖ�e�ȄƨǹǿǛdzƨŚęƘƼ�
2. Ħƪȅ¾ĶRřzƩƎƺȅřzşƨĥÿƔƽƝǪLJǽǵǿŜƥřzşƩ
ņŧȃĐªƌ"ƄƲƝƪn®q�ƩĔn(ȄƩ�²ǀĝ�Ɩƣȅřzǀŕ
ƭüAƘƼ�
3. /řzǀĶRƘƼƲƤ2ǀĘƻŋƘ�
4. uSŇƩøƔȃĐªƌ"ƩQƄƲƝƪƄn®q�ƩĔn(ƩQȄƨ�Ɨƣȅ
ōőƖƝǻȁǡƨǪLJǽǵǿǀ5ÖƘƼŶȃµşƥƥƵƨǪLJǽǵǿƪ¬ðƘ
ƼȄ�
5. 1Ȇ4ǀĘƻŋƖȅºŔŊ�IJǀ�Ƽ�
Ŗ�e�A&
Ŗ�e�B&
Ŗ�e�C& Ŗ�e�E&
Ŗ�e�F&Ŗ�e�şƩņŧƪ&Đªƌ"ƲƝƪn®q�ƩĔn(&
źŷ&
ĦƩüAµƩǻȁǻ�
• µ;tƨƍƒƼĦkƏřziƎƺřzjƨüAƘƼùç�
.][)]([
][)]([)( k
iNl ilil
ijijkij Nj
dtdt
tpki
∈∀=∑∈
−
−
βα
βα
ττ
ijd)(tijτ ȋŶǪLJǽǵǿŜ� ȋŶřzşņŧ�
ņŧƏœƑƣƵƄđŇƨĥÿƔƽƝǪLJǽǵǿŜƏaƒƽƫƄƹƻŕ�ƔƽƶƘƑƧƼ&ȃƝƞƖƄĔnƨŕ�ƔƽƼƥƊƋƿƒƤƪƧƑƄƇùçƏ�ƏƼƈƥƊƋƓƥƨØ�ƘƼȄ&&ƓƽƨƹƠƣƄƇpƖœSƻǀƖƝƱƋƏƄjƪŊŒƈƥƊƋƹƋƧđŇƏŕƫƽƼƹƋƨƧƼ&
źŸ&
⎩⎨⎧
∉∈
=Δ)(),(if0)(),(if)(
)(tTjitTjitLQ
t k
kkkijτ
ǪLJǽǵǿŜƩ·±�
∑ =Δ+=+ m
kkijijij ttt
1)()()1( τρττ
uSņŧƏøƊĦƪƄƛƩđŇƨƹƻaƑƩǪLJǽǵǿǀ5ÖƘƼȃŜƪƄņŧƨGÓ Ȅ&
uSņŧƏŞƊĦƪƄǪLJǽǵǿ5ÖŜƏpƧƊ&
kíôƩĦƏ5ÖƘƼŜ&
ƝƞƖƄ&
ρƪ1ƹƻoƔƊƅƘƧƿƟ&ǪLJǽǵǿƪµşƥƥƵƨ¬ðƘƼ&
źŹ�
ĦƩħŐƏĘƻŋƔƽƼƥ…&
º7ƪƄƊƾƊƾƧđŇƨ&ĥÿƔƽƊƝǪLJǽǵǿƏ…&
ƞǁƞǁƄºøđŇƯƩĥÿŜ&ƏaƑƧƠƣƊƑ&
Đªƌ"ƩQƄƲƝƪƄn®q�ƩĔn(ƩQǀºoBƘƼ&Ŗ�e�ƩŪ�ƏÔƴƺƽƼ&
źź�
ŏŝWU��ǙǪǡdžLJǃAntMap�• AntMapƤƪųŎƨŏŝWUǀ��ƤƐƼ�
• ŏŝWU��Ʃ�ÄƏ/ĠABƔƽƣƍƻƄǷȁǓȁƪǯǚǿ�ơƤWU
ǀ�īƘƼƓƥƵƤƐƼ�
• AntMapƤƪƄųŎǃǻǑǺǗdzƩ9ßǀêƎƖƄ��ƔƽƝWUƩĊ�ǀ
ǫȁǡǖǡǹǝǬÃhƨƹƠƣùļƤƐ
Ƽ�
• ÇƆƧOSȃWindows,&Mac,&Linux,&SolarisƧƦȄ�ƤjħKĞƅǬǽǏǹDzǿǏijĽƪJava&
h:p://lbm.ab.a.u-tokyo.ac.jp/~iwata/antmap/& -2���� źŻ�
ƭƏ�ùƎƧŖ�e�ǀƄǫȁǡǖǡǹǝǬ(ƥƖƣùļƤƐƼ�
Ŗ�e�ƭ¹ƌƩųŎBƨƹƻƄ�ƭƩùƎƺƖƔƩĺ"ƵKĞƥƧƠƝ�
ǫȁǡǖǡǹǝǬȃbootstrapȄ×ȋƃ/ǒǿǬǻƎƺśĬǀĸƖO®ǒǿǬǻǀ2¡4ƘƼƅ2¡4ǠȁǚƨơƊƣôóƩēĴŜǀ2ĴąƘƼƅƓƩǬǽǘǖǀa®SĘƻŋƖƣƄ§hƔƽƼēĴŜƩǦǹǞǍǀĺ"ƘƼƅ�
Ŗ�eƩƭƩùƎƺƖƔƩĺ"�
źż�
jǠȁǚƯƩŔë&
• RGPƩWebǮȁǕȃh:p://rgp.dna.affrc.go.jp/Ȅƹƻ.�ƖƝDžǤƩȊmċƩ5ŧǠȁǚǀëƊƣȅAntMapƨƹƼŏŝWU��ǀĻƳƝȇ&&mċ①ǃǍǨNjǺ/ǐǔǨNjǺŶDHLsȃ169DZȁNjȁȄ&&mċ②³¿¶/Kasalath//³¿¶ŶBILsȃ245DZȁNjȁȄ&&mċ③ƉƛƳƩƻ/IR24&RILsȃ375DZȁNjȁȄ&
źŽ�
ŏŝWU&
• Ŗ�eƩŪ��ƒƨĭƖƝµş&– mċ①ƃ2.5û&– mċ②ƃ3û&– mċ③ƃ7û&
ƃƃȃCPUȋIntel&Mobile&Pen^um&1.6GHȏȄ&
• 100G�ƩǫȁǡǖǡǹǝǬÃh&– mċ③ƃ385&
• §hƔƽƝŪ�ƪȅ�ŘDZȁNjȁȃǠȁǚ①ȅ②ȅ③ƤƛƽƜƽ5,&2,&5DZȁNjȁȄǀţƊƣȅ0ĩƔƽƣƊƼŪ�ƨ�ġƖƝȇ&
mċ①&
mċ②&
mċ③&źž�
ǔDzǶǼȁǔǸǿǠȁǚƯƩŔë&
• 1,200DZȁNjȁ&• 10ŏxĚ&
Č2.55&
źſ�
Ě÷Ğ&(Swarm&intelligence)�
• EčƧǻȁǻƨ[ƢƊƣħAƘƼLjȁǕLJǿǡȃ ƌƫƄĦƄŵƄŴȄƏƄ�ƊƨȃaƑƩ]NƄs�óƨȄ�ëƖƉƋƓƥƨƹƻƄcťÊóƧ:�ǔǖǟdzǀƵƝƧƑƣƵƄĚƽƥƖƣƩħAƏ>ðƔƽƼƓƥƏƉƼȃ ƌƫƄŴƩĚƽƩAƐȄ&
• ƓƩƹƋƧêãƩĚƽƥƖƣƩħAǀǐǿǩǶȁǚ�ƤÉ'ƘƼƓƥƨƹƻƄĬŦƧºŔBRŬƩjëóƧIJǀ¥ĎƖƹƋƥƘƼ�Ĩ�
źƀ�
ĉeĚºŔB&(par^cle&swarm&op^mizer:&PSO)&
• ǧǹǴȁǚĀş�ǀŮƭSƼĉeĚǀëƊƣºŔBǀħƋǃǻǑǺǗdz&
• ƉƼµßƤƩMĉeƩŔ��ƪƄƛƩµßƨĉeƏ�ęƘƼĀş�Èƨ[ƢƊƣĴąƔƽƼƅ&
• MĉeƪŔ��ƩĠyǭǖǡƏ�ƺƽƝ�Èǀĵ�ƤƐƼƅƲƝƄĉeşƤ�\�ªƘƼƓƥƤƄĚ1ǭǖǡ�ÈƵĵ�ƖƣƊƼƅ&
• ĉeƏĀş1ǀŮƭSƼƥƐƩŎ�ǭǎǡǻƪƄȈơ<ƩµßƤƩŎ�ǭǎǡǻƄĠyǭǖǡ�ÈƄĚ1ǭǖǡ�Èƨ[ƢƊƣÕhƔƽƼ&
PSO&UIJ&
Ě1ǭǖǡ&
ĉeA&
ĉeB&
ĉeC&
ĉeAƩ&Ġyǭǖǡ&
€
c1r1(x p − x)€
c2r2(xg − x)
€
xg
€
x p
€
x€
wvt−1
€
vt
€
x t ← x t−1 + vtvt ← wvt−1 + c1r1(x p − x) + c2r2(xg − x)$ % &
��3� 6!>?B@CA!
€
f (x)
w,!c1,&c2ƪƉƺƎƗƴķhƔƽƝ(ȃ1ƨŊƊ((w&<&1)ƏƹƊƥƔƽƼȄƄ&r1,&r2ƪǃǝǬǠȁǡÒƨê�ƔƽƼ[0,&1]Ʃ�Ç�®ƤƉƼƅ&
PSOƨƹƼºŔBƩ �
ŻŹ�
f (x) =10n+ xi2 −10cos2π xi( )
i=1
n
∑
Global&minimum&(bZóºo(Ȅȋ�
x = 0 f (x) = 0ƩƥƐ�
RastriginŠ®�
n&=&2�
R&code�
Żź�
# load required packagesrequire(rgl)require(pso)# set x and y arraysx <- (-50:50)/10y <- x# set a objective function (surface)z <- matrix(NA, length(x), length(y))for(i in 1:length(x)) {
for(j in 1:length(y)) { # Rastrigin function z[i,j] <- 20 + (x[i]^2 - 10 * cos(2 * pi * x[i]))
+ (y[j]^2 - 10 * cos(2 * pi * y[j]))}
}# show the objective surfaceopen3d()persp3d(x, y, z, col = "green")# optimizationo1 <- psoptim(rep(NA,2),function(x) 20+sum(x^2-10*cos(2*pi*x)), lower=-5,upper=5,control=list(abstol=1e-8))show(o1)�
Ʋƥƴ&• ŏŝWUƪƄŖ�eşƩĐªƌū�ǀ�«ƘƼƝƴƨ¼ëƤƉƼƅQTLIJÀƨƵŏŝWUƏëƊƺƽƼ&
• ŏŝWUƪƄŖ�eşƤĴąƔƽƼĐªƌ"ƨ[ƢƊƣƄŖ�e�ƩĚ5ƒƄƭƎƌƏħƿƽƼƓƥƨƹƻÆćƔƽƼ&
• ƭƎƌƨƪ]NƩ®ƏÅƴƣbƐƊƅŖ�e�®ƏaƊ]NƨƪäƨƄ@çóƧǃǻǑǺǗdzƨƹƼºŔŪ�Ʃ¥ĎƏ�ĭƥƧƼ&
ŻŻ�
Żż�
ac!!1.0�
AC&0.5&(1-r)�
AaCc&&&&&&&89�
Ac&0.5r�
Aacc&&&&&&&&2�
aC&0.5r�
aaCc&&&&&&&&12�
ac&&0.5&(1-r)�
aacc&&&&&&&&97�
ab!!1.0�
AB&0.5&(1-r)�
AaBb&&&&&&&80�
Ab&0.5r�
Aabb&&&&&&&&11�
aB&0.5r�
aaBb&&&&&&&&21�
ab&&0.5&(1-r)�
aabb&&&&&&&&88�
bc!!1.0�
BC&0.5&(1-r)�
BbCc&&&&&&&89�
Bc&0.5r�
Bbcc&&&&&&&&12�
bC&0.5r�
bbCc&&&&&&&&12�
bc&&0.5&(1-r)�
bbcc&&&&&&&&87�
ŻŽ�
ǦDžNJǴǡǺǝǎǖă6SoǼǰȁǡĿŬ�
1. <ƩǖǹDžǢƩĩƨ[ƢƊƣƄŖ�e�A-BşƨơƊƣƄn®q�ǀĐªƌ"rƩŠ®lnL(r)ƥƖƣ¸Ɛ�ƖƄƛƽǀƵƥƨrƩºq§hŜǀÔƴƣ�ƔƊƅƲƝOƗ�ƨ[ƢƐƄŖ�e�A-C,&B-CşƨơƊƣƵĐªƌ"ǀÔƴƣ�ƔƊƅ&&
2. 1ƤÔƴƝĐªƌ"ƨn�ƘƼHaldaneƩWUņŧƥƄKosambiƩWUņŧǀǘǿǜǵǻnjǿE�ƤÔƴƣ�ƔƊƅ&ȃƧƍƄlogƩĴąƪgoogleƤƵƤƐƲƘƅ ƌƫƄln(0.3)ƥƖƣÃĎƘƼƥ0.3ƩĠán®ƏĴąƤƐƲƘƅȄ&
3. Ŗ�e�AƄBƄCƪƦƩŪƨǁƤƊƼƥĝƌƼƩƏd�ƤƖƸƋƎȎ�
Żž�
ǼǰȁǡƩ©4×�
Ȍ©4Ʃ�²ȍ&ǴȁǻƩ¿¯ƄƲƝƪƄPDFƶMS−WordǪljȁDZǝǡȃƤƐƽƫPDFȄƩǪǂDžǻƥƖƣ��ƖƄe-mailƨÚ�ƖƣƄ[email protected]ƨŌƠƣ�ƔƊȃĠAƤIéǴȁǻƏŌƻŋƔƽƲƘȄƅ&
Ȍ5Ŝȍ&äƨ¢hƪƖƲƚǁƏƄaƑƣƵA4Ƥ1ǮȁǕ�1ƨƲƥƴƣ�ƔƊƅ&Ȍ½Ţȍ&2014�5»29³ƲƤ&
Żſ�