From cottrell@wfu.edu Sat Mar 12 16:50:37 2011
From: Allin Cottrell <cottrell@wfu.edu>
To: gretl-users@gretlml.univpm.it
Subject: Re: [Gretl-users] Threshold VAR: Coding Help
Date: Sat, 12 Mar 2011 16:50:36 -0500
Message-ID: <Pine.A41.4.58.1103121636260.217446@f1n11.sp2net.wfu.edu>
In-Reply-To: =?utf-8?q?21700=5F1299964457=5Fp2CLEG5c021904=5FAANLkTi=3DP?=
 =?utf-8?q?oeaGiRZfod=3D=3Dkypht-2RMvT1Gfen-q3zFMDi=40mail=2Egmail=2Ecom?=
MIME-Version: 1.0
Content-Type: multipart/mixed; boundary="===============7703799114538640060=="

--===============7703799114538640060==
Content-Type: text/plain; charset="utf-8"
Content-Transfer-Encoding: 7bit

On Sat, 12 Mar 2011, Henrique Andrade wrote:

> *Y(t) = X(t-1) + e1(t); if z(t-d)>=tau*
> *Y(t) = X(t-1) + e2(t); if z(t-d)<tau*

If it's a TVAR I _think_ it should look more like

Y_t = D1 + B1(L)Y_t + I_t*(D2 + B2(L)Y_t) + U_t

where I_t = 1 if z(t-d) >= tau, otherwise 0.

In that case something like the following should do it:

<script>
open australia.gdt
list Y = PAU IAU E
# VAR without any regime-changing
var 1 Y
# make a fake threshold
series z = normal()
scalar tau = 1.5
series r2 = z(-1) >= tau
loop foreach i Y
  series r2$i = r2 * $i(-1)
endloop
# VAR with regime changing
var 1 Y ; r2*
</script>

Allin Cottrell

--===============7703799114538640060==--


From henrique.coelho@gmail.com Sat Mar 12 19:48:45 2011
From: Henrique Andrade <henrique.coelho@gmail.com>
To: gretl-users@gretlml.univpm.it
Subject: Re: [Gretl-users] Threshold VAR: Coding Help
Date: Sat, 12 Mar 2011 21:48:14 -0300
Message-ID: <AANLkTimjS4BFqw6p1vokbKXqpLkmUaYXCQZcOiJtnNLm@mail.gmail.com>
In-Reply-To: Pine.A41.4.58.1103121636260.217446@f1n11.sp2net.wfu.edu
Content-Type: multipart/mixed; boundary="===============1870581064347757528=="

--===============1870581064347757528==
Content-Type: text/html
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="attachment.html"
MIME-Version: 1.0

