¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å...

12
060310391 ȃʼngŘȄ Ȁ 0560565 ȃégŘȄ ǦDžNJǴǡǺǝǎǖ ă6S ŏŝWU 2014/5/22 8:40-10:10 @1-ă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ǝǗ WUŠ® bZIJƄsIJ ĐNƚºŔBRŬƄuSǘȁǻǖDZǿRŬ Ě÷ĞƄĦǐǽǣȁºŔB ź ŏŝWU vs. ãéWU ãéWU DNA^[n®ƧƦƄãéóƧņŧƨ[ƢƑWU ŏŝWU ŖeşƩŏŝƩƔȃĐªƌ"ȄǀņŧƥƖƝ WU ĵ2þŭƩņŧƪƊƨÓƖƧƊ ŖeşƩĐªƌƩūǀ«ƘƼƨƪƄŏŝ WUƩ\Əśĭ Ż

Transcript of ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å...

Page 1: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

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Ʃ�\Əśĭ&

Ż�

Page 2: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

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�

ſ�

Page 3: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

Ȏ�ĆƩƩòèƩ=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ŜǀÔƴƹƅ�

Ȏ�

ŸŹ�

Page 4: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

&Ü®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

ŸŽ�

Page 5: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

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.%�'������,���

Źŷ�

Page 6: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

ŏŝ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 θ

ŹŻ�

Page 7: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

ƌƫƄ&Đªƌ"Ʃ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ƔƽƼȄ&

Źſ�

Page 8: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

ĦǐǽǣȁºŔ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ǽǵǿƪµşƥƥƵƨ¬ðƘƼ&

źŹ�

Page 9: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

ĦƩħŐƏĘƻŋƔƽƼƥ…&

º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ȁȄ&

źŽ�

Page 10: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

ŏŝ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ƔƽƼ&

Page 11: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

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ƨƹƼºŔŪ�Ʃ¥ĎƏ�ĭƥƧƼ&

ŻŻ�

Page 12: ¿³ © Í þ â 060310391& - lecture.ecc.u-tokyo.ac.jp060310391& Ig X v &0560565& ég X & && æ Å Ê ô á ú Ý Î Ö& 6S& O ]WU & 2014/5/22&&&8:40-10:10& @1L p- 4 A k & x V e O

Żż�

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³ƲƤ&

Żſ