kiểm định giả thuyết \(H_0: p=0,5\),
với đối thuyết \(H_0: p≠0,5\) mức ý
nghĩa 5%. Vì cỡ mẫu n = 15 là bé, nên tiêu chuẩn bác bỏ \(H_0\) là:
\(f≥\frac{\frac{x \alpha}{2}}{n}\),
hoặc \(f≤ \frac{\frac{\bar{x}
\alpha}{2}}{n}\), tức là \(f∉
(\frac{\frac{\bar{x} \alpha}{2}}{n};\frac{\frac{x
\alpha}{2}}{n})\)
Từ mẫu ta có: \(n=15,
f=\frac{8}{15}\). Với \(\alpha =
0,5\), từ bảng phụ lục 1, điều kiện:
\[
\sum_{x=x_{0,025}}^{15} C_{15}^x.(0,5)^{15}≤ 0,025
\]
cho ta \(x_{0,025}=12\), nên \(\frac{x_{\frac{\alpha}{2}}}{n}=\frac{12}{15}\),
và từ điều kiện:
\[
\sum_{x=0}^{\bar{x}_{0,025}} C_{15}^x.(0,5)^{15}≤ 0,025
\]
cho ta \(\bar{x}_{0,025}=3\), nên
\(\frac{\bar{x}_{\frac{\alpha}{2}}}{n}=\frac{3}{15}\).
Ta có: \(f=\frac{8}{15}∈(\frac{3}{15};\frac{12}{15})\)
ta chấp nhận \(H_0:p=0,5\) và như vậy
là chưa có cơ sở để cho rằng trong số khách hàng của hệ dịch vụ A, tỷ lệ
chọn hai mức phí là khác nhau.
Khoảng tin cậy cho một tỷ lệ
Trong hầu hết các ứng dụng, người ta sử dụng khoảng tin cậy để xác
định phạm vi của các giá trị hợp lý cho tham số. Chúng ta biết rằng ước
lượng ML cho tỷ lệ \(\mathrm{p}\), tỷ
lệ tính chất A trong tổng thể là \(f=f(A)\), tỷ lệ tính chất A trên mẫu và khi
cỡ mẫu \(n \rightarrow+\infty\), phân
phối của tỷ lệ mẫu \(f\) dần tới phân
phối chuẩn, hay phân phối của \(\frac{(f-p)
\sqrt{n}}{\sqrt{p(1-p)}}\) dần tới phân phối chuẩn chính tắc. Vậy
nên khi cỡ mẫu n “đủ lớn”, người ta có thể xấp xỉ phân phối của \(\frac{(f-p) \sqrt{n}}{\sqrt{p(1-p)}}\) với
phân phối chuẩn chính tắc. Tuy nhiên n thế nào là “đủ lớn” để có xấp xỉ
nói trên? Trở lại phân phối Nhị thức \(\mathrm{B}(\mathrm{n}, \mathrm{p})\), ta
biết rằng khi trung bình n.p của phân phối Nhị thức càng lớn thì phân
phối này càng gần với phân phối chuẩn. Mặt khác ta biết rằng \(n \cdot f(A)=\) \(m(A)\) và \(n-n .
f(A)\) là số lần xuất hiện và không xuất hiện A trong n lần thử
có phân phối \(\mathrm{B}(\mathrm{n},
\mathrm{p})\) và \(\mathrm{B}(\mathrm{n}, 1-\mathrm{p})\). Vậy
điều kiện để xấp xỉ phân phối của \(f\)
và \(1-\) \(f\) với phân phối chuẩn là \(\mathrm{n}\) đủ lớn sao cho n.p và \(\mathrm{n}(1-\mathrm{p})\) cũng khá lớn.
Trong thực hành vì p chưa biết, người ta thay nó bởi ước lượng ML là
\(f\), vì vậy điều kiện n.p và \(\mathrm{n} .(1-\mathrm{p})\) khá lớn được
thay bằng \(n\). \(f\) và \(n\). \((1-f)\) khá lớn. Vì vậy đối với bài toán
ước lượng khoảng tin cậy cho tỷ lệ, người ta thường yêu cầu cỡ mẫu khá
lớn và chia làm hai trường hợp:
Khi \(f\)
không quá gần 0 hoặc 1:
Khoảng tin cậy, với độ tin cậy \(\gamma=1-\alpha\) cho tỷ lệ \(p\) là:
\[
\left(f-u\left(\frac{\alpha}{2}\right) \cdot \sqrt{\frac{f(1-f)}{n}},
f+u\left(\frac{\alpha}{2}\right) \cdot \sqrt{\frac{f(1-f)}{n}}\right)
\]
Chú ý: Thông thường độ tin cậy là \(95
\%\), khi đó khoảng tin cậy cho \(p\) là: \[
\left(f-1,96 . \sqrt{\frac{f(1-f)}{n}}, f+1,96 .
\sqrt{\frac{f(1-f)}{n}}\right)
\] và khi không nói gì tới độ tin cậy thì mặc định độ tin cậy là
\(95 \%\).
Ví du 4: Trong mẫu điều tra ở ví dụ 2, hãy ước lượng khoảng
tin cậy \(95 \%\) cho tỷ lệ người Mỹ
ủng hộ việc phá thai hợp pháp cho phụ nữ có thai mà không muốn có
con.
Giải: Từ mẫu, ta có tỷ lệ ủng hộ việc phá thai hợp pháp là \(f=0,446\).
Khoảng tin cậy cần tìm có dạng: \(\left(f-1,96 . \sqrt{\frac{f(1-f)}{n}}, f+1,96 .
\sqrt{\frac{f(1-f)}{n}}\right)\) thay số liệu vào ta có: \(f-1,96 . \sqrt{\frac{f(1-f)}{n}}=0,414 ; f+1,96 .
\sqrt{\frac{f(1-f)}{n}}=0,478\)
Vậy có thể cho rằng với mức độ tin cậy \(95
\%\) thì tỷ lệ dân số Mỹ ủng hộ việc phá thai hợp pháp cho phụ nữ
có thai mà không muốn có con nhiểu hơn là vào khoảng từ \(41,4 \%\) đến \(47,8 \%\).
Khi \(f \leq
0,1\), hoặc \(f \geq
0,9\) :
Trong trường hợp này, việc xấp xỉ phân phối nhị thức bởi phân phối
chuẩn nói chung là thô, vì mặc dù cỡ mẫu n khá lớn nhưng có thể \(n . f\) và \(n
.(1-f)\) lại không đủ lớn,vậy nên xấp xỉ bởi phân phối Poisson là
hợp lý hơn. Vì thế, bài toán khoảng tin cậy cho tỷ lệ \(\pi\) chia làm 2 bước:
Bước 1. Tìm khoảng tin cậy cho \(n .
p\) là khoảng \(\left(n p_1, n
p_2\right)\), trong đó \(n
p_1\), \(n p_2\) tìm được qua
tra bảng phụ lục IV, căn cứ vào giá trị tần số mẫu: \(m=\) n.f. (Trong bảng này: cột đầu tiên là
giá trị chẵn chục của \(m\) và dòng
trên cùng bổ sung chữ số hàng đơn vị của \(m\), giao của hàng (chẵn chục của \(m\) ) và cột (hàng đơn vị của \(m\) ) là ô chứa 2 cận: số trên là \(n p_1\), số dưới là \(\left.n p_2\right)\).
Bước 2. Chia 2 cận \(n p_1, n p_2\)
cho \(n\) (cỡ mẫu) nhận được khoảng
\(\left(p_1, p_2\right)\) là khoảng tin
cậy \(95 \%\) cho tỷ lệ \(p\).
Ví du 5: Điều tra mức thu nhập của 250 hộ dân ở một khu vực, thấy có
7 hộ thu nhập ở mức nghèo khổ. Hãy ước lượng khoảng tin cậy cho tỷ lệ hộ
thu nhập ở mức nghèo khổ ở khu vực này.
Giải: Gọi p là tỷ lệ hộ có thu nhập ở mức nghèo khổ của khu vực. Từ
mẫu điều tra có tỷ lệ hộ có thu nhập ở mức nghèo khổ là: \(f=\frac{7}{250}=0,014<0,1\). Vậy ta tìm
khoảng tin cậy cho p theo các bước sau:
- Ước lượng khoảng tin cậy cho n.p: Với \(m=n . f=7\), tra bảng phụ lục IV, có khoảng
tin cậy cho n.p là: \((2,31 ;
14,42)\).
-Suy ra khoảng tin cậy cho p là: \((0,00924
; 0,05768)\).
Vậy với độ tin cậy \(95 \%\), có thể
cho rằng tỷ lệ hộ có mức thu nhập nghèo khổ ở khu vực này vào khoảng từ
\(0,0924 \%\) đến \(5,768 \%\).
Nhận xét: Với trường hợp mẫu này, nếu ta cố xấp xỉ phân phối của
\(f\) với phân phối chuẩn, tức là dùng
công thức 1.5.1a, thì ước lượng sẽ quá thô, thậm chí có cận dưới: \[
f-1,96 . \sqrt{\frac{f(1-f)}{n}}=0,014-1,96 .
\sqrt{\frac{0,014(1-0,014)}{250}}=-0,0321<0 .
\]
Chú ý: Trong thực hành, kiểm định \(H_0:
p=0,5\) chỉ yêu cầu \(n>10\)
trong khi kiểm định \(H_0: p=0.1\) hoặc
\(H_0: p=0.9\) yêu cầu \(n>50\). Yêu cầu về kích cỡ mẫu phản ánh
sự thay đổi ngày càng lệch của việc phân phối mẫu của \(f\) khi \(p\) đạt đến 0 hoặc 1.
Tóm lại, chương này đã giới thiệu hai phân phối chính cho phân tích
dữ liệu định tính: nhị phân và Poisson. Đồng thời chúng tôi cũng đã giới
thiệu phương pháp ước lượng hợp lý cực đại và minh hoạ việc sử dụng nó
cho dữ liệu tỷ lệ. Phần còn lại của nội dung sử dụng suy luận ML cho các
tham số nhị phân và Poisson trong nhiều trường hợp khác nhau.
Bài tập chương 1
1.1. Trong các ví dụ sau, xác định biến đáp ứng và
các biến giải thích.
a. Thái độ đối với kiểm soát súng (ủng hộ, phản đối), giới tính (nữ,
nam), Giáo dục của mẹ (trung học, cao đẳng).
b. Bệnh tim (có, không), Huyết áp, Mức cholesterol (Cholesterol).
c. Chủng tộc (trắng, không trắng), Tôn giáo (Công giáo, Do Thái, Tin
lành), Bình chọn Tổng thống (Dân chủ, Cộng hòa, Khác), Thu nhập hàng
năm.
1.2. Thang đo nào phù hợp nhất cho các biến sau đây-
danh nghĩa, hoặc thứ tự?
a. Liên kết đảng chính trị (Dân chủ, Cộng hòa, không liên kết).
b. Có bằng cao nhất (không có, trung học, cử nhân, thạc sĩ, tiến
sĩ).
c. Tình trạng bệnh nhân (tốt, công bằng, nghiêm trọng, phê bình).
d. Vị trí của bệnh viện (London, Boston, Madison, Rochester,
Toronto).
e. Nước giải khát yêu thích (bia, nước trái cây, sữa, nước ngọt, rượu
vang, khác).
1.3. Các tấm silicon cho chip máy tính được sản xuất
bởi một công ty công nghệ cao có một số khuyết tật trung bình là 1.0 mỗi
tấm. Nếu số lượng khuyết tật có phân phối Poisson, tìm xác suất mà một
tấm có (a) đúng 0 khiếm tật, (b) đúng 1 khiếm tật, (c) ít nhất 2 khiếm
khuyết.
1.4. Mỗi trong 100 câu hỏi trắc nghiệm trong một kỳ
thi có bốn đáp án có thể nhưng chỉ có một đáp án đúng. Đối với mỗi câu
hỏi, học sinh chọn ngẫu nhiên một đáp án là câu trả lời. Xác định phân
phối số câu trả lời đúng của học sinh trong kỳ thi. Dựa trên trung bình
và độ lệch chuẩn của bản phân phối đó, liệu có đáng ngạc nhiên nếu học
sinh có ít nhất 50 câu trả lời đúng?
1.5. Một đồng xu cân đối được tung hai lần. Gọi là
số lượng mặt sấp thu được.
a. Xác định xác suất cho các giá trị có thể cho , giá trị trung bình
và phương sai của Y.
b. Tính các xác suất cho phân phối Poisson có cùng giá trị trung bình
của Y. Sự khác biệt của nó so với phương sai trong (a) như thế nào?
c. Đối với mỗi lần tung một đồng xu không đồng chất, gọi 𝜋 là biểu
thị xác suất của mặt sấp. Giả sử có 0 mặt sấp trong 2 lần tung. Tìm ước
lượng ML của 𝜋.
1.6. Quay lại các bài tập trước đó.
a. Tính các xác suất nhị thức cho \(n=2\) khi xác suất của mặt sấp cho mỗi lần
tung bằng (i) \(\pi=0.6\), (ii) \(\pi=0.4\).
b. Giả sử chúng ta quan sát \(Y=1\).
Tính toán và phác hoạ hàm hợp lý.
c. Sử dụng hàm hợp lý vẽ từ (b), chỉ ra ước lượng ML của \(\pi=0.5\).
1.7. Trong cuốn tự truyện A Sort of Life, tác giả
người Anh Graham Greene mô tả giai đoạn trầm cảm nghiêm trọng trong thời
gian đó ông chơi Russian Roulette. “Trò chơi” này bao gồm đặt một viên
đạn vào một trong sáu khoang súng, quay các buồng để chọn ngẫu nhiên, và
sau đó nhằm một đầu người và bóp cò.
a. Greene chơi game này sáu lần, và may mắn là không ai trong số họ
đã bị bắn một viên đạn. Tìm xác suất của kết quả này.
b. Giả sử một người tiếp tục chơi trò chơi này cho đến khi đạn bắn
trúng đích. Gọi Y là số viên đạn đã bắn. Chứng tỏ rằng xác suất của kết
quả y bằng với \((5/6)^{y-1}(1/6)\) ,
cho y=1,2,3,… (được gọi là phân phối hình học).
1.8. Một mẫu gồm những phụ nữ bị chứng khó thở do
thuốc nhuộm đã dùng thuốc giảm đau được chỉ định để làm giảm tác hại.
Một thuốc giảm đau mới được tuyên bố hiệu quả hơn. Sau khi thử thuốc
giảm đau mới, 40 phụ nữ cho thấy giảm đau hơn với thuốc đã chỉ định, và
60 người cho biết có giảm đau tốt hơn với loại thuốc mới.a. Hãy kiểm
định giả thuyết cho rằng khả năng làm giảm đau của loại thuốcđã chỉ định
và của loại thuốc mới là như nhau.b. Xây dựng và giải thích khoảng tin
cậy 95% cho xác suất giảm đau tốt hơn của thuốc giảm đau mới.
1.9. Trở lại các bài toán trước. Các nhà nghiên cứu
muốn có một mẫu đủ lớn để có thể ước lượng xác suất lựa chọn thuốc giảm
đau mới với độ chính xác 0.08, với độ tin cậy 0.95.
1.10. Tạp chí Newsweek (27/3/1989) báo cáo kết quả
cuộc thăm dò ý kiến tôn giáo, do Tổ chức Gallup tiến hành. Trong số 750
người Mỹ trưởng thành, 24% tin vào luân hồi. Xem đây như là một mẫu ngẫu
nhiên, xây dựng và giải thích một khoảng tin cậy 95% cho tỷ lệ thực sự
của người lớn Mỹ tin tưởng vào luân hồi.
1.11. Một nhà tội phạm học muốn ước tính tỷ lệ công
dân Mỹ sống trong một căn nhà có súng ống. Cuộc Khảo sát Xã hội năm 1991
đã hỏi những người trả lời, “Trong nhà bạn ở có bất kỳ súng hoặc súng
lục ổ quay nào?” Trong số những người được hỏi, 393 trả lời “có” và 583
trả lời “không”. Xây dựng khoảng tin cậy 90% cho tỷ lệ thực của “có”.
Giải thích.
1.12. Giả sử \(Y\)
là một biến ngẫu nhiên có phân phối nhị thức B(n, p), và \(f=Y/n\) . Chỉ ra rằng \(σ(f)=\sqrt{p(1-p)/n}\) .
1.13. Một biến ngẫu nhiên có phân phối Poisson, với
tham số chưa biết \(μ\) . Quan sát duy
nhất bằng 0.
a. Tìm và vẽ các hàm hợp lý trên không gian các giá trị có thể có của
các giá trị cho \(μ\) .
b. Ước lượng ML của \(μ\) là gì?
1.14. Tìm ước lượng ML cho tham số 𝜆 của phân phối
mũ.
1.15. Chỉ ra một giá trị \(p_0\) với thống kê \(z=(f-p_0)/\sqrt{p_0(1-p_0)/n}\) nhận một
vài giá trị cố định là lời giải của phương trình \((1+z_0^2 /n)\pi_0^2+(-2f-z_0^2/n)p_0+f^2\)
. Từ đó, sử dụng công thức \(x=-b±\sqrt{b^2-4ac}/2a\) để giải phương
trình bậc hai \(ax^2+bx+c=0\), lấy các
giới hạn cho khoảng tin cậy 95% trong ví dụ 3 cho tỷ lệ Người Mỹ ủng hộ
phá thai hợp pháp.
Chương 2: BẢNG NGẪU NHIÊN VÀ SUY DIỄN THỐNG KÊ
Phân tích các mối quan hệ trong một véc tơ quan sát là trọng tâm của
hầu hết các phân tích thống kê đa biến. Chương này đề cập đến mối quan
hệ giữa hai biến định tính và phần cuối đề cập đến quan hệ của nhiều
biến định tính. Trong đó giới thiệu các tham số mô tả sự kết hợp và đưa
ra các phương pháp suy diễn cho các tham số đó.
Nhiều ứng dụng liên quan đến việc so sánh giữa hai nhóm. Phần 2.1 sẽ
trình bày cấu trúc xác suất cho bảng ngẫu nhiên. Phần 2.2 sẽ trình bày
các phương pháp phân tích sự khác biệt về tỷ lệ. Phần 2.3 trình bày về
độ đo odd ratio, đóng vai trò rất quan trọng đối với một số phương pháp
được đề cập trong nội dung học phần này. Phần 2.4 và 2.5 giới thiệu một
số phương pháp phân tích để trả lời câu hỏi có tồn tại hay không sự liên
kết giữa hai biến định tính. Phần 2.6 là mở rộng cho bảng ngẫu nhiên
nhiều chiều.
Cấu trúc xác suất cho các bảng ngẫu nhiên
Dữ liệu định tính là tần số suất hiện các biểu hiện của các biến. Cho
X và Y là hai biến định tính, X có \(k\) biểu hiện: \(A_1, A_2,...,A_k\) và Y có \(m\) biểu hiện: \(B_1, B_2, ...,B_m\). Chúng ta có thể sử
dụng một bảng gồm \(k\) hàng và \(m\) cột để thể hiện kết quả có thể xảy ra
từ việc khảo sát:
Bảng 2.1
| \(A_1\) |
\(n_{11}\) |
\(n_{12}\) |
… |
\(n_{1m}\) |
| \(A_2\) |
\(n_{21}\) |
\(n_{22}\) |
… |
\(n_{2m}\) |
| … |
… |
… |
… |
… |
| \(A_k\) |
\(n_{k1}\) |
\(n_{k2}\) |
… |
\(n_{km}\) |
Bảng này được gọi là bảng ngẫu nhiên hai chiều \(k×m\) trong đó \(n_{ij}\) là số lần quan sát được cặp thuộc
tính \((A_i, B_j)\), còn gọi là tần số
của \((A_i, B_j)\). Một bảng ngẫu nhiên
hai biến được gọi là bảng hai chiều; một bảng ngẫu nhiên 3 biến gọi là
bảng ngẫu nhiên 3 chiều. Thường thì khi trình bày bảng dữ liệu 2 chiều,
người ta không kẻ đường phân chia cột mà chỉ kẻ một số đường ngang cần
thiết.
Ví dụ 1: Khảo sát 1.091 người Mỹ về niềm tin của họ đối với
thế giới bên kia (sau khi chết): 435 nữ tin và 147 nữ không tin hoặc
không đưa ra quyết định; 375 tin. Vấn đề là có tồn tại mối liên hệ giữa
giữa giới tính và niềm tin vào thế giới bên kia không?
Trong mối quan hệ này, biến Giới tính (có hai thuộc tính: “Nam” và
“Nữ”) là biến giải thích, biến “Niềm tin vào thế giới bên kia” (có 2
thuộc tính: “Tin” và “Không tin”) là biến đáp ứng (hay biến phụ thuộc).
Dữ liệu nói trên được trình bày bởi bảng ngẫu nhiên 2 chiều 2×2 theo một
trong hai dạng sau:
Nguồn: Data from 1991 General Social Survery
Bảng 2.2. Kết quả điều tra về niềm tin vào thế giới bên kia theo giới
tính
| Giới tính |
Niềm tin vào thế giới bên kia |
|
|
Tin |
Không tin hoặc không quyết định |
Nữ
Nam |
435
375 |
147
134 |
Hoặc
Bảng 2.3
| Nữ |
435 |
147 |
| Nam |
375 |
134 |
Phân phối đồng thời, phân phối biên duyên và xác suất có
điều kiện
Trước tiên chúng ta trình bày một số khái niện cơ bản của xác suất
cho các bảng ngẫu nhiên hai chiều. Giả sử rằng mỗi quan sát trong mẫu
được chọn ngẫu nhiên từ một tổng thể và mỗi quan sát gắn liền với 2 biến
ngẫu nhiên nhị phân X và Y. Đặt \(\pi_{ij}=P(X=A_i,Y=B_j)\) là xác suất để
(X,Y) nhận giá trị \((A_i, B_j)\).
\(\pi_{ij}\) gọi là phân phối xác suất
đồng thời của X và Y và thỏa \(\sum_{ij}\pi_{ij}=1\).
Phân phối xác suất biên duyên là tổng theo từng hàng và từng cột của
bảng phân phối xác suất đồng thời của X và Y và được ký hiệu là {\({\pi_{i+}}\)} và {\({\pi_{+j}}\)}, dấu + kí hiệu cho việc tổng
theo hàng hoặc cột tùy vào vị trí của nó nằm ở vị trí hàng hay cột,
nghĩa là: \(\pi_{i+}=\pi_{i1}+\pi_{i2}\) và \(\pi_{+j}=\pi_{1j}+\pi_{2j}\)
Phân phối đồng thời và phân phối biên duyên của các biến nhị phân X,
Y được thể hiện qua bảng sau:
Bảng 2.4
| \(A_1\) |
\(\pi_{11}\) |
\(\pi_{12}\) |
\(\pi_{1+}\) |
| \(A_1\) |
\(\pi_{21}\) |
\(\pi_{22}\) |
\(\pi_{2+}\) |
| \(\pi_{+j}\) |
\(\pi_{+1}\) |
\(\pi_{+2}\) |
1 |
Tương tự chúng ta ký hiệu \(n_{ij}\)
là số quan sát có biểu hiện \(A_i\) đối
với X và \(B_j\) đối với Y, ký hiệu
\(f_{ij} = n_{ij}/n\) với \(n =\sum n_{ij}\) . Bảng sau đây gọi là bảng
phân phối mẫu đồng thời và phân phối mẫu biên duyên của các biến nhị
phân X, Y:
Bảng 2.5
| \(A_1\) |
\(f_{11}\) |
\(f_{12}\) |
\(f_{1+}\) |
| \(A_1\) |
\(f_{21}\) |
\(f_{22}\) |
\(f_{2+}\) |
| \(f_{+j}\) |
\(f_{+1}\) |
\(f_{+2}\) |
1 |
Trong nhiều bảng ngẫu nhiên, một biến (thường là biến cột, Y) là biến
đáp ứng và biến khác (biến hàng, X) là biến giải thích, với mỗi giá trị
của X là thông tin để xây dựng một phân phối xác suất riêng cho Y, phân
phối như vậy gọi là phân phối xác suất có điều kiện của Y tương ứng với
từng giá trị của X.
Với hai biến nhị phân của \(X, Y\)
nói trên:
- Phân phối có điều kiện của Y với điều kiện \(X = A_1\) là:
| \(P\) |
\(\pi_{11}/\pi_{1+}\) |
\(\pi_{12}/\pi_{1+}\) |
1 |
-Phân phối mẫu có điều kiện của \(Y\) với điều kiện \(X=A_1\) là:
| \(P\) |
\(f_{11}/f_{1+}\) |
\(f_{12}/f_{1+}\) |
1 |
-Phân phối có điều kiện của Y với điều kiện \(X=A_2\) là:
| \(P\) |
\(\pi_{21}/\pi_{2+}\) |
\(\pi_{22}/\pi_{2+}\) |
1 |
-Phân phối mẫu có điều kiện của \(Y\) với điều kiện \(X=A_2\) là:
| \(P\) |
\(f_{21}/f_{2+}\) |
\(f_{22}/f_{2+}\) |
1 |
Ví dụ 2: Từ bảng 2.2, là một bảng ngẫu nhiên 2 chiều với
tổng số quan sát là 1091 được định tính theo giới tính X và niềm tin Y
của họ sau khi chết, ta có bảng:
Bảng 2.6
| Nữ |
435 |
147 |
\(n_{1+}=582\) |
| Nam |
375 |
134 |
\(n_{2+}=509\) |
| \(n_{+j}\) |
\(n_{+1}=810\) |
\(n_{+2}=281\) |
\(n=1091\) |
Tính theo công thức: \(f_{ij}=\frac{n_{ij}}{n}\), ta nhận được một
thể hiện của phân phối mẫu đồng thời và phân phối mẫu biên duyên của các
biến X (giới tính: Nam, Nữ) và Y (niềm tin: Có, Không):
Bảng 2.7
| Nữ |
\(f_{11}=0,3987\) |
\(f_{12}=0,1348\) |
\(f_{1+}=0,5335\) |
| Nam |
\(f_{21}=0,3437\) |
\(f_{22}=0,1228\) |
\(f_{2+}=0,4665\) |
| \(f_{+j}\) |
\(f_{+1}=0,7424\) |
\(f_{+2}=0,2576\) |
1 |
Chúng ta sẽ tìm phân phối mẫu có điều kiện của niềm tin vào thế giới
bên kia với giới tính.
- Đối với phụ nữ, tỷ lệ trả lời “có tin vào thế giới bên kia” là
\(435/582 = 0,7473\); và trả lời “không
tin vào thế giới bên kia” là \(147/582 =
0,2527\).
- Đối với đàn ông, tỷ lệ trả lời “có tin vào thế giới bên kia” là
\(375/509 = 0,7367\); và trả lời “không
tin vào thế giới bên kia” là \(147/582 =
0,2633\).
Sự độc lập
Hai biến được cho là độc lập về mặt thống kê nếu các phân phối có
điều kiện của \(Y\) là không thay đổi
khi \(X\) thay đổi giá trị. Khi hai
biến độc lập, xác suất của bất kỳ phản ứng nào của cột \(j\) giống nhau ở mỗi hàng. Ví dụ, niềm tin
vào một thế giới bên kia là không phụ thuộc vào giới tính nếu xác suất
vào thế giới bên kia của cả nam và nữ là 0,7367.
Khi cả hai biến là các biến đáp ứng, người ta có thể mô tả mối quan
hệ của chúng bằng cách sử dụng phân phối đồng thời của chúng, hoặc phân
phối có điều kiện của \(Y\) theo \(X\) hoặc phân phối có điều kiện của \(X\) theo \(Y\). Tính độc lập thống kê sẽ tương đương
với tính chất xác suất đồng thời bằng với tích các xác suất biên nghĩa
là
\[
\pi_{ij}=\pi_{i+}.\pi_{+j},∀i,j
\]
Phân phối đa thức
Mô hình phân phối Nhị thức và mô hình phân phối Poisson đã được trình
bày trong chương 1. Trong phần này sẽ giới thiệu về phân phối đa thức.
Theo một nghĩa nào đó thì phân phối đa thức là sự mở rộng của phân phối
Nhị thức.
Mô hình phân phối đa thức
Xét một dãy \(n\) phép thử độc lập ,
trong mỗi phép thử có một và chỉ một trong \(k\) sự kiện \(A_1,A_2,...,A_k\) xảy ra với xác suất tương
ứng \(p_1, p_2, ..., p_k,
(p_1+p_2+...+p_k)\) (tức là trong mỗi phép thử, \(k\) biến cố này lập thành một hệ đầy đủ).
Gọi \(X_i\) là số lần xuất hiện sự kiện
\(A_i (i=1, 2, ..., k)\) trong \(n\) lần thử. Khi đó luật phân phối của véc
tơ ngẫu nhiên \(X=(X_1, X_2, ...,
X_k)\) được gọi là phân phối đa thức \(k\) chiều với các tham số \(n, p_1, p_2,..., p_{k-1}\)
Phân phối xác suất của véc tơ ngẫu nhiên có phân phối đa thức
k chiều.
Giả sử \(X = (X_1, X_2,… , X_k)\) là
véc tơ ngẫu nhiên có phân phối đa thức k chiều. Ta biết rằng mỗi thành
phần \(X_i ~ B(n,p_i) (i=1,2,...,k)\)
và: \(X_1+X_2+...+X_k =n\), vì thế mà
véc tơ \(X = (X_1, X_2,… , X_k)\) là
véc tơ ngẫu nhiên rời rạc k chiều, mỗi giá trị của nó là một điểm k
chiều có tọa độ nguyên không âm: \((r_1, r_2,
..., r_k)\), sao cho: \(r_1+r_2+...+r_k
=n\). Mỗi kết cục thuận lợi cho biến cố (\(X_1=r_1, X_2=r_2,...,X_k =r_k)\) là một dãy
gồm n sự kiện liên kết với nhau bởi phép giao, trong đó: sự kiện \(A_1\) xuất hiện \(r_1\) lần, sự kiện \(A_2\) xuất hiện \(r_2\) lần, ,…, sự kiện \(A_k\) xuất hiện \(r_k\) lần. Như vậy số kết cục thuận lợi cho
biến cố \((X_1=r_1, X_2=r_2,...,X_k
=r_k)\) chính là số cách chọn \(r_1\) vị trí cho \(A_1, r_2\) vị trí cho \(A_2, \ldots, r_k\) vị trí cho \(A_k\) trong dãy \(\mathrm{n}\) vị trí nói trên, tức là bằng:
\[
C_n^{r_1} \cdot C_{n-r_1}^{r_2} \cdot C_{n-r_1-r_2}^{r_3} \ldots .
C_{n-r_1-r_2-\cdots r_{k-2}}^{r_{k-1}}=\frac{n!}{r_{1}!r_{2}!\ldots
r_{k}!}
\]
Mặt khác do các lần thử độc lập nhau nên mỗi kết cục này đều có xác
suất là:
\(p_1^{r_1} \cdot p_2^{r_2} \ldots . .
p_k^{r_k}\)
Từ đó ta nhận được phân phối xác suất của véc tơ ngẫu nhiên có phân
phối đa thức \(k\) chiều:
\(P\left(X_1=r_1, X_2=r_2, \ldots,
X_k=r_k\right)=\frac{n!}{r_{1}!r_{2}!\ldots r_{k}!} \cdot p_1^{r_1}
\cdot p_2^{r_2} \ldots . . p_k^{r_k}\) (2.1.1)
với mọi điểm \(\left(r_1, r_2, \ldots,
r_k\right)\) có tọa độ nguyên, không âm mà: \(r_1+r_2+\cdots+r_k=n\).
Nhận xét:
- Đối với véc tơ ngẫu nhiên \(X=\left(X_1,
X_2, \ldots, X_k\right)\) có phân phối đa thức \(k\) chiều, với các tham số \(n, p_1, p_2, \ldots, p_{k-1}\), ta có véc
tơ kỳ vọng:
\(E X=\left(E X_1, E X_2, \ldots, E
X_k\right)=\left(n p_1, n p_2, \ldots, n p_k\right)\)
- Phân phối đa thức với \(k=2\), là
phân phối đồng thời của 2 biến nhị thức \(X_1\), \(n-X_1\), trong đó \(X_1 \sim \mathrm{B}\left(\mathrm{n},
p_1\right)\), \(\left(\mathrm{n}-X_1\right) \sim
\mathrm{B}\left(\mathrm{n}, 1-\mathrm{p}_1\right)\).
Ví dụ 3. Mỗi khách hàng vào hệ dịch vụ A có thể chọn một
trong 3 mức phí phục vụ: 100 ngàn đồng, 150 ngàn đồng và 200 ngàn đồng.
Được biết lượng khách chọn các mức phí này tương ứng theo tỷ lệ: 5: 3:
2. Có 3 khách hàng vào hệ dịch vụ A và họ độc lập nhau trong việc chọn
phí dịch vụ. Tìm xác suất để trong số đó có ít nhất 2 khách chọn mức phí
200 ngàn đồng.
Giải: Gọi \(X_1, X_2, X_3\) lần lượt
là số khách chọn mức phí 100 ngàn đồng, 150 ngàn đồng, 200 ngàn đồng
trong số 3 khách nói trên. Khi đó véc tơ \(X=(X_1, X_2, X_3)\) là véc tơ ngẫu nhiên có
phân phối đa thức 3 chiều. Ký hiệu S là biến cố trong số 3 khách hàng có
ít nhất 2 khách chọn mức phí 200 ngàn đồng, ta có biểu diễn:
\(S=(X_1=0, X_2=1, X_3=2) \cup(X_1=1,
X_2=0, X_3=2) \cup(X_1=0, X_2=0, X_3=3)\)
Từ đó suy ra xác suất cần tính:
\(P(S)=\frac{3!}{0!1!2!} \cdot 0,3 \cdot
0,2^2+\frac{3!}{1!0!2!} \cdot 0,5 \cdot 0,2^2+\frac{3!}{0!0!3!} \cdot
0,2^3=0,\)
Lấy mẫu cho phân phối Poisson, Nhị thức và Đa
thức
Các mô hình lấy mẫu được giới thiệu trong phần 1.3 được mở rộng ra
cho tần số trong các ô của bảng ngẫu nhiên. Ví dụ, mô hình lấy mẫu
Poisson cho bảng hai chiều \(2 \times
2\), các ô trong bảng được xem như là các biến Poisson độc
lập.
Khi các hàng của một bảng ngẫu nhiên đại diện cho các nhóm khác nhau,
cỡ mẫu cho các nhóm này thường đã được cố định trước trong giai đoạn
thiết kế mẫu. Để đơn giản, chúng ta sẽ trình bày về điều này cho bảng
\(2 \times 2\) cho hai biến \(X\) (biến giải thích) và \(Y\) (biến đáp ứng). Khi tổng số cận biên
cho các biểu hiện của \(X\) được cố
định chứ không phải ngẫu nhiên, một phân phối chung cho \(X\) và \(Y\) không còn ý nghĩa nữa, nhưng phân phối
có điều kiện cho \(Y\) ở mỗi biểu hiện
của \(X\) thì có ý nghĩa.
Khi tổng cỡ mẫu trong bảng được cố định nhưng tổng của hàng hoặc cột
không cố định, mô hình lấy mẫu đa thức là phù hợp, trong đó các ô là
những kết quả có thể xảy ra. Chẳng hạn, bảng 2.1 định tính một mẫu ngẫu
nhiên của 1.091 đối tượng theo giới tính và niềm tin vào thế giới bên
kia. Giá trị trong 4 ô là tần số của phân phối đa thức bốn chiều.
Đối với các mẫu của phân phối đa thức trên các ô của một bảng ngẫu
nhiên, các cột là biến đáp ứng và các hàng là biến giải thích. Sau đó,
để mô tả dữ liệu, một cách hợp lý, ta phân chia các ô đếm theo tổng số
hàng để tạo thành phân phối có điều kiện về đáp ứng. Khi làm như vậy,
chúng ta vốn đã xử lý các tổng số hàng là cố định và phân tích dữ liệu
giống như khi chúng hình thành các mẫu độc lập riêng biệt. Ví dụ như
trong Bảng 2.1, chúng ta có thể coi kết quả cho nữ giới như một mẫu nhị
thức với các loại kết quả là “có” và “không hoặc không quyết định” cho
niềm tin vào một thế giới bên kia, và kết quả cho nam giới như một mẫu
nhị thức độc lập về cùng một phản ứng. Nếu có nhiều hơn hai loại phản
hồi, chẳng hạn như (“có”, “không”, “không quyết định”), chúng ta sẽ coi
các mẫu là mẫu đa thức độc lập.
Đối với hầu hết các phân tích, người ta không cần quan tâm về mô hình
lấy mẫu nào có ý nghĩa nhất. Đối với các phương pháp suy diễn chính
trong tài liệu này, các mô hình lấy mẫu là Poisson, đa thức.
So sánh tỷ lệ trong một bảng hai chiều
Biến đáp ứng có hai biểu hiện được gọi là biến nhị phân. Ví dụ, “niềm
tin vào thế giới bên kia” là nhị phân khi được đo bằng các biểu hiện
(có, không). Nhiều nghiên cứu cần so sánh giữa hai nhóm đối với đáp ứng
nhị phân của biến \(Y\). Dữ liệu có thể
được hiển thị trong bảng ngẫu nhiên \(2 \times
2\), trong đó các hàng là hai nhóm và các cột là mức phản ứng của
\(Y\). Phần này trình bày các phương
pháp để so sánh mức độ phản ứng trên biển nhị phân.
Sự khác biệt giữa hai tỷ lệ
Với biến đáp ứng nhị phân Y, chúng ta sử dụng thuật ngữ chung “thành
công” cho một đáp ứng và “thất bại” đối với đáp ứng còn lại. Đối với các
đối tượng trong hàng i, ký hiệu \(\pi_i\) là xác suất thành công thì \(1-\pi_i\) là xác suất của một thất bại. Cặp
xác suất \(\left(\pi_i,
1-\pi_i\right)\) là phân phối xác suất có điều kiện của Y trong
hàng i.
Định nghĩa 1: Khác biệt của 2 tỷ lệ: \(\pi_1-\pi_2\) so sánh xác suất thành công
của hai hàng.
Ta có: \(-1 \leq \pi_1-\pi_2 \leq
1\) và \(\pi_1-\pi_2=0 \Leftrightarrow
\pi_1=\pi_2\).
TRANG 28
Ký hiệu \(f_i\) là tỷ lệ mẫu của
biểu hiện “thành công” ở hàng thứ \(i(i=1,2)\), sự khác biệt \(\pi_1-\pi_2\) giữa hai tỷ lệ sẽ được ước
lượng bởi \(f_1-f_2\).
Và khi cỡ mẫu n khá lớn, chúng ta có công thức ước lượng khoảng tin
cậy với độ tin cậy \(\gamma=1-\alpha\),
cho sự sai khác về tỷ lệ: \(\pi_1-\pi_2\) là:
\(((f_1-f_2)-u_{\frac{\alpha}{2}}
\hat{\sigma}(f_1-f_2),(f_1-f_2)+u_{\frac{\alpha}{2}}
\hat{\sigma}(f_1-f_2))\) (2.2.1)
với: \(\quad
\hat{\sigma}\left(f_1-f_2\right)=\sqrt{\frac{f_1\left(1-f_1\right)}{n_{1+}}+\frac{f_2\left(1-f_2\right)}{n_{2+}}}\)
(2.2.2)
và \(u_{\alpha / 2}\) tra trong bảng
phân phối chuẩn chính tắc theo hệ thức: \(\Phi(u_\lambda)=1-\lambda\).(\(u_\lambda\) được gọi là giá trị tới hạn mức
\(\lambda\) của phân phối chuẩn chính
tắc).
Ví du 4: Dựa vào điều tra ở bảng 2.2, cần ước lượng khoảng
tin cậy 95% cho sai khác về tỷ lệ tin vào thế giới bên kia giữa phụ nữ
và đàn ông Mỹ.
Với \(\gamma=1-\alpha=0,95\), có
\(\alpha=0,05,
u_{\frac{\alpha}{2}}=u_{0,025}=1,96\). Từ bảng 2.2, có:
\(f_1=\frac{n_{11}}{n_{1+}}=0,7473 ;
f_2=\frac{n_{21}}{n_{2+}}=0,7367 ;
\hat{\sigma}(f_1-f_2)=0,0266;\)
\((f_1-f_2)-u_{\alpha/2}\hat{\sigma}(f_1-f_2)=-0,0415
;(f_1-f_2)+u_{\alpha / 2} \hat{\sigma}(f_1-f_2)=0,0627\)
Vậy khoảng tin cậy \(95 \%\) cho sai
lệch này là: \((-0,0415 ;
0,0627)\).
Ví du 5: Bảng 2.8 là số liệu báo cáo về mối quan hệ giữa sử
dụng thuốc aspirin và bệnh nhồi máu cơ tim (myocardial
infarction) do nhóm nghiên cứu sức khoẻ của các bác sĩ tại Trường y
Harvard. Nghiên cứu được thực hiện trong thời gian 5 năm trên những
người bị bệnh tim được chọn một cách ngẫu nhiên và chia làm hai nhóm:
Nhóm dùng aspirin và nhóm dùng giả dược (placebo). Mỗi ngày,
các bác sĩ sẽ cho bệnh nhân uống một viên aspirin hoặc một viên giả
dược.
|
Myocardial Infarction |
|
|
| Group |
Yes |
No |
Total |
| Placebo |
189 |
10845 |
11034 |
| Aspirin |
104 |
10933 |
11037 |
Nguồn: Preliminary Report: Findings from the Aspirin Component of
the Ongoing Physicians’ Health Study. N. Engl. J. Med., 318: 262-264
(1988)
Bảng 2.8. Cross Classification of Aspirin Use and Myocardial
Infarction (MI)
Chúng ta xem hai hàng trong Bảng 2.8 như là các mẫu nhị phân độc lập,
trong hàng 1 : \(\mathrm{n}_{1+}=11034\) người và có 189
người bị nhồi máu cơ tim tỷ lệ là \(f_1= 189 /
11034=0,171\), trong hàng 2 : \(\mathrm{n}_{2+}=11037\) người và có 104 bị
nhồi máu cơ tim tỷ lệ là \(f_2=104 /
11037=0,0094\), và: \(u_{\frac{\alpha}{2}}
\hat{\sigma}\left(f_1-f_2\right)=1,96 .(0,015)=0,0294\). Áp dụng
công thức trên với độ tin cậy 95% chúng ta khoảng tin cậy 95% cho \(\pi_1-\pi_2\) là: \((0,005;0,011)\), hay \(0,005<\pi_1-\pi_2<0,011\left(0,5
\%<\pi_1-\pi_2<1,1 \%\right)\) nghĩa là giữa người uống và
không uống thuốc asparin thì tỷ lệ bị nhồi máu cơ tim có sự khác biệt từ
0,5 đến 1,1%.
Rủi ro tương đối (relative risk)
Sự khác biệt giữa hai tỷ lệ \(\pi_1\) và \(\pi_2\) được đánh giá qua sai số tuyệt đối
\(\left(\pi_1-\right.\) \(\pi_2\) ) chỉ dựa vào khoảng cách giữa hai
tỷ lệ này mà bỏ qua độ lớn và ý nghĩa của chúng, nghĩa là khoảng cách
như nhau thì đánh giá sai lệch như nhau. Chẳng hạn khi \(\pi_1=0,01\) và \(\pi_2=0,001\) và khi \(\pi_1=0,410\) và \(\pi_2=0,401\) đều được đánh giá là có sai
lệch như nhau (cùng mức sai lệch là 0,009 ). Tuy nhiên nếu chỉ dựa vào
khoảng cách này để đánh giá, so sánh hai tương quan giữa hai tỷ lệ thì
không đủ, thậm chí sẽ bỏ qua những ý nghĩa quan trọng của sự tương quan
của chúng, đặc biệt là khi hai tỷ lệ này cùng gần 0 hoặc cùng gần 1 ,
hoặc cùng gần 0,5 .
Ví dụ khi chúng ta so sánh \(\pi_1\)
và \(\pi_2\) là tỷ lệ bị tác dụng phụ
(phản ứng bất lợi) của hai loại thuốc, trong hai trường hợp:
- Trường hợp 1: Sự khác biệt nằm trong khoảng từ 0,001 đến 0,01 .
- Trường hợp 2: Sự khác biệt nằm trong khoảng từ 0,401 đến 0,410
.
Cả hai trường hợp này, sự khác biệt về tỷ lệ bị tác dụng phụ là như
nhau (đều là 0,009\()\). Nhưng nếu ta
để ý đến tỷ số \(\frac{\pi_1}{\pi_2}\),
thì ở trường hợp \(1:
\frac{\pi_1}{\pi_2}=10\), cho thấy khi bệnh nhân điều trị, phản
ứng bất lợi của loại thuốc thứ nhất gấp 10 lần loại thuốc thứ hai, còn ở
trường hợp 2: \(\frac{\pi_1}{\pi_2}=1,0224\), cho thấy khi
bệnh nhân điều trị, phản ứng bất lợi của loại thuốc thứ nhất gấp 1,0224
lần loại thuốc thứ hai, tức là gần như nhau. Rõ ràng sự khác biệt ở
trường hợp 1 là đáng chú ý hơn.
Trong những trường hợp như vậy, tỉ số giữa các tỷ lệ cũng là một
thước đo hữu ích.
Định nghĩa 2: Giả sử \(\pi_1\) và \(\pi_2\) là xác suât “thành công” của hai
nhóm. Khi đó rủi ro tương đối (relative risk) giữa hai nhóm này là tỷ lệ
xác suất “thành công” của hai nhóm: \[
\frac{\pi_1}{\pi_2}
\](2.2.3)
Chẳng hạn: Trong ví dụ trên tỷ lệ về phản ứng bất lợi của loại thuốc
thứ nhất và loại thuốc thứ hai, trong trường hợp 1 có rủi ro tương đối
là 10,00 và trong trường hợp 2 có rủi ro tương đối là 1,02 .
Nhận xét:
- Rủi ro tương đối là số không âm.
- Rủi ro tương đối bằng 1 khi và chỉ khi biến đáp ứng độc lập theo
nhóm.
- Các tỷ lệ tổng thể \(\pi_1\) và
\(\pi_2\) nói chung là chưa biết, mà
qua điều tra ta nhận được các tỷ lệ mẫu tương ứng là \(f_1\) và \(f_2\). Khi đó ta có tỷ lệ: \(\frac{f_1}{f_2}\) gọi là rủi ro tương đối
mẫu, ta dùng nó để xấp xỉ cho rủi ro tương đối \(\frac{\pi_1}{\pi_2}\) trên tổng thể.
Người ta chỉ ra rằng phân phối mẫu rủi ro tương đối có thể bị lệch
rất lớn trừ khi kích thước mẫu đủ lớn. Do đó công thức ước lượng khoảng
tin cậy của nó khá phức tạp, người ta phải thông qua ước lượng khoảng
tin cậy cho \(\log
\left(\frac{\pi_1}{\pi_2}\right)\) là:
\[
\log \left(\frac{f_1}{f_2}\right) \pm u(\alpha / 2)
\sqrt{\frac{1-f_1}{n_{1+} \cdot f_1}+\frac{1-f_2}{n_{2+} \cdot f_2}}
\](2.2.4)
Ví du 7: Trong Bảng 2.3, có \(f_1=0,0171, f_2=0,0094, n_{1+}=\) \(11034, n_{2+}=11037\)
- Rủi ro tương đối mẫu là \(\mathrm{f}_1 /
\mathrm{f}_2=0,0171 / 0,0094=1,82\), cho thấy: tỷ lệ bị bệnh nhồi
máu cơ tim đối với người sử dụng thuốc asporin cao hơn \(82 \%\) so với nhóm sử dụng giả dược.
- Sử dụng công thức (2.2.4), ta có khoảng tin cậy \(95 \%\) cho \(\log \left(\frac{\pi_1}{\pi_2}\right)\) là:
\((0,3607 ; 0,8329)\). Suy ra khoảng
tin cậy cho rủi ro tương đối \(\frac{\pi_1}{\pi_2}\) là: \((1,43 ; 2,30)\)
Như vậy, chúng ta có thể tin tuởng \(95
\%\) rằng, sau năm năm, nếu sử dụng thuốc asparin thì nguy cơ bị
nhồi máu cơ tim cao gấp từ 1,43 đến 2,30 lần so với người không sử dụng
thuốc asparin.
Qua đó thấy rằng: khoảng tin cậy của rủi ro tương đối cho thấy nguy
cơ bị nhồi máu cơ tim cao hơn ít nhất \(43
\%\) đối với nhóm dùng giả dược. Khoảng tin cậy cho sự khác biệt
về tỷ lệ là \((0,005 ; 0,011)\) có vẻ
như hai nhóm khác rất nhỏ, nhưng nguy cơ tương đối cho thấy sự khác biệt
có thể có ý nghĩa quan trọng đối với sức khoẻ cộng đồng. Nếu chúng ta
chỉ sử dụng sự chênh lệch tỷ lệ để so sánh hai nhóm có thể phần nào gây
nhầm lẫn khi tỷ lệ của cả hai nhóm đều rất gần không.
Nhiều khi, thay vì khảo sát xác suất “thành công”, người ta cần khảo
sát thông tin để tính toán cho tỷ lệ của xác suất “thất bại” \(\left(1-\pi_1\right)
/\left(1-\pi_2\right)\).
Tỷ lệ chênh (Odds ratio)
Trong phần này chúng ta trình bày một công cụ khác để đánh giá cho
bảng ngẫu nhiên \(2 \times 2\), được
gọi là tỷ lệ chênh (odds ratio). Đây là một tham số cơ bản được sử dụng
trong các mô hình sẽ được trình bày ở các phần sau của tài liệu bài
giảng này.
Định nghĩa 3: Trong hàng \(i\) của bảng phân phối đồng thời \(2 \times 2\) của các biến quan sát \(X\) và \(Y\), tỷ lệ cược (odds) của “thành công”
được định nghĩa là tỷ số giữa xác suất “thành công” và xác suất “không
thành công”, ký hiệu odds \(s_i\) :
\[
o d d s_i=\frac{\pi_i}{1-\pi_i}
\](2.3.1)
Nhận xét:
- odds \(s_i\) là tỷ lệ giữa khả
năng “thành công” và khả năng “không thành công” trong cùng một điều
kiện tác động của biến giải thích.
- Odds luôn không âm: Odds = 1 có nghĩa là khả năng “thành công” và
“không thành công” là ngang nhau; Odds > 1 có nghĩa là khả năng
“thành công” cao hơn khả năng “không thành công”
- Khi biết odds, ta có thể tính được xác suất thành công, vì ta
có:
LS0tDQp0aXRsZTogIkdpYcyBbyB0cmnMgG5oIHRyYW5nIDE2LTMwIg0KZGF0ZTogImByIGZvcm1hdChTeXMudGltZSgpLCAnJUg6JU06JVMsICVkIC0gJW0gLSAlWScpYCINCm91dHB1dDoNCiAgaHRtbF9kb2N1bWVudDoNCiAgICB0b2M6IFRSVUUNCiAgICB0b2NfZmxvYXQ6IFRSVUUNCiAgICBkZl9wcmludDogcGFnZWQNCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlDQogICAgY29kZV9mb2xkaW5nOiBoaWRlDQogIHBkZl9kb2N1bWVudDoNCiAgICBleHRyYV9kZXBlbmRlbmNpZXM6DQogICAgICB2aWV0bmFtOiB1dGY4DQogICAgdG9jOiB5ZXMNCiAgICBudW1iZXJfc2VjdGlvbnM6IHllcw0KICB3b3JkX2RvY3VtZW50Og0KICAgIHRvYzogeWVzDQogICAgbnVtYmVyX3NlY3Rpb25zOiB5ZXMNCmdlb21ldHJ5Og0KICAgICAgLSBpbm5lcj0zY20NCiAgICAgIC0gb3V0ZXI9NGNtDQogICAgICAtIHRvcD0zY20NCiAgICAgIC0gYm90dG9tPTRjbQ0KICAgICAgLSBoZWFkc2VwPTIycHQNCiAgICAgIC0gaGVhZGhlaWdodD0xMXB0DQogICAgICAtIGZvb3Rza2lwPTMzcHQNCiAgICAgIC0gaWdub3JlaGVhZA0KICAgICAgLSBpZ25vcmVmb290DQogICAgICAtIGhlaWdodHJvdW5kZWQNCi0tLQ0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0NCmtuaXRyOjpvcHRzX2NodW5rJHNldChlY2hvID0gVFJVRSwgb3B0aW9ucyhkaWdpdHMgPSA0KSwgIGF0dHIuc291cmNlPScubnVtYmVyTGluZXMnKQ0KYGBgDQoNCmtp4buDbSDEkeG7i25oIGdp4bqjIHRodXnhur90ICRIXzA6IHA9MCw1JCwgduG7m2kgxJHhu5FpIHRodXnhur90ICRIXzA6IHDiiaAwLDUkIG3hu6ljIMO9IG5naMSpYSA1JS4gVsOsIGPhu6EgbeG6q3UgbiA9IDE1IGzDoCBiw6ksIG7Dqm4gdGnDqnUgY2h14bqpbiBiw6FjIGLhu48gJEhfMCQgbMOgOg0KDQokZuKJpVxmcmFje1xmcmFje3ggXGFscGhhfXsyfX17bn0kLCBob8SDzKNjICRm4omkIFxmcmFje1xmcmFje1xiYXJ7eH0gXGFscGhhfXsyfX17bn0kLCB0xrDMgWMgbGHMgCAkZuKIiSAoXGZyYWN7XGZyYWN7XGJhcnt4fSBcYWxwaGF9ezJ9fXtufTtcZnJhY3tcZnJhY3t4IFxhbHBoYX17Mn19e259KSQNCg0KVOG7qyBt4bqrdSB0YSBjw7M6ICRuPTE1LCBmPVxmcmFjezh9ezE1fSQuIFbhu5tpICRcYWxwaGEgPSAwLDUkLCB04burIGLhuqNuZyBwaOG7pSBs4bulYyAxLCDEkWnhu4F1IGtp4buHbjoNCg0KJCQNClxzdW1fe3g9eF97MCwwMjV9fV57MTV9IENfezE1fV54LigwLDUpXnsxNX3iiaQgMCwwMjUNCiQkDQoNCmNobyB0YSAkeF97MCwwMjV9PTEyJCwgbsOqbiAkXGZyYWN7eF97XGZyYWN7XGFscGhhfXsyfX19e259PVxmcmFjezEyfXsxNX0kLCB2YcyAIHTGsMyAIMSRacOqzIB1IGtpw6rMo246DQoNCiQkDQpcc3VtX3t4PTB9XntcYmFye3h9X3swLDAyNX19IENfezE1fV54LigwLDUpXnsxNX3iiaQgMCwwMjUNCiQkDQoNCmNobyB0YSAkXGJhcnt4fV97MCwwMjV9PTMkLCBuw6puICRcZnJhY3tcYmFye3h9X3tcZnJhY3tcYWxwaGF9ezJ9fX17bn09XGZyYWN7M317MTV9JC4gVGEgY2/MgTogJGY9XGZyYWN7OH17MTV94oiIKFxmcmFjezN9ezE1fTtcZnJhY3sxMn17MTV9KSQgdGEgY2jhuqVwIG5o4bqtbiAkSF8wOnA9MCw1JCB2w6AgbmjGsCB24bqteSBsw6AgY2jGsGEgY8OzIGPGoSBz4bufIMSR4buDIGNobyBy4bqxbmcgdHJvbmcgc+G7kSBraMOhY2ggaMOgbmcgY+G7p2EgaOG7hyBk4buLY2ggduG7pSBBLCB04bu3IGzhu4cgY2jhu41uIGhhaSBt4bupYyBwaMOtIGzDoCBraMOhYyBuaGF1Lg0KDQojIyMgS2hv4bqjbmcgdGluIGPhuq15IGNobyBt4buZdCB04bu3IGzhu4cNCg0KVHJvbmcgaOG6p3UgaOG6v3QgY8OhYyDhu6luZyBk4bulbmcsIG5nxrDhu51pIHRhIHPhu60gZOG7pW5nIGtob+G6o25nIHRpbiBj4bqteSDEkeG7gyB4w6FjIMSR4buLbmggcGjhuqFtIHZpIGPhu6dhIGPDoWMgZ2nDoSB0cuG7iyBo4bujcCBsw70gY2hvIHRoYW0gc+G7kS4gQ2jDum5nIHRhIGJp4bq/dCBy4bqxbmcgxrDhu5tjIGzGsOG7o25nIE1MIGNobyB04bu3IGzhu4cgJFxtYXRocm17cH0kLCB04bu3IGzhu4cgdMOtbmggY2jhuqV0IEEgdHJvbmcgdOG7lW5nIHRo4buDIGzDoCAkZj1mKEEpJCwgdOG7tyBs4buHIHTDrW5oIGNo4bqldCBBIHRyw6puIG3huqt1IHbDoCBraGkgY+G7oSBt4bqrdSAkbiBccmlnaHRhcnJvdytcaW5mdHkkLCBwaMOibiBwaOG7kWkgY+G7p2EgdOG7tyBs4buHIG3huqt1ICRmJCBk4bqnbiB04bubaSBwaMOibiBwaOG7kWkgY2h14bqpbiwgaGF5IHBow6JuIHBo4buRaSBj4bunYSAkXGZyYWN7KGYtcCkgXHNxcnR7bn19e1xzcXJ0e3AoMS1wKX19JCBk4bqnbiB04bubaSBwaMOibiBwaOG7kWkgY2h14bqpbiBjaMOtbmggdOG6r2MuIFbhuq15IG7Dqm4ga2hpIGPhu6EgbeG6q3UgbiAixJHhu6cgbOG7m24iLCBuZ8aw4budaSB0YSBjw7MgdGjhu4MgeOG6pXAgeOG7iSBwaMOibiBwaOG7kWkgY+G7p2EgJFxmcmFjeyhmLXApIFxzcXJ0e259fXtcc3FydHtwKDEtcCl9fSQgduG7m2kgcGjDom4gcGjhu5FpIGNodeG6qW4gY2jDrW5oIHThuq9jLiBUdXkgbmhpw6puIG4gdGjhur8gbsOgbyBsw6AgIsSR4bunIGzhu5tuIiDEkeG7gyBjw7MgeOG6pXAgeOG7iSBuw7NpIHRyw6puPyBUcuG7nyBs4bqhaSBwaMOibiBwaOG7kWkgTmjhu4sgdGjhu6ljICRcbWF0aHJte0J9KFxtYXRocm17bn0sIFxtYXRocm17cH0pJCwgdGEgYmnhur90IHLhurFuZyBraGkgdHJ1bmcgYsOsbmggbi5wIGPhu6dhIHBow6JuIHBo4buRaSBOaOG7iyB0aOG7qWMgY8OgbmcgbOG7m24gdGjDrCBwaMOibiBwaOG7kWkgbsOgeSBjw6BuZyBn4bqnbiB24bubaSBwaMOibiBwaOG7kWkgY2h14bqpbi4gTeG6t3Qga2jDoWMgdGEgYmnhur90IHLhurFuZyAkbiBcY2RvdCBmKEEpPSQgJG0oQSkkIHbDoCAkbi1uIC4gZihBKSQgbMOgIHPhu5EgbOG6p24geHXhuqV0IGhp4buHbiB2w6Aga2jDtG5nIHh14bqldCBoaeG7h24gQSB0cm9uZyBuIGzhuqduIHRo4butIGPDsyBwaMOibiBwaOG7kWkgJFxtYXRocm17Qn0oXG1hdGhybXtufSwgXG1hdGhybXtwfSkkIHbDoCAkXG1hdGhybXtCfShcbWF0aHJte259LCAxLVxtYXRocm17cH0pJC4gVuG6rXkgxJFp4buBdSBraeG7h24gxJHhu4MgeOG6pXAgeOG7iSBwaMOibiBwaOG7kWkgY+G7p2EgJGYkIHbDoCAkMS0kICRmJCB24bubaSBwaMOibiBwaOG7kWkgY2h14bqpbiBsw6AgJFxtYXRocm17bn0kIMSR4bunIGzhu5tuIHNhbyBjaG8gbi5wIHbDoCAkXG1hdGhybXtufSgxLVxtYXRocm17cH0pJCBjxaluZyBraMOhIGzhu5tuLiBUcm9uZyB0aOG7sWMgaMOgbmggdsOsIHAgY2jGsGEgYmnhur90LCBuZ8aw4budaSB0YSB0aGF5IG7DsyBi4bufaSDGsOG7m2MgbMaw4bujbmcgTUwgbMOgICRmJCwgdsOsIHbhuq15IMSRaeG7gXUga2nhu4duIG4ucCB2w6AgJFxtYXRocm17bn0gLigxLVxtYXRocm17cH0pJCBraMOhIGzhu5tuIMSRxrDhu6NjIHRoYXkgYuG6sW5nICRuJC4gJGYkIHbDoCAkbiQuICQoMS1mKSQga2jDoSBs4bubbi4gVsOsIHbhuq15IMSR4buRaSB24bubaSBiw6BpIHRvw6FuIMaw4bubYyBsxrDhu6NuZyBraG/huqNuZyB0aW4gY+G6rXkgY2hvIHThu7cgbOG7hywgbmfGsOG7nWkgdGEgdGjGsOG7nW5nIHnDqnUgY+G6p3UgY+G7oSBt4bqrdSBraMOhIGzhu5tuIHbDoCBjaGlhIGzDoG0gaGFpIHRyxrDhu51uZyBo4bujcDoNCg0KKipLaGkqKiAkZiQgKipraMO0bmcgcXXDoSBn4bqnbiAwIGhv4bq3YyAxOiAqKg0KDQpLaG/huqNuZyB0aW4gY+G6rXksIHbhu5tpIMSR4buZIHRpbiBj4bqteSAkXGdhbW1hPTEtXGFscGhhJCBjaG8gdOG7tyBs4buHICRwJCBsw6A6DQoNCiQkDQpcbGVmdChmLXVcbGVmdChcZnJhY3tcYWxwaGF9ezJ9XHJpZ2h0KSBcY2RvdCBcc3FydHtcZnJhY3tmKDEtZil9e259fSwgZit1XGxlZnQoXGZyYWN7XGFscGhhfXsyfVxyaWdodCkgXGNkb3QgXHNxcnR7XGZyYWN7ZigxLWYpfXtufX1ccmlnaHQpDQokJA0KDQpDaMO6IMO9OiBUaMO0bmcgdGjGsOG7nW5nIMSR4buZIHRpbiBj4bqteSBsw6AgJDk1IFwlJCwga2hpIMSRw7Mga2hv4bqjbmcgdGluIGPhuq15IGNobyAkcCQgbMOgOiAkJA0KXGxlZnQoZi0xLDk2IC4gXHNxcnR7XGZyYWN7ZigxLWYpfXtufX0sIGYrMSw5NiAuIFxzcXJ0e1xmcmFje2YoMS1mKX17bn19XHJpZ2h0KQ0KJCQgdsOgIGtoaSBraMO0bmcgbsOzaSBnw6wgdOG7m2kgxJHhu5kgdGluIGPhuq15IHRow6wgbeG6t2MgxJHhu4tuaCDEkeG7mSB0aW4gY+G6rXkgbMOgICQ5NSBcJSQuDQoNCipWw60gZHUgNDoqIFRyb25nIG3huqt1IMSRaeG7gXUgdHJhIOG7nyB2w60gZOG7pSAyLCBow6N5IMaw4bubYyBsxrDhu6NuZyBraG/huqNuZyB0aW4gY+G6rXkgJDk1IFwlJCBjaG8gdOG7tyBs4buHIG5nxrDhu51pIE3hu7kg4bunbmcgaOG7mSB2aeG7h2MgcGjDoSB0aGFpIGjhu6NwIHBow6FwIGNobyBwaOG7pSBu4buvIGPDsyB0aGFpIG3DoCBraMO0bmcgbXXhu5FuIGPDsyBjb24uDQoNCkdp4bqjaTogVOG7qyBt4bqrdSwgdGEgY8OzIHThu7cgbOG7hyDhu6duZyBo4buZIHZp4buHYyBwaMOhIHRoYWkgaOG7o3AgcGjDoXAgbMOgICRmPTAsNDQ2JC4NCg0KS2hv4bqjbmcgdGluIGPhuq15IGPhuqduIHTDrG0gY8OzIGThuqFuZzogJFxsZWZ0KGYtMSw5NiAuIFxzcXJ0e1xmcmFje2YoMS1mKX17bn19LCBmKzEsOTYgLiBcc3FydHtcZnJhY3tmKDEtZil9e259fVxyaWdodCkkIHRoYXkgc+G7kSBsaeG7h3UgdsOgbyB0YSBjw7M6ICRmLTEsOTYgLiBcc3FydHtcZnJhY3tmKDEtZil9e259fT0wLDQxNCA7IGYrMSw5NiAuIFxzcXJ0e1xmcmFje2YoMS1mKX17bn19PTAsNDc4JA0KDQpW4bqteSBjw7MgdGjhu4MgY2hvIHLhurFuZyB24bubaSBt4bupYyDEkeG7mSB0aW4gY+G6rXkgJDk1IFwlJCB0aMOsIHThu7cgbOG7hyBkw6JuIHPhu5EgTeG7uSDhu6duZyBo4buZIHZp4buHYyBwaMOhIHRoYWkgaOG7o3AgcGjDoXAgY2hvIHBo4bulIG7hu68gY8OzIHRoYWkgbcOgIGtow7RuZyBtdeG7kW4gY8OzIGNvbiBuaGnhu4N1IGjGoW4gbMOgIHbDoG8ga2hv4bqjbmcgdOG7qyAkNDEsNCBcJSQgxJHhur9uICQ0Nyw4IFwlJC4NCg0KKipLaGkqKiAkZiBcbGVxIDAsMSQqKiwgaG/hurdjKiogJGYgXGdlcSAwLDkkICoqOioqDQoNClRyb25nIHRyxrDhu51uZyBo4bujcCBuw6B5LCB2aeG7h2MgeOG6pXAgeOG7iSBwaMOibiBwaOG7kWkgbmjhu4sgdGjhu6ljIGLhu59pIHBow6JuIHBo4buRaSBjaHXhuqluIG7Ds2kgY2h1bmcgbMOgIHRow7QsIHbDrCBt4bq3YyBkw7kgY+G7oSBt4bqrdSBuIGtow6EgbOG7m24gbmjGsG5nIGPDsyB0aOG7gyAkbiAuIGYkIHbDoCAkbiAuKDEtZikkIGzhuqFpIGtow7RuZyDEkeG7pyBs4bubbix24bqteSBuw6puIHjhuqVwIHjhu4kgYuG7n2kgcGjDom4gcGjhu5FpIFBvaXNzb24gbMOgIGjhu6NwIGzDvSBoxqFuLiBWw6wgdGjhur8sIGLDoGkgdG/DoW4ga2hv4bqjbmcgdGluIGPhuq15IGNobyB04bu3IGzhu4cgJFxwaSQgY2hpYSBsw6BtIDIgYsaw4bubYzoNCg0KQsawxqHMgWMgMS4gVMOsbSBraG/huqNuZyB0aW4gY+G6rXkgY2hvICRuIC4gcCQgbMOgIGtob+G6o25nICRcbGVmdChuIHBfMSwgbiBwXzJccmlnaHQpJCwgdHJvbmcgxJHDsyAkbiBwXzEkLCAkbiBwXzIkIHTDrG0gxJHGsOG7o2MgcXVhIHRyYSBi4bqjbmcgcGjhu6UgbOG7pWMgSVYsIGPEg24gY+G7qSB2w6BvIGdpw6EgdHLhu4sgdOG6p24gc+G7kSBt4bqrdTogJG09JCBuLmYuIChUcm9uZyBi4bqjbmcgbsOgeTogY+G7mXQgxJHhuqd1IHRpw6puIGzDoCBnacOhIHRy4buLIGNo4bq1biBjaOG7pWMgY+G7p2EgJG0kIHbDoCBkw7JuZyB0csOqbiBjw7luZyBi4buVIHN1bmcgY2jhu68gc+G7kSBow6BuZyDEkcahbiB24buLIGPhu6dhICRtJCwgZ2lhbyBj4bunYSBow6BuZyAoY2jhurVuIGNo4bulYyBj4bunYSAkbSQgKSB2w6AgY+G7mXQgKGjDoG5nIMSRxqFuIHbhu4sgY+G7p2EgJG0kICkgbMOgIMO0IGNo4bupYSAyIGPhuq1uOiBz4buRIHRyw6puIGzDoCAkbiBwXzEkLCBz4buRIGTGsOG7m2kgbMOgICRcbGVmdC5uIHBfMlxyaWdodCkkLg0KDQpCxrDGocyBYyAyLiBDaGlhIDIgY+G6rW4gJG4gcF8xLCBuIHBfMiQgY2hvICRuJCAoY+G7oSBt4bqrdSkgbmjhuq1uIMSRxrDhu6NjIGtob+G6o25nICRcbGVmdChwXzEsIHBfMlxyaWdodCkkIGzDoCBraG/huqNuZyB0aW4gY+G6rXkgJDk1IFwlJCBjaG8gdOG7tyBs4buHICRwJC4NCg0KVsOtIGR1IDU6IMSQaeG7gXUgdHJhIG3hu6ljIHRodSBuaOG6rXAgY+G7p2EgMjUwIGjhu5kgZMOibiDhu58gbeG7mXQga2h1IHbhu7FjLCB0aOG6pXkgY8OzIDcgaOG7mSB0aHUgbmjhuq1wIOG7nyBt4bupYyBuZ2jDqG8ga2jhu5UuIEjDo3kgxrDhu5tjIGzGsOG7o25nIGtob+G6o25nIHRpbiBj4bqteSBjaG8gdOG7tyBs4buHIGjhu5kgdGh1IG5o4bqtcCDhu58gbeG7qWMgbmdow6hvIGto4buVIOG7nyBraHUgduG7sWMgbsOgeS4NCg0KR2nhuqNpOiBH4buNaSBwIGzDoCB04bu3IGzhu4cgaOG7mSBjw7MgdGh1IG5o4bqtcCDhu58gbeG7qWMgbmdow6hvIGto4buVIGPhu6dhIGtodSB24buxYy4gVOG7qyBt4bqrdSDEkWnhu4F1IHRyYSBjw7MgdOG7tyBs4buHIGjhu5kgY8OzIHRodSBuaOG6rXAg4bufIG3hu6ljIG5naMOobyBraOG7lSBsw6A6ICRmPVxmcmFjezd9ezI1MH09MCwwMTQ8MCwxJC4gVuG6rXkgdGEgdMOsbSBraG/huqNuZyB0aW4gY+G6rXkgY2hvIHAgdGhlbyBjw6FjIGLGsOG7m2Mgc2F1Og0KDQpcLSDGr+G7m2MgbMaw4bujbmcga2hv4bqjbmcgdGluIGPhuq15IGNobyBuLnA6IFbhu5tpICRtPW4gLiBmPTckLCB0cmEgYuG6o25nIHBo4bulIGzhu6VjIElWLCBjw7Mga2hv4bqjbmcgdGluIGPhuq15IGNobyBuLnAgbMOgOiAkKDIsMzEgOyAxNCw0MikkLg0KDQotU3V5IHJhIGtob+G6o25nIHRpbiBj4bqteSBjaG8gcCBsw6A6ICQoMCwwMDkyNCA7IDAsMDU3NjgpJC4NCg0KVuG6rXkgduG7m2kgxJHhu5kgdGluIGPhuq15ICQ5NSBcJSQsIGPDsyB0aOG7gyBjaG8gcuG6sW5nIHThu7cgbOG7hyBo4buZIGPDsyBt4bupYyB0aHUgbmjhuq1wIG5naMOobyBraOG7lSDhu58ga2h1IHbhu7FjIG7DoHkgdsOgbyBraG/huqNuZyB04burICQwLDA5MjQgXCUkIMSR4bq/biAkNSw3NjggXCUkLg0KDQpOaOG6rW4geMOpdDogVuG7m2kgdHLGsOG7nW5nIGjhu6NwIG3huqt1IG7DoHksIG7hur91IHRhIGPhu5EgeOG6pXAgeOG7iSBwaMOibiBwaOG7kWkgY+G7p2EgJGYkIHbhu5tpIHBow6JuIHBo4buRaSBjaHXhuqluLCB04bupYyBsw6AgZMO5bmcgY8O0bmcgdGjhu6ljIDEuNS4xYSwgdGjDrCDGsOG7m2MgbMaw4bujbmcgc+G6vSBxdcOhIHRow7QsIHRo4bqtbSBjaMOtIGPDsyBj4bqtbiBkxrDhu5tpOiAkJA0KZi0xLDk2IC4gXHNxcnR7XGZyYWN7ZigxLWYpfXtufX09MCwwMTQtMSw5NiAuIFxzcXJ0e1xmcmFjezAsMDE0KDEtMCwwMTQpfXsyNTB9fT0tMCwwMzIxPDAgLg0KJCQNCg0KQ2jDuiDDvTogVHJvbmcgdGjhu7FjIGjDoG5oLCBraeG7g20gxJHhu4tuaCAkSF8wOiBwPTAsNSQgY2jhu4kgecOqdSBj4bqndSAkbj4xMCQgdHJvbmcga2hpIGtp4buDbSDEkeG7i25oICRIXzA6IHA9MC4xJCBob+G6t2MgJEhfMDogcD0wLjkkIHnDqnUgY+G6p3UgJG4+NTAkLiBZw6p1IGPhuqd1IHbhu4Ega8OtY2ggY+G7oSBt4bqrdSBwaOG6o24gw6FuaCBz4buxIHRoYXkgxJHhu5VpIG5nw6B5IGPDoG5nIGzhu4djaCBj4bunYSB2aeG7h2MgcGjDom4gcGjhu5FpIG3huqt1IGPhu6dhICRmJCBraGkgJHAkIMSR4bqhdCDEkeG6v24gMCBob+G6t2MgMS4NCg0KVMOzbSBs4bqhaSwgY2jGsMahbmcgbsOgeSDEkcOjIGdp4bubaSB0aGnhu4d1IGhhaSBwaMOibiBwaOG7kWkgY2jDrW5oIGNobyBwaMOibiB0w61jaCBk4buvIGxp4buHdSDEkeG7i25oIHTDrW5oOiBuaOG7iyBwaMOibiB2w6AgUG9pc3Nvbi4gxJDhu5NuZyB0aOG7nWkgY2jDum5nIHTDtGkgY8WpbmcgxJHDoyBnaeG7m2kgdGhp4buHdSBwaMawxqFuZyBwaMOhcCDGsOG7m2MgbMaw4bujbmcgaOG7o3AgbMO9IGPhu7FjIMSR4bqhaSB2w6AgbWluaCBob+G6oSB2aeG7h2Mgc+G7rSBk4bulbmcgbsOzIGNobyBk4buvIGxp4buHdSB04bu3IGzhu4cuIFBo4bqnbiBjw7JuIGzhuqFpIGPhu6dhIG7hu5lpIGR1bmcgc+G7rSBk4bulbmcgc3V5IGx14bqtbiBNTCBjaG8gY8OhYyB0aGFtIHPhu5Egbmjhu4sgcGjDom4gdsOgIFBvaXNzb24gdHJvbmcgbmhp4buBdSB0csaw4budbmcgaOG7o3Aga2jDoWMgbmhhdS4NCg0KIyBCw6BpIHThuq1wIGNoxrDGoW5nIDENCg0KKioxLjEuKiogVHJvbmcgY8OhYyB2w60gZOG7pSBzYXUsIHjDoWMgxJHhu4tuaCBiaeG6v24gxJHDoXAg4bupbmcgdsOgIGPDoWMgYmnhur9uIGdp4bqjaSB0aMOtY2guDQoNCmFcLiBUaMOhaSDEkeG7mSDEkeG7kWkgduG7m2kga2nhu4NtIHNvw6F0IHPDum5nICjhu6duZyBo4buZLCBwaOG6o24gxJHhu5FpKSwgZ2nhu5tpIHTDrW5oIChu4buvLCBuYW0pLCBHacOhbyBk4bulYyBj4bunYSBt4bq5ICh0cnVuZyBo4buNYywgY2FvIMSR4bqzbmcpLg0KDQpiXC4gQuG7h25oIHRpbSAoY8OzLCBraMO0bmcpLCBIdXnhur90IMOhcCwgTeG7qWMgY2hvbGVzdGVyb2wgKENob2xlc3Rlcm9sKS4NCg0KY1wuIENo4bunbmcgdOG7mWMgKHRy4bqvbmcsIGtow7RuZyB0cuG6r25nKSwgVMO0biBnacOhbyAoQ8O0bmcgZ2nDoW8sIERvIFRow6FpLCBUaW4gbMOgbmgpLCBCw6xuaCBjaOG7jW4gVOG7lW5nIHRo4buRbmcgKETDom4gY2jhu6csIEPhu5luZyBow7JhLCBLaMOhYyksIFRodSBuaOG6rXAgaMOgbmcgbsSDbS4NCg0KKioxLjIqKi4gVGhhbmcgxJFvIG7DoG8gcGjDuSBo4bujcCBuaOG6pXQgY2hvIGPDoWMgYmnhur9uIHNhdSDEkcOieS0gZGFuaCBuZ2jEqWEsIGhv4bq3YyB0aOG7qSB04buxPw0KDQphXC4gTGnDqm4ga+G6v3QgxJHhuqNuZyBjaMOtbmggdHLhu4sgKETDom4gY2jhu6csIEPhu5luZyBow7JhLCBraMO0bmcgbGnDqm4ga+G6v3QpLg0KDQpiXC4gQ8OzIGLhurFuZyBjYW8gbmjhuqV0IChraMO0bmcgY8OzLCB0cnVuZyBo4buNYywgY+G7rSBuaMOibiwgdGjhuqFjIHPEqSwgdGnhur9uIHPEqSkuDQoNCmNcLiBUw6xuaCB0cuG6oW5nIGLhu4duaCBuaMOibiAodOG7kXQsIGPDtG5nIGLhurFuZywgbmdoacOqbSB0cuG7jW5nLCBwaMOqIGLDrG5oKS4NCg0KZFwuIFbhu4sgdHLDrSBj4bunYSBi4buHbmggdmnhu4duIChMb25kb24sIEJvc3RvbiwgTWFkaXNvbiwgUm9jaGVzdGVyLCBUb3JvbnRvKS4NCg0KZVwuIE7GsOG7m2MgZ2nhuqNpIGtow6F0IHnDqnUgdGjDrWNoIChiaWEsIG7GsOG7m2MgdHLDoWkgY8OieSwgc+G7r2EsIG7GsOG7m2Mgbmfhu410LCByxrDhu6N1IHZhbmcsIGtow6FjKS4NCg0KKioxLjMuKiogQ8OhYyB04bqlbSBzaWxpY29uIGNobyBjaGlwIG3DoXkgdMOtbmggxJHGsOG7o2Mgc+G6o24geHXhuqV0IGLhu59pIG3hu5l0IGPDtG5nIHR5IGPDtG5nIG5naOG7hyBjYW8gY8OzIG3hu5l0IHPhu5Ega2h1eeG6v3QgdOG6rXQgdHJ1bmcgYsOsbmggbMOgIDEuMCBt4buXaSB04bqlbS4gTuG6v3Ugc+G7kSBsxrDhu6NuZyBraHV54bq/dCB04bqtdCBjw7MgcGjDom4gcGjhu5FpIFBvaXNzb24sIHTDrG0geMOhYyBzdeG6pXQgbcOgIG3hu5l0IHThuqVtIGPDsyAoYSkgxJHDum5nIDAga2hp4bq/bSB04bqtdCwgKGIpIMSRw7puZyAxIGtoaeG6v20gdOG6rXQsIChjKSDDrXQgbmjhuqV0IDIga2hp4bq/bSBraHV54bq/dC4NCg0KKioxLjQuKiogTeG7l2kgdHJvbmcgMTAwIGPDonUgaOG7j2kgdHLhuq9jIG5naGnhu4dtIHRyb25nIG3hu5l0IGvhu7MgdGhpIGPDsyBi4buRbiDEkcOhcCDDoW4gY8OzIHRo4buDIG5oxrBuZyBjaOG7iSBjw7MgbeG7mXQgxJHDoXAgw6FuIMSRw7puZy4gxJDhu5FpIHbhu5tpIG3hu5dpIGPDonUgaOG7j2ksIGjhu41jIHNpbmggY2jhu41uIG5n4bqrdSBuaGnDqm4gbeG7mXQgxJHDoXAgw6FuIGzDoCBjw6J1IHRy4bqjIGzhu51pLiBYw6FjIMSR4buLbmggcGjDom4gcGjhu5FpIHPhu5EgY8OidSB0cuG6oyBs4budaSDEkcO6bmcgY+G7p2EgaOG7jWMgc2luaCB0cm9uZyBr4buzIHRoaS4gROG7sWEgdHLDqm4gdHJ1bmcgYsOsbmggdsOgIMSR4buZIGzhu4djaCBjaHXhuqluIGPhu6dhIGLhuqNuIHBow6JuIHBo4buRaSDEkcOzLCBsaeG7h3UgY8OzIMSRw6FuZyBuZ+G6oWMgbmhpw6puIG7hur91IGjhu41jIHNpbmggY8OzIMOtdCBuaOG6pXQgNTAgY8OidSB0cuG6oyBs4budaSDEkcO6bmc/DQoNCioqMS41LioqIE3hu5l0IMSR4buTbmcgeHUgY8OibiDEkeG7kWkgxJHGsOG7o2MgdHVuZyBoYWkgbOG6p24uIEfhu41pIGzDoCBz4buRIGzGsOG7o25nIG3hurd0IHPhuqVwIHRodSDEkcaw4bujYy4NCg0KYVwuIFjDoWMgxJHhu4tuaCB4w6FjIHN14bqldCBjaG8gY8OhYyBnacOhIHRy4buLIGPDsyB0aOG7gyBjaG8gLCBnacOhIHRy4buLIHRydW5nIGLDrG5oIHbDoCBwaMawxqFuZyBzYWkgY+G7p2EgWS4NCg0KYlwuIFTDrW5oIGPDoWMgeMOhYyBzdeG6pXQgY2hvIHBow6JuIHBo4buRaSBQb2lzc29uIGPDsyBjw7luZyBnacOhIHRy4buLIHRydW5nIGLDrG5oIGPhu6dhIFkuIFPhu7Ega2jDoWMgYmnhu4d0IGPhu6dhIG7DsyBzbyB24bubaSBwaMawxqFuZyBzYWkgdHJvbmcgKGEpIG5oxrAgdGjhur8gbsOgbz8NCg0KY1wuIMSQ4buRaSB24bubaSBt4buXaSBs4bqnbiB0dW5nIG3hu5l0IMSR4buTbmcgeHUga2jDtG5nIMSR4buTbmcgY2jhuqV0LCBn4buNaSDwnZyLIGzDoCBiaeG7g3UgdGjhu4sgeMOhYyBzdeG6pXQgY+G7p2EgbeG6t3Qgc+G6pXAuIEdp4bqjIHPhu60gY8OzIDAgbeG6t3Qgc+G6pXAgdHJvbmcgMiBs4bqnbiB0dW5nLiBUw6xtIMaw4bubYyBsxrDhu6NuZyBNTCBj4bunYSDwnZyLLg0KDQoqKjEuNi4qKiBRdWF5IGzhuqFpIGPDoWMgYsOgaSB04bqtcCB0csaw4bubYyDEkcOzLg0KDQphXC4gVMOtbmggY8OhYyB4w6FjIHN14bqldCBuaOG7iyB0aOG7qWMgY2hvICRuPTIkIGtoaSB4w6FjIHN14bqldCBj4bunYSBt4bq3dCBz4bqlcCBjaG8gbeG7l2kgbOG6p24gdHVuZyBi4bqxbmcgKGkpICRccGk9MC42JCwgKGlpKSAkXHBpPTAuNCQuDQoNCmJcLiBHaeG6oyBz4butIGNow7puZyB0YSBxdWFuIHPDoXQgJFk9MSQuIFTDrW5oIHRvw6FuIHbDoCBwaMOhYyBob+G6oSBow6BtIGjhu6NwIGzDvS4NCg0KY1wuIFPhu60gZOG7pW5nIGjDoG0gaOG7o3AgbMO9IHbhur0gdOG7qyAoYiksIGNo4buJIHJhIMaw4bubYyBsxrDhu6NuZyBNTCBj4bunYSAkXHBpPTAuNSQuDQoNCioqMS43LioqIFRyb25nIGN14buRbiB04buxIHRydXnhu4duIEEgU29ydCBvZiBMaWZlLCB0w6FjIGdp4bqjIG5nxrDhu51pIEFuaCBHcmFoYW0gR3JlZW5lIG3DtCB04bqjIGdpYWkgxJFv4bqhbiB0cuG6p20gY+G6o20gbmdoacOqbSB0cuG7jW5nIHRyb25nIHRo4budaSBnaWFuIMSRw7Mgw7RuZyBjaMahaSBSdXNzaWFuIFJvdWxldHRlLiAiVHLDsiBjaMahaSIgbsOgeSBiYW8gZ+G7k20gxJHhurd0IG3hu5l0IHZpw6puIMSR4bqhbiB2w6BvIG3hu5l0IHRyb25nIHPDoXUga2hvYW5nIHPDum5nLCBxdWF5IGPDoWMgYnXhu5NuZyDEkeG7gyBjaOG7jW4gbmfhuqt1IG5oacOqbiwgdsOgIHNhdSDEkcOzIG5o4bqxbSBt4buZdCDEkeG6p3UgbmfGsOG7nWkgdsOgIGLDs3AgY8OyLg0KDQphXC4gR3JlZW5lIGNoxqFpIGdhbWUgbsOgeSBzw6F1IGzhuqduLCB2w6AgbWF5IG3huq9uIGzDoCBraMO0bmcgYWkgdHJvbmcgc+G7kSBo4buNIMSRw6MgYuG7iyBi4bqvbiBt4buZdCB2acOqbiDEkeG6oW4uIFTDrG0geMOhYyBzdeG6pXQgY+G7p2Ega+G6v3QgcXXhuqMgbsOgeS4NCg0KYlwuIEdp4bqjIHPhu60gbeG7mXQgbmfGsOG7nWkgdGnhur9wIHThu6VjIGNoxqFpIHRyw7IgY2jGoWkgbsOgeSBjaG8gxJHhur9uIGtoaSDEkeG6oW4gYuG6r24gdHLDum5nIMSRw61jaC4gR+G7jWkgWSBsw6Agc+G7kSB2acOqbiDEkeG6oW4gxJHDoyBi4bqvbi4gQ2jhu6luZyB04buPIHLhurFuZyB4w6FjIHN14bqldCBj4bunYSBr4bq/dCBxdeG6oyB5IGLhurFuZyB24bubaSAkKDUvNilee3ktMX0oMS82KSQgLCBjaG8geT0xLDIsMyzigKYgKMSRxrDhu6NjIGfhu41pIGzDoCBwaMOibiBwaOG7kWkgaMOsbmggaOG7jWMpLg0KDQoqKjEuOC4qKiBN4buZdCBt4bqrdSBn4buTbSBuaOG7r25nIHBo4bulIG7hu68gYuG7iyBjaOG7qW5nIGtow7MgdGjhu58gZG8gdGh14buRYyBuaHXhu5ltIMSRw6MgZMO5bmcgdGh14buRYyBnaeG6o20gxJFhdSDEkcaw4bujYyBjaOG7iSDEkeG7i25oIMSR4buDIGzDoG0gZ2nhuqNtIHTDoWMgaOG6oWkuIE3hu5l0IHRodeG7kWMgZ2nhuqNtIMSRYXUgbeG7m2kgxJHGsOG7o2MgdHV5w6puIGLhu5EgaGnhu4d1IHF14bqjIGjGoW4uIFNhdSBraGkgdGjhu60gdGh14buRYyBnaeG6o20gxJFhdSBt4bubaSwgNDAgcGjhu6UgbuG7ryBjaG8gdGjhuqV5IGdp4bqjbSDEkWF1IGjGoW4gduG7m2kgdGh14buRYyDEkcOjIGNo4buJIMSR4buLbmgsIHbDoCA2MCBuZ8aw4budaSBjaG8gYmnhur90IGPDsyBnaeG6o20gxJFhdSB04buRdCBoxqFuIHbhu5tpIGxv4bqhaSB0aHXhu5FjIG3hu5tpLmEuIEjDo3kga2nhu4NtIMSR4buLbmggZ2nhuqMgdGh1eeG6v3QgY2hvIHLhurFuZyBraOG6oyBuxINuZyBsw6BtIGdp4bqjbSDEkWF1IGPhu6dhIGxv4bqhaSB0aHXhu5FjxJHDoyBjaOG7iSDEkeG7i25oIHbDoCBj4bunYSBsb+G6oWkgdGh14buRYyBt4bubaSBsw6AgbmjGsCBuaGF1LmIuIFjDonkgZOG7sW5nIHbDoCBnaeG6o2kgdGjDrWNoIGtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIHjDoWMgc3XhuqV0IGdp4bqjbSDEkWF1IHThu5F0IGjGoW4gY+G7p2EgdGh14buRYyBnaeG6o20gxJFhdSBt4bubaS4NCg0KKioxLjkuKiogVHLhu58gbOG6oWkgY8OhYyBiw6BpIHRvw6FuIHRyxrDhu5tjLiBDw6FjIG5ow6AgbmdoacOqbiBj4bupdSBtdeG7kW4gY8OzIG3hu5l0IG3huqt1IMSR4bunIGzhu5tuIMSR4buDIGPDsyB0aOG7gyDGsOG7m2MgbMaw4bujbmcgeMOhYyBzdeG6pXQgbOG7sWEgY2jhu41uIHRodeG7kWMgZ2nhuqNtIMSRYXUgbeG7m2kgduG7m2kgxJHhu5kgY2jDrW5oIHjDoWMgMC4wOCwgduG7m2kgxJHhu5kgdGluIGPhuq15IDAuOTUuDQoNCioqMS4xMC4qKiBU4bqhcCBjaMOtIE5ld3N3ZWVrICgyNy8zLzE5ODkpIGLDoW8gY8OhbyBr4bq/dCBxdeG6oyBjdeG7mWMgdGjEg20gZMOyIMO9IGtp4bq/biB0w7RuIGdpw6FvLCBkbyBU4buVIGNo4bupYyBHYWxsdXAgdGnhur9uIGjDoG5oLiBUcm9uZyBz4buRIDc1MCBuZ8aw4budaSBN4bu5IHRyxrDhu59uZyB0aMOgbmgsIDI0JSB0aW4gdsOgbyBsdcOibiBo4buTaS4gWGVtIMSRw6J5IG5oxrAgbMOgIG3hu5l0IG3huqt1IG5n4bqrdSBuaGnDqm4sIHjDonkgZOG7sW5nIHbDoCBnaeG6o2kgdGjDrWNoIG3hu5l0IGtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIHThu7cgbOG7hyB0aOG7sWMgc+G7sSBj4bunYSBuZ8aw4budaSBs4bubbiBN4bu5IHRpbiB0xrDhu59uZyB2w6BvIGx1w6JuIGjhu5NpLg0KDQoqKjEuMTEuKiogTeG7mXQgbmjDoCB04buZaSBwaOG6oW0gaOG7jWMgbXXhu5FuIMaw4bubYyB0w61uaCB04bu3IGzhu4cgY8O0bmcgZMOibiBN4bu5IHPhu5FuZyB0cm9uZyBt4buZdCBjxINuIG5ow6AgY8OzIHPDum5nIOG7kW5nLiBDdeG7mWMgS2jhuqNvIHPDoXQgWMOjIGjhu5lpIG7Eg20gMTk5MSDEkcOjIGjhu49pIG5o4buvbmcgbmfGsOG7nWkgdHLhuqMgbOG7nWksICJUcm9uZyBuaMOgIGLhuqFuIOG7nyBjw7MgYuG6pXQga+G7syBzw7puZyBob+G6t2Mgc8O6bmcgbOG7pWMg4buVIHF1YXkgbsOgbz8iIFRyb25nIHPhu5Egbmjhu69uZyBuZ8aw4budaSDEkcaw4bujYyBo4buPaSwgMzkzIHRy4bqjIGzhu51pICJjw7MiIHbDoCA1ODMgdHLhuqMgbOG7nWkgImtow7RuZyIuIFjDonkgZOG7sW5nIGtob+G6o25nIHRpbiBj4bqteSA5MCUgY2hvIHThu7cgbOG7hyB0aOG7sWMgY+G7p2EgImPDsyIuIEdp4bqjaSB0aMOtY2guDQoNCioqMS4xMi4qKiBHaeG6oyBz4butICRZJCBsw6AgbeG7mXQgYmnhur9uIG5n4bqrdSBuaGnDqm4gY8OzIHBow6JuIHBo4buRaSBuaOG7iyB0aOG7qWMgQihuLCBwKSwgdsOgICRmPVkvbiQgLiBDaOG7iSByYSBy4bqxbmcgJM+DKGYpPVxzcXJ0e3AoMS1wKS9ufSQgLg0KDQoqKjEuMTMuKiogTeG7mXQgYmnhur9uIG5n4bqrdSBuaGnDqm4gY8OzIHBow6JuIHBo4buRaSBQb2lzc29uLCB24bubaSB0aGFtIHPhu5EgY2jGsGEgYmnhur90ICTOvCQgLiBRdWFuIHPDoXQgZHV5IG5o4bqldCBi4bqxbmcgMC4NCg0KYVwuIFTDrG0gdsOgIHbhur0gY8OhYyBow6BtIGjhu6NwIGzDvSB0csOqbiBraMO0bmcgZ2lhbiBjw6FjIGdpw6EgdHLhu4sgY8OzIHRo4buDIGPDsyBj4bunYSBjw6FjIGdpw6EgdHLhu4sgY2hvICTOvCQgLg0KDQpiXC4gxq/hu5tjIGzGsOG7o25nIE1MIGPhu6dhICTOvCQgbMOgIGfDrD8NCg0KKioxLjE0LioqIFTDrG0gxrDhu5tjIGzGsOG7o25nIE1MIGNobyB0aGFtIHPhu5Eg8J2chiBj4bunYSBwaMOibiBwaOG7kWkgbcWpLg0KDQoqKjEuMTUuKiogQ2jhu4kgcmEgbeG7mXQgZ2nDoSB0cuG7iyAkcF8wJCB24bubaSB0aOG7kW5nIGvDqiAkej0oZi1wXzApL1xzcXJ0e3BfMCgxLXBfMCkvbn0kIG5o4bqtbiBt4buZdCB2w6BpIGdpw6EgdHLhu4sgY+G7kSDEkeG7i25oIGzDoCBs4budaSBnaeG6o2kgY+G7p2EgcGjGsMahbmcgdHLDrG5oICQoMSt6XzBeMiAvbilccGlfMF4yKygtMmYtel8wXjIvbilwXzArZl4yJCAuIFThu6sgxJHDsywgc+G7rSBk4bulbmcgY8O0bmcgdGjhu6ljICR4PS1iwrFcc3FydHtiXjItNGFjfS8yYSQgxJHhu4MgZ2nhuqNpIHBoxrDGoW5nIHRyw6xuaCBi4bqtYyBoYWkgJGF4XjIrYngrYz0wJCwgbOG6pXkgY8OhYyBnaeG7m2kgaOG6oW4gY2hvIGtob+G6o25nIHRpbiBj4bqteSA5NSUgdHJvbmcgdsOtIGThu6UgMyBjaG8gdOG7tyBs4buHIE5nxrDhu51pIE3hu7kg4bunbmcgaOG7mSBwaMOhIHRoYWkgaOG7o3AgcGjDoXAuDQoNCiMgQ2jGsMahbmcgMjogQuG6ok5HIE5H4bqqVSBOSEnDik4gVsOAIFNVWSBESeG7hE4gVEjhu5BORyBLw4oNCg0KUGjDom4gdMOtY2ggY8OhYyBt4buRaSBxdWFuIGjhu4cgdHJvbmcgbeG7mXQgdsOpYyB0xqEgcXVhbiBzw6F0IGzDoCB0cuG7jW5nIHTDom0gY+G7p2EgaOG6p3UgaOG6v3QgY8OhYyBwaMOibiB0w61jaCB0aOG7kW5nIGvDqiDEkWEgYmnhur9uLiBDaMawxqFuZyBuw6B5IMSR4buBIGPhuq1wIMSR4bq/biBt4buRaSBxdWFuIGjhu4cgZ2nhu69hIGhhaSBiaeG6v24gxJHhu4tuaCB0w61uaCB2w6AgcGjhuqduIGN14buRaSDEkeG7gSBj4bqtcCDEkeG6v24gcXVhbiBo4buHIGPhu6dhIG5oaeG7gXUgYmnhur9uIMSR4buLbmggdMOtbmguIFRyb25nIMSRw7MgZ2nhu5tpIHRoaeG7h3UgY8OhYyB0aGFtIHPhu5EgbcO0IHThuqMgc+G7sSBr4bq/dCBo4bujcCB2w6AgxJHGsGEgcmEgY8OhYyBwaMawxqFuZyBwaMOhcCBzdXkgZGnhu4VuIGNobyBjw6FjIHRoYW0gc+G7kSDEkcOzLg0KDQpOaGnhu4F1IOG7qW5nIGThu6VuZyBsacOqbiBxdWFuIMSR4bq/biB2aeG7h2Mgc28gc8OhbmggZ2nhu69hIGhhaSBuaMOzbS4gUGjhuqduIDIuMSBz4bq9IHRyw6xuaCBiw6B5IGPhuqV1IHRyw7pjIHjDoWMgc3XhuqV0IGNobyBi4bqjbmcgbmfhuqt1IG5oacOqbi4gUGjhuqduIDIuMiBz4bq9IHRyw6xuaCBiw6B5IGPDoWMgcGjGsMahbmcgcGjDoXAgcGjDom4gdMOtY2ggc+G7sSBraMOhYyBiaeG7h3QgduG7gSB04bu3IGzhu4cuIFBo4bqnbiAyLjMgdHLDrG5oIGLDoHkgduG7gSDEkeG7mSDEkW8gb2RkIHJhdGlvLCDEkcOzbmcgdmFpIHRyw7IgcuG6pXQgcXVhbiB0cuG7jW5nIMSR4buRaSB24bubaSBt4buZdCBz4buRIHBoxrDGoW5nIHBow6FwIMSRxrDhu6NjIMSR4buBIGPhuq1wIHRyb25nIG7hu5lpIGR1bmcgaOG7jWMgcGjhuqduIG7DoHkuIFBo4bqnbiAyLjQgdsOgIDIuNSBnaeG7m2kgdGhp4buHdSBt4buZdCBz4buRIHBoxrDGoW5nIHBow6FwIHBow6JuIHTDrWNoIMSR4buDIHRy4bqjIGzhu51pIGPDonUgaOG7j2kgY8OzIHThu5NuIHThuqFpIGhheSBraMO0bmcgc+G7sSBsacOqbiBr4bq/dCBnaeG7r2EgaGFpIGJp4bq/biDEkeG7i25oIHTDrW5oLiBQaOG6p24gMi42IGzDoCBt4bufIHLhu5luZyBjaG8gYuG6o25nIG5n4bqrdSBuaGnDqm4gbmhp4buBdSBjaGnhu4F1Lg0KDQojIyAqKkPhuqV1IHRyw7pjIHjDoWMgc3XhuqV0IGNobyBjw6FjIGLhuqNuZyBuZ+G6q3Ugbmhpw6puKioNCg0KROG7ryBsaeG7h3UgxJHhu4tuaCB0w61uaCBsw6AgdOG6p24gc+G7kSBzdeG6pXQgaGnhu4duIGPDoWMgYmnhu4N1IGhp4buHbiBj4bunYSBjw6FjIGJp4bq/bi4gQ2hvIFggdsOgIFkgbMOgIGhhaSBiaeG6v24gxJHhu4tuaCB0w61uaCwgWCBjw7MgJGskIGJp4buDdSBoaeG7h246ICRBXzEsIEFfMiwuLi4sQV9rJCB2w6AgWSBjw7MgJG0kIGJp4buDdSBoaeG7h246ICRCXzEsIEJfMiwgLi4uLEJfbSQuIENow7puZyB0YSBjw7MgdGjhu4Mgc+G7rSBk4bulbmcgbeG7mXQgYuG6o25nIGfhu5NtICRrJCBow6BuZyB2w6AgJG0kIGPhu5l0IMSR4buDIHRo4buDIGhp4buHbiBr4bq/dCBxdeG6oyBjw7MgdGjhu4MgeOG6o3kgcmEgdOG7qyB2aeG7h2Mga2jhuqNvIHPDoXQ6DQoNCnwgKipYIFxcIFkqKiB8ICRCXzEkICAgIHwgJEJfMiQgICAgfCAuLi4gfCAkQl9tJCAgICB8DQp8LS0tLS0tLS0tLS0tfC0tLS0tLS0tLS18LS0tLS0tLS0tLXwtLS0tLXwtLS0tLS0tLS0tfA0KfCAkQV8xJCAgICAgIHwgJG5fezExfSQgfCAkbl97MTJ9JCB8IC4uLiB8ICRuX3sxbX0kIHwNCnwgJEFfMiQgICAgICB8ICRuX3syMX0kIHwgJG5fezIyfSQgfCAuLi4gfCAkbl97Mm19JCB8DQp8IC4uLiAgICAgICAgfCAuLi4gICAgICB8IC4uLiAgICAgIHwgLi4uIHwgLi4uICAgICAgfA0KfCAkQV9rJCAgICAgIHwgJG5fe2sxfSQgfCAkbl97azJ9JCB8IC4uLiB8ICRuX3trbX0kIHwNCg0KOiBCYcyJbmcgMi4xDQoNCkLhuqNuZyBuw6B5IMSRxrDhu6NjIGfhu41pIGzDoCBi4bqjbmcgbmfhuqt1IG5oacOqbiBoYWkgY2hp4buBdSAka8OXbSQgdHJvbmcgxJHDsyAkbl97aWp9JCBsw6Agc+G7kSBs4bqnbiBxdWFuIHPDoXQgxJHGsOG7o2MgY+G6t3AgdGh14buZYyB0w61uaCAkKEFfaSwgQl9qKSQsIGPDsm4gZ+G7jWkgbMOgIHThuqduIHPhu5EgY+G7p2EgJChBX2ksIEJfaikkLiBN4buZdCBi4bqjbmcgbmfhuqt1IG5oacOqbiBoYWkgYmnhur9uIMSRxrDhu6NjIGfhu41pIGzDoCBi4bqjbmcgaGFpIGNoaeG7gXU7IG3hu5l0IGLhuqNuZyBuZ+G6q3Ugbmhpw6puIDMgYmnhur9uIGfhu41pIGzDoCBi4bqjbmcgbmfhuqt1IG5oacOqbiAzIGNoaeG7gXUuIFRoxrDhu51uZyB0aMOsIGtoaSB0csOsbmggYsOgeSBi4bqjbmcgZOG7ryBsaeG7h3UgMiBjaGnhu4F1LCBuZ8aw4budaSB0YSBraMO0bmcga+G6uyDEkcaw4budbmcgcGjDom4gY2hpYSBj4buZdCBtw6AgY2jhu4kga+G6uyBt4buZdCBz4buRIMSRxrDhu51uZyBuZ2FuZyBj4bqnbiB0aGnhur90Lg0KDQoqVsOtIGThu6UgMToqIEto4bqjbyBzw6F0IDEuMDkxIG5nxrDhu51pIE3hu7kgduG7gSBuaeG7gW0gdGluIGPhu6dhIGjhu40gxJHhu5FpIHbhu5tpIHRo4bq/IGdp4bubaSBiw6puIGtpYSAoc2F1IGtoaSBjaOG6v3QpOiA0MzUgbuG7ryB0aW4gdsOgIDE0NyBu4buvIGtow7RuZyB0aW4gaG/hurdjIGtow7RuZyDEkcawYSByYSBxdXnhur90IMSR4buLbmg7IDM3NSB0aW4uIFbhuqVuIMSR4buBIGzDoCBjw7MgdOG7k24gdOG6oWkgbeG7kWkgbGnDqm4gaOG7hyBnaeG7r2EgZ2nhu69hIGdp4bubaSB0w61uaCB2w6Agbmnhu4FtIHRpbiB2w6BvIHRo4bq/IGdp4bubaSBiw6puIGtpYSBraMO0bmc/DQoNClRyb25nIG3hu5FpIHF1YW4gaOG7hyBuw6B5LCBiaeG6v24gR2nhu5tpIHTDrW5oIChjw7MgaGFpIHRodeG7mWMgdMOtbmg6IOKAnE5hbeKAnSB2w6Ag4oCcTuG7r+KAnSkgbMOgIGJp4bq/biBnaeG6o2kgdGjDrWNoLCBiaeG6v24g4oCcTmnhu4FtIHRpbiB2w6BvIHRo4bq/IGdp4bubaSBiw6puIGtpYeKAnSAoY8OzIDIgdGh14buZYyB0w61uaDog4oCcVGlu4oCdIHbDoCDigJxLaMO0bmcgdGlu4oCdKSBsw6AgYmnhur9uIMSRw6FwIOG7qW5nIChoYXkgYmnhur9uIHBo4bulIHRodeG7mWMpLiBE4buvIGxp4buHdSBuw7NpIHRyw6puIMSRxrDhu6NjIHRyw6xuaCBiw6B5IGLhu59pIGLhuqNuZyBuZ+G6q3Ugbmhpw6puIDIgY2hp4buBdSAyw5cyIHRoZW8gbeG7mXQgdHJvbmcgaGFpIGThuqFuZyBzYXU6DQoNCistLS0tLS0tLS0tLS0tLS0rLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLSstLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0rDQp8IEdpxqHMgWkgdGnMgW5oICAgICB8IE5pw6rMgG0gdGluIHZhzIBvIHRow6rMgSBnacahzIFpIGLDqm4ga2lhIHwgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB8DQorLS0tLS0tLS0tLS0tLS0tKy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0rLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tKw0KfCAgICAgICAgICAgICAgIHwgVGluICAgICAgICAgICAgICAgICAgICAgICAgICAgfCBLaMO0bmcgdGluIGhvxIPMo2Mga2jDtG5nIHF1ecOqzIF0IMSRacyjbmggfA0KKy0tLS0tLS0tLS0tLS0tLSstLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tKy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLSsNCnwgTsawzIMgICAgICAgICAgICB8IDQzNSAgICAgICAgICAgICAgICAgICAgICAgICAgIHwgMTQ3ICAgICAgICAgICAgICAgICAgICAgICAgICAgICB8DQp8ICAgICAgICAgICAgICAgfCAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB8ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgfA0KfCBOYW0gICAgICAgICAgIHwgMzc1ICAgICAgICAgICAgICAgICAgICAgICAgICAgfCAxMzQgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHwNCistLS0tLS0tLS0tLS0tLS0rLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLSstLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0rDQoNCjogTmd1w7TMgG46IERhdGEgZnJvbSAxOTkxIEdlbmVyYWwgU29jaWFsIFN1cnZlcnlcDQpCYcyJbmcgMi4yLiBLw6rMgXQgcXVhzIkgxJFpw6rMgHUgdHJhIHbDqsyAIG5pw6rMgG0gdGluIHZhzIBvIHRow6rMgSBnacahzIFpIGLDqm4ga2lhIHRoZW8gZ2nGocyBaSB0acyBbmgNCg0KSG/Eg8yjYw0KDQp8IEdpxqHMgWkgdGnMgW5oIFxcIE5pw6rMgG0gdGluIHwgQ2/MgSAgfCBLaMO0bmcgfA0KfC0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tfC0tLS0tfC0tLS0tLS18DQp8IE7GsMyDICAgICAgICAgICAgICAgICAgICB8IDQzNSB8IDE0NyAgIHwNCnwgTmFtICAgICAgICAgICAgICAgICAgIHwgMzc1IHwgMTM0ICAgfA0KDQo6IEJhzIluZyAyLjMNCg0KIyMjICoqUGjDom4gcGjhu5FpIMSR4buTbmcgdGjhu51pLCBwaMOibiBwaOG7kWkgYmnDqm4gZHV5w6puIHbDoCB4w6FjIHN14bqldCBjw7MgxJFp4buBdSBraeG7h24qKg0KDQpUcsaw4bubYyB0acOqbiBjaMO6bmcgdGEgdHLDrG5oIGLDoHkgbeG7mXQgc+G7kSBraMOhaSBuaeG7h24gY8ahIGLhuqNuIGPhu6dhIHjDoWMgc3XhuqV0IGNobyBjw6FjIGLhuqNuZyBuZ+G6q3Ugbmhpw6puIGhhaSBjaGnhu4F1LiBHaeG6oyBz4butIHLhurFuZyBt4buXaSBxdWFuIHPDoXQgdHJvbmcgbeG6q3UgxJHGsOG7o2MgY2jhu41uIG5n4bqrdSBuaGnDqm4gdOG7qyBt4buZdCB04buVbmcgdGjhu4MgdsOgIG3hu5dpIHF1YW4gc8OhdCBn4bqvbiBsaeG7gW4gduG7m2kgMiBiaeG6v24gbmfhuqt1IG5oacOqbiBuaOG7iyBwaMOibiBYIHbDoCBZLiDEkOG6t3QgJFxwaV97aWp9PVAoWD1BX2ksWT1CX2opJCBsYcyAIHjDoWMgc3XhuqV0IMSR4buDIChYLFkpIG5o4bqtbiBnacOhIHRy4buLICQoQV9pLCBCX2opJC4gJFxwaV97aWp9JCBn4buNaSBsw6AgcGjDom4gcGjhu5FpIHjDoWMgc3XhuqV0IMSR4buTbmcgdGjhu51pIGPhu6dhIFggdsOgIFkgdsOgIHRo4buPYSAkXHN1bV97aWp9XHBpX3tpan09MSQuDQoNClBow6JuIHBo4buRaSB4w6FjIHN14bqldCBiacOqbiBkdXnDqm4gbMOgIHThu5VuZyB0aGVvIHThu6tuZyBow6BuZyB2w6AgdOG7q25nIGPhu5l0IGPhu6dhIGLhuqNuZyBwaMOibiBwaOG7kWkgeMOhYyBzdeG6pXQgxJHhu5NuZyB0aOG7nWkgY+G7p2EgWCB2w6AgWSB2w6AgxJHGsOG7o2Mga8O9IGhp4buHdSBsw6AgeyR7XHBpX3tpK319JH0gdsOgIHske1xwaV97K2p9fSR9LCBk4bqldSArIGvDrSBoaeG7h3UgY2hvIHZp4buHYyB04buVbmcgdGhlbyBow6BuZyBob+G6t2MgY+G7mXQgdMO5eSB2w6BvIHbhu4sgdHLDrSBj4bunYSBuw7MgbuG6sW0g4bufIHbhu4sgdHLDrSBow6BuZyBoYXkgY+G7mXQsIG5naMSpYSBsw6A6ICRccGlfe2krfT1ccGlfe2kxfStccGlfe2kyfSQgdmHMgCAkXHBpX3sran09XHBpX3sxan0rXHBpX3syan0kDQoNClBow6JuIHBo4buRaSDEkeG7k25nIHRo4budaSB2w6AgcGjDom4gcGjhu5FpIGJpw6puIGR1ecOqbiBj4bunYSBjw6FjIGJp4bq/biBuaOG7iyBwaMOibiBYLCBZIMSRxrDhu6NjIHRo4buDIGhp4buHbiBxdWEgYuG6o25nIHNhdToNCg0KfCBYXFxZICAgICAgIHwgJEJfMSQgICAgICB8ICRCXzIkICAgICAgfCAkXHBpX3tpK30kIHwNCnwtLS0tLS0tLS0tLS18LS0tLS0tLS0tLS0tfC0tLS0tLS0tLS0tLXwtLS0tLS0tLS0tLS18DQp8ICRBXzEkICAgICAgfCAkXHBpX3sxMX0kIHwgJFxwaV97MTJ9JCB8ICRccGlfezErfSQgfA0KfCAkQV8xJCAgICAgIHwgJFxwaV97MjF9JCB8ICRccGlfezIyfSQgfCAkXHBpX3syK30kIHwNCnwgJFxwaV97K2p9JCB8ICRccGlfeysxfSQgfCAkXHBpX3srMn0kIHwgMSAgICAgICAgICB8DQoNCjogQmHMiW5nIDIuNA0KDQpUxrDGoW5nIHThu7EgY2jDum5nIHRhIGvDvSBoaeG7h3UgJG5fe2lqfSQgbMOgIHPhu5EgcXVhbiBzw6F0IGPDsyBiaeG7g3UgaGnhu4duICRBX2kkIMSR4buRaSB24bubaSBYIHbDoCAkQl9qJCDEkeG7kWkgduG7m2kgWSwga8O9IGhp4buHdSAkZl97aWp9ID0gbl97aWp9L24kIHbhu5tpICRuID1cc3VtIG5fe2lqfSQgLiBC4bqjbmcgc2F1IMSRw6J5IGfhu41pIGzDoCBi4bqjbmcgcGjDom4gcGjhu5FpIG3huqt1IMSR4buTbmcgdGjhu51pIHbDoCBwaMOibiBwaOG7kWkgbeG6q3UgYmnDqm4gZHV5w6puIGPhu6dhIGPDoWMgYmnhur9uIG5o4buLIHBow6JuIFgsIFk6DQoNCnwgWFxcWSAgICAgfCAkQl8xJCAgICB8ICRCXzIkICAgIHwgJGZfe2krfSQgfA0KfC0tLS0tLS0tLS18LS0tLS0tLS0tLXwtLS0tLS0tLS0tfC0tLS0tLS0tLS18DQp8ICRBXzEkICAgIHwgJGZfezExfSQgfCAkZl97MTJ9JCB8ICRmX3sxK30kIHwNCnwgJEFfMSQgICAgfCAkZl97MjF9JCB8ICRmX3syMn0kIHwgJGZfezIrfSQgfA0KfCAkZl97K2p9JCB8ICRmX3srMX0kIHwgJGZfeysyfSQgfCAxICAgICAgICB8DQoNCjogQmHMiW5nIDIuNQ0KDQpUcm9uZyBuaGnhu4F1IGLhuqNuZyBuZ+G6q3Ugbmhpw6puLCBt4buZdCBiaeG6v24gKHRoxrDhu51uZyBsw6AgYmnhur9uIGPhu5l0LCBZKSBsw6AgYmnhur9uIMSRw6FwIOG7qW5nIHbDoCBiaeG6v24ga2jDoWMgKGJp4bq/biBow6BuZywgWCkgbMOgIGJp4bq/biBnaeG6o2kgdGjDrWNoLCB24bubaSBt4buXaSBnacOhIHRy4buLIGPhu6dhIFggbMOgIHRow7RuZyB0aW4gxJHhu4MgeMOieSBk4buxbmcgbeG7mXQgcGjDom4gcGjhu5FpIHjDoWMgc3XhuqV0IHJpw6puZyBjaG8gWSwgcGjDom4gcGjhu5FpIG5oxrAgduG6rXkgZ+G7jWkgbMOgIHBow6JuIHBo4buRaSB4w6FjIHN14bqldCBjw7MgxJFp4buBdSBraeG7h24gY+G7p2EgWSB0xrDGoW5nIOG7qW5nIHbhu5tpIHThu6tuZyBnacOhIHRy4buLIGPhu6dhIFguDQoNClbGocyBaSBoYWkgYmnDqsyBbiBuaGnMoyBwaMOibiBjdcyJYSAkWCwgWSQgbm/MgWkgdHLDqm46DQoNClwtIFBow6JuIHBo4buRaSBjw7MgxJFp4buBdSBraeG7h24gY+G7p2EgWSB24bubaSDEkWnhu4F1IGtp4buHbiAkWCA9IEFfMSQgbMOgOg0KDQp8IFkgICB8ICRCXzEkICAgICAgICAgICAgICAgfCAkQl8yJCAgICAgICAgICAgICAgIHwgJFxzdW0kIHwNCnwtLS0tLXwtLS0tLS0tLS0tLS0tLS0tLS0tLS18LS0tLS0tLS0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tfA0KfCAkUCQgfCAkXHBpX3sxMX0vXHBpX3sxK30kIHwgJFxwaV97MTJ9L1xwaV97MSt9JCB8IDEgICAgICB8DQoNCi1QaMOibiBwaMO0zIFpIG3DosyDdSBjb8yBIMSRacOqzIB1IGtpw6rMo24gY3XMiWEgJFkkIHbGocyBaSDEkWnDqsyAdSBracOqzKNuICRYPUFfMSQgbGHMgDoNCg0KfCBZICAgfCAkQl8xJCAgICAgICAgICAgfCAkQl8yJCAgICAgICAgICAgfCAkXHN1bSQgfA0KfC0tLS0tfC0tLS0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tfA0KfCAkUCQgfCAkZl97MTF9L2ZfezErfSQgfCAkZl97MTJ9L2ZfezErfSQgfCAxICAgICAgfA0KDQotUGjDom4gcGjDtMyBaSBjb8yBIMSRacOqzIB1IGtpw6rMo24gY3XMiWEgWSB2xqHMgWkgxJFpw6rMgHUga2nDqsyjbiAkWD1BXzIkIGxhzIA6DQoNCnwgWSAgIHwgJEJfMSQgICAgICAgICAgICAgICB8ICRCXzIkICAgICAgICAgICAgICAgfCAkXHN1bSQgfA0KfC0tLS0tfC0tLS0tLS0tLS0tLS0tLS0tLS0tLXwtLS0tLS0tLS0tLS0tLS0tLS0tLS18LS0tLS0tLS18DQp8ICRQJCB8ICRccGlfezIxfS9ccGlfezIrfSQgfCAkXHBpX3syMn0vXHBpX3syK30kIHwgMSAgICAgIHwNCg0KLVBow6JuIHBow7TMgWkgbcOizIN1IGNvzIEgxJFpw6rMgHUga2nDqsyjbiBjdcyJYSAkWSQgdsahzIFpIMSRacOqzIB1IGtpw6rMo24gJFg9QV8yJCBsYcyAOg0KDQp8IFkgICB8ICRCXzEkICAgICAgICAgICB8ICRCXzIkICAgICAgICAgICB8ICRcc3VtJCB8DQp8LS0tLS18LS0tLS0tLS0tLS0tLS0tLS18LS0tLS0tLS0tLS0tLS0tLS18LS0tLS0tLS18DQp8ICRQJCB8ICRmX3syMX0vZl97Mit9JCB8ICRmX3syMn0vZl97Mit9JCB8IDEgICAgICB8DQoNCipWw60gZOG7pSAyOiogVOG7qyBi4bqjbmcgMi4yLCBsw6AgbeG7mXQgYuG6o25nIG5n4bqrdSBuaGnDqm4gMiBjaGnhu4F1IHbhu5tpIHThu5VuZyBz4buRIHF1YW4gc8OhdCBsw6AgMTA5MSDEkcaw4bujYyDEkeG7i25oIHTDrW5oIHRoZW8gZ2nhu5tpIHTDrW5oIFggdsOgIG5p4buBbSB0aW4gWSBj4bunYSBo4buNIHNhdSBraGkgY2jhur90LCB0YSBjw7MgYuG6o25nOg0KDQp8IFhcXFkgICAgIHwgQ2/MgSAgICAgICAgICAgfCBLaMO0bmcgICAgICAgIHwgJG5fe2krfSQgICAgIHwNCnwtLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tfA0KfCBOxrDMgyAgICAgICB8IDQzNSAgICAgICAgICB8IDE0NyAgICAgICAgICB8ICRuX3sxK309NTgyJCB8DQp8IE5hbSAgICAgIHwgMzc1ICAgICAgICAgIHwgMTM0ICAgICAgICAgIHwgJG5fezIrfT01MDkkIHwNCnwgJG5feytqfSQgfCAkbl97KzF9PTgxMCQgfCAkbl97KzJ9PTI4MSQgfCAkbj0xMDkxJCAgICAgfA0KDQo6IEJhzIluZyAyLjYNCg0KVMOtbmggdGhlbyBjw7RuZyB0aOG7qWM6ICRmX3tpan09XGZyYWN7bl97aWp9fXtufSQsIHRhIG5o4bqtbiDEkcaw4bujYyBt4buZdCB0aOG7gyBoaeG7h24gY+G7p2EgcGjDom4gcGjhu5FpIG3huqt1IMSR4buTbmcgdGjhu51pIHbDoCBwaMOibiBwaOG7kWkgbeG6q3UgYmnDqm4gZHV5w6puIGPhu6dhIGPDoWMgYmnhur9uIFggKGdp4bubaSB0w61uaDogTmFtLCBO4buvKSB2w6AgWSAobmnhu4FtIHRpbjogQ8OzLCBLaMO0bmcpOg0KDQp8IFhcXFkgICAgIHwgQ2/MgSAgICAgICAgICAgICAgfCBLaMO0bmcgICAgICAgICAgIHwgJGZfe2krfSQgICAgICAgIHwNCnwtLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tLS0tfA0KfCBOxrDMgyAgICAgICB8ICRmX3sxMX09MCwzOTg3JCB8ICRmX3sxMn09MCwxMzQ4JCB8ICRmX3sxK309MCw1MzM1JCB8DQp8IE5hbSAgICAgIHwgJGZfezIxfT0wLDM0MzckIHwgJGZfezIyfT0wLDEyMjgkIHwgJGZfezIrfT0wLDQ2NjUkIHwNCnwgJGZfeytqfSQgfCAkZl97KzF9PTAsNzQyNCQgfCAkZl97KzJ9PTAsMjU3NiQgfCAxICAgICAgICAgICAgICAgfA0KDQo6IEJhzIluZyAyLjcNCg0KQ2jDum5nIHRhIHPhur0gdMOsbSBwaMOibiBwaOG7kWkgbeG6q3UgY8OzIMSRaeG7gXUga2nhu4duIGPhu6dhIG5p4buBbSB0aW4gdsOgbyB0aOG6vyBnaeG7m2kgYsOqbiBraWEgduG7m2kgZ2nhu5tpIHTDrW5oLg0KDQpcLSDEkOG7kWkgduG7m2kgcGjhu6UgbuG7rywgdOG7tyBs4buHIHRy4bqjIGzhu51pIOKAnGPDsyB0aW4gdsOgbyB0aOG6vyBnaeG7m2kgYsOqbiBraWHigJ0gbMOgICQ0MzUvNTgyID0gMCw3NDczJDsgdsOgIHRy4bqjIGzhu51pIOKAnGtow7RuZyB0aW4gdsOgbyB0aOG6vyBnaeG7m2kgYsOqbiBraWHigJ0gbMOgICQxNDcvNTgyID0gMCwyNTI3JC4NCg0KXC0gxJDhu5FpIHbhu5tpIMSRw6BuIMO0bmcsIHThu7cgbOG7hyB0cuG6oyBs4budaSDigJxjw7MgdGluIHbDoG8gdGjhur8gZ2nhu5tpIGLDqm4ga2lh4oCdIGzDoCAkMzc1LzUwOSA9IDAsNzM2NyQ7IHbDoCB0cuG6oyBs4budaSDigJxraMO0bmcgdGluIHbDoG8gdGjhur8gZ2nhu5tpIGLDqm4ga2lh4oCdIGzDoCAkMTQ3LzU4MiA9IDAsMjYzMyQuDQoNCiMjIyAqKlPhu7EgxJHhu5ljIGzhuq1wKioNCg0KSGFpIGJp4bq/biDEkcaw4bujYyBjaG8gbMOgIMSR4buZYyBs4bqtcCB24buBIG3hurd0IHRo4buRbmcga8OqIG7hur91IGPDoWMgcGjDom4gcGjhu5FpIGPDsyDEkWnhu4F1IGtp4buHbiBj4bunYSAkWSQgbMOgIGtow7RuZyB0aGF5IMSR4buVaSBraGkgJFgkIHRoYXkgxJHhu5VpIGdpw6EgdHLhu4suIEtoaSBoYWkgYmnhur9uIMSR4buZYyBs4bqtcCwgeMOhYyBzdeG6pXQgY+G7p2EgYuG6pXQga+G7syBwaOG6o24g4bupbmcgbsOgbyBj4bunYSBj4buZdCAkaiQgZ2nhu5FuZyBuaGF1IOG7nyBt4buXaSBow6BuZy4gVsOtIGThu6UsIG5p4buBbSB0aW4gdsOgbyBt4buZdCB0aOG6vyBnaeG7m2kgYsOqbiBraWEgbMOgIGtow7RuZyBwaOG7pSB0aHXhu5ljIHbDoG8gZ2nhu5tpIHTDrW5oIG7hur91IHjDoWMgc3XhuqV0IHbDoG8gdGjhur8gZ2nhu5tpIGLDqm4ga2lhIGPhu6dhIGPhuqMgbmFtIHbDoCBu4buvIGzDoCAwLDczNjcuDQoNCktoaSBj4bqjIGhhaSBiaeG6v24gbMOgIGPDoWMgYmnhur9uIMSRw6FwIOG7qW5nLCBuZ8aw4budaSB0YSBjw7MgdGjhu4MgbcO0IHThuqMgbeG7kWkgcXVhbiBo4buHIGPhu6dhIGNow7puZyBi4bqxbmcgY8OhY2ggc+G7rSBk4bulbmcgcGjDom4gcGjhu5FpIMSR4buTbmcgdGjhu51pIGPhu6dhIGNow7puZywgaG/hurdjIHBow6JuIHBo4buRaSBjw7MgxJFp4buBdSBraeG7h24gY+G7p2EgJFkkIHRoZW8gJFgkIGhv4bq3YyBwaMOibiBwaOG7kWkgY8OzIMSRaeG7gXUga2nhu4duIGPhu6dhICRYJCB0aGVvICRZJC4gVMOtbmggxJHhu5ljIGzhuq1wIHRo4buRbmcga8OqIHPhur0gdMawxqFuZyDEkcawxqFuZyB24bubaSB0w61uaCBjaOG6pXQgeMOhYyBzdeG6pXQgxJHhu5NuZyB0aOG7nWkgYuG6sW5nIHbhu5tpIHTDrWNoIGPDoWMgeMOhYyBzdeG6pXQgYmnDqm4gbmdoxKlhIGzDoA0KDQokJA0KXHBpX3tpan09XHBpX3tpK30uXHBpX3sran0s4oiAaSxqDQokJA0KDQojIyMgKipQaMOibiBwaOG7kWkgxJFhIHRo4bupYyoqDQoNCk3DtCBow6xuaCBwaMOibiBwaOG7kWkgTmjhu4sgdGjhu6ljIHbDoCBtw7QgaMOsbmggcGjDom4gcGjhu5FpIFBvaXNzb24gxJHDoyDEkcaw4bujYyB0csOsbmggYsOgeSB0cm9uZyBjaMawxqFuZyAxLiBUcm9uZyBwaOG6p24gbsOgeSBz4bq9IGdp4bubaSB0aGnhu4d1IHbhu4EgcGjDom4gcGjhu5FpIMSRYSB0aOG7qWMuIFRoZW8gbeG7mXQgbmdoxKlhIG7DoG8gxJHDsyB0aMOsIHBow6JuIHBo4buRaSDEkWEgdGjhu6ljIGzDoCBz4buxIG3hu58gcuG7mW5nIGPhu6dhIHBow6JuIHBo4buRaSBOaOG7iyB0aOG7qWMuDQoNCioqTcO0IGjDrG5oIHBow6JuIHBo4buRaSDEkWEgdGjhu6ljKioNCg0KWMOpdCBt4buZdCBkw6N5ICRuJCBwaMOpcCB0aOG7rSDEkeG7mWMgbOG6rXAgLCB0cm9uZyBt4buXaSBwaMOpcCB0aOG7rSBjw7MgbeG7mXQgdsOgIGNo4buJIG3hu5l0IHRyb25nICRrJCBz4buxIGtp4buHbiAkQV8xLEFfMiwuLi4sQV9rJCB44bqjeSByYSB24bubaSB4w6FjIHN14bqldCB0xrDGoW5nIOG7qW5nICRwXzEsIHBfMiwgLi4uLCBwX2ssIChwXzErcF8yKy4uLitwX2spJCAodOG7qWMgbMOgIHRyb25nIG3hu5dpIHBow6lwIHRo4butLCAkayQgYmnhur9uIGPhu5EgbsOgeSBs4bqtcCB0aMOgbmggbeG7mXQgaOG7hyDEkeG6p3kgxJHhu6cpLiBH4buNaSAkWF9pJCBsw6Agc+G7kSBs4bqnbiB4deG6pXQgaGnhu4duIHPhu7Ega2nhu4duICRBX2kgIChpPTEsIDIsIC4uLiwgaykkIHRyb25nICRuJCBs4bqnbiB0aOG7rS4gS2hpIMSRw7MgbHXhuq10IHBow6JuIHBo4buRaSBj4bunYSB2w6ljIHTGoSBuZ+G6q3Ugbmhpw6puICRYPShYXzEsIFhfMiwgLi4uLCBYX2spJCDEkcaw4bujYyBn4buNaSBsw6AgcGjDom4gcGjhu5FpIMSRYSB0aOG7qWMgJGskIGNoaeG7gXUgduG7m2kgY8OhYyB0aGFtIHPhu5EgJG4sIHBfMSwgcF8yLC4uLiwgcF97ay0xfSRcDQoNCioqUGjDom4gcGjhu5FpIHjDoWMgc3XhuqV0IGPhu6dhIHbDqWMgdMahIG5n4bqrdSBuaGnDqm4gY8OzIHBow6JuIHBo4buRaSDEkWEgdGjhu6ljIGsgY2hp4buBdS4qKg0KDQpHaeG6oyBz4butICRYID0gKFhfMSwgWF8yLOKApiAsIFhfaykkIGzDoCB2w6ljIHTGoSBuZ+G6q3Ugbmhpw6puIGPDsyBwaMOibiBwaOG7kWkgxJFhIHRo4bupYyBrIGNoaeG7gXUuIFRhIGJp4bq/dCBy4bqxbmcgbeG7l2kgdGjDoG5oIHBo4bqnbiAkWF9pIH4gQihuLHBfaSkgKGk9MSwyLC4uLixrKSQgdmHMgDogJFhfMStYXzIrLi4uK1hfayA9biQsIHbDrCB0aOG6vyBtw6AgdsOpYyB0xqEgJFggPSAoWF8xLCBYXzIs4oCmICwgWF9rKSQgbMOgIHbDqWMgdMahIG5n4bqrdSBuaGnDqm4gcuG7nWkgcuG6oWMgayBjaGnhu4F1LCBt4buXaSBnacOhIHRy4buLIGPhu6dhIG7DsyBsw6AgbeG7mXQgxJFp4buDbSBrIGNoaeG7gXUgY8OzIHThu41hIMSR4buZIG5ndXnDqm4ga2jDtG5nIMOibTogJChyXzEsIHJfMiwgLi4uLCByX2spJCwgc2FvIGNobzogJHJfMStyXzIrLi4uK3JfayA9biQuIE3hu5dpIGvhur90IGPhu6VjIHRodeG6rW4gbOG7o2kgY2hvIGJp4bq/biBj4buRICgkWF8xPXJfMSwgWF8yPXJfMiwuLi4sWF9rID1yX2spJCBsw6AgbeG7mXQgZMOjeSBn4buTbSBuIHPhu7Ega2nhu4duIGxpw6puIGvhur90IHbhu5tpIG5oYXUgYuG7n2kgcGjDqXAgZ2lhbywgdHJvbmcgxJHDszogc+G7sSBraeG7h24gJEFfMSQgeHXhuqV0IGhp4buHbiAkcl8xJCBs4bqnbiwgc+G7sSBraeG7h24gJEFfMiQgeHXhuqV0IGhp4buHbiAkcl8yJCBs4bqnbiwgLOKApiwgc+G7sSBraeG7h24gJEFfayQgeHXhuqV0IGhp4buHbiAkcl9rJCBs4bqnbi4gTmjGsCB24bqteSBz4buRIGvhur90IGPhu6VjIHRodeG6rW4gbOG7o2kgY2hvIGJp4bq/biBj4buRICQoWF8xPXJfMSwgWF8yPXJfMiwuLi4sWF9rID1yX2spJCBjaMOtbmggbMOgIHPhu5EgY8OhY2ggY2jhu41uICRyXzEkIHbhu4sgdHLDrSBjaG8gJEFfMSwgcl8yJCB24buLIHRyw60gY2hvICRBXzIsIFxsZG90cywgcl9rJCB24buLIHRyw60gY2hvICRBX2skIHRyb25nIGTDo3kgJFxtYXRocm17bn0kIHbhu4sgdHLDrSBuw7NpIHRyw6puLCB04bupYyBsw6AgYuG6sW5nOiAkJA0KQ19uXntyXzF9IFxjZG90IENfe24tcl8xfV57cl8yfSBcY2RvdCBDX3tuLXJfMS1yXzJ9XntyXzN9IFxsZG90cyAuIENfe24tcl8xLXJfMi1cY2RvdHMgcl97ay0yfX1ee3Jfe2stMX19PVxmcmFje24hfXtyX3sxfSFyX3syfSFcbGRvdHMgcl97a30hfQ0KJCQNCg0KTeG6t3Qga2jDoWMgZG8gY8OhYyBs4bqnbiB0aOG7rSDEkeG7mWMgbOG6rXAgbmhhdSBuw6puIG3hu5dpIGvhur90IGPhu6VjIG7DoHkgxJHhu4F1IGPDsyB4w6FjIHN14bqldCBsw6A6DQoNCiRwXzFee3JfMX0gXGNkb3QgcF8yXntyXzJ9IFxsZG90cyAuIC4gcF9rXntyX2t9JA0KDQpU4burIMSRw7MgdGEgbmjhuq1uIMSRxrDhu6NjIHBow6JuIHBo4buRaSB4w6FjIHN14bqldCBj4bunYSB2w6ljIHTGoSBuZ+G6q3Ugbmhpw6puIGPDsyBwaMOibiBwaOG7kWkgxJFhIHRo4bupYyAkayQgY2hp4buBdToNCg0KJFBcbGVmdChYXzE9cl8xLCBYXzI9cl8yLCBcbGRvdHMsIFhfaz1yX2tccmlnaHQpPVxmcmFje24hfXtyX3sxfSFyX3syfSFcbGRvdHMgcl97a30hfSBcY2RvdCBwXzFee3JfMX0gXGNkb3QgcF8yXntyXzJ9IFxsZG90cyAuIC4gcF9rXntyX2t9JCAoMi4xLjEpDQoNCnbhu5tpIG3hu41pIMSRaeG7g20gJFxsZWZ0KHJfMSwgcl8yLCBcbGRvdHMsIHJfa1xyaWdodCkkIGPDsyB04buNYSDEkeG7mSBuZ3V5w6puLCBraMO0bmcgw6JtIG3DoDogJHJfMStyXzIrXGNkb3RzK3Jfaz1uJC4NCg0KKioqTmjhuq1uIHjDqXQ6KioqDQoNClwtIMSQ4buRaSB24bubaSB2w6ljIHTGoSBuZ+G6q3Ugbmhpw6puICRYPVxsZWZ0KFhfMSwgWF8yLCBcbGRvdHMsIFhfa1xyaWdodCkkIGPDsyBwaMOibiBwaOG7kWkgxJFhIHRo4bupYyAkayQgY2hp4buBdSwgduG7m2kgY8OhYyB0aGFtIHPhu5EgJG4sIHBfMSwgcF8yLCBcbGRvdHMsIHBfe2stMX0kLCB0YSBjw7MgdsOpYyB0xqEga+G7syB24buNbmc6DQoNCiRFIFg9XGxlZnQoRSBYXzEsIEUgWF8yLCBcbGRvdHMsIEUgWF9rXHJpZ2h0KT1cbGVmdChuIHBfMSwgbiBwXzIsIFxsZG90cywgbiBwX2tccmlnaHQpJA0KDQpcLSBQaMOibiBwaOG7kWkgxJFhIHRo4bupYyB24bubaSAkaz0yJCwgbMOgIHBow6JuIHBo4buRaSDEkeG7k25nIHRo4budaSBj4bunYSAyIGJp4bq/biBuaOG7iyB0aOG7qWMgJFhfMSQsICRuLVhfMSQsIHRyb25nIMSRw7MgJFhfMSBcc2ltIFxtYXRocm17Qn1cbGVmdChcbWF0aHJte259LCBwXzFccmlnaHQpJCwgJFxsZWZ0KFxtYXRocm17bn0tWF8xXHJpZ2h0KSBcc2ltIFxtYXRocm17Qn1cbGVmdChcbWF0aHJte259LCAxLVxtYXRocm17cH1fMVxyaWdodCkkLg0KDQoqVsOtIGThu6UgMyouIE3hu5dpIGtow6FjaCBow6BuZyB2w6BvIGjhu4cgZOG7i2NoIHbhu6UgQSBjw7MgdGjhu4MgY2jhu41uIG3hu5l0IHRyb25nIDMgbeG7qWMgcGjDrSBwaOG7pWMgduG7pTogMTAwIG5nw6BuIMSR4buTbmcsIDE1MCBuZ8OgbiDEkeG7k25nIHbDoCAyMDAgbmfDoG4gxJHhu5NuZy4gxJDGsOG7o2MgYmnhur90IGzGsOG7o25nIGtow6FjaCBjaOG7jW4gY8OhYyBt4bupYyBwaMOtIG7DoHkgdMawxqFuZyDhu6luZyB0aGVvIHThu7cgbOG7hzogNTogMzogMi4gQ8OzIDMga2jDoWNoIGjDoG5nIHbDoG8gaOG7hyBk4buLY2ggduG7pSBBIHbDoCBo4buNIMSR4buZYyBs4bqtcCBuaGF1IHRyb25nIHZp4buHYyBjaOG7jW4gcGjDrSBk4buLY2ggduG7pS4gVMOsbSB4w6FjIHN14bqldCDEkeG7gyB0cm9uZyBz4buRIMSRw7MgY8OzIMOtdCBuaOG6pXQgMiBraMOhY2ggY2jhu41uIG3hu6ljIHBow60gMjAwIG5nw6BuIMSR4buTbmcuDQoNCkdp4bqjaTogR+G7jWkgJFhfMSwgWF8yLCBYXzMkIGzhuqduIGzGsOG7o3QgbMOgIHPhu5Ega2jDoWNoIGNo4buNbiBt4bupYyBwaMOtIDEwMCBuZ8OgbiDEkeG7k25nLCAxNTAgbmfDoG4gxJHhu5NuZywgMjAwIG5nw6BuIMSR4buTbmcgdHJvbmcgc+G7kSAzIGtow6FjaCBuw7NpIHRyw6puLiBLaGkgxJHDsyB2w6ljIHTGoSAkWD0oWF8xLCBYXzIsIFhfMykkIGzDoCB2w6ljIHTGoSBuZ+G6q3Ugbmhpw6puIGPDsyBwaMOibiBwaOG7kWkgxJFhIHRo4bupYyAzIGNoaeG7gXUuIEvDvSBoaeG7h3UgUyBsw6AgYmnhur9uIGPhu5EgdHJvbmcgc+G7kSAzIGtow6FjaCBow6BuZyBjw7Mgw610IG5o4bqldCAyIGtow6FjaCBjaOG7jW4gbeG7qWMgcGjDrSAyMDAgbmfDoG4gxJHhu5NuZywgdGEgY8OzIGJp4buDdSBkaeG7hW46DQoNCiRTPShYXzE9MCwgWF8yPTEsIFhfMz0yKSBcY3VwKFhfMT0xLCBYXzI9MCwgWF8zPTIpIFxjdXAoWF8xPTAsIFhfMj0wLCBYXzM9MykkDQoNClThu6sgxJHDsyBzdXkgcmEgeMOhYyBzdeG6pXQgY+G6p24gdMOtbmg6DQoNCiRQKFMpPVxmcmFjezMhfXswITEhMiF9IFxjZG90IDAsMyBcY2RvdCAwLDJeMitcZnJhY3szIX17MSEwITIhfSBcY2RvdCAwLDUgXGNkb3QgMCwyXjIrXGZyYWN7MyF9ezAhMCEzIX0gXGNkb3QgMCwyXjM9MCwkDQoNCiMjIyAqKkzhuqV5IG3huqt1IGNobyBwaMOibiBwaOG7kWkgUG9pc3NvbiwgTmjhu4sgdGjhu6ljIHbDoCDEkGEgdGjhu6ljKioNCg0KQ8OhYyBtw7QgaMOsbmggbOG6pXkgbeG6q3UgxJHGsOG7o2MgZ2nhu5tpIHRoaeG7h3UgdHJvbmcgcGjhuqduIDEuMyDEkcaw4bujYyBt4bufIHLhu5luZyByYSBjaG8gdOG6p24gc+G7kSB0cm9uZyBjw6FjIMO0IGPhu6dhIGLhuqNuZyBuZ+G6q3Ugbmhpw6puLiBWw60gZOG7pSwgbcO0IGjDrG5oIGzhuqV5IG3huqt1IFBvaXNzb24gY2hvIGLhuqNuZyBoYWkgY2hp4buBdSAkMiBcdGltZXMgMiQsIGPDoWMgw7QgdHJvbmcgYuG6o25nIMSRxrDhu6NjIHhlbSBuaMawIGzDoCBjw6FjIGJp4bq/biBQb2lzc29uIMSR4buZYyBs4bqtcC4NCg0KS2hpIGPDoWMgaMOgbmcgY+G7p2EgbeG7mXQgYuG6o25nIG5n4bqrdSBuaGnDqm4gxJHhuqFpIGRp4buHbiBjaG8gY8OhYyBuaMOzbSBraMOhYyBuaGF1LCBj4buhIG3huqt1IGNobyBjw6FjIG5ow7NtIG7DoHkgdGjGsOG7nW5nIMSRw6MgxJHGsOG7o2MgY+G7kSDEkeG7i25oIHRyxrDhu5tjIHRyb25nIGdpYWkgxJFv4bqhbiB0aGnhur90IGvhur8gbeG6q3UuIMSQ4buDIMSRxqFuIGdp4bqjbiwgY2jDum5nIHRhIHPhur0gdHLDrG5oIGLDoHkgduG7gSDEkWnhu4F1IG7DoHkgY2hvIGLhuqNuZyAkMiBcdGltZXMgMiQgY2hvIGhhaSBiaeG6v24gJFgkIChiaeG6v24gZ2nhuqNpIHRow61jaCkgdsOgICRZJCAoYmnhur9uIMSRw6FwIOG7qW5nKS4gS2hpIHThu5VuZyBz4buRIGPhuq1uIGJpw6puIGNobyBjw6FjIGJp4buDdSBoaeG7h24gY+G7p2EgJFgkIMSRxrDhu6NjIGPhu5EgxJHhu4tuaCBjaOG7qSBraMO0bmcgcGjhuqNpIG5n4bqrdSBuaGnDqm4sIG3hu5l0IHBow6JuIHBo4buRaSBjaHVuZyBjaG8gJFgkIHbDoCAkWSQga2jDtG5nIGPDsm4gw70gbmdoxKlhIG7hu69hLCBuaMawbmcgcGjDom4gcGjhu5FpIGPDsyDEkWnhu4F1IGtp4buHbiBjaG8gJFkkIOG7nyBt4buXaSBiaeG7g3UgaGnhu4duIGPhu6dhICRYJCB0aMOsIGPDsyDDvSBuZ2jEqWEuDQoNCktoaSB04buVbmcgY+G7oSBt4bqrdSB0cm9uZyBi4bqjbmcgxJHGsOG7o2MgY+G7kSDEkeG7i25oIG5oxrBuZyB04buVbmcgY+G7p2EgaMOgbmcgaG/hurdjIGPhu5l0IGtow7RuZyBj4buRIMSR4buLbmgsIG3DtCBow6xuaCBs4bqleSBt4bqrdSDEkWEgdGjhu6ljIGzDoCBwaMO5IGjhu6NwLCB0cm9uZyDEkcOzIGPDoWMgw7QgbMOgIG5o4buvbmcga+G6v3QgcXXhuqMgY8OzIHRo4buDIHjhuqN5IHJhLiBDaOG6s25nIGjhuqFuLCBi4bqjbmcgMi4xIMSR4buLbmggdMOtbmggbeG7mXQgbeG6q3Ugbmfhuqt1IG5oacOqbiBj4bunYSAxLjA5MSDEkeG7kWkgdMaw4bujbmcgdGhlbyBnaeG7m2kgdMOtbmggdsOgIG5p4buBbSB0aW4gdsOgbyB0aOG6vyBnaeG7m2kgYsOqbiBraWEuIEdpw6EgdHLhu4sgdHJvbmcgNCDDtCBsw6AgdOG6p24gc+G7kSBj4bunYSBwaMOibiBwaOG7kWkgxJFhIHRo4bupYyBi4buRbiBjaGnhu4F1Lg0KDQrEkOG7kWkgduG7m2kgY8OhYyBt4bqrdSBj4bunYSBwaMOibiBwaOG7kWkgxJFhIHRo4bupYyB0csOqbiBjw6FjIMO0IGPhu6dhIG3hu5l0IGLhuqNuZyBuZ+G6q3Ugbmhpw6puLCBjw6FjIGPhu5l0IGzDoCBiaeG6v24gxJHDoXAg4bupbmcgdsOgIGPDoWMgaMOgbmcgbMOgIGJp4bq/biBnaeG6o2kgdGjDrWNoLiBTYXUgxJHDsywgxJHhu4MgbcO0IHThuqMgZOG7ryBsaeG7h3UsIG3hu5l0IGPDoWNoIGjhu6NwIGzDvSwgdGEgcGjDom4gY2hpYSBjw6FjIMO0IMSR4bq/bSB0aGVvIHThu5VuZyBz4buRIGjDoG5nIMSR4buDIHThuqFvIHRow6BuaCBwaMOibiBwaOG7kWkgY8OzIMSRaeG7gXUga2nhu4duIHbhu4EgxJHDoXAg4bupbmcuIEtoaSBsw6BtIG5oxrAgduG6rXksIGNow7puZyB0YSB24buRbiDEkcOjIHjhu60gbMO9IGPDoWMgdOG7lW5nIHPhu5EgaMOgbmcgbMOgIGPhu5EgxJHhu4tuaCB2w6AgcGjDom4gdMOtY2ggZOG7ryBsaeG7h3UgZ2nhu5FuZyBuaMawIGtoaSBjaMO6bmcgaMOsbmggdGjDoG5oIGPDoWMgbeG6q3UgxJHhu5ljIGzhuq1wIHJpw6puZyBiaeG7h3QuIFbDrSBk4bulIG5oxrAgdHJvbmcgQuG6o25nIDIuMSwgY2jDum5nIHRhIGPDsyB0aOG7gyBjb2kga+G6v3QgcXXhuqMgY2hvIG7hu68gZ2nhu5tpIG5oxrAgbeG7mXQgbeG6q3Ugbmjhu4sgdGjhu6ljIHbhu5tpIGPDoWMgbG/huqFpIGvhur90IHF14bqjIGzDoCAiY8OzIiB2w6AgImtow7RuZyBob+G6t2Mga2jDtG5nIHF1eeG6v3QgxJHhu4tuaCIgY2hvIG5p4buBbSB0aW4gdsOgbyBt4buZdCB0aOG6vyBnaeG7m2kgYsOqbiBraWEsIHbDoCBr4bq/dCBxdeG6oyBjaG8gbmFtIGdp4bubaSBuaMawIG3hu5l0IG3huqt1IG5o4buLIHRo4bupYyDEkeG7mWMgbOG6rXAgduG7gSBjw7luZyBt4buZdCBwaOG6o24g4bupbmcuIE7hur91IGPDsyBuaGnhu4F1IGjGoW4gaGFpIGxv4bqhaSBwaOG6o24gaOG7k2ksIGNo4bqzbmcgaOG6oW4gbmjGsCAoImPDsyIsICJraMO0bmciLCAia2jDtG5nIHF1eeG6v3QgxJHhu4tuaCIpLCBjaMO6bmcgdGEgc+G6vSBjb2kgY8OhYyBt4bqrdSBsw6AgbeG6q3UgxJFhIHRo4bupYyDEkeG7mWMgbOG6rXAuDQoNCsSQ4buRaSB24bubaSBo4bqndSBo4bq/dCBjw6FjIHBow6JuIHTDrWNoLCBuZ8aw4budaSB0YSBraMO0bmcgY+G6p24gcXVhbiB0w6JtIHbhu4EgbcO0IGjDrG5oIGzhuqV5IG3huqt1IG7DoG8gY8OzIMO9IG5naMSpYSBuaOG6pXQuIMSQ4buRaSB24bubaSBjw6FjIHBoxrDGoW5nIHBow6FwIHN1eSBkaeG7hW4gY2jDrW5oIHRyb25nIHTDoGkgbGnhu4d1IG7DoHksIGPDoWMgbcO0IGjDrG5oIGzhuqV5IG3huqt1IGzDoCBQb2lzc29uLCDEkWEgdGjhu6ljLg0KDQojIyAqKlNvIHPDoW5oIHThu7cgbOG7hyB0cm9uZyBt4buZdCBi4bqjbmcgaGFpIGNoaeG7gXUqKg0KDQpCaeG6v24gxJHDoXAg4bupbmcgY8OzIGhhaSBiaeG7g3UgaGnhu4duIMSRxrDhu6NjIGfhu41pIGzDoCBiaeG6v24gbmjhu4sgcGjDom4uIFbDrSBk4bulLCAibmnhu4FtIHRpbiB2w6BvIHRo4bq/IGdp4bubaSBiw6puIGtpYSIgbMOgIG5o4buLIHBow6JuIGtoaSDEkcaw4bujYyDEkW8gYuG6sW5nIGPDoWMgYmnhu4N1IGhp4buHbiAoY8OzLCBraMO0bmcpLiBOaGnhu4F1IG5naGnDqm4gY+G7qXUgY+G6p24gc28gc8OhbmggZ2nhu69hIGhhaSBuaMOzbSDEkeG7kWkgduG7m2kgxJHDoXAg4bupbmcgbmjhu4sgcGjDom4gY+G7p2EgYmnhur9uICRZJC4gROG7ryBsaeG7h3UgY8OzIHRo4buDIMSRxrDhu6NjIGhp4buDbiB0aOG7iyB0cm9uZyBi4bqjbmcgbmfhuqt1IG5oacOqbiAkMiBcdGltZXMgMiQsIHRyb25nIMSRw7MgY8OhYyBow6BuZyBsw6AgaGFpIG5ow7NtIHbDoCBjw6FjIGPhu5l0IGzDoCBt4bupYyBwaOG6o24g4bupbmcgY+G7p2EgJFkkLiBQaOG6p24gbsOgeSB0csOsbmggYsOgeSBjw6FjIHBoxrDGoW5nIHBow6FwIMSR4buDIHNvIHPDoW5oIG3hu6ljIMSR4buZIHBo4bqjbiDhu6luZyB0csOqbiBiaeG7g24gbmjhu4sgcGjDom4uDQoNCiMjIyAqKlPhu7Ega2jDoWMgYmnhu4d0IGdp4buvYSBoYWkgdOG7tyBs4buHKioNCg0KVuG7m2kgYmnhur9uIMSRw6FwIOG7qW5nIG5o4buLIHBow6JuIFksIGNow7puZyB0YSBz4butIGThu6VuZyB0aHXhuq10IG5n4buvIGNodW5nICJ0aMOgbmggY8O0bmciIGNobyBt4buZdCDEkcOhcCDhu6luZyB2w6AgInRo4bqldCBi4bqhaSIgxJHhu5FpIHbhu5tpIMSRw6FwIOG7qW5nIGPDsm4gbOG6oWkuIMSQ4buRaSB24bubaSBjw6FjIMSR4buRaSB0xrDhu6NuZyB0cm9uZyBow6BuZyBpLCBrw70gaGnhu4d1ICRccGlfaSQgbMOgIHjDoWMgc3XhuqV0IHRow6BuaCBjw7RuZyB0aMOsICQxLVxwaV9pJCBsw6AgeMOhYyBzdeG6pXQgY+G7p2EgbeG7mXQgdGjhuqV0IGLhuqFpLiBD4bq3cCB4w6FjIHN14bqldCAkXGxlZnQoXHBpX2ksIDEtXHBpX2lccmlnaHQpJCBsw6AgcGjDom4gcGjhu5FpIHjDoWMgc3XhuqV0IGPDsyDEkWnhu4F1IGtp4buHbiBj4bunYSBZIHRyb25nIGjDoG5nIGkuDQoNCioqxJDhu4tuaCBuZ2jEqWEgMToqKiBLaMOhYyBiaeG7h3QgY+G7p2EgMiB04bu3IGzhu4c6ICRccGlfMS1ccGlfMiQgc28gc8OhbmggeMOhYyBzdeG6pXQgdGjDoG5oIGPDtG5nIGPhu6dhIGhhaSBow6BuZy4NCg0KVGEgY8OzOiAkLTEgXGxlcSBccGlfMS1ccGlfMiBcbGVxIDEkIHbDoCAkXHBpXzEtXHBpXzI9MCBcTGVmdHJpZ2h0YXJyb3cgXHBpXzE9XHBpXzIkLg0KDQoqKlRSQU5HIDI4KioNCg0KS8O9IGhp4buHdSAkZl9pJCBsw6AgdOG7tyBs4buHIG3huqt1IGPhu6dhIGJp4buDdSBoaeG7h24gInRow6BuaCBjw7RuZyIg4bufIGjDoG5nIHRo4bupICRpKGk9MSwyKSQsIHPhu7Ega2jDoWMgYmnhu4d0ICRccGlfMS1ccGlfMiQgZ2nhu69hIGhhaSB04bu3IGzhu4cgc+G6vSDEkcaw4bujYyDGsOG7m2MgbMaw4bujbmcgYuG7n2kgJGZfMS1mXzIkLg0KDQpWw6Aga2hpIGPhu6EgbeG6q3UgbiBraMOhIGzhu5tuLCBjaMO6bmcgdGEgY8OzIGPDtG5nIHRo4bupYyDGsOG7m2MgbMaw4bujbmcga2hv4bqjbmcgdGluIGPhuq15IHbhu5tpIMSR4buZIHRpbiBj4bqteSAkXGdhbW1hPTEtXGFscGhhJCwgY2hvIHPhu7Egc2FpIGtow6FjIHbhu4EgdOG7tyBs4buHOiAkXHBpXzEtXHBpXzIkIGzDoDoNCg0KJCgoZl8xLWZfMiktdV97XGZyYWN7XGFscGhhfXsyfX0gXGhhdHtcc2lnbWF9KGZfMS1mXzIpLChmXzEtZl8yKSt1X3tcZnJhY3tcYWxwaGF9ezJ9fSBcaGF0e1xzaWdtYX0oZl8xLWZfMikpJCAoMi4yLjEpDQoNCnbhu5tpOiAkXHF1YWQgXGhhdHtcc2lnbWF9XGxlZnQoZl8xLWZfMlxyaWdodCk9XHNxcnR7XGZyYWN7Zl8xXGxlZnQoMS1mXzFccmlnaHQpfXtuX3sxK319K1xmcmFje2ZfMlxsZWZ0KDEtZl8yXHJpZ2h0KX17bl97Mit9fX0kICgyLjIuMikNCg0KdsOgICR1X3tcYWxwaGEgLyAyfSQgdHJhIHRyb25nIGLhuqNuZyBwaMOibiBwaOG7kWkgY2h14bqpbiBjaMOtbmggdOG6r2MgdGhlbyBo4buHIHRo4bupYzogJFxQaGkodV9cbGFtYmRhKT0xLVxsYW1iZGEkLigkdV9cbGFtYmRhJCDEkcaw4bujYyBn4buNaSBsw6AgZ2nDoSB0cuG7iyB04bubaSBo4bqhbiBt4bupYyAkXGxhbWJkYSQgY+G7p2EgcGjDom4gcGjhu5FpIGNodeG6qW4gY2jDrW5oIHThuq9jKS4NCg0KKlbDrSBkdSA0OiogROG7sWEgdsOgbyDEkWnhu4F1IHRyYSDhu58gYuG6o25nIDIuMiwgY+G6p24gxrDhu5tjIGzGsOG7o25nIGtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIHNhaSBraMOhYyB24buBIHThu7cgbOG7hyB0aW4gdsOgbyB0aOG6vyBnaeG7m2kgYsOqbiBraWEgZ2nhu69hIHBo4bulIG7hu68gdsOgIMSRw6BuIMO0bmcgTeG7uS4NCg0KVuG7m2kgJFxnYW1tYT0xLVxhbHBoYT0wLDk1JCwgY8OzICRcYWxwaGE9MCwwNSwgdV97XGZyYWN7XGFscGhhfXsyfX09dV97MCwwMjV9PTEsOTYkLiBU4burIGLhuqNuZyAyLjIsIGPDszoNCg0KJGZfMT1cZnJhY3tuX3sxMX19e25fezErfX09MCw3NDczIDsgZl8yPVxmcmFje25fezIxfX17bl97Mit9fT0wLDczNjcgOyBcaGF0e1xzaWdtYX0oZl8xLWZfMik9MCwwMjY2OyQNCg0KJChmXzEtZl8yKS11X3tcYWxwaGEvMn1caGF0e1xzaWdtYX0oZl8xLWZfMik9LTAsMDQxNSA7KGZfMS1mXzIpK3Vfe1xhbHBoYSAvIDJ9IFxoYXR7XHNpZ21hfShmXzEtZl8yKT0wLDA2MjckDQoNClbhuq15IGtob+G6o25nIHRpbiBj4bqteSAkOTUgXCUkIGNobyBzYWkgbOG7h2NoIG7DoHkgbMOgOiAkKC0wLDA0MTUgOyAwLDA2MjcpJC4NCg0KKlbDrSBkdSA1OiogQuG6o25nIDIuOCBsw6Agc+G7kSBsaeG7h3UgYsOhbyBjw6FvIHbhu4EgbeG7kWkgcXVhbiBo4buHIGdp4buvYSBz4butIGThu6VuZyB0aHXhu5FjIGFzcGlyaW4gdsOgIGLhu4duaCBuaOG7k2kgbcOhdSBjxqEgdGltICoobXlvY2FyZGlhbCBpbmZhcmN0aW9uKSogZG8gbmjDs20gbmdoacOqbiBj4bupdSBz4bupYyBraG/hursgY+G7p2EgY8OhYyBiw6FjIHPEqSB04bqhaSBUcsaw4budbmcgeSBIYXJ2YXJkLiBOZ2hpw6puIGPhu6l1IMSRxrDhu6NjIHRo4buxYyBoaeG7h24gdHJvbmcgdGjhu51pIGdpYW4gNSBuxINtIHRyw6puIG5o4buvbmcgbmfGsOG7nWkgYuG7iyBi4buHbmggdGltIMSRxrDhu6NjIGNo4buNbiBt4buZdCBjw6FjaCBuZ+G6q3Ugbmhpw6puIHbDoCBjaGlhIGzDoG0gaGFpIG5ow7NtOiBOaMOzbSBkw7luZyBhc3BpcmluIHbDoCBuaMOzbSBkw7luZyBnaeG6oyBkxrDhu6NjICoocGxhY2VibykqLiBN4buXaSBuZ8OgeSwgY8OhYyBiw6FjIHPEqSBz4bq9IGNobyBi4buHbmggbmjDom4gdeG7kW5nIG3hu5l0IHZpw6puIGFzcGlyaW4gaG/hurdjIG3hu5l0IHZpw6puIGdp4bqjIGTGsOG7o2MuDQoNCnwgICAgICAgICB8ICAgICAgICAgICAgICAgICAgICAgICB8ICAgICAgIHwgICAgICAgfA0KfDotLS0tLS0tOnw6LS0tLS0tLS0tLS0tLS0tLS0tLS0tOnw6LS0tLS06fDotLS0tLTp8DQp8ICAgICAgICAgfCBNeW9jYXJkaWFsIEluZmFyY3Rpb24gfCAgICAgICB8ICAgICAgIHwNCnwgIEdyb3VwICB8ICAgICAgICAgIFllcyAgICAgICAgICB8ICBObyAgIHwgVG90YWwgfA0KfCBQbGFjZWJvIHwgICAgICAgICAgMTg5ICAgICAgICAgIHwgMTA4NDUgfCAxMTAzNCB8DQp8IEFzcGlyaW4gfCAgICAgICAgICAxMDQgICAgICAgICAgfCAxMDkzMyB8IDExMDM3IHwNCg0KKk5ndeG7k246IFByZWxpbWluYXJ5IFJlcG9ydDogRmluZGluZ3MgZnJvbSB0aGUgQXNwaXJpbiBDb21wb25lbnQgb2YgdGhlIE9uZ29pbmcgUGh5c2ljaWFucycgSGVhbHRoIFN0dWR5LiBOLiBFbmdsLiBKLiBNZWQuLCAzMTg6IDI2Mi0yNjQgKDE5ODgpKg0KDQoqQuG6o25nIDIuOC4gQ3Jvc3MgQ2xhc3NpZmljYXRpb24gb2YgQXNwaXJpbiBVc2UgYW5kIE15b2NhcmRpYWwgSW5mYXJjdGlvbiAoTUkpKg0KDQpDaMO6bmcgdGEgeGVtIGhhaSBow6BuZyB0cm9uZyBC4bqjbmcgMi44IG5oxrAgbMOgIGPDoWMgbeG6q3Ugbmjhu4sgcGjDom4gxJHhu5ljIGzhuq1wLCB0cm9uZyBow6BuZyAxIDogJFxtYXRocm17bn1fezErfT0xMTAzNCQgbmfGsOG7nWkgdsOgIGPDsyAxODkgbmfGsOG7nWkgYuG7iyBuaOG7k2kgbcOhdSBjxqEgdGltIHThu7cgbOG7hyBsw6AgJGZfMT0gMTg5IC8gMTEwMzQ9MCwxNzEkLCB0cm9uZyBow6BuZyAyIDogJFxtYXRocm17bn1fezIrfT0xMTAzNyQgbmfGsOG7nWkgdsOgIGPDsyAxMDQgYuG7iyBuaOG7k2kgbcOhdSBjxqEgdGltIHThu7cgbOG7hyBsw6AgJGZfMj0xMDQgLyAxMTAzNz0wLDAwOTQkLCB2w6A6ICR1X3tcZnJhY3tcYWxwaGF9ezJ9fSBcaGF0e1xzaWdtYX1cbGVmdChmXzEtZl8yXHJpZ2h0KT0xLDk2IC4oMCwwMTUpPTAsMDI5NCQuIMOBcCBk4bulbmcgY8O0bmcgdGjhu6ljIHRyw6puIHbhu5tpIMSR4buZIHRpbiBj4bqteSA5NSUgY2jDum5nIHRhIGtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvICRccGlfMS1ccGlfMiQgbMOgOiAkKDAsMDA1OzAsMDExKSQsIGhheSAkMCwwMDU8XHBpXzEtXHBpXzI8MCwwMTFcbGVmdCgwLDUgXCU8XHBpXzEtXHBpXzI8MSwxIFwlXHJpZ2h0KSQgbmdoxKlhIGzDoCBnaeG7r2EgbmfGsOG7nWkgdeG7kW5nIHbDoCBraMO0bmcgdeG7kW5nIHRodeG7kWMgYXNwYXJpbiB0aMOsIHThu7cgbOG7hyBi4buLIG5o4buTaSBtw6F1IGPGoSB0aW0gY8OzIHPhu7Ega2jDoWMgYmnhu4d0IHThu6sgMCw1IMSR4bq/biAxLDElLg0KDQojIyMgKipS4bunaSBybyB0xrDGoW5nIMSR4buRaSAqKHJlbGF0aXZlIHJpc2spKioqDQoNClPhu7Ega2jDoWMgYmnhu4d0IGdp4buvYSBoYWkgdOG7tyBs4buHICRccGlfMSQgdsOgICRccGlfMiQgxJHGsOG7o2MgxJHDoW5oIGdpw6EgcXVhIHNhaSBz4buRIHR1eeG7h3QgxJHhu5FpICRcbGVmdChccGlfMS1ccmlnaHQuJCAkXHBpXzIkICkgY2jhu4kgZOG7sWEgdsOgbyBraG/huqNuZyBjw6FjaCBnaeG7r2EgaGFpIHThu7cgbOG7hyBuw6B5IG3DoCBi4buPIHF1YSDEkeG7mSBs4bubbiB2w6Agw70gbmdoxKlhIGPhu6dhIGNow7puZywgbmdoxKlhIGzDoCBraG/huqNuZyBjw6FjaCBuaMawIG5oYXUgdGjDrCDEkcOhbmggZ2nDoSBzYWkgbOG7h2NoIG5oxrAgbmhhdS4gQ2jhurNuZyBo4bqhbiBraGkgJFxwaV8xPTAsMDEkIHbDoCAkXHBpXzI9MCwwMDEkIHbDoCBraGkgJFxwaV8xPTAsNDEwJCB2w6AgJFxwaV8yPTAsNDAxJCDEkeG7gXUgxJHGsOG7o2MgxJHDoW5oIGdpw6EgbMOgIGPDsyBzYWkgbOG7h2NoIG5oxrAgbmhhdSAoY8O5bmcgbeG7qWMgc2FpIGzhu4djaCBsw6AgMCwwMDkgKS4gVHV5IG5oacOqbiBu4bq/dSBjaOG7iSBk4buxYSB2w6BvIGtob+G6o25nIGPDoWNoIG7DoHkgxJHhu4MgxJHDoW5oIGdpw6EsIHNvIHPDoW5oIGhhaSB0xrDGoW5nIHF1YW4gZ2nhu69hIGhhaSB04bu3IGzhu4cgdGjDrCBraMO0bmcgxJHhu6csIHRo4bqtbSBjaMOtIHPhur0gYuG7jyBxdWEgbmjhu69uZyDDvSBuZ2jEqWEgcXVhbiB0cuG7jW5nIGPhu6dhIHPhu7EgdMawxqFuZyBxdWFuIGPhu6dhIGNow7puZywgxJHhurdjIGJp4buHdCBsw6Aga2hpIGhhaSB04bu3IGzhu4cgbsOgeSBjw7luZyBn4bqnbiAwIGhv4bq3YyBjw7luZyBn4bqnbiAxICwgaG/hurdjIGPDuW5nIGfhuqduIDAsNSAuDQoNClbDrSBk4bulIGtoaSBjaMO6bmcgdGEgc28gc8OhbmggJFxwaV8xJCB2w6AgJFxwaV8yJCBsw6AgdOG7tyBs4buHIGLhu4sgdMOhYyBk4bulbmcgcGjhu6UgKHBo4bqjbiDhu6luZyBi4bqldCBs4bujaSkgY+G7p2EgaGFpIGxv4bqhaSB0aHXhu5FjLCB0cm9uZyBoYWkgdHLGsOG7nW5nIGjhu6NwOg0KDQpcLSBUcsaw4budbmcgaOG7o3AgMTogU+G7sSBraMOhYyBiaeG7h3QgbuG6sW0gdHJvbmcga2hv4bqjbmcgdOG7qyAwLDAwMSDEkeG6v24gMCwwMSAuDQoNClwtIFRyxrDhu51uZyBo4bujcCAyOiBT4buxIGtow6FjIGJp4buHdCBu4bqxbSB0cm9uZyBraG/huqNuZyB04burIDAsNDAxIMSR4bq/biAwLDQxMCAuDQoNCkPhuqMgaGFpIHRyxrDhu51uZyBo4bujcCBuw6B5LCBz4buxIGtow6FjIGJp4buHdCB24buBIHThu7cgbOG7hyBi4buLIHTDoWMgZOG7pW5nIHBo4bulIGzDoCBuaMawIG5oYXUgKMSR4buBdSBsw6AgMCwwMDkkKSQuIE5oxrBuZyBu4bq/dSB0YSDEkeG7gyDDvSDEkeG6v24gdOG7tyBz4buRICRcZnJhY3tccGlfMX17XHBpXzJ9JCwgdGjDrCDhu58gdHLGsOG7nW5nIGjhu6NwICQxOiBcZnJhY3tccGlfMX17XHBpXzJ9PTEwJCwgY2hvIHRo4bqleSBraGkgYuG7h25oIG5ow6JuIMSRaeG7gXUgdHLhu4ssIHBo4bqjbiDhu6luZyBi4bqldCBs4bujaSBj4bunYSBsb+G6oWkgdGh14buRYyB0aOG7qSBuaOG6pXQgZ+G6pXAgMTAgbOG6p24gbG/huqFpIHRodeG7kWMgdGjhu6kgaGFpLCBjw7JuIOG7nyB0csaw4budbmcgaOG7o3AgMjogJFxmcmFje1xwaV8xfXtccGlfMn09MSwwMjI0JCwgY2hvIHRo4bqleSBraGkgYuG7h25oIG5ow6JuIMSRaeG7gXUgdHLhu4ssIHBo4bqjbiDhu6luZyBi4bqldCBs4bujaSBj4bunYSBsb+G6oWkgdGh14buRYyB0aOG7qSBuaOG6pXQgZ+G6pXAgMSwwMjI0IGzhuqduIGxv4bqhaSB0aHXhu5FjIHRo4bupIGhhaSwgdOG7qWMgbMOgIGfhuqduIG5oxrAgbmhhdS4gUsO1IHLDoG5nIHPhu7Ega2jDoWMgYmnhu4d0IOG7nyB0csaw4budbmcgaOG7o3AgMSBsw6AgxJHDoW5nIGNow7ogw70gaMahbi4NCg0KVHJvbmcgbmjGsMyDbmcgdHLGsOG7nW5nIGjhu6NwIG5oxrAgduG6rXksIHThu4kgc+G7kSBnaeG7r2EgY8OhYyB04bu3IGzhu4cgY8WpbmcgbMOgIG3hu5l0IHRoxrDhu5tjIMSRbyBo4buvdSDDrWNoLg0KDQoqKsSQ4buLbmggbmdoxKlhIDI6KiogR2nhuqMgc+G7rSAkXHBpXzEkIHbDoCAkXHBpXzIkIGzDoCB4w6FjIHN1w6J0ICJ0aMOgbmggY8O0bmciIGPhu6dhIGhhaSBuaMOzbS4gS2hpIMSRw7MgcuG7p2kgcm8gdMawxqFuZyDEkeG7kWkgKHJlbGF0aXZlIHJpc2spIGdp4buvYSBoYWkgbmjDs20gbsOgeSBsw6AgdOG7tyBs4buHIHjDoWMgc3XhuqV0ICJ0aMOgbmggY8O0bmciIGPhu6dhIGhhaSBuaMOzbTogJCQNClxmcmFje1xwaV8xfXtccGlfMn0NCiQkKDIuMi4zKQ0KDQpDaOG6s25nIGjhuqFuOiBUcm9uZyB2w60gZOG7pSB0csOqbiB04bu3IGzhu4cgduG7gSBwaOG6o24g4bupbmcgYuG6pXQgbOG7o2kgY+G7p2EgbG/huqFpIHRodeG7kWMgdGjhu6kgbmjhuqV0IHbDoCBsb+G6oWkgdGh14buRYyB0aOG7qSBoYWksIHRyb25nIHRyxrDhu51uZyBo4bujcCAxIGPDsyBy4bunaSBybyB0xrDGoW5nIMSR4buRaSBsw6AgMTAsMDAgdsOgIHRyb25nIHRyxrDhu51uZyBo4bujcCAyIGPDsyBy4bunaSBybyB0xrDGoW5nIMSR4buRaSBsw6AgMSwwMiAuDQoNCioqTmjhuq1uIHjDqXQ6KioNCg0KXC0gUuG7p2kgcm8gdMawxqFuZyDEkeG7kWkgbMOgIHPhu5Ega2jDtG5nIMOibS4NCg0KXC0gUuG7p2kgcm8gdMawxqFuZyDEkeG7kWkgYuG6sW5nIDEga2hpIHbDoCBjaOG7iSBraGkgYmnhur9uIMSRw6FwIOG7qW5nIMSR4buZYyBs4bqtcCB0aGVvIG5ow7NtLg0KDQpcLSBDw6FjIHThu7cgbOG7hyB04buVbmcgdGjhu4MgJFxwaV8xJCB2w6AgJFxwaV8yJCBuw7NpIGNodW5nIGzDoCBjaMawYSBiaeG6v3QsIG3DoCBxdWEgxJFp4buBdSB0cmEgdGEgbmjhuq1uIMSRxrDhu6NjIGPDoWMgdOG7tyBs4buHIG3huqt1IHTGsMahbmcg4bupbmcgbMOgICRmXzEkIHbDoCAkZl8yJC4gS2hpIMSRw7MgdGEgY8OzIHThu7cgbOG7hzogJFxmcmFje2ZfMX17Zl8yfSQgZ+G7jWkgbMOgIHLhu6dpIHJvIHTGsMahbmcgxJHhu5FpIG3huqt1LCB0YSBkw7luZyBuw7MgxJHhu4MgeOG6pXAgeOG7iSBjaG8gcuG7p2kgcm8gdMawxqFuZyDEkeG7kWkgJFxmcmFje1xwaV8xfXtccGlfMn0kIHRyw6puIHThu5VuZyB0aOG7gy4NCg0KTmfGsOG7nWkgdGEgY2jhu4kgcmEgcuG6sW5nIHBow6JuIHBo4buRaSBt4bqrdSBy4bunaSBybyB0xrDGoW5nIMSR4buRaSBjw7MgdGjhu4MgYuG7iyBs4buHY2ggcuG6pXQgbOG7m24gdHLhu6sga2hpIGvDrWNoIHRoxrDhu5tjIG3huqt1IMSR4bunIGzhu5tuLiBEbyDEkcOzIGPDtG5nIHRo4bupYyDGsOG7m2MgbMaw4bujbmcga2hv4bqjbmcgdGluIGPhuq15IGPhu6dhIG7DsyBraMOhIHBo4bupYyB04bqhcCwgbmfGsOG7nWkgdGEgcGjhuqNpIHRow7RuZyBxdWEgxrDhu5tjIGzGsOG7o25nIGtob+G6o25nIHRpbiBj4bqteSBjaG8gJFxsb2cgXGxlZnQoXGZyYWN7XHBpXzF9e1xwaV8yfVxyaWdodCkkIGzDoDoNCg0KJCQNClxsb2cgXGxlZnQoXGZyYWN7Zl8xfXtmXzJ9XHJpZ2h0KSBccG0gdShcYWxwaGEgLyAyKSBcc3FydHtcZnJhY3sxLWZfMX17bl97MSt9IFxjZG90IGZfMX0rXGZyYWN7MS1mXzJ9e25fezIrfSBcY2RvdCBmXzJ9fQ0KJCQoMi4yLjQpDQoNCipWw60gZHUgNzoqIFRyb25nIELhuqNuZyAyLjMsIGPDsyAkZl8xPTAsMDE3MSwgZl8yPTAsMDA5NCwgbl97MSt9PSQgJDExMDM0LCBuX3syK309MTEwMzckDQoNClwtIFLhu6dpIHJvIHTGsMahbmcgxJHhu5FpIG3huqt1IGzDoCAkXG1hdGhybXtmfV8xIC8gXG1hdGhybXtmfV8yPTAsMDE3MSAvIDAsMDA5ND0xLDgyJCwgY2hvIHRo4bqleTogdOG7tyBs4buHIGLhu4sgYuG7h25oIG5o4buTaSBtw6F1IGPGoSB0aW0gxJHhu5FpIHbhu5tpIG5nxrDhu51pIHPhu60gZOG7pW5nIHRodeG7kWMgYXNwb3JpbiBjYW8gaMahbiAkODIgXCUkIHNvIHbhu5tpIG5ow7NtIHPhu60gZOG7pW5nIGdp4bqjIGTGsOG7o2MuDQoNClwtIFPhu60gZOG7pW5nIGPDtG5nIHRo4bupYyAoMi4yLjQpLCB0YSBjw7Mga2hv4bqjbmcgdGluIGPhuq15ICQ5NSBcJSQgY2hvICRcbG9nIFxsZWZ0KFxmcmFje1xwaV8xfXtccGlfMn1ccmlnaHQpJCBsw6A6ICQoMCwzNjA3IDsgMCw4MzI5KSQuIFN1eSByYSBraG/huqNuZyB0aW4gY+G6rXkgY2hvIHLhu6dpIHJvIHTGsMahbmcgxJHhu5FpICRcZnJhY3tccGlfMX17XHBpXzJ9JCBsw6A6ICQoMSw0MyA7IDIsMzApJA0KDQpOaMawIHbhuq15LCBjaMO6bmcgdGEgY8OzIHRo4buDIHRpbiB0deG7n25nICQ5NSBcJSQgcuG6sW5nLCBzYXUgbsSDbSBuxINtLCBu4bq/dSBz4butIGThu6VuZyB0aHXhu5FjIGFzcGFyaW4gdGjDrCBuZ3V5IGPGoSBi4buLIG5o4buTaSBtw6F1IGPGoSB0aW0gY2FvIGfhuqVwIHThu6sgMSw0MyDEkeG6v24gMiwzMCBs4bqnbiBzbyB24bubaSBuZ8aw4budaSBraMO0bmcgc+G7rSBk4bulbmcgdGh14buRYyBhc3BhcmluLg0KDQpRdWEgxJHDsyB0aOG6pXkgcuG6sW5nOiBraG/huqNuZyB0aW4gY+G6rXkgY+G7p2EgcuG7p2kgcm8gdMawxqFuZyDEkeG7kWkgY2hvIHRo4bqleSBuZ3V5IGPGoSBi4buLIG5o4buTaSBtw6F1IGPGoSB0aW0gY2FvIGjGoW4gw610IG5o4bqldCAkNDMgXCUkIMSR4buRaSB24bubaSBuaMOzbSBkw7luZyBnaeG6oyBkxrDhu6NjLiBLaG/huqNuZyB0aW4gY+G6rXkgY2hvIHPhu7Ega2jDoWMgYmnhu4d0IHbhu4EgdOG7tyBs4buHIGzDoCAkKDAsMDA1IDsgMCwwMTEpJCBjw7MgduG6uyBuaMawIGhhaSBuaMOzbSBraMOhYyBy4bqldCBuaOG7jywgbmjGsG5nIG5ndXkgY8ahIHTGsMahbmcgxJHhu5FpIGNobyB0aOG6pXkgc+G7sSBraMOhYyBiaeG7h3QgY8OzIHRo4buDIGPDsyDDvSBuZ2jEqWEgcXVhbiB0cuG7jW5nIMSR4buRaSB24bubaSBz4bupYyBraG/hursgY+G7mW5nIMSR4buTbmcuIE7hur91IGNow7puZyB0YSBjaOG7iSBz4butIGThu6VuZyBz4buxIGNow6puaCBs4buHY2ggdOG7tyBs4buHIMSR4buDIHNvIHPDoW5oIGhhaSBuaMOzbSBjw7MgdGjhu4MgcGjhuqduIG7DoG8gZ8OieSBuaOG6p20gbOG6q24ga2hpIHThu7cgbOG7hyBj4bunYSBj4bqjIGhhaSBuaMOzbSDEkeG7gXUgcuG6pXQgZ+G6p24ga2jDtG5nLg0KDQpOaGnhu4F1IGtoaSwgdGhheSB2w6wga2jhuqNvIHPDoXQgeMOhYyBzdeG6pXQgInRow6BuaCBjw7RuZyIsIG5nxrDhu51pIHRhIGPhuqduIGto4bqjbyBzw6F0IHRow7RuZyB0aW4gxJHhu4MgdMOtbmggdG/DoW4gY2hvIHThu7cgbOG7hyBj4bunYSB4w6FjIHN14bqldCAidGjhuqV0IGLhuqFpIiAkXGxlZnQoMS1ccGlfMVxyaWdodCkgL1xsZWZ0KDEtXHBpXzJccmlnaHQpJC4NCg0KIyMjICoqVOG7tyBs4buHIGNow6puaCAoT2RkcyByYXRpbykqKg0KDQpUcm9uZyBwaOG6p24gbsOgeSBjaMO6bmcgdGEgdHLDrG5oIGLDoHkgbeG7mXQgY8O0bmcgY+G7pSBraMOhYyDEkeG7gyDEkcOhbmggZ2nDoSBjaG8gYuG6o25nIG5n4bqrdSBuaGnDqm4gJDIgXHRpbWVzIDIkLCDEkcaw4bujYyBn4buNaSBsw6AgdOG7tyBs4buHIGNow6puaCAob2RkcyByYXRpbykuIMSQw6J5IGzDoCBt4buZdCB0aGFtIHPhu5EgY8ahIGLhuqNuIMSRxrDhu6NjIHPhu60gZOG7pW5nIHRyb25nIGPDoWMgbcO0IGjDrG5oIHPhur0gxJHGsOG7o2MgdHLDrG5oIGLDoHkg4bufIGPDoWMgcGjhuqduIHNhdSBj4bunYSB0w6BpIGxp4buHdSBiw6BpIGdp4bqjbmcgbsOgeS4NCg0KKirEkOG7i25oIG5naMSpYSAzOioqIFRyb25nIGjDoG5nICRpJCBj4bunYSBi4bqjbmcgcGjDom4gcGjhu5FpIMSR4buTbmcgdGjhu51pICQyIFx0aW1lcyAyJCBj4bunYSBjw6FjIGJp4bq/biBxdWFuIHPDoXQgJFgkIHbDoCAkWSQsIHThu7cgbOG7hyBjxrDhu6NjIChvZGRzKSBj4bunYSAidGjDoG5oIGPDtG5nIiDEkcaw4bujYyDEkeG7i25oIG5naMSpYSBsw6AgdOG7tyBz4buRIGdp4buvYSB4w6FjIHN14bqldCAidGjDoG5oIGPDtG5nIiB2w6AgeMOhYyBzdeG6pXQgImtow7RuZyB0aMOgbmggY8O0bmciLCBrw70gaGnhu4d1IG9kZHMgJHNfaSQgOiAkJA0KbyBkIGQgc19pPVxmcmFje1xwaV9pfXsxLVxwaV9pfQ0KJCQoMi4zLjEpDQoNCioqTmjhuq1uIHjDqXQ6KioNCg0KXC0gb2RkcyAkc19pJCBsw6AgdOG7tyBs4buHIGdp4buvYSBraOG6oyBuxINuZyAidGjDoG5oIGPDtG5nIiB2w6Aga2jhuqMgbsSDbmcgImtow7RuZyB0aMOgbmggY8O0bmciIHRyb25nIGPDuW5nIG3hu5l0IMSRaeG7gXUga2nhu4duIHTDoWMgxJHhu5luZyBj4bunYSBiaeG6v24gZ2nhuqNpIHRow61jaC4NCg0KXC0gT2RkcyBsdcO0biBraMO0bmcgw6JtOiBPZGRzID0gMSBjw7MgbmdoxKlhIGzDoCBraOG6oyBuxINuZyAidGjDoG5oIGPDtG5nIiB2w6AgImtow7RuZyB0aMOgbmggY8O0bmciIGzDoCBuZ2FuZyBuaGF1OyBPZGRzIFw+IDEgY8OzIG5naMSpYSBsw6Aga2jhuqMgbsSDbmcgInRow6BuaCBjw7RuZyIgY2FvIGjGoW4ga2jhuqMgbsSDbmcgImtow7RuZyB0aMOgbmggY8O0bmciDQoNClwtIEtoaSBiaeG6v3Qgb2RkcywgdGEgY8OzIHRo4buDIHTDrW5oIMSRxrDhu6NjIHjDoWMgc3XhuqV0IHRow6BuaCBjw7RuZywgdsOsIHRhIGPDszoNCg==