RW0gMTIgZGUgbWFy528gZGUgMjAxMaBBbGxpbiBlc2NyZXZldTo8ZGl2Pjxicj48ZGl2IGNsYXNz
PSJnbWFpbF9xdW90ZSI+PGJsb2NrcXVvdGUgY2xhc3M9ImdtYWlsX3F1b3RlIiBzdHlsZT0ibWFy
Z2luOjAgMCAwIC44ZXg7Ym9yZGVyLWxlZnQ6MXB4ICNjY2Mgc29saWQ7cGFkZGluZy1sZWZ0OjFl
eCI+T24gU2F0LCAxMiBNYXIgMjAxMSwgSGVucmlxdWUgQW5kcmFkZSB3cm90ZTo8YnI+CgoKCjxi
cj4KJmd0OyAqWSh0KSA9IFgodC0xKSArIGUxKHQpOyBpZiB6KHQtZCkmZ3Q7PXRhdSo8YnI+CiZn
dDsgKlkodCkgPSBYKHQtMSkgKyBlMih0KTsgaWYgeih0LWQpJmx0O3RhdSo8YnI+Cjxicj4KSWYg
aXQmIzM5O3MgYSBUVkFSIEkgX3RoaW5rXyBpdCBzaG91bGQgbG9vayBtb3JlIGxpa2U8YnI+Cjxi
cj4KWV90ID0gRDEgKyBCMShMKVlfdCArIElfdCooRDIgKyBCMihMKVlfdCkgKyBVX3Q8YnI+Cjxi
cj4Kd2hlcmUgSV90ID0gMSBpZiB6KHQtZCkgJmd0Oz0gdGF1LCBvdGhlcndpc2UgMC48YnI+PC9i
bG9ja3F1b3RlPjxkaXY+PGJyPjwvZGl2PjxkaXY+RGVhciBBbGxpbiwgSSBuZWVkIHRoYXQgdGhl
IGVycm9ycyAoYW5kIHZhcmlhbmNlcykgZGlmZmVyIGluIGVhY2ggcmVnaW1lOjwvZGl2PjxkaXY+
PGJyPjwvZGl2PjxkaXY+WV90ID0gRDEgKyBCMShMKVlfdCArIFUxX3QsIKBpZqB6KHQtZCkgJmd0
Oz0gdGF1PC9kaXY+CgoKPGRpdj5ZX3QgPSBEMiArIEIyKEwpWV90ICsgVTJfdCwgoG90aGVyd2lz
ZTwvZGl2PjxkaXY+PGJyPjwvZGl2PjxkaXY+U28gSSB0aGluayBJIG5lZWQgdG8gdXNlIHRoZSB0
aHJlc2hvbGQgdmFsdWVzIHRvIGJ1aWxkIHR3byBkYXRhIHNlZ21lbnRzLiBUbyBkbyB0aGlzIEkm
IzM5O20gZG9pbmcgc29tZXRoaW5nIGNsb3NlIHRvIHlvdXIgc2NyaXB0OjwvZGl2PjxkaXY+PGJy
PjwvZGl2PgoKPGRpdj48Yj4KJmx0O3NjcmlwdCZndDs8L2I+PC9kaXY+PGRpdj48ZGl2PjxiPjxm
b250IGNvbG9yPSIjMDA5OTAwIj5zZXJpZXM8L2ZvbnQ+PC9iPiB6ID0gPGZvbnQgY29sb3I9IiMz
Mzk5OTkiPm5vcm1hbDwvZm9udD4oKTwvZGl2PjxkaXY+PGI+PGZvbnQgY29sb3I9IiMwMDk5MDAi
PnNjYWxhcjwvZm9udD48L2I+IHRhdSA9IDEuNTwvZGl2Pgo8ZGl2Pjxmb250IGNvbG9yPSIjMDA5
OTAwIj48Yj5zZXJpZXM8L2I+PC9mb250PiByMSA9IHooLTEpICZndDs9IHRhdTwvZGl2PjxkaXY+
PGZvbnQgY29sb3I9IiMwMDk5MDAiPjxiPnNlcmllczwvYj48L2ZvbnQ+IHIyID0gMS1yMTwvZGl2
PjxkaXY+PGRpdj48Yj48Zm9udCBjb2xvcj0iIzY2MzMwMCI+bG9vcCBmb3JlYWNoPC9mb250Pjwv
Yj4gaSBZPC9kaXY+CjxkaXY+oKAgoDxmb250IGNvbG9yPSIjMDA5OTAwIj48Yj5zZXJpZXM8L2I+
PC9mb250PiByMSRpID0gPGZvbnQgY29sb3I9IiMzMzk5OTkiPnplcm9taXNzPC9mb250PihyMSAq
ICRpKTwvZGl2PjxkaXY+oKAgoDxmb250IGNvbG9yPSIjMDA5OTAwIj48Yj5zZXJpZXM8L2I+PC9m
b250PiByMiRpID0gPGZvbnQgY29sb3I9IiMzMzk5OTkiPnplcm9taXNzPC9mb250PihyMiAqICRp
KTwvZGl2PgoKCjxkaXY+PGI+PGZvbnQgY29sb3I9IiM2NjMzMDAiPmVuZGxvb3A8L2ZvbnQ+PC9i
PjwvZGl2PjwvZGl2PjwvZGl2PjxkaXY+PGI+PGZvbnQgY29sb3I9IiM2NjMzMDAiPjxkaXY+PHNw
YW4gY2xhc3M9IkFwcGxlLXN0eWxlLXNwYW4iIHN0eWxlPSJmb250LXdlaWdodDogbm9ybWFsOyI+
PGZvbnQgY2xhc3M9IkFwcGxlLXN0eWxlLXNwYW4iIGNvbG9yPSIjMDAwMDAwIj48c3BhbiBjbGFz
cz0iQXBwbGUtc3R5bGUtc3BhbiIgc3R5bGU9ImNvbG9yOiByZ2IoMCwgMTUzLCAwKTsgIj48Yj5s
aXN0PC9iPjwvc3Bhbj6gcmVnaW1lMSA9IHIxUEFVIHIxSUFVIHIxRTwvZm9udD48L3NwYW4+PC9k
aXY+Cgo8ZGl2PjxzcGFuIGNsYXNzPSJBcHBsZS1zdHlsZS1zcGFuIiBzdHlsZT0iZm9udC13ZWln
aHQ6IG5vcm1hbDsiPjxmb250IGNsYXNzPSJBcHBsZS1zdHlsZS1zcGFuIiBjb2xvcj0iIzAwMDAw
MCI+PHNwYW4gY2xhc3M9IkFwcGxlLXN0eWxlLXNwYW4iIHN0eWxlPSJjb2xvcjogcmdiKDAsIDE1
MywgMCk7ICI+PGI+bGlzdDwvYj48L3NwYW4+oHJlZ2ltZTIgPSByMlBBVSByMklBVSByMkWgPC9m
b250Pjwvc3Bhbj48L2Rpdj4KCjwvZm9udD48L2I+PC9kaXY+PGRpdj48Yj4mbHQ7L3NjcmlwJmd0
OzwvYj48L2Rpdj48ZGl2Pjxicj48L2Rpdj48ZGl2Pk5vdyBJIGp1c3QgbmVlZCB0byBlc3RpbWF0
ZSB0d28gc2VwYXJhdGUgVkFSczo8L2Rpdj48ZGl2PjxkaXY+PGRpdiBjbGFzcz0iZ21haWxfcXVv
dGUiPjxkaXY+PGJyPjwvZGl2PjxkaXY+PGI+Jmx0O3NjcmlwdCZndDs8L2I+PC9kaXY+PGRpdj48
ZGl2PjxiPjxmb250IGNvbG9yPSIjNjYzMzAwIj48c3BhbiBjbGFzcz0iQXBwbGUtc3R5bGUtc3Bh
biIgc3R5bGU9ImNvbG9yOiByZ2IoMCwgMCwgMCk7IGZvbnQtd2VpZ2h0OiBub3JtYWw7ICI+oDwv
c3Bhbj48c3BhbiBjbGFzcz0iQXBwbGUtc3R5bGUtc3BhbiIgc3R5bGU9ImZvbnQtd2VpZ2h0OiBu
b3JtYWw7Ij48Zm9udCBjbGFzcz0iQXBwbGUtc3R5bGUtc3BhbiIgY29sb3I9IiNDQzMzQ0MiPiZx
dW90O1ZBUiByZWdpbWUgMSZxdW90OzwvZm9udD48L3NwYW4+PC9mb250PjwvYj48Zm9udCBjbGFz
cz0iQXBwbGUtc3R5bGUtc3BhbiIgY29sb3I9IiNDQzMzQ0MiPqA8L2ZvbnQ+Jmx0Oy0gPGZvbnQg
Y2xhc3M9IkFwcGxlLXN0eWxlLXNwYW4iIGNvbG9yPSIjNjYzMzAwIj48Yj52YXI8L2I+PC9mb250
PqAxIHJlZ2ltZTE8L2Rpdj4KCjxkaXY+PGRpdj48Zm9udCBjbGFzcz0iQXBwbGUtc3R5bGUtc3Bh
biIgY29sb3I9IiM2NjMzMDAiPjxiPjxzcGFuIGNsYXNzPSJBcHBsZS1zdHlsZS1zcGFuIiBzdHls
ZT0iY29sb3I6IHJnYigwLCAwLCAwKTsgZm9udC13ZWlnaHQ6IG5vcm1hbDsgIj48Yj48Zm9udCBj
b2xvcj0iIzY2MzMwMCI+PHNwYW4gY2xhc3M9IkFwcGxlLXN0eWxlLXNwYW4iIHN0eWxlPSJjb2xv
cjogcmdiKDAsIDAsIDApOyBmb250LXdlaWdodDogbm9ybWFsOyAiPqA8L3NwYW4+PHNwYW4gY2xh
c3M9IkFwcGxlLXN0eWxlLXNwYW4iIHN0eWxlPSJmb250LXdlaWdodDogbm9ybWFsOyAiPjxmb250
IGNsYXNzPSJBcHBsZS1zdHlsZS1zcGFuIiBjb2xvcj0iI0NDMzNDQyI+JnF1b3Q7VkFSIHJlZ2lt
ZSAyJnF1b3Q7PC9mb250Pjwvc3Bhbj48L2ZvbnQ+PC9iPjxmb250IGNsYXNzPSJBcHBsZS1zdHls
ZS1zcGFuIiBjb2xvcj0iI0NDMzNDQyI+oDwvZm9udD4mbHQ7LaA8Zm9udCBjbGFzcz0iQXBwbGUt
c3R5bGUtc3BhbiIgY29sb3I9IiM2NjMzMDAiPjxiPnZhcjwvYj48L2ZvbnQ+oDEgcmVnaW1lMjwv
c3Bhbj48L2I+PC9mb250PjwvZGl2PgoKPC9kaXY+PC9kaXY+PGRpdj48Yj48Zm9udCBjb2xvcj0i
IzY2MzMwMCI+PGRpdj48c3BhbiBjbGFzcz0iQXBwbGUtc3R5bGUtc3BhbiIgc3R5bGU9ImNvbG9y
OiByZ2IoMCwgMCwgMCk7IGZvbnQtd2VpZ2h0OiBub3JtYWw7ICI+PGI+Jmx0Oy9zY3JpcCZndDs8
L2I+PC9zcGFuPjwvZGl2PjwvZm9udD48L2I+PC9kaXY+PC9kaXY+PC9kaXY+PC9kaXY+PGRpdj48
YnI+PC9kaXY+PGRpdj5CdXQgd2UgaGF2ZSBtaXNzaW5nIG9ic2VydmF0aW9ucy4uLiBIb3cgY2Fu
IHdlIGhhbmRsZSB0aGlzPzwvZGl2PgoKPGRpdj48YnI+PC9kaXY+PGRpdj5GaW5hbGx5LCBJIHdv
dWxkIGxpa2UgdG8gbWFrZSBvbmUgbW9yZSBxdWVzdGlvbi4gTG9va2luZyBhdCB5b3VyIG9yaWdp
bmFsIGNvZGUsIHdlIGhhdmUsIGluIHRoZSBsYXN0IGxpbmUsIHRoaXMgZXhwcmVzc2lvbjqgPGk+
dmFyIDEgWSA7IHIyKjwvaT4uIFdoYXQgdGhlIHN5bWJvbCAmcXVvdDsqJnF1b3Q7IChpbiByMiop
IG1lYW5zPzwvZGl2PjxkaXY+Cgo8YnI+PC9kaXY+PGRpdj5CZXN0IHJlZ2FyZHMsPC9kaXY+PGRp
dj4tLaA8L2Rpdj48L2Rpdj48Yj48Zm9udCBmYWNlPSJ0YWhvbWEsIHNhbnMtc2VyaWYiPkhlbnJp
cXVlIEMuIGRlIEFuZHJhZGU8L2ZvbnQ+PC9iPjxmb250IGZhY2U9InRhaG9tYSwgc2Fucy1zZXJp
ZiI+PGJyPkRvdXRvcmFuZG8gZW0gRWNvbm9taWEgQXBsaWNhZGE8YnI+VW5pdmVyc2lkYWRlIEZl
ZGVyYWwgZG8gUmlvIEdyYW5kZSBkbyBTdWw8YnI+CgoKPGEgaHJlZj0iaHR0cDovL3d3dy51ZnJn
cy5ici9wcGdlIiB0YXJnZXQ9Il9ibGFuayI+d3d3LnVmcmdzLmJyL3BwZ2U8L2E+PC9mb250Pjxi
cj4KPC9kaXY+Cg==

