>512 :@>;0A/@5;:@; 5-3:;?5:3 #1>2;>9-:/1treethinkers.org/wp-content/uploads/2013/03/MCMC.pdf · '...
Transcript of >512 :@>;0A/@5;:@; 5-3:;?5:3 #1>2;>9-:/1treethinkers.org/wp-content/uploads/2013/03/MCMC.pdf · '...
���>512��:@>;0A/@5;:�@;��5-3:;?5:3� � ��#1>2;>9-:/1
�=4,9�$��::=0�1<->@91:@�;2��B;8A@5;:����/;8;3E
'����-B5?
;018�.-?10�5:21>1:/1�5?�.-?10�;:�@41�9;018
�:/07��,>0/�"3D7:2090?4.��910=09.0
���:/07�>;0.414.,?4:9
8:/07�>070.?4:9
8:/07�,/0<@,.D
8:/07�@9.0=?,49?D,A0=,2492
��>?48,?492�@9/0=�?30�8:/07
7460743::/�:;?484E,?4:9
�����>48@7,?4:9
� ��5:�@41;>E�-:0�<>-/@5/1
������;;=:C48,?4:9�:1�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D
�����49�?30:=D
,9�,;;=:;=4,?07D�.:9>?=@.?0/�,9/�,/0<@,?07D�=@9�.3,49�4>�2@,=,9?00/�?:�;=:A4/0�,9�,=-4?=,=47D�;=0.4>0�/0>.=4;?4:9�:1�?30�5:49?�>?,?4:9,=D�/09>4?D�
�����49�;=,.?4.0
,7?3:@23�,�24A09�>,8;70=�8,D�B:=6�B077�49�8:>?�.,>0>��,77�>,8;70=>�B477�1,47�49�>:80�.,>0>��,9/�4>�9:?�2@,=,9?00/�?:�B:=6�1:=�,9D�24A09�.,>0�
#�(309�/:�B0�69:B�?3,?�?30������;=:A4/0>�,9�,..@=,?0�,;;=:C48,?4:9����1:=�,�24A09�08;4=4.,7�,9,7D>4>�
!�(�$��
� ��<1>2;>9-:/1�-:0�"��
������;;=:C48,?4:9�:1�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D
�?�4>�9:?�>@114.409?�?:�80=07D�-0�/00;7D�.:9.0=90/�,-:@?������;0=1:=8,9.0
1:=�,9D��,D0>4,9�4910=09.0�-,>0/�:9������D:@�900/�?:�-0�.:8;70?07D�:->0>>0/�,-:@?�4?�
;,=?4.@7,=7D�1:=�.:8;70C�8:/07>4910=09.0�;=:-708>�
!@?7490
����5-3:;?5:3� � ��<1>2;>9-:/1
����>1B51C�;2�C41>1�C1’B1�.11:�-:0�C4E
�����5-3:;?@5/?�.-?10�;:�?5:381�/4-5:?
�(��5-3:;?@5/?�.-?10�;:�@41�<>5;>�
(��5-3:;?@5/?�.-?10�;:�9A8@5<81�/4-5:?�
#>18595:->E��:-8E?1?
$-@1?�;2��5:1-31��5B1>?525/-@5;:
#4E8;31:E��?@59-@5;:
� >��-E1?�
�5?/>1@1�&>-5@��B;8A@5;:
�-E1?&>-5@?
�5B1>31:/1�&591��?@59-@5;:
����%&�∆∆ ∆∆ ∆∆ ∆∆ ∆∆
∆∆ ∆∆
�;:@5:A;A?�&>-5@�B;8A@5;:
�-E1?&>-5@?
?E991&$�� 95?/�$
95?/�$
&"#"�"�* ��$"!"�$�
>�-E1?��<-5>10� � �/4-5:?<->@5@5;:�?<1/525/�9;018
���%&��/4-5:?
'��!�>18-D10�/8;/7
95?/�$
%?,?4>?4.,7��>?48,?4:9�:1�"3D7:209D���9�!@?7490
$-@1?�;2�&>-5@��B;8A@5;:
95?/�$95?/�$
#>18595:->E��:-8E?1?
$-@1?�;2��5:1-31��5B1>?525/-@5;:
#4E8;31:E��?@59-@5;:
� >��-E1?�
�5?/>1@1�&>-5@��B;8A@5;:
�-E1?&>-5@?
�5B1>31:/1�&591��?@59-@5;:
����%&�∆∆ ∆∆ ∆∆ ∆∆ ∆∆
∆∆ ∆∆
�;:@5:A;A?�&>-5@�B;8A@5;:
�-E1?&>-5@?
?E991&$�� 95?/�$
95?/�$
&"#"�"�* ��$"!"�$�
>�-E1?��<-5>10� � �/4-5:?<->@5@5;:�?<1/525/�9;018
���%&��/4-5:?
'��!�>18-D10�/8;/7
95?/�$
%?,?4>?4.,7��>?48,?4:9�:1�"3D7:209D���9�!@?7490
$-@1?�;2�&>-5@��B;8A@5;:
95?/�$95?/�$
�,D0>4,9��910=09.0�:1�"3D7:209D�49�!90�%74/0
��&=00
Pr(Parameter!Data) = Pr(Data!Parameter) x Pr(Parameter) Pr(Data)
<;?@1>5;>�<>;.-.585@E 8571854;;0�2A:/@5;: <>5;>�<>;.-.585@E
9->35:-8�8571854;;0
<;?@1>5;>�<>;.-.585@E
���#4E8;31:1@5/�9;018�<->-91@1>?�>>@80�,9�,7429809?,�X,�:1�N�>4?0>�1:=�S�>;0.40>�
���-@-
X = x1,x2,...,xN( )!
?:;:7:2D τ = τ1,τ2,...,τ(2s−5)!!( ) !-=,9.3�7092?3> ν = ν1,ν2,...,ν(2S−3)( ) !
>@->?4?@?4:9�=,?0> θ = θAC ,θAG ,θAT ,θCG,θCT ,θGT( )! ��:/07�:1�.3,=,.?0=�.3,920 Φ = θ,π,α,T( ) !
>?,?4:9,=D�1=0<@09.40> π = π A ,πC ,πG,πT( )!,8:92�>4?0�=,?0�A,=4,?4:9 α !
����#4E8;31:1@5/�8571854;;0�2A:/@5;:
�(�#>5;>?�;:�<->-91@1>?~Uniform~Dirichlet (1,...,1)
~Dirichlet (1,1,1,1,1,1)~Dirichlet (1,1,1,1)~Uniform (0.1,100)
L(τ,ν,Φ) = f (X |τ,ν,Φ) = f (xi |τ j ,ν k,Φ) f (ν k,Φ)dν ki=1
N
∏ν k ,Φ
2S−3
∫j=1
(2S−5)!!
∑ dΦ !
(�#;?@1>5;>�#>;.-.585@E
f (τ,ν,Φ | X) = f (X |τ,ν,Φ) f (τ,ν,Φ)f (X)
%-9<81?�2>;9�@41�/4-5:�-<<>;D59-@1�@41�6;5:@�<;?@1>5;>�
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
&30�1=0<@09.D�:1�>,8;70/�;,=,80?0=�A,7@0>�;=:A4/0>�,�A,74/�0>?48,?0�:1�����?30�;:>?0=4:=�;=:-,-474?D�:1�?3,?�;,=,80?0=
(0�.,9�<@0=D�?30�5:49?�;:>?0=4:=�B4?3�=0>;0.?�?:�,9D�49/4A4/@,7�;,=,80?0=����:1�49?0=0>?��?30�8,=249,7�;:>?0=4:=�;=:-,-474?D
#>;3>-995:3�;A>� � ��>;.;@!@=�=:-:?�;,=,.3@?0>�49?:�,�=,9/:8�7:.,?4:9�49�?30�5:49?�;:>?0=4:=�/09>4?D
,9/�B477�0C;7:=0�;,=,80?0=�>;,.0�-D�1:77:B492�?30>0�>48;70�=@70>�
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
1@>;<;85?�1@�-8����� ����-?@5:3?������
#>;3>-995:3�;A>� � ��>;.;@!@=�=:-:?�;,=,.3@?0>�49?:�,�=,9/:8�7:.,?4:9�49�?30�5:49?�;:>?0=4:=�/09>4?D
,9/�B477�0C;7:=0�;,=,80?0=�>;,.0�-D�1:77:B492�?30>0�>48;70�=@70>�
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
1@>;<;85?�1@�-8����� ����-?@5:3?������
#>;3>-995:3�;A>� � ��>;.;@!@=�=:-:?�;,=,.3@?0>�49?:�,�=,9/:8�7:.,?4:9�49�?30�5:49?�;:>?0=4:=�/09>4?D
,9/�B477�0C;7:=0�;,=,80?0=�>;,.0�-D�1:77:B492�?30>0�>48;70�=@70>�
���1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�@;3477��4?�,@?:8,?4.,77D�?,60>�?30�>?0;
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
1@>;<;85?�1@�-8����� ����-?@5:3?������
Pr[Accept] = 1�
#>;3>-995:3�;A>� � ��>;.;@!@=�=:-:?�;,=,.3@?0>�49?:�,�=,9/:8�7:.,?4:9�49�?30�5:49?�;:>?0=4:=�/09>4?D
,9/�B477�0C;7:=0�;,=,80?0=�>;,.0�-D�1:77:B492�?30>0�>48;70�=@70>�
���1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�@;3477��4?�,@?:8,?4.,77D�?,60>�?30�>?0;
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
1@>;<;85?�1@�-8����� ����-?@5:3?������
#>;3>-995:3�;A>� � ��>;.;@!@=�=:-:?�;,=,.3@?0>�49?:�,�=,9/:8�7:.,?4:9�49�?30�5:49?�;:>?0=4:=�/09>4?D
,9/�B477�0C;7:=0�;,=,80?0=�>;,.0�-D�1:77:B492�?30>0�>48;70�=@70>�
���1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�@;3477��4?�,@?:8,?4.,77D�?,60>�?30�>?0;
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
��1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�/:B93477��4?�/4A4/0>�?30�070A,?4:9�:1?30�;=:;:>0/�7:.,?4:9�-D�?30�.@==09?�7:.,?4:9��,9/�4?�:97D�?,60>�?30�>?0;�41�?30�<@:?409?�4>�70>>�?3,9�,�@941:=8�=,9/:8�A,=4,-70��U[0,1]�
1@>;<;85?�1@�-8����� ����-?@5:3?������
Pr[Accept] = new height old height�
#>;3>-995:3�;A>� � ��>;.;@!@=�=:-:?�;,=,.3@?0>�49?:�,�=,9/:8�7:.,?4:9�49�?30�5:49?�;:>?0=4:=�/09>4?D
,9/�B477�0C;7:=0�;,=,80?0=�>;,.0�-D�1:77:B492�?30>0�>48;70�=@70>�
���1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�@;3477��4?�,@?:8,?4.,77D�?,60>�?30�>?0;
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
��1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�/:B93477��4?�/4A4/0>�?30�070A,?4:9�:1?30�;=:;:>0/�7:.,?4:9�-D�?30�.@==09?�7:.,?4:9��,9/�4?�:97D�?,60>�?30�>?0;�41�?30�<@:?409?�4>�70>>�?3,9�,�@941:=8�=,9/:8�A,=4,-70��U[0,1]�
��&30�;=:;:>,7�/4>?=4-@?4:9�4>�>D880?=4.,7��>:�Pr[A B] = Pr[B A]�
1@>;<;85?�1@�-8����� ����-?@5:3?������
#>;3>-995:3�;A>� � ��>;.;@!@=�=:-:?�;,=,.3@?0>�49?:�,�=,9/:8�7:.,?4:9�49�?30�5:49?�;:>?0=4:=�/09>4?D
,9/�B477�0C;7:=0�;,=,80?0=�>;,.0�-D�1:77:B492�?30>0�>48;70�=@70>�
���1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�@;3477��4?�,@?:8,?4.,77D�?,60>�?30�>?0;
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
��1�?30�;=:;:>0/�>?0;�B477�?,60�?30�=:-:?�/:B93477��4?�/4A4/0>�?30�070A,?4:9�:1?30�;=:;:>0/�7:.,?4:9�-D�?30�.@==09?�7:.,?4:9��,9/�4?�:97D�?,60>�?30�>?0;�41�?30�<@:?409?�4>�70>>�?3,9�,�@941:=8�=,9/:8�A,=4,-70��U[0,1]�
��&30�;=:;:>,7�/4>?=4-@?4:9�4>�>D880?=4.,7��>:�Pr[A B] = Pr[B A]�
1@>;<;85?�1@�-8����� ����-?@5:3?������
��?8534@8E�9;>1�2;>9-8�01?/>5<@5;:�;2�@41� 1@>;<;85?��-?@5:3?�-83;>5@49
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
���94?4,74E0�?30�.3,49�B4?3�>:80�=,9/:8�A,7@0>�1:=�,77�;,=,80?0=>��49.7@/492?30�?=00�B4?3�-=,9.3�7092?3>��Θ = τ,ν( )
�%070.?�,�;,=,80?0=�?:�.3,920�,..:=/492�?:�4?’>�;=:;:>,7�;=:-,-474?D��02� The MCMC sampler will use the following moves: With prob. Chain will change 2.86 % param. 1 (tratio) with Dirichlet proposal 2.86 % param. 2 (revmat) with Dirichlet proposal 2.86 % param. 3 (revmat) with Dirichlet proposal 2.86 % param. 4 (state frequencies) with Dirichlet proposal 2.86 % param. 5 (state frequencies) with Dirichlet proposal 2.86 % param. 6 (state frequencies) with Dirichlet proposal 2.86 % param. 7 (state frequencies) with Dirichlet proposal 2.86 % param. 8 (state frequencies) with Dirichlet proposal 2.86 % param. 9 (gamma shape) with multiplier 2.86 % param. 10 (gamma shape) with multiplier 2.86 % param. 11 (gamma shape) with multiplier 2.86 % param. 12 (gamma shape) with multiplier 2.86 % param. 13 (gamma shape) with multiplier 2.86 % param. 14 (prop. invar. sites) with beta proposal 2.86 % param. 15 (rate multiplier) with Dirichlet proposal 14.29 % param. 16 (topology and branch lengths) with LOCAL 42.86 % param. 16 (topology and branch lengths) with extending TBR
1@>;<;85?�1@�-8����� ����-?@5:3?������
��?8534@8E�9;>1�2;>9-8�01?/>5<@5;:�;2�@41� 1@>;<;85?��-?@5:3?�-83;>5@49
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
���94?4,74E0�?30�.3,49�B4?3�>:80�=,9/:8�A,7@0>�1:=�,77�;,=,80?0=>��49.7@/492?30�?=00�B4?3�-=,9.3�7092?3>��
�%070.?�,�;,=,80?0=�?:�.3,920�,..:=/492�?:�4?’>�;=:;:>,7�;=:-,-474?D
��"=:;:>0�,�.3,920�?:�?30�>070.?0/�;,=,80?0=�@>492�?30�;,=,80?0=�>;0.414.���;=:;:>,7�80.3,94>8�?3,?�4>�
Θ = τ,ν( )
(i) >?:.3,>?4.(ii) 4==0/@.4-70(iii) ,;0=4:/4.
1@>;<;85?�1@�-8����� ����-?@5:3?������
��?8534@8E�9;>1�2;>9-8�01?/>5<@5;:�;2�@41� 1@>;<;85?��-?@5:3?�-83;>5@49
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
1@>;<;85?�1@�-8����� ����-?@5:3?������
���94?4,74E0�?30�.3,49�B4?3�>:80�=,9/:8�A,7@0>�1:=�,77�;,=,80?0=>��49.7@/492?30�?=00�B4?3�-=,9.3�7092?3>��
�%070.?�,�;,=,80?0=�?:�.3,920�,..:=/492�?:�4?’>�;=:;:>,7�;=:-,-474?D
��"=:;:>0�,�.3,920�?:�?30�>070.?0/�;,=,80?0=�@>492�?30�;,=,80?0=�>;0.414.���;=:;:>,7�80.3,94>8
Θ = τ,ν( )
���,7.@7,?0�?30�;=:-,-474?D�:1�,..0;?492�?30�;=:;:>0/�.3,920�
R =min 1, f X |Θ'( )f X |Θ( )
⋅f Θ'( )f Θ( )
⋅f Θ |Θ'( )f Θ' |Θ( )
$
% &
'
( ) !
8571854;;0�>-@5;� <>5;>�>-@5;� <>;<;?-8�>-@5;�
��?8534@8E�9;>1�2;>9-8�01?/>5<@5;:�;2�@41� 1@>;<;85?��-?@5:3?�-83;>5@49
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
1@>;<;85?�1@�-8����� ����-?@5:3?������
���94?4,74E0�?30�.3,49�B4?3�>:80�=,9/:8�A,7@0>�1:=�,77�;,=,80?0=>��49.7@/492?30�?=00�B4?3�-=,9.3�7092?3>��
�%070.?�,�;,=,80?0=�?:�.3,920�,..:=/492�?:�4?’>�;=:;:>,7�;=:-,-474?D
��"=:;:>0�,�.3,920�?:�?30�>070.?0/�;,=,80?0=�@>492�?30�;,=,80?0=�>;0.414.���;=:;:>,7�80.3,94>8
Θ = τ,ν( )
���,7.@7,?0�?30�;=:-,-474?D�:1�,..0;?492�?30�;=:;:>0/�.3,920
���090=,?0�,�@941:=8�=,9/:8�A,=4,-70��U[0,1]��,..0;?�41�R > U
��?8534@8E�9;>1�2;>9-8�01?/>5<@5;:�;2�@41� 1@>;<;85?��-?@5:3?�-83;>5@49
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3�����
1@>;<;85?�1@�-8����� ����-?@5:3?������
���94?4,74E0�?30�.3,49�B4?3�>:80�=,9/:8�A,7@0>�1:=�,77�;,=,80?0=>��49.7@/492?30�?=00�B4?3�-=,9.3�7092?3>��
�%070.?�,�;,=,80?0=�?:�.3,920�,..:=/492�?:�4?’>�;=:;:>,7�;=:-,-474?D
��"=:;:>0�,�.3,920�?:�?30�>070.?0/�;,=,80?0=�@>492�?30�;,=,80?0=�>;0.414.���;=:;:>,7�80.3,94>8
Θ = τ,ν( )
���,7.@7,?0�?30�;=:-,-474?D�:1�,..0;?492�?30�;=:;:>0/�.3,920
���090=,?0�,�@941:=8�=,9/:8�A,=4,-70��U[0,1]��,..0;?�41�R > U
��$0;0,?�>?0;>� �����-:�‘-01=A-@1’�:A9.1>�;2�@591?
!@?7490
����5-3:;?5:3� � ��<1>2;>9-:/1
����>1B51C�;2�C41>1�C1’B1�.11:�-:0�C4E
�����5-3:;?@5/?�.-?10�;:�?5:381�/4-5:?
�(��5-3:;?@5/?�.-?10�;:�@41�<>5;>�
(��5-3:;?@5/?�.-?10�;:�9A8@5<81�/4-5:?�
���:9A0=209.0
��4C492
��%,8;7492�49?09>4?D
�,>�?30�=:-:?�.:770.?0/�09:@23�>,8;70>�?:�,/0<@,?07D�/0>.=4-0�?30�;:>?0=4:=����;=:-,-474?D�/4>?=4-@?4:9�
�>>0>>492������"0=1:=8,9.0�&3=00��,49��>>@0>
�>�?30�.3,49��=:-:?��>@..0>>1@77D�49?02=,?492�:A0=�?30�5:49?�;:>?0=4:=�;=:-,-474?D�
�,>�?30�.3,49��=:-:?��>@..0>>1@77D�?,=20?0/�?30�>?,?4:9,=D�/4>?=4-@?4:9��
!@?7490
����5-3:;?5:3� � ��<1>2;>9-:/1
����>1B51C�;2�C41>1�C1’B1�.11:�-:0�C4E
�����5-3:;?@5/?�.-?10�;:�?5:381�/4-5:?
�(��5-3:;?@5/?�.-?10�;:�@41�<>5;>�
(��5-3:;?@5/?�.-?10�;:�9A8@5<81�/4-5:?�
���:9A0=209.0�/4,29:>?4.>
(i) &480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�>:80�;,=,80?0=>�,=0�8:=0�=074,-70�?3,9�:?30=>
�>?0;>�8,D�:..@=�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
-,/�.:9A0=209.0 -0??0=�.:9A0=209.0
�C,8;70��&=,.0=�;7:?>�:1�?=00�7092?3�,?�?B:�>?,20>�:1�,�>49270��=�,D0>�=@9
Cornus_3.8_1.tre
elength.out
State0 1000000 2000000 3000000 4000000 5000000 6000000 7000000 8000000 9000000
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
1.1
1.2
79�
�?;91C4-@�0-@-�?1@�01<1:01:@
1,>?� >7:B�
-,>0�1=0< >@-�=,?0> �%$' &� ?:;:7:2D
Cornus_3.8_1.tre
elength_short.out
State0 500000 1000000 1500000 2000000 2500000 3000000
0.6
0.7
0.8
0.9
1
1.1
1.2
�����;,?3:7:240>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
%7:B�.:9A0=209.0�:1�?480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�@9/0=�;,=,80?0=4E0/�8:/07
�>.,70�:1�?@9492�;,=,80?0=>�?::�>8,77��,..0;?,9.0�=,?0>�?::�3423�
�49,;;=:;=4,?0�;=4:=>��02��?=00�7092?3�
�����;,?3:7:240>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
�0.=0,>0�49�?480�>0=40>�;7:?�:1�79��0>?48,?0>��‘-@=9:@?’�
�:A0=�;,=,80?0=4E0/�8:/07
�;::=7D�.3:>09@9=0,74>?4.�;=4:=>
�:90�:=�8:=0�B0,6�;,=,80?0=>��8,=249,7�;:>?0=4:=>���;:>?0=4:=>�>4847,=�
���:9A0=209.0�/4,29:>?4.>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
(i) &480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�>:80�;,=,80?0=>�,=0�8:=0�=074,-70�?3,9�:?30=>
�>?0;>�8,D�:..@=�
�/4>.=0?0�;,=,80?0=>��02��.@8@7,?4A0�-4�;,=?4?4:9�1=0<@09.D����(&)
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
�C,8;70���(&)�;7:?>�:1�.@8@7,?4A0�-4�;,=?4?4:9�1=0<@09.D�:1���9:/0>
-,/�.:9A0=209.0 -0??0=�.:9A0=209.0
!E8-:01>�1@�-8�����
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��4C492�/4,29:>?4.>
(i) �:=8�:1�?30�?480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=B,=8�,9/�1@EED�.,?0=;477,=>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
-,/�84C492 -0??0=�84C492
�C,8;70��&=,.0=�;7:?>�:1�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
cynm
ix_m
b_rate_m
ultiplier.run2.p
State0 2000000 4000000 6000000 8000000 10000000
0.1
0.125
0.15
0.175
0.2
0.225
0.25
0.275
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
-,/�84C492 -0??0=�84C492
�C,8;70��&=,.0=�;7:?>�:1�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
cynm
ix_m
b_rate_m
ultiplier.run2.p
State0 2000000 4000000 6000000 8000000 10000000
0.1
0.125
0.15
0.175
0.2
0.225
0.25
0.275
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��4C492�/4,29:>?4.>
(i) �:=8�:1�?30�?480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�/4>.=0?0�;,=,80?0=>�
B,=8�,9/�1@EED�.,?0=;477,=>
�/4>?,9.0>�,8:92�>,8;70/�?:;:7:240>��&=00%0?'4E
A,G,',-1&#'(?,/D.5-./&.(U(A5..!."I1Y(
TreeSetViz: http://comet.lehman.cuny.edu/treeviz/ Mesquite: http://www.mesquiteproject.org/mesquite/mesquite.html
Hillis et al., Analysis and Visualization of Tree Space, Syst. Biol., 54(3): 471-482.
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
&=00%0?'4E
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
&=00%.,;0=
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��4C492�/4,29:>?4.>
(ii) �..0;?,9.0�=,?0>�:1�;,=,80?0=�@;/,?0>
�.:9?49@:@>���/4>.=0?0�;,=,80?0=>���=�,D0>�����%&��0?.=,?0>�>3:@7/�4/0,77D�1,77�49�?30�~20-70%�=,920
(i) �:=8�:1�?30�?480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�/4>.=0?0�;,=,80?0=>�
B,=8�,9/�1@EED�.,?0=;477,=>
�/4>?,9.0>�,8:92�>,8;70/�?:;:7:240>��&=00%0?'4E
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
Acceptance rates for the moves in the "cold" chain of run 1: With prob. Chain accepted changes to 13.61 % param. 1 (revmat) with Dirichlet proposal
.
.
. 0.04 % param. 34 (rate multiplier) Dirichlet proposal 6.59 % param. 35 (topology and branch lengths) TBR 14.06 % param. 35 (topology and branch lengths) LOCAL
Acceptance rates for the moves in the "cold" chain of run 1: With prob. Chain accepted changes to 33.30 % param. 1 (revmat) with Dirichlet proposal
.
.
. 19.13 % param. 34 (rate multiplier) Dirichlet proposal 17.40 % param. 35 (topology and branch lengths) TBR 29.76 % param. 35 (topology and branch lengths) LOCAL
-,/�84C492 -0??0=�84C492
�C,8;70��&=,.0=�;7:?>�:1�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��4C492�/4,29:>?4.>
�,..0;?,9.0�=,?0>�.,9�-0�.:9?=:770/�-D�A,=D492�?30�>.,70�:1�?30���?@9492�;,=,80?0=>�1:=�?30�=070A,9?�;=:;:>,7�80.3,94>8>
?:�49.=0,>0�=,?0>��/0.=0,>0�>.,70���A4.0�A0=>,
(ii) �..0;?,9.0�=,?0>�:1�;,=,80?0=�@;/,?0>
�.:9?49@:@>���/4>.=0?0�;,=,80?0=>���=�,D0>�����%&��0?.=,?0>�>3:@7/�4/0,77D�1,77�49�?30�~20-70%�=,920
(i) �:=8�:1�?30�?480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�/4>.=0?0�;,=,80?0=>�
B,=8�,9/�1@EED�.,?0=;477,=>
�/4>?,9.0>�,8:92�>,8;70/�?:;:7:240>��&=00%0?'4E
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��4C492�/4,29:>?4.>
(iii) �:=8�:1�?30�8,=249,7�;:>?0=4:=�;=:-,-474?D�/09>4?40>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
-0B,=0�:1�;:=.@;490�=:,/6477
�,..0;?,9.0�=,?0>�.,9�-0�.:9?=:770/�-D�A,=D492�?30�>.,70�:1�?30���?@9492�;,=,80?0=>�1:=�?30�=070A,9?�;=:;:>,7�80.3,94>8>
?:�49.=0,>0�=,?0>��/0.=0,>0�>.,70���A4.0�A0=>,
(ii) �..0;?,9.0�=,?0>�:1�;,=,80?0=�@;/,?0>
�.:9?49@:@>���/4>.=0?0�;,=,80?0=>���=�,D0>�����%&��0?.=,?0>�>3:@7/�4/0,77D�1,77�49�?30�~20-70%�=,920
(i) �:=8�:1�?30�?480�>0=40>�;7:?>�:1�;,=,80?0=�0>?48,?0>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�/4>.=0?0�;,=,80?0=>�
B,=8�,9/�1@EED�.,?0=;477,=>
�/4>?,9.0>�,8:92�>,8;70/�?:;:7:240>��&=00%0?'4E
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
Acceptance rates for the moves in the "cold" chain of run 1: With prob. Chain accepted changes to 13.61 % param. 1 (revmat) with Dirichlet proposal
.
.
. 0.04 % param. 34 (rate multiplier) Dirichlet proposal 6.59 % param. 35 (topology and branch lengths) TBR 14.06 % param. 35 (topology and branch lengths) LOCAL
Acceptance rates for the moves in the "cold" chain of run 1: With prob. Chain accepted changes to 33.30 % param. 1 (revmat) with Dirichlet proposal
.
.
. 19.13 % param. 34 (rate multiplier) Dirichlet proposal 17.40 % param. 35 (topology and branch lengths) TBR 29.76 % param. 35 (topology and branch lengths) LOCAL
Frequency
m{2}0.2 0.25 0.3 0.35 0.4 0.450
100
200
300
400
500
600
700
800
Frequency
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-4 9E-4 1E-30
100
200
300
400
500
-,/�84C492
�C,8;70��",=,80?0=�0>?48,?0>�1:=�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
-0??0=�84C492
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
Acceptance rates for the moves in the "cold" chain of run 1: With prob. Chain accepted changes to 13.61 % param. 1 (revmat) with Dirichlet proposal
.
.
. 0.04 % param. 34 (rate multiplier) Dirichlet proposal 6.59 % param. 35 (topology and branch lengths) TBR 14.06 % param. 35 (topology and branch lengths) LOCAL
Acceptance rates for the moves in the "cold" chain of run 1: With prob. Chain accepted changes to 33.30 % param. 1 (revmat) with Dirichlet proposal
.
.
. 19.13 % param. 34 (rate multiplier) Dirichlet proposal 17.40 % param. 35 (topology and branch lengths) TBR 29.76 % param. 35 (topology and branch lengths) LOCAL
Frequency
m{2}0.2 0.25 0.3 0.35 0.4 0.450
100
200
300
400
500
600
700
800
Frequency
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-4 9E-4 1E-30
100
200
300
400
500
-,/�84C492
�C,8;70��",=,80?0=�0>?48,?0>�1:=�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
-0??0=�84C492
(iv)��A@;/;>>18-@5;:�@591����&��;2�<->-91@1>�?-9<81?
⇢k =
Pn�ki=1 (xi � x̄)(xi+k � x̄)Pn
i=1(xi � x̄)2
)1�C;A80�1D<1/@�@41�kth�8-3�-A@;/;>>18-@5;:�@;�.1�?9-881>�-?�k�5:/>1-?1?��;A>��?@�-:0��@4�0>-C?�?4;A80�.1�81??�/;>>18-@10�@4-:�;A>��?@�-:0��:0�0>-C?�
�2�-A@;/;>>18-@5;:�5?�?@588�>18-@5B18E�4534�2;>�45341>�B-8A1?�;2�k��@45?�5:05/-@1?�4534�013>11�;2�/;>>18-@5;:�.1@C11:�;A>�0>-C?�-:0�?8;C�95D5:3
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��4C492�/4,29:>?4.>
�&41�8-3�k�-A@;/;>>18-@5;:�+>4;,7�5?�@41�/;>>18-@5;:�1B1>E�0>-C�-:0�5@?�kth�8-3�
0 5 10 15 20 25 30
0.0
0.2
0.4
0.6
0.8
1.0
Lag
ACF
Good Mixing
0 10 20 30 40
0.0
0.2
0.4
0.6
0.8
1.0
Lag
ACF
Bad Mixing
(iv)��A@;/;>>18-@5;:�@591����&��;2�<->-91@1>�?-9<81?
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��4C492�/4,29:>?4.>
-0??0=�84C492 -,/�84C492
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��%,8;70�>4E0�/4,29:>?4.>
(i) �110.?4A0�%,8;70�%4E0���%%��/4,29:>?4.
�9@8-0=�:1�>,8;70>,@?:.:==07,?4:9�?480����&�
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
�C,8;70���%%�A,7@0>�1:=�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
7:B�49?09>4?D
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
��%,8;70�>4E0�/4,29:>?4.>
(i) �110.?4A0�%,8;70�%4E0���%%��/4,29:>?4.
(ii) �:=8�:1�?30�8,=249,7�;:>?0=4:=�;=:-,-474?D�/09>4?40>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
-=:?30=�:1�;:=.@;490�=:,/6477
�9@8-0=�:1�>,8;70>,@?:.:==07,?4:9�?480����&�
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
09>@=0�%���.:8;74,9.0�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
�C,8;70��",=,80?0=�0>?48,?0>�1:=�80,9�=,?0�8@7?4;740=>�1=:8����%&�=@9>
7:B�49?09>4?D
Frequency
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-40
2.5
5
7.5
10
12.5
���.D.70>
-0??0=�49?09>4?D
Frequency
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-4 9E-40
10
20
30
40
50
60
70
���.D.70>
�49,/0<@,?0�.3,49�7092?3
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
�C,8;70��",=,80?0=�0>?48,?0>�1:=�80,9�=,?0�8@7?4;740=>�1=:8����%&�=@9>
7:B�49?09>4?D
Frequency
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-40
2.5
5
7.5
10
12.5
���.D.70>
-0??0=�49?09>4?D
Frequency
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-4 9E-40
25
50
75
100
125
����.D.70>
�49,/0<@,?0�.3,49�7092?3
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
Frequency
m{2}0.2 0.25 0.3 0.35 0.4 0.450
100
200
300
400
500
600
700
800
Frequency
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-4 9E-4 1E-30
100
200
300
400
500
�C,8;70��",=,80?0=�0>?48,?0>�1:=�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
7:B�49?09>4?D -0??0=�49?09>4?D
���.D.70> ����.D.70>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
alpha p
�����;,?3:7:240>
",=,80?0=�49?0=,.?4:9�-0?B009�I+G�84C?@=0�1:=�,8:92�>4?0�=,?0�A,=4,?4:9
Frequency
mtDNA2.alpha0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.60
100
200
300
400
500
Frequency
mtDNA2.pInv0.1 0.15 0.2 0.25 0.3 0.35 0.40
100
200
300
400
500
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�%49270�$@9>
�8@7?4�8:/,7�8,=249,7�/09>4?40>�>?08�1=:8�9:9�4/09?414,-474?D�
�@>0�G B4?3�,//4?4:9,7�/4>.=0?0�=,?0�.,?02:=40>�
alpha p
�����;,?3:7:240>
",=,80?0=�49?0=,.?4:9�-0?B009�I+G�84C?@=0�1:=�,8:92�>4?0�=,?0�A,=4,?4:9
Frequency
mtDNA2.alpha0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.60
100
200
300
400
500
Frequency
mtDNA2.pInv0.1 0.15 0.2 0.25 0.3 0.35 0.40
100
200
300
400
500
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3��0?=:;:74>��:@;70>�����
$;.;@�%=A-0>;:��
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3��0?=:;:74>��:@;70>�����
$;.;@�%=A-0>;:��
�!���������������������������������������������$�������� ������������������ �������
���������������#�����!���������������������������������$�������$��� �����������������������$�����������������$����� �
N����������������������������������������������������������������������������������$�������$�
����������������������N–1 ����������������$��������
���������������"����������������������
��������������������������������������������������$�������$��� �����������������������������$��������!�����������������������������$�������� �����
����������$�����"�������������������"�����������������������������������$���������������������������"�������� ������������������������"�������������������������������������������$�
������������� i = 1/(1 + iT)����� 0.25 0.20 0.15 0.101 1.00 1.00 1.00 1.002 0.80 0.83 0.87 0.913 0.66 0.71 0.77 0.834 0.57 0.62 0.69 0.77
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3��0?=:;:74>��:@;70>�����
��?8534@8E�9;>1�2;>9-8�01?/>5<@5;:
�5-3:;?5:3� ���<1>2;>9-:/1
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3��0?=:;:74>��:@;70>�����
&41�<>59->E�05-3:;?@5/�5?�@41�-//1<@-:/1�>-@1?�2;>�<>;<;?10�/4-5:�?C-<?�
�52�-//1<@-:/1�>-@1?�->1�@;;�8;C��01/>1-?1�@41�B-8A1�;2�@41�@19<1>-@A>1�<->-91@1>�52�-//1<@-:/1�>-@1?�->1�@;;�4534��5:/>1-?1�@41�B-8A1�;2�@41�@19<1>-@A>1�<->-91@1>
�?�-�>A81�;2�@4A9.��-//1<@-:/1�>-@1?�2;>�<>;<;?-8?�?4;A80�2-88�5:�~20–70%�>-:31
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3��0?=:;:74>��:@;70>�����
-,/�84C492 -0??0=�84C492
�C,8;70��&=,.0=�;7:?>�:1�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
�5-3:;?5:3� ���<1>2;>9-:/1
�;;=:C48,?492�?30��:49?�":>?0=4:=�"=:-,-474?D��09>4?D�B4?3��0?=:;:74>��:@;70>�����
&41�<>59->E�05-3:;?@5/�5?�@41�-//1<@-:/1�>-@1?�2;>�<>;<;?10�/4-5:�?C-<?�
�52�-//1<@-:/1�>-@1?�->1�@;;�8;C��01/>1-?1�@41�B-8A1�;2�@41�@19<1>-@A>1�<->-91@1>�52�-//1<@-:/1�>-@1?�->1�@;;�4534��5:/>1-?1�@41�B-8A1�;2�@41�@19<1>-@A>1�<->-91@1>
"@41>�-?<1/@?�/;:@>;885:3�@41�.14-B5;>�;2�@41� 1@>;<;85?�/;A<85:3�/-:�.1�9;052510�@;����59<>;B1� � ��<1>2;>9-:/1���
�5:/>1-?1�@41�:A9.1>�;2�5:/>191:@-88E�41-@10�/4-5:?��13��nchains�<->-91@1>��5:/>1-?1�@41�2>1=A1:/E�;2�-@@19<@10�/4-5:�?C-<�1B1:@?��13��swapfreq�<->-91@1>��5:/>1-?1�@41�:A9.1>�;2�?C-<?�-@@19<@10�<1>�1B1:@��13��nswaps�<->-91@1>�
�?�-�>A81�;2�@4A9.��-//1<@-:/1�>-@1?�2;>�<>;<;?-8?�?4;A80�2-88�5:�~20–70%�>-:31
!@?7490
����5-3:;?5:3� � ��<1>2;>9-:/1
����>1B51C�;2�C41>1�C1’B1�.11:�-:0�C4E
�����5-3:;?@5/?�.-?10�;:�?5:381�/4-5:?
�(��5-3:;?@5/?�.-?10�;:�@41�<>5;>�
(��5-3:;?@5/?�.-?10�;:�9A8@5<81�/4-5:?�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�?30�"=4:=
�>?48,?492�@9/0=�?30�;=4:=
�,=249,7�;:>?0=4:=�/09>4?40>�1:=�;,=,80?0=>�,=0�@;/,?0/�A0=>4:9>�:1�?30����.:==0>;:9/492�;=4:=�;=:-,-474?D�/09>4?40>��?30D�,=0�@;/,?0/�-D�?30�����491:=8,?4:9�49�?30�/,?,�A4,�?30�7460743::/�1@9.?4:9
Pr[τ i | X] =Pr[X |τ i] × Pr[τ i]Pr[X |τ j ] × Pr[τ j ]j=1
B( s )∑!
;:>?0=4:=�;=:-,-474?D7460743::/�1@9.?4:9
8,=249,7�7460743::/
;=4:=�;=:-,-474?D
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9�?30�"=4:=
������.,9�-0�=@9�?:�?,=20?�?30�5:49?�;=4:=�04?30=�-D�0>?48,?492�B4?3�9:����/,?,�:=�-D�1:=.492�?30�7460743::/�1@9.?4:9�=0?@=9��
�,=249,7�;:>?0=4:=�/09>4?40>�1:=�;,=,80?0=>�,=0�@;/,?0/�A0=>4:9>�:1�?30����.:==0>;:9/492�;=4:=�;=:-,-474?D�/09>4?40>��?30D�,=0�@;/,?0/�-D�?30�����491:=8,?4:9�49�?30�/,?,�A4,�?30�7460743::/�1@9.?4:9
R =min 1, f X |Θ'( )f X |Θ( )
⋅f Θ'( )f Θ( )
⋅f Θ |Θ'( )f Θ' |Θ( )
$
% &
'
( ) !
8571854;;0�>-@5;� <>5;>�>-@5;� <>;<;?-8�>-@5;�
(0�.,9�.:8;,=0�?30�8,=249,7�;=4:=�/09>4?40>�?:�?304=�;:>?0=4:=�.:@9?0=;,=?>����?:�307;�4/09?41D�B0,6�;,=,80?0=>
�>?48,?492�@9/0=�?30�;=4:=
!@?7490
����5-3:;?5:3� � ��<1>2;>9-:/1
����>1B51C�;2�C41>1�C1’B1�.11:�-:0�C4E
�����5-3:;?@5/?�.-?10�;:�?5:381�/4-5:?
�(��5-3:;?@5/?�.-?10�;:�@41�<>5;>�
(��5-3:;?@5/?�.-?10�;:�9A8@5<81�/4-5:?�
&30�2090=,7�4/0,�4>�?:�.:8;,=0�0>?48,?0>�1=:8�8@7?4;70�49/0;09/09?�.3,49>����494?4,?0/�1=:8�=,9/:8�;,=,80?0=�A,7@0>�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�:=8�:1�?30�8,=249,7�;:>?0=4:=�/09>4?40>�1:=�,77�;,=,80?0=>
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�C,8;70��&=,.0=�;7:?>�:1�8,=249,7�/09>4?40>�1=:8�8@7?4;70��=�,D0>�=@9>
-,/�.:9A0=209.0 -0??0=�.:9A0=209.0
Alpha[7]
0.5 1 1.5 2 2.5 3 3.50
0.5
1
1.5
2
2.5
3
meanRate2E-4 3E-4 4E-4 5E-4 6E-4 7E-4 8E-4 9E-4 1E-30
1000
2000
3000
4000
5000
6000
�&=,.0=�/08:
&30�2090=,7�4/0,�4>�?:�.:8;,=0�0>?48,?0>�1=:8�8@7?4;70�49/0;09/09?�.3,49>����494?4,?0/�1=:8�=,9/:8�;,=,80?0=�A,7@0>�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�:=8�:1�?30�8,=249,7�;:>?0=4:=�/09>4?40>�1:=�,77�;,=,80?0=
�.:9?49@:@>�;,=,80?0=>�
�"%$��/4,29:>?4.���=�,D0>
1. Run m ≥ 2 chains of length 2n from overdispersed starting values.
3. Calculate the within-chain and between-chain variance.
4. Calculate the estimated variance of the parameter as a weighted sum of the within-chain and between-chain variance.
2. Discard the first n draws of each chain.
5. Calculate the PSRF.
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
-,/�.:9A0=209.0
-0??0=�.:9A0=209.0
�C,8;70��"%$��A,7@0>�1:=�=07,?4A0�=,?0�8@7?4;740=>�1=:8�?B:��=�,D0>�=@9>
95% Cred. Interval ---------------------- Parameter Mean Variance Lower Upper Median PSRF * ------------------------------------------------------------------------------------------- TL{all} 4.921609 2.998138 2.836000 7.295000 5.056000 9.084 kappa{4,5} 3.095696 0.054125 2.667623 3.587024 3.085271 1.000 alpha{5} 1.006544 0.087721 0.606472 1.738482 0.950093 1.000 pinvar{1} 0.307396 0.009357 0.095913 0.471070 0.316173 1.000 m{1} 0.264226 0.009315 0.146502 0.421870 0.244468 5.507 m{2} 0.040919 0.000227 0.022205 0.065884 0.037425 5.279 m{3} 2.721453 7.157157 0.039001 5.544253 5.030560 69.564 m{4} 2.125810 3.568002 0.199137 4.044249 3.917338 150.012 m{5} 0.188768 0.004373 0.109303 0.295129 0.170624 5.749 -------------------------------------------------------------------------------------------
95% Cred. Interval ---------------------- Parameter Mean Variance Lower Upper Median PSRF * ------------------------------------------------------------------------------------------- TL{all} 0.073893 0.000034 0.063000 0.086000 0.074000 1.000 kappa{2,3} 3.236308 0.366904 2.199024 4.587719 3.190195 1.000 m{1} 1.285838 0.028345 0.980634 1.630387 1.278161 1.000 m{2} 1.423906 0.015507 1.182596 1.664627 1.423610 1.000 m{3} 0.589346 0.005341 0.453175 0.736459 0.587617 1.001 -------------------------------------------------------------------------------------------
&30�2090=,7�4/0,�4>�?:�.:8;,=0�0>?48,?0>�1=:8�8@7?4;70�49/0;09/09?�.3,49>����494?4,?0/�1=:8�=,9/:8�;,=,80?0=�A,7@0>�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�:=8�:1�?30�8,=249,7�;:>?0=4:=�/09>4?40>�1:=�,77�;,=,80?0=
�.:9?49@:@>�;,=,80?0=>�
�/4>.=0?0�;,=,80?0=>�
�&:;:7:2D
�>4847,=4?D�:1�;,4=0/�.3,49>��02���%�%��/4,29:>?4.�49��=�,D0>�
�"%$��/4,29:>?4.���=�,D0>
�>4847,=4?D�:1�8,=249,7�/09>4?40>��&=,.0=
�����;,?3:7:240>
&480�>0=40>�;7:?�:1��%�%��/4,29:>?4.
�>?:;�>,8;7492�B309��%�%��< 0.01
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
0.01
ASDSF
�@9/0=�;,=,80?0=4E0/�8:/07
�����;,?3:7:240>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
0.01
ASDSF
�:A0=�;,=,80?0=4E0/�8:/07
&480�>0=40>�;7:?�:1��%�%��/4,29:>?4.
��%�%��:>.477,?0�,-:@?�?30 0.01�?3=0>3:7/�A,7@0
�;::=7D�.3:>09�;=4:=>
&30�2090=,7�4/0,�4>�?:�.:8;,=0�0>?48,?0>�1=:8�8@7?4;70�49/0;09/09?�.3,49>����494?4,?0/�1=:8�=,9/:8�;,=,80?0=�A,7@0>�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�:=8�:1�?30�8,=249,7�;:>?0=4:=�/09>4?40>�1:=�,77�;,=,80?0=
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=�/4>.=0?0�;,=,80?0=>�
�&:;:7:2D
�>4847,=4?D�:1�;,4=0/�.3,49>��02���%�%��/4,29:>?4.�49��=�,D0>�
�/4>?,9.0>�,8:92�>,8;70/�?:;:7:240>��&=00%0?'4E&=00%.,;0=
�>;74?�1=0<@09.40>���;=0>09.0,->09.0���(&)
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�C,8;70��>;74?�1=0<@09.40>���;=0>09.0,->09.0�49��(&)
-0??0=�.:9A0=209.0-,/�.:9A0=209.0
&30�2090=,7�4/0,�4>�?:�.:8;,=0�0>?48,?0>�1=:8�8@7?4;70�49/0;09/09?�.3,49>����494?4,?0/�1=:8�=,9/:8�;,=,80?0=�A,7@0>�
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�:=8�:1�?30�8,=249,7�;:>?0=4:=�/09>4?40>�1:=�,77�;,=,80?0=
�.:9?49@:@>�;,=,80?0=>��02��>@->?4?@?4:9�=,?0>���&=,.0=�/4>.=0?0�;,=,80?0=>�
�&:;:7:2D
�>4847,=4?D�:1�;,4=0/�.3,49>��02���%�%��/4,29:>?4.�49��=�,D0>�
�/4>?,9.0>�,8:92�>,8;70/�?:;:7:240>��&=00%0?'4E&=00%.,;0=
�>;74?�1=0<@09.40>���;=0>09.0,->09.0���(&)�"�
�9:/,7�>@;;:=?���(&)�=�,D0>
�>>0>>492������"0=1:=8,9.0��4,29:>?4.>��,>0/�:9��@7?4;70�$@9>
�C,8;70��‘comparetrees’�;7:?�:1�?=00>�>,8;70/�-D�?B:��=�,D0>�=@9>
!E8-:01>�1@�-8�����
-0??0=�.:9A0=209.0-,/�.:9A0=209.0
You can never be absolutely certain that the MCMC is reliable, you can only identify when something has gone wrong. Gelman
%@88,=D��%:80��090=,7�%?=,?0240>�1:=��>>0>>492������"0=1:=8,9.0
��)41:�0;�E;A�:110�@;�-??1??� � ��<1>2;>9-:/1���(�)%
��)41:�?4;A80�E;A�-??1??�@41�<1>2;>9-:/1�;2�5:05B50A-8�>A:?����(�)%
�)45/4�05-3:;?@5/?�?4;A80�E;A�A?1�@;�-??1??�5:05B50A-8�>A:?�����?3,?�,=0�=070A,9?�1:=�?30�8:/07>;,=,80?0=>�D:@�,=0�0>?48,?492�@9/0=
��)41:�5?�-�?5:381�>A:�?A225/51:@�@;�-??1??� � ��<1>2;>9-:/1� �'�$
��)41:�?4;A80�E;A�1?@59-@1�A:01>�@41�<>5;>�(�� �'�$�"!%%������,9/�-0�B,=D�:1�;=:2=,8>�B30=0�4?�4>�9:?�;:>>4-70�
You can never be absolutely certain that the MCMC is reliable, you can only identify when something has gone wrong. Gelman
%@88,=D��%:80��090=,7�%?=,?0240>�1:=��>>0>>492������"0=1:=8,9.0
��)41:�?4;A80�E;A�<1>2;>9�9A8@5<81�5:01<1:01:@� � ��>A:?���(�)%��,9/�-0�B,=D�:1�;>0@/:�49/0;09/09.0�
��)45/4�05-3:;?@5/?�?4;A80�E;A�A?1�@;�-??1??�5:05B50A-8�>A:?�����?3,?�,=0�=070A,9?�1:=�?30�8:/07>;,=,80?0=>�D:@�,=0�0>?48,?492�@9/0=
���;C�9-:E�5:01<1:01:@� � ��>A:?�->1�?A225/51:@��%��� )��%�"!%%������40��,>�8,9D�,>�D:@�?3496�D:@=�/,?,;=:-708�/0>0=A0�
���;C�8;:3�?4;A80�E;A�>A:�1-/4� � ��-:-8E?5?��%��! ���%�"!%%������40��,>�7:92�,>�D:@�?3496�D:@=�/,?,;=:-708�/0>0=A0�
��)41:�?4;A80�E;A�A?1� 1@>;<;85?��;A<85:3�(3090A0=�D:@�.,99:?�-0�.0=?,49�?3,?�>?,9/,=/������4>�,/0<@,?0
40����(�)%��,9/�-0�B,=D�:1�;=:2=,8>�B30=0�4?�4>�9:?�;:>>4-70�
&=,.0=
�>>0>>492������"0=1:=8,9.0�%:1?B,=0�&::7>
�(&)���"��
&=00%0?'4E
&=00%.,;0=
�!�
.:/,