lambda ~ gamma(response*d,d) abundance ~ zip(pi,lambda)response temp

abundance

l…+.P\?;%("EF n! "#$"&%('*),+.-0/1*2("43ƒpj?;8 89':3;"&<,F

! "y8T+VUD2?;%("–%("46W'7-6@"–6+.2?1*2+='7-A>$1*- e "–>@

1*%@17>&2("&%+.#$"43Xo1*-1*<='7PR'7?6(<.‡‚176]+.-BD"4>&2+='7-hEFHG0o e ‡T+.-d2%('D30?>&+.-;P91–<J1*2("&-:2Z%@1*-3;'78^\1*%+J1 e <=" —,-7’4`4`4`b’/.106?>@

2 1*2s2 "n%("46'7-6("89'D3;"&<32w+=6S6("&<="4>&2("43y+VK L{42F5g"n3;"&-'72(" e ‡ % L76–{"2 ” {

98;:6;%(' e 1 e

+.<.+.2z‡‚'*Ks2 "t%("46'7-6("t89':3;"&<2KM'7%!2 "w2@1bU;'7- <QKY'7%1*<.<)=L{>-7’4`4`4`b’@?no

<Q{1-7’4`4`4` ’@A 1*-3B2€{>-7’4`4`4` ’/.F„O>4>4'7%(30+.-;P\<.‡Ro

C

ED76{ Ž\76R’L76–{"2W6:’;F  :6:’4`4`4`$’;F a :6 ” { C 8 ED76{ Ž\760L76{42\’;FG8;:6 ” 8/:6R`

NO "&%("HFG8/:6C3;"&-'72("46 1*%@1*89"&2("&%(6 '*K€2 "5%("46'7-6(" 89'D3;"&<I2k'*K€2 " 2@1bU;'7-<o

C 8 3;"&-'72("46C2 " <.+.[R"&<.+. ':'D3 KY?;->&2+='7-'*Kh2 " %("46'7-6("~89'D3;"&<J2 1*-3 6{

 8;: ’4`4`4`$’ a :K ” r†%("4>$1*<.<ML

a

NPO

Q

N:6SR>- vuF…l;?;%2 "&%$o72 "Z2('72@1*<d8T+VU02?;%("m<.+.[R"&<.+. 'D'D3

'*K…2 "–1 e ?;-31*->4"–'*KS2 "2@1bU;'7-T<1*2O6+.2("U=n+=6!2 "&-™' e 2@1*+.-"43 e ‡I8Q1*%P\+.-1*<.+.#4+.-;P

C

EDV6{ Ž\760W  ’ Q

6:’;F  :6d’4`4`4`$’;F a :6 ” { r†EF$ŠRv

a

X8 O  C 8

ED76{ Ž\76DL76–{"2R’@W



’;FG8/:6 ” Q

8;:6d` r†EFRv

' ^R"&%92 "h<J1*2("&-d29%@1*-3;'78 ^71*%+J1 e <=" V6 ! d?6$oZ2 "h%("&P\%("46(6+='7- K†?;->&2+='7- '*Kt2 "

8T+VU02?;%("%("46'7-6("89'D3;"&<x+=6

YZED769{ Ž\76DW  ’ Q

6d’;F  :6R’4`4`4`b’;F a :6 ” { a X8 O

[Z

8 \W ^] F_8/:6 ” Q 8/:6R`

n%(' e 1 e +.<.+.2+="46

Q

61*%("89'D3;"&<="43‚?6+.-;PT2 "a`w+.%+=>@ ;<="&2w30+=62%+e

?;2+='7-Q  :6d’4`4`4` ’ Q

a bX :60c ’4`4`4` ’/ca ”ed `t+.%+=>( ;<="&27fc$’4`4`4`b’/ca ” ’ r†EF v

NO "&%("gc

8

+=6L1 :‡ "&%1*%@1*89"&2("&%L%(" ;%("46("&-d2+.-;PBhijlkfmnj_kR"&^:+=3;"&->4"ZKM'7%S2 "Z%("46'7-6("

89'D3;"&<n2Foc  >$1*- e "m+.-d2("&%;%("&2("439176phij_kfmnjlkR6 1*8;<="Z6+.#$"nKY'7%…2 "Z%("46W'7-6@"n89'D3;"&<

2Cˆrqw"4>$Š\‹bŒiF„ -'7-0),+.-0KM'7%8Q1*2+.^R"Q>( '7+=>4"€+=62('™1 ;<.‡tsR"_u%("&‡v‰6;%+='7%ac xw{y.

bX

KM'7%™1*<.<z-|{}=~{. ˆ‰pB0Š*Eox€‚qUƒwŠ\‹bŒiF…„w- 2 "ƒ'72 "&%‚ 1*-3 176@6?;89"ƒ2 1*2

2 "&%("–1*%("2 ;%("4">4'78"&2+.-;PT%("46'7-6("t89':3;"&<=6$o02 "t<.+.-"$1*%O%("46'7-6("t89'D3;"&<mr

Z

uvuo

2 "†€1*?6(6+J1*-%("46'7-6("I89':3;"&<r

Z(‡

vuo1*-32 "Q#$"&%(' %("46'7-6("Q89'D3;"&<,on1*-3 2 "

1*<="4'D"4>4'7<='7P\+=62 176xhgijlkfmnj_kd"&^D+=3;"&->4"Z2 1*2

Z ‡

+=6L89'7%(" ;%(' e 1 e <="!2

1*-Z

$o7NO ;+=>( xo

+.-2?;%-xoD+=6Z89'7%(" ;%(' e 1 e <="t2

1*-Z(ˆ

F! "&-6( "w6 '7?;<=3‚173;' ;2O6('789"t>( '7+=>4"6?>@

2 1*2c

‡Š‰

c

‰ c ˆ F

‹aŒŽ 

]] o‘t

@’

(w w}(

;“

z-–2 "mpj?;8T89"&%s89':3;"&<:89':3;"&<="43%("46;”W'7-6("46

Z

\W Ž] F•6 ” KY'7%S2 "m17>&2?1*<D1 e ?;-31*->4"46

—–

KY'7%1 e ?;-31*->4"46–'*Km2@1bU1h1*26(+.2("J=ˆ€‚qUƒwŠ\‹0o)qt"4>$Š\‹bŒio)”;%('$^D+=3;"4352 1*2–1Q89'D3;"&<

–w{

Z

\W ^] F•6 ” +=6…173;'n”;2("43XFs']2 ;+=6L"&-3Xo7":?;+.^71*<="&-:2L6 1*8 ”;<="Z6+.#$"Z 176 e "4"&-?6("43

176–1hp!1$‡R"46+J1*- h ijlkfmnj_kS%("&P\?;<J1*%+.# 1*2+='7- 2(':'7<L+.-H`t+.%+=>( ;<="&2–8?;<.2+.-'78T+J1*<m89'D3;"&<=6

KY'7%>4'7?;-d2–31*2@10oX"\FPF.oW-1*+.^R"€pO14‡R"46>&<J176(6+"&%(6Tˆ w'7-ŠoWl€‚€Š *Œio1*-3ƒ+.-53;"&2("&% )

8T+.-;+.-;P1pO14‡R"46+J1*-Q-"&2iNZ'7%[I62%?>&2?;%("\oR+,F‰"\F.oD16("&2n'*KX>4'7-30+.2+='7-1*<+.-3;"G”"&-3;"&->4"

176(6?;8 ”;2+='7-6$oKM'7%>4'7?;-d2y31*2@1?;-3;"&%y2 "9176(6?;8 ”;2+='7-g'*K`t+.%+=>@ ;<="&2–8T+VUD2?;%("Q'*K

8y?;<.2+.-'78T+J1*<=6ˆrqw"4>$Š\‹;oB9s *ŒiF

„]-™":?;+.^\1*<="&-d2]6@1*8 ”;<="y6+.#$"

Z

\W ] F  ” • Z

\W ] F ‡ ” •*•

Z

\W ] F K ”

+=6€1 %@1*-3;'78 ^\1*%+J1 e <="+.-2 "hpj?;8T89"&%T89'D3;"&<,on6+.->4"‚1*2€<="$1762 FR4’4`4`4` ’;F K 1*%("

%@1*-3;'78 ^\1*%+J1 e <="46w1*-3‚+.2]>$1*- e "+.-:2("&%@”;%("&2("43ƒ2 1*22 "&‡‚3;"&2("&%8T+.-"y ' N 8y?>@

NZ"T2%?621;U0"43 6("&2–'*KO%("&P\%("46(6+='7- KY?;->&2+='7-63;"&2("&%8T+.-"43 e ‡ 2 "z”1*%@1*89"&2("&%(6

FR4’4`4`4`$’;F K 1*-3C1*-g"&-:^:+.%('7-;89"&-:2JWnFT! :?6$os1*-*hij_kfmnjlk…3;"&2("&%8T+.-1*2+='7-'*K!2 "

d?1*-:2+.2i‡

Z ” w {

: , Z

\W

] F ” w

‰ +.2] '7<=3;6O2 1*2

C

‚EF;’

u””

‰ 0R’ r†EFRv

NO "&%("‚EF;’ u” +=61*- 'n”"&- te 1*<.<LNO+.2

~F'7%+.P\+.-xoX+=6"46(6("&-d2+J1*<,oX6(+.->4"+.2tP\+.^R"46–1*-?M”M”W"&% e '7?;-3CKM'7%T1*-~":?;+.^\1*<="&-d296@1*8 ”;<="‚6+.#$"\F l'7%T+.-62@1*->4"\onNO "&-~+.2T>4'789"46

2(' 1‚Pd1*?6@6+J1*- %("46;”W'7-6("\o



Z "!"#%$&$

” +=6"d?1*<Z2('A2 "I6?M”;%("&8y?;8 '*K2 "6(>$1*<.+.-;P

KM17>&2('7%g”;%(' ^:+=3;"43‚2 1*2O2 ">4'7%%("46;”W'7-30+.-;PQ6(?M”;%("&8y?;8 176ahŠij_kEm+jlkW30+=62%+e

?;2+='7-6?>( ‚2 1*2OKM'7%1*<.<x+.2(6 )i"&-:^:+.%('7-;89"&-:2(6]2 "&%("–+=61-'7-0),#$"&%(' hŠijlkfmnj_k”;%(' e 1 e +.<.+.2z‡Ro

NO ;+=>( xo;+.-™2?;%-xo8y?62

e

"6@”"4>&+"43 hŠijlkfmnj_k†F

i-'7%(3;"&%€2(' 6("4"I '$N2 "':?1*-d2+.2i‡

Z

”–e '7?;-3;6^\1*%+J1*->4"h'*K2 " ”;%('n”'7% )

2+='7-6 6?M”M”'\6("\oKM'7% 6(+.8 ”;<.+=>&+.2i‡Ro2 1*2 2 " %("46@”'7-6(" KY?;->&2+='7- +=6ƒ'*K92 " KM'7%8

Z

\W ] F•6R’

” { Z “

\W ] F ” onNO "&%("

+=6T1™P\<=' e 1*<6@>$1*<.+.-;P K,17>&2('7%€NO '\6("-'7-0),#$"&%('

”;%+='7%”;%(' e 1 e +.<.+.2+="46–1*%("y+.- 2 "6("&2 ’)( ” Fw!

P\+.^R"461*- ?M”M”"&% e '7?;-3‚KM'7%

2 "^71*%+J1*->4"–'*KL89'D3;"&<="43†”;%('n”W'7%2+='7-6176!KM'7<.<='$N6+*

,

C

—– Fd$’4`4`4`$’;F K ’@W i” { Z “

\W  ] F36 ”  Z

“.-/

Z “

\W  ] F•6 ””

‡ Z “‡

0- Z

“0-• - ” { r†EF1Rv

Z

“0- Z

“.-/ Z “

\W Ž] F•6 ””

‡ Z “‡

0- Z

“.-• - ” { r†EF*E:v

‡ Z “‡

‡ Z “‡

0- Z

“0-• - ” {

Z

“0-32 176

2 ( ’ r†EF\‹Rv

4&576%8981:<;>=@?.ACBEDGFHFGFGDIA1J KCBMLON PRQSIQTMU9V=HW"?.X1BGDGFGFHFGD>X7JYL&Z [\U%=G]_^)`"S?.AYaMbX1BEDGFGFHFGD>X7JYLOc

dGe>fdg

K

dEeih

dGjg fdglk

B

h7m7n U9=HSI=oX

k

cqp

a

X%a\rsotvuRwyx+z{iZ

NO "&%("

Z

“0-{ L KO  Z “

\W ] F i” 1‚6(?;8 '*KO2 "T%("46@”'7-6("46y1*2y6+.2("z= e "uKY'7%("96@>$1*<.+.-;P

?6+.-;P

FmB:+.8T+.<J1*%<.‡Ro2 "Š”W'\6(2("&%+='7%^71*%+J1*->4"–+=6 e '7?;-3;"43 e ‡

,

C

—– Fd$’4`4`4`$’;F K ’@W  ’ i” {

Z

“.-•~Ž

.-’ r†EF\GRv

NO "&%("  —

K

+=6m1^R"4>&2('7%m'*Kx1 e ?;-31*->4"46!1*2m6+.2("p=…1*-39Ž .- { L K6 O  Ž —– o:6+.->4"

2 " `w+.%+=>@ ;<="&2…30+=6(2%+e ?;2+='7-y+=6s1>4'7- ?;Pd1*2(" ”;%+='7%s'*K;2 "n8?;<.2+.-'78T+J1*<D30+=62%+e

?;2+='7-ˆ€‚qUƒwŠ\‹bŒiF

iK%(":?;+.%("43XoR1*-":?;+.^71*<="&-:2L6 1*8 ”;<="Z6+.#$"Z>$1*- e "m>4'7-:2%('7<.<="4391*-3y89'D3;"&<="436("G”0)

1*%@1*2("&<.‡KM'7%"$17>( 6(+.2(" e ‡+.-d2%('D30?>&+.-;P"&+.2 "&%w1€6(2(':>@ 1762+=>–'7%>4'7-62@1*-:2]^\1*%+J1 e <="

 2 1*2+=66+.2("u)i6;”W"4>&+>–1*-3™3;"&2("&%8T+.-"461*-™":?;+.^71*<="&-:2w6@1*8 ”;<="6+.#$"176!KY'7<.<=' N6

C B

’4`4`4`$’

C



KCB

W  ’;Fd$’4`4`4`$’;F K”ed `t+.%+=>@ ;<="&2b 7–

 ’4`4`4`$’ 7–

i” ’ r†EF \v

NO "&%("

—–

{ Z

\W  ] F•6 ”

L

6

Z

\W Ž] F•6” { YZ

C

—–

FR$’4`4`4` ’;F36:’@W

i”

+=62 ""uU9”"4>&2("43 ”;%('n”'7%2+='7-f'*K2@1bU0'7-4< 1*26+.2(" =@Fƒ!

 +=6€1*-~":?;+.^71*<="&-:2

6@1*8 ”;<="6+.#$"]KM'7%O16+.2("=@F

 8Q1$‡ e "t1%@1*-3;'78^71*%+J1 e <="1*-3+.2!8Q14‡ e "t?6("43h176

189'D3;"&<.+.-;P‚2(':'7<L2('h89'D3;"&<^71*%+J1*2+='7- e "&2iNj"4"&-C30+VuX"&%("&-d2y6("430+.89"&-:2–<J14‡R"&%(6–2 1*2

8T+.P\ :2 e " ”;%("46("&-d2$o!"\FPF.oj30?"h2('g30+Vu"&%+.-;P 6("430+.89"&-d2@1*2+='7- %@1*2("46$F4qt'$NZ"&^R"&%$oj+VK

2 "&%("]+=6m-'%("$176@'7-QNO d‡



6 '7?;<=3 e "]30+Vu"&%("&-:2mKY%('78 2 1*2j'*K

8

r\={"2Dvn1–P\<=' e 1*<

":?;+.^71*<="&-:2j6 1*8 ”;<="w6+.#$"w6 '7?;<=3 e "w173;'n”;2("43XFo`t?"]2('–2 "w+.-"d?1*<.+.2z‡ WFEF\‹0o



e '7?;-3;62 "^\1*%+J1*->4"46]'*KL2 "–'D>4>&?;%("&->4"Š”;%(' e 1 e +.<.+.2+="46t176

,

C

—– FR$’4`4`4`$’;F K ’@W i” {

 `

j'7-6(":?"&-d2<.‡Ro…2 "T^\1*%+J1*->4"9'*K

C

—– %("430?>4"462('‚#$"&%('A176–2 " ”;%+='7%6@1*8 ”;<="Q6+.#$"

 1*2]6+.2("S=n+=6!+.->&%("$176("43™NO+.2 '7?;2<.+.8T+.2(6$o;+,F‰"\F.o

 2 ( -'€8Q1*22("&%NO 1*2t1*%("2 "

”;%+='7%€6@”"4>&+>$1*2+='7-6yKM'7%y2 "Q2@1bU0'7-~6@”"4>&+>z”1*%@1*89"&2("&%(6Fd&’4`4`4`b’;F

K

1*-3"&-d^D+V)

%('7-;89"&-:2xW  F '7%("t+.8 ”'7%2@1*-d2<.‡Ro1761*-A1•”M”;<.+=>$1*2+='7- '*KeZ " e ‡D6("&^)v‰6!+.-":?1*<.+.2i‡

Nj"' e 6("&%^R"2 1*2t19%@1*-3;'78 ^\1*%+J1 e <="

C

—– P\+.^R"&-F•6€1*-3TW  %("430?>4"46t2('9+.2(6”;%+='7%

"uU[”W"4>&2@1*2+='7-9176…2 "”;%+='7%6@1*8 ”;<="!6+.#$"m+=6…+.->&%("$176("43NO+.2 '7?;2<.+.8T+.2(6$F z2…6 '7?;<=3 e "

-'72("43XoS '$NZ"&^R"&%$os2 1*2–2 ":?1*<.+.2@1*2+.^R"962%?>&2?;%("9'*Kj2 "T89'D30+>$1*2+='7- WFEF

2(' 2 "Apj?;8T89"&%h89':3;"&<w%("46?;<.2(6+.- 1C30+Vu"&%("&-:2-"&2zNZ'7%[ 62%?>&2?;%("A2 1*2 ”'\6 )

6("46@6€89'7%("‚>4'78 ”;<="uUf3;"G”"&-3;"&->&+="46 e "&2iNj"4"&-k2 "h^\1*%+J1 e <="46T2 1*-f?6+.-;P 89'D3;"&<

WFEF$G0F

B:+.2("+=3 B”W"4>&+="46]„ B”"4>&+="46p „]+.%O2("&8 ”"&%@1*2?;%(" S%+J1*<=6

B:+.2(" Š

B:+.2("

B:+.2(" Š

l…+.P\?;%("EF n„w-™"uU;1*8 ”;<="y31*2@1

x"&2!?6O>4'7-6+=3;"&%]16+.8 ”;<=""uU;1*8 ”;<="2 1*2>4'7-6+=6(2(6j'*K…2 ;%("4"6+.2("46$oD2iNj'T6@”"4>&+="46O„

1*-3Qpo:1*-3I6?M”M”W'\6@"2 "t"&-d^D+.%('7-;89"&-d2!+=6m+=3;"&-:2+=>$1*<+.-‚1*<.<2 ;%("4"t6+.2("46$o:"\FPF.od2 "&‡

14^R"–2 "6@1*89"1*+.%w2("&8 ”W"&%@1*2?;%("46 F i-™2 "31*2@196("&2]+.<.<.?62%@1*2("43A+.- ls+.P\?;%("EF0o

2 "T"&8 ”;+.%+=>$1*<e”;%('n”W'7%2+='7-g'*Kj2@1bU0'7-+=6 0FHŠ+.- 6+.2("46 T1*-3 ‚1*-3 0F€+.- 6+.2("

0o1*-3™2 ""&-d^D+.%('7-;89"&-d2+=6+.-:^71*%+J1*-:2 e "&2iNj"4"&- 6+.2("46$F!l'7%w6(+.8 ”;<.+=>&+.2i‡RoNZ"?6("y1

>4'7-62@1*-d2%("46;”W'7-6("tKY?;->&2+='7-Z

\W ] F•6 ” {*F36:’!KM'7%]1*<.<S"&-d^D+.%('7-;89"&-d2(6UW’

NO "&%("UF•6y—™˜ -

+=6O162%+=>&2<.‡z”'\6+.2+.^R"6(>$1*<J1*%j%("46?;<.2+.-;PT+.-h1%("&P\%("46(6(+='7-QKY?;->&2+='7-YZ

C –WFd4’;F

‡

’@W i” { F•6

Fd…• F ‡

<—~,-7’@|0

KY'7%–"uU[”W"4>&2("43H”;%('n”W'7%2+='7-6'*K!6;”W"4>&+="46$FB5g"Q176(6?;89"92 1*2aF36I+=6Bhzij_kEm+jlk1h?;-;+V)

KY'7%8T<.‡‚30+=62%+e ?;2("43 +.-A6("&2€

’iFIj'7-6(":?"&-d2<.‡Ro

 Z ” { rY6("4" FEFRvuFOl…+.P7)

?;%("AEF C+.<.<.?62%@1*2("46 '$N 1*- 173 ?6289"&-:2‚'*K1*- ?M”M”"&% e '7?;-3 13u"4>&2(62 "

”W'\62("&%+='7%30+=62%+e ?;2+='7-g'*K”;%(' e 1 e +.<.+.2i‡ '*KZ6@”"4>&+="46! 'D>4>&?;%("&->4"9176"g{ - 1*-3

{}-11

Fkl;?;%2 "&%$oZl…+.P\?;%("‚EF ƒ%("G”;%("46@"&-d2(692 "h<.+.8T+.2+.-;P 30+=62%+e ?;2+='7- 2 1*2

176w-' `t+.%+=>( ;<="&2<J14‡R"&%$o+,F‰"\F.o#

2 ( Ft! "a”'\62("&%+='7%3;"&-6+.2+="46t+.-ƒl…+.P\?;%("EF

1*%("t1•”M”;%('$UD+.8Q1*2("43?6+.-;P  ~ %@1*-3;'786 1*8 ”;<.+.-;PNO+.2 \oH‹0o+.2("&%@1*2+='7-6

NO+.2 1 e ?;%-0),+.- ”1*%@1*89"&2("&%^\1*<.?" 0oƒ1*-3 1‚2 ;+.-;-;+.-;P ^\1*<.?" \o0F‚l…+.P7)

?;%("wEF %("&^R"$1*<=62 1*2ŠhŠij_kEm+jlk6;”W"4>&+>$1*2+='7-A'*K…2 " d?1*-:2+.2i‡

 Z ” +=6$oD+.-3;"4"43Xo1

6+.P\-;+>$1*-d2t3;"&2@1*+.<x+.-†”;%@17>&2+=>4"\F

x"&2T?6€-"uU02€>4'7-6(+=3;"&%Q1*- "uU9”;<.+=>&+.29>4'7-:2%('7<'*Kw2 "‚":?;+.^71*<="&-:296@1*8 ”;<="‚6+.#$"

rZB;B;vuoX+,F‰"\F.oS1h6("&22+.-;P‚+.- NO ;+=>@ g2 " `t+.%+=>( ;<="&2–<J1$‡R"&%–+=6'*KjKY'7%8 WFEF DF€! ;+=6

%("46?;<.2(6O+.-†”;%+='7%]30+=62%+e

?;2+='7-C B F  ’;F

‡

’@W ” d p!"&2@1

%$  ’

&$

‡ ” KM'7%]1*<.<S6+.2("46x=…{>-7’@|Q1*-3';’

KY'7%9'D>4>&?;%("&->4" ”;%(' e 1 e +.2+="465r

C

)( { - / C B

vuoZNO "&%(" $ { YZ

C – W  ’;F  ’;F

‡ ” {

F  EFRj• F ‡ ” bX

1*-3

 +=691*-f"d?;+.^\1*<="&-d296 1*8 ”;<="6(+.#$"h1*2T6+.2(" =@F „ e 6("&%^R"h2 1*2

*,+l=HW-.1= `)] `S/.9QW>SM`"S%0 SM`"]21%:43 ^<`SIQ`4.9V=o`"]216587:97Z [\U%=G]!;S?Ib-=<?> ?@- LGbBA:5&L#CD

fE<h

Fj

0 0.2 0.4 0.6 0.8 1

Posterior density

In document Computational methods and models for paleoecology (sivua 51-56)