--===============1870581064347757528==--


From peter.summers@ttu.edu Sat Mar 12 20:08:25 2011
From: Summers, Peter <peter.summers@ttu.edu>
To: gretl-users@gretlml.univpm.it
Subject: Re: [Gretl-users] Threshold VAR: Coding Help
Date: Sat, 12 Mar 2011 19:04:27 -0600
Message-ID: <8885DA157278484489B1EBC8B8C043E3FDC8A13364@COTTUS.ttu.edu>
In-Reply-To: AANLkTimjS4BFqw6p1vokbKXqpLkmUaYXCQZcOiJtnNLm@mail.gmail.com
MIME-Version: 1.0
Content-Type: multipart/mixed; boundary="===============8943951163478508042=="

--===============8943951163478508042==
Content-Type: text/plain; charset="utf-8"
Content-Transfer-Encoding: quoted-printable

FWIW, in a univariate TAR model the data gets re-ordered by (say) ascending v=
alues of the threshold variable, then separate regressions get run on the 2 s=
ub-samples. Could this be done with an "msortby" command followed by smpl com=
mands?

PS

________________________________________
From: gretl-users-bounces(a)lists.wfu.edu [gretl-users-bounces(a)lists.wfu.ed=
u] On Behalf Of Henrique Andrade [henrique.coelho(a)gmail.com]
Sent: Saturday, March 12, 2011 6:48 PM
To: Gretl list
Subject: Re: [Gretl-users] Threshold VAR: Coding Help

Em 12 de mar=C3=A7o de 2011 Allin escreveu:

On Sat, 12 Mar 2011, Henrique Andrade wrote:

> *Y(t) =3D X(t-1) + e1(t); if z(t-d)>=3Dtau*
> *Y(t) =3D X(t-1) + e2(t); if z(t-d)<tau*

If it's a TVAR I _think_ it should look more like

Y_t =3D D1 + B1(L)Y_t + I_t*(D2 + B2(L)Y_t) + U_t

where I_t =3D 1 if z(t-d) >=3D tau, otherwise 0.

Dear Allin, I need that the errors (and variances) differ in each regime:

Y_t =3D D1 + B1(L)Y_t + U1_t,  if z(t-d) >=3D tau
Y_t =3D D2 + B2(L)Y_t + U2_t,  otherwise

So I think I need to use the threshold values to build two data segments. To =
do this I'm doing something close to your script:

<script>
series z =3D normal()
scalar tau =3D 1.5
series r1 =3D z(-1) >=3D tau
series r2 =3D 1-r1
loop foreach i Y
    series r1$i =3D zeromiss(r1 * $i)
    series r2$i =3D zeromiss(r2 * $i)
endloop
list regime1 =3D r1PAU r1IAU r1E
list regime2 =3D r2PAU r2IAU r2E
</scrip>

Now I just need to estimate two separate VARs:

<script>
 "VAR regime 1" <- var 1 regime1
 "VAR regime 2" <- var 1 regime2
</scrip>

But we have missing observations... How can we handle this?

Finally, I would like to make one more question. Looking at your original cod=
e, we have, in the last line, this expression: var 1 Y ; r2*. What the symbol=
 "*" (in r2*) means?

Best regards,
--
Henrique C. de Andrade
Doutorando em Economia Aplicada
Universidade Federal do Rio Grande do Sul
www.ufrgs.br/ppge<http://www.ufrgs.br/ppge>


--===============8943951163478508042==--


From henrique.coelho@gmail.com Sat Mar 12 23:08:16 2011
From: Henrique Andrade <henrique.coelho@gmail.com>
To: gretl-users@gretlml.univpm.it
Subject: Re: [Gretl-users] Threshold VAR: Coding Help
Date: Sun, 13 Mar 2011 01:07:41 -0300
Message-ID: <AANLkTim=AJaoPJg8aux3nsTYCN9VF=sB1Sao6n1+pNjc@mail.gmail.com>
In-Reply-To: 8885DA157278484489B1EBC8B8C043E3FDC8A13364@COTTUS.ttu.edu
Content-Type: multipart/mixed; boundary="===============0174297264483648822=="

--===============0174297264483648822==
Content-Type: text/html
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="attachment.html"
MIME-Version: 1.0

RW0gMTIgZGUgbWFy528gZGUgMjAxMaBQZXRlciA8c3BhbiBkaXI9Imx0ciI+Jmx0OzxhIGhyZWY9
Im1haWx0bzpwZXRlci5zdW1tZXJzQHR0dS5lZHUiPnBldGVyLnN1bW1lcnNAdHR1LmVkdTwvYT4m
Z3Q7IGVzY3JldmV1Ojwvc3Bhbj48ZGl2PjxzcGFuIGRpcj0ibHRyIj48L3NwYW4+PGJyPjxkaXYg
Y2xhc3M9ImdtYWlsX3F1b3RlIj48YmxvY2txdW90ZSBjbGFzcz0iZ21haWxfcXVvdGUiIHN0eWxl
PSJtYXJnaW46MCAwIDAgLjhleDtib3JkZXItbGVmdDoxcHggI2NjYyBzb2xpZDtwYWRkaW5nLWxl
ZnQ6MWV4OyI+CgpGV0lXLCBpbiBhIHVuaXZhcmlhdGUgVEFSIG1vZGVsIHRoZSBkYXRhIGdldHMg
cmUtb3JkZXJlZCBieSAoc2F5KSBhc2NlbmRpbmcgdmFsdWVzIG9mIHRoZSB0aHJlc2hvbGQgdmFy
aWFibGUsIHRoZW4gc2VwYXJhdGUgcmVncmVzc2lvbnMgZ2V0IHJ1biBvbiB0aGUgMiBzdWItc2Ft
cGxlcy4gQ291bGQgdGhpcyBiZSBkb25lIHdpdGggYW4gJnF1b3Q7bXNvcnRieSZxdW90OyBjb21t
YW5kIGZvbGxvd2VkIGJ5IHNtcGwgY29tbWFuZHM/PGJyPgoKPC9ibG9ja3F1b3RlPjxkaXY+PGJy
PjwvZGl2PjxkaXY+RGVhciBQZXRlciwgYmFzZWQgdXBvbiBFeGFtcGxlIDQuMyBmcm9tIHRoZSBS
QVRTIFByb2dyYW1taW5nIE1hbm51YWwgYnkgV2FsdGVyIEVuZGVycywgSSB3cm90ZSBhIHNjcmlw
dCB0aGF0IGNhbiBiZSB1c2VkIHRvIGVzdGltYXRlIGEgc2ltcGxlIFRBUiBtb2RlbCB3aXRoIDIg
cmVnaW1lcy4gTWF5YmUgeW91JiMzOTtkIGxpa2UgdG8gdGFrZSBhIGxvb2sgaW4gaXQuPC9kaXY+
Cgo8L2Rpdj48ZGl2Pjxicj48L2Rpdj5CZXN0LDxicj4tLSA8YnI+PGI+PGZvbnQgZmFjZT0idGFo
b21hLCBzYW5zLXNlcmlmIj5IZW5yaXF1ZSBDLiBkZSBBbmRyYWRlPC9mb250PjwvYj48Zm9udCBm
YWNlPSJ0YWhvbWEsIHNhbnMtc2VyaWYiPjxicj5Eb3V0b3JhbmRvIGVtIEVjb25vbWlhIEFwbGlj
YWRhPGJyPlVuaXZlcnNpZGFkZSBGZWRlcmFsIGRvIFJpbyBHcmFuZGUgZG8gU3VsPGJyPgoKPGEg
aHJlZj0iaHR0cDovL3d3dy51ZnJncy5ici9wcGdlIiB0YXJnZXQ9Il9ibGFuayI+d3d3LnVmcmdz
LmJyL3BwZ2U8L2E+PC9mb250Pjxicj4KPC9kaXY+Cg==

--===============0174297264483648822==
Content-Type: application/octet-stream
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="TAR.inp"
MIME-Version: 1.0

IyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIwoj
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAjCiMg
ICAgVGhyZXNob2xkIEF1dG9yZWdyZXNpb24gRXN0aW1hdGlvbiB3aXRoIEdyZXRsICAgICMKIyAg
ICAgICAgICAgICAgICAgICAgICAgICAgSGVucmlxdWUgQW5kcmFkZSAoMjAxMSkgICAgIwojICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAjCiMgICAg
T2JzOiBCYXNlZCB1cG9uIEV4YW1wbGUgNC4zIGZyb20gdGhlIFJBVFMgUHJvLSAgICMKIyAgICAg
ICAgIGdyYW1taW5nIE1hbm51YWwgYnkgV2FsdGVyIEVuZGVycyAgICAgICAgICAgIwojICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAjCiMjIyMjIyMj
IyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMKCnNldCBlY2hv
IG9mZgpzZXQgbWVzc2FnZXMgb2ZmCm9wZW4gYXVzdHJhbGlhLmdkdAoKc2VyaWVzIGxkX1BBVSA9
IGxkaWZmKFBBVSkKCnNjYWxhciBsb3cgPSAxICsgaW50KDAuMTUqJG5vYnMpCnNjYWxhciBoaWdo
ID0gJG5vYnMgLSBpbnQoMC4xNSokbm9icykKc2NhbGFyIGVzc190ZXN0ID0gMTAwMDAwMApzZXJp
ZXMgdGhyZXNoX3Rlc3QgPSBzb3J0KGxkX1BBVSkKc2NhbGFyIHRocmVzaG9sZCA9IHRocmVzaF90
ZXN0W2xvd10KCm1hdHJpeCBFU1MgPSB6ZXJvcygkbm9icywxKQoKbG9vcCBpID0gbG93Li5oaWdo
IC0tcXVpZXQKICAgIHNlcmllcyBwbHVzID0gKGxkX1BBVSgtMSk+PXRocmVzaF90ZXN0WyRpXSkK
ICAgIHNlcmllcyBtaW51cyA9IDEgLSBwbHVzCiAgICBzZXJpZXMgeTFfcGx1cyA9IHBsdXMqbGRf
UEFVKC0xKQogICAgc2VyaWVzIHkyX3BsdXMgPSBwbHVzKmxkX1BBVSgtMikKICAgIHNlcmllcyB5
MV9taW51cyA9IG1pbnVzKmxkX1BBVSgtMSkKICAgIHNlcmllcyB5Ml9taW51cyA9IG1pbnVzKmxk
X1BBVSgtMikKICAgIG9scyBsZF9QQVUgcGx1cyB5MV9wbHVzIHkyX3BsdXMgbWludXMgeTFfbWlu
dXMgeTJfbWludXMgLS1xdWlldAogICAgc2NhbGFyIGVzc18kaSA9ICRlc3MKICAgIEVTU1skaSxd
ID0gZXNzXyRpCiAgICBpZiAkZXNzIDwgZXNzX3Rlc3QKICAgICAgICBzY2FsYXIgZXNzX3Rlc3Qg
PSAkZXNzCiAgICAgICAgc2NhbGFyIHRocmVzaG9sZCA9IHRocmVzaF90ZXN0WyRpXQogICAgZW5k
aWYKZW5kbG9vcApwcmludCB0aHJlc2hvbGQKCnNlcmllcyBwbHVzID0gKGxkX1BBVSgtMSk+PXRo
cmVzaG9sZCkKc2VyaWVzIG1pbnVzID0gMSAtIHBsdXMKc2VyaWVzIHkxX3BsdXMgPSBwbHVzKmxk
X1BBVSgtMSkKc2VyaWVzIHkyX3BsdXMgPSBwbHVzKmxkX1BBVSgtMikKc2VyaWVzIHkxX21pbnVz
ID0gbWludXMqbGRfUEFVKC0xKQpzZXJpZXMgeTJfbWludXMgPSBtaW51cypsZF9QQVUoLTIpCm9s
cyBsZF9QQVUgcGx1cyB5MV9wbHVzIHkyX3BsdXMgbWludXMgeTFfbWludXMgeTJfbWludXMKCm1h
dHJpeCBFU1MgPSBFU1NbbG93OiRub2JzLWxvdysxLF0KCnNtcGwgbG93IGhpZ2gKc2VyaWVzIGVz
cyA9IEVTU1ssMV0KZ251cGxvdCBlc3MgdGhyZXNoX3Rlc3Qge3NldCB0aXRsZSAiU29tYSBkb3Mg
UXVhZHJhZG9zIGRvcyBSZXPDrWR1b3MiOyBzZXQgeGxhYmVsICJUaHJlc2hvbGQiOyBzZXQgeWxh
YmVsICJlc3MiOyBzZXQgc3R5bGUgbGluZSAxIGx3IDIgbGMgcmdiICIjMGYwMDAwIjt9IC0tc3Vw
cHJlc3MtZml0dGVkIC0td2l0aC1saW5lcyAtLW91dHB1dD1kaXNwbGF5Cg==

--===============0174297264483648822==--