Chương 1. Tổng quan về đề tài

1.1. Giới thiệu dữ liệu

NMES1988 là một bộ dữ liệu chéo, trích từ kết quả của cuộc khảo sát Dịch vụ chăm sóc sức khỏe tại Hoa Kỳ (NMES). NMES là cuộc khảo sát những người đã không được khám, chữa bệnh tại các cơ sở chăm sóc sức khỏe trong một thời gian dài trong suốt năm 1987 và 1988.

Bộ dữ liệu NMES1988 chỉ lấy kết quả khảo sát của những người trên 66 tuổi và vào năm 1988.

Tập dữ liệu NMES1988 bao gồm 4406 quan sát với 19 biến:

visit : số lần khám bệnh với bác sĩ tại bệnh viện.

nvisit : số lần khám bệnh với nhân viên y tế như y tá, điều dưỡng (không phải bác sĩ) tại bệnh viện.

ovisits : số lần khám bệnh với bác sĩ tại phòng khám ngoại trú.

novisits : số lần thăm khám với nhân viên y tế tại phòng khám ngoại trú.

emergency : số lần cấp cứu.

hospital : số lần nhập viện.

health: trạng thái sức khỏe. “poor” là không khỏe, “average” là bình thường, “excellent” là rất khỏe.

chronic : số bệnh mãn tính mà người được khảo sát mắc phải.

adl : người được khảo sát có gặp khó khăn trong việc thực hiện các hoạt động hàng ngày như vệ sinh cá nhân, di chuyển, ăn uống… hay không? “limited” là có gặp khó khăn trong việc thực hiện các hoạt động hàng ngày, “normal” là không có.

region : nói về vùng miền của người thực hiện khảo sát. “northeast” là miền đông bắc, “midwest” là miền trung tây Hoa Kỳ, “west” là miền tây, “orther” là vùng miền khác.

age : độ tuổi (độ tuổi chia 10).

afam : người được khảo sát có phải người gốc Mỹ gôc Phi hay không? “yes” là phải, “no” là không.

gender : giới tính, “male” là nam, “female” là nữ.

married : đã kết hôn chưa? “yes” là kết hôn rồi, “no” là chưa kết hôn.

school : là số năm đi học.

income : thu nhập gia đình (đơn vị : 10ngàn đô).

employed : người được khảo sát có đang đi làm không? “yes” là có, “no” là không.

insurance : người được khảo sát có tự mua bảo hiểm cá nhân hay không? “yes” là có, “no” là không. Bảo hiểm cá nhân hay còn gọi là bảo hiểm tư nhân là cách gọi chung cho các loại sản phẩm bảo hiểm khác nhau như bảo hiểm nhân thọ, bảo hiểm sức khỏe, bảo hiểm y tế,… Bảo hiểm cá nhân bắt buộc người được bảo hiểm cần phải bỏ tiền ra mua bảo hiểm.

medicaid : người được khảo sát có được bảo hiểm bởi tổ chức bảo hiểm Medicaid không? “yes” là có, “no” là không. Medicaid là chương trình bảo hiểm chăm sóc sức khỏe, mà người được bảo hiểm không cần phải chi tiền để mua bảo hiểm. Medicaid bảo hiểm cho công dân Mỹ, thỏa mãn một hoặc nhiều điều kiện về thu nhập, độ tuổi cũng như giá trị tài sản mà người đó sở hữu.

1.2. Giới thiệu đề tài

Trong bài phân tích này, nhóm tác giả sẽ phân tích các yếu tố ảnh hưởng đến việc mua bảo hiểm cá nhân.

Nhóm đã chọn 6 biến để nghiên cứu bao gồm:

Biến định tính : insurance, employed, medicaid
Biến định lượng : income, hospital, chronic

Nghiên cứu này sẽ phân tích các yếu tố như thu nhập, tìm trạng việc làm, số lần nhập viện, bệnh mãn tính, được Medicaid bảo hiểm ảnh hưởng như thế nào đến việc một cá nhân ra quyết định mua bảo hiểm cá nhân.

Biến phụ thuộc : insurance
Biến độc lập : income, hospital, chronic, employed, medicaid

library(tidyverse)
library(ggplot2)
library(AER)
library(DT)
library(dplyr)
library(DescTools)
library(epitools)
library(scales)
library(VGAM)
data("NMES1988")
t <- NMES1988
datatable(t)

Chương 2. Thống kê mô tả 1 biến

2.1. Biến định tính

2.1.1. Đổi tên biểu hiện của các biến định tính

Để tránh nhầm sự nhầm lẫn trong quá trình phân tích, ta đổi tên các biểu hiện “yes”, “no” của 3 biến định tính insurance, employed, medicaid

Đối với biến insurance : “yes” đổi thành “coBHCN”, nghĩa là Người khảo sát có mua bảo hiểm cá nhân. “no” đổi thành “khongBHCN”, nghĩa là Người khảo sát không mua bảo hiểm cá nhân.
Đối với biến employed : “yes” đổi thành “covieclam”, nghĩa là Người khảo sát đang có việc làm. “no” đổi thành “khongvieclam”, nghĩa là Người khảo sát đang không có việc làm.
Đối với biến medicaid : “yes” đổi thành “coBHMedi”, nghĩa là Người khảo sát được bảo hiểm bởi tổ chức Medicaid. “no” đổi thành “khongBHMedi”, nghĩa là Người khảo sát không được bảo hiểm bởi tổ chức Medicaid.

t$insurance <- factor(t$insurance, levels = c('yes', 'no'), labels = c('coBHCN','khongBHCN'))

t$employed <- factor(t$employed, levels = c('yes', 'no'), labels = c('covieclam', 'khongvieclam' ))

t$medicaid <- factor(t$medicaid, levels = c('yes','no'), labels = c('coBHMedi','khongBHMedi' ))

2.1.2. Biến insurance

table(t$insurance)

## 
##    coBHCN khongBHCN 
##      3421       985

prop.table(table(t$insurance))

## 
##    coBHCN khongBHCN 
## 0.7764412 0.2235588

ggplot(t, aes(x = insurance, y= after_stat(count)))+
  geom_bar(fill = 'pink')+
  geom_text(aes(label = scales::percent(after_stat(count/sum(count)))), stat = 'count', color = 'black', vjust = -.1)+
  labs(title=  'Số lượng người mua bảo hiểm cá nhân')

Có thể thấy có sự chênh lệch lớn giữa số lượng người có bảo hiểm cá nhân và số lượng người không có. Trong tổng số người tham gia khảo sát có 3421 người có mua bảo hiểm cá nhân (chiếm 77.64%), còn lại là 985 người không sở hữu bảo hiểm cá nhân (chiếm 22.35%).

2.1.3. Biến employed

addmargins(table(t$employed))

## 
##    covieclam khongvieclam          Sum 
##          455         3951         4406

ggplot(t, aes(x = employed, y= after_stat(count)))+
  geom_bar(fill = 'black')+
  geom_text(aes(label = scales::percent(after_stat(count/sum(count)))), stat = 'count', color = 'black', vjust = -.5)+
  labs(title=  'Tình trạng việc làm')

Dựa vào biểu đồ, ta thấy được trong số những người tham gia khảo sát thì phần lớn là người không có việc. Biểu hiện không có việc làm có sự chênh lệch rõ ràng so với biểu hiện có việc làm. Cụ thể, trong tổng số 4406 người, thì có 455 người đang có việc (chỉ chiếm 10%) và 3951 người còn lại là không có việc làm (chiếm tận 90%).

2.1.4. Biến medicaid

addmargins(table(t$medicaid))

## 
##    coBHMedi khongBHMedi         Sum 
##         402        4004        4406

ggplot(t, aes(x = medicaid, y= after_stat(count)))+
  geom_bar(fill = 'lightblue')+
  geom_text(aes(label = scales::percent(after_stat(count/sum(count)))), stat = 'count', color = 'black', vjust = -.5)+
  labs(title=  'Có được tổ chức Medicaid bảo hiểm không?')

Biến medicaid cũng có sự chênh lệch lớn giữa hai biểu hiện. Tại thời điểm khảo sát, chỉ có 402 người là được tổ chức Medicaid bảo hiểm (chiếm phần nhỏ 9%), 4004 người còn lại không được tổ chức Medicaid tài trợ bảo hiểm (chiếm 91%).

2.2. Biến định lượng

2.2.1. Biến income

Chuyển đổi từ biến định lượng sang biến định tính bằng cách chia khoảng trong dữ liệu.

Đầu tiên, ta nhìn tổng quan các biến income.

summary(t$income)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
## -1.0125  0.9122  1.6982  2.5271  3.1728 54.8351

Trung bình thu nhập của nhúng người được khảo sát nằm ở 2.5271.
Mức thu nhập thấp nhất là -1.0125, cao nhất là 54.8351.
Trong đó ở mức thu nhập 0.9122 thì có 25% dữ liệu thấp hơn lượng thu nhập này Ở mức 1.6982 thì có 50% dữ liệu có thu nhập thấp hơn và 50% dữ liệu cao hơn mức này. Còn ở mức 3.1728 thì có 75% dữ liệu có lượng thu nhập thấp hơn mức này.

Tiếp theo, ta tạo một biến mới tên là income1. Biến này sẽ chứa dữ liệu về thu nhập của người tham gia khảo sát và cho biết ở mức thu nhập nào là thấp, mức nào là cao.

Khoảng từ -1.1 đến 2.5 sẽ có tên là ‘thap’ , được xem mức thu nhập thấp.
Khoảng từ 2.5 đến 55 tên là ‘cao’, được xem như mức thu nhập cao.

t$incomed <- cut(t$income, breaks = c(-1.1 ,2.5,55), labels = c("thap","cao"))

table(t$incomed)

## 
## thap  cao 
## 2932 1474

ggplot(t, aes(x = incomed, y= after_stat(count)))+
  geom_bar(fill = 'pink')+
  geom_text(aes(label = scales::percent(after_stat(count/sum(count)))), stat = 'count', color = 'black', vjust = -.5)+
  labs(title=  'Thu nhập của những người khảo sát')

Dựa trên biểu đồ, ta có thể thấy những người có thu nhập thấp chiếm số lượng nhiều hơn rất nhiều so với người có thu nhập cao. Cụ thể, có 2943 người trong tổng số 4406 người có thu thập thấp (chiếm 67%), còn lại là 1463 người có thu nhập cao (chiếm 33%).

2.2.2. Biến hospital

Tương tự cho các biến hospital

summary (t$hospital)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   0.000   0.000   0.000   0.296   0.000   8.000

Trung bình những người khảo sát không nhập viện lần nào, (giá trị trung bình là 0.296)
Số lần thấp nhất là 0, cao nhất là 8 lần.
Có thể thấy ở tứ phân vị thứ nhất, thứ hai và thứ ba, thì số lần nhập viện đều là 0 lần. Điều này cho thấy số lần nhập viện của những người được khảo sát rất ít thậm chí là không có.

Tiếp theo, ta tạo một biến mới tên là hospital1. Biến này sẽ chứa dữ liệu về số lần nhập viện của người tham gia khảo sát và cho biết ở mấy lần nào là thấp, mấy lần là cao.

Khoảng từ -0.5 đến 2.5 sẽ có tên là ‘it’ , được xem là không nhập viện.
Khoảng từ 2.5 đến 9 tên là ‘nhieu’, được xem như nhập viện nhiều.

t$hospital1 <- cut(t$hospital, breaks = c(-0.5,2.5,9), labels = c("it","nhieu"))

addmargins(table(t$hospital1))

## 
##    it nhieu   Sum 
##  4316    90  4406

ggplot(t, aes(x = hospital1, y= after_stat(count)))+
  geom_bar(fill = 'blue')+
  geom_text(aes(label = scales::percent(after_stat(count/sum(count)))), stat = 'count', color = 'black', vjust = -.5)+
  labs(title=  'Số lần nhập viện của những người khảo sát')

Dựa trên biểu đồ, ta có thể thấy có khoảng chênh lệch rất lớn giữa giữa người nhập viện ít và người nhập viện nhiều. Cụ thể, số người không nhập viện chiếm tận 3541 người (80%), trong khi nhóm người có nhập viện chiếm số lượng rất ít, chỉ có 865 người(20%).

2.2.3. Biến chronic

Đối với biến chronic:

summary (t$chronic)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   0.000   1.000   1.000   1.542   2.000   8.000

Trung bình số bệnh mãn tính mà những người được khảo sát mắc phải là 1.542
Số bệnh thấp nhất là 0, cao nhất là 8.
Trong đó ở mức bị 1 bệnh mãn tính thì có 25% dữ liệu thấp hơn mức này. Ở mức 1 bệnh thì có 50% dữ liệu có số lượng bệnh thấp hơn và 50% dữ liệu cao hơn mức này. Còn ở mức bị 2 bệnh mãn tính thì có 75% dữ liệu thấp hơn mức này.

Tiếp theo, ta tạo một biến mới tên là chronic1. Biến này sẽ chứa dữ liệu về số lượng bệnh mãn tính mà người tham gia khảo sát mắc phải.

Khoảng từ -0.5 đến 0.5 sẽ có tên là ‘khongbenh’ , được hiểu là người khảo sát không mắc bệnh mãn tính.
Khoảng từ 0.5 đến 9 tên là ‘cobenh’, được hiểu là người khảo sát mắc 1 hoặc nhiều hơn 1 bệnh mãn tính.

t$chronic1 <- cut(t$chronic, breaks = c(-0.5,0.5,9), labels = c("khongbenh","cobenh"))

addmargins(table(t$chronic1))

## 
## khongbenh    cobenh       Sum 
##      1025      3381      4406

ggplot(t, aes(x = chronic1, y= after_stat(count)))+
  geom_bar(fill = 'pink')+
  geom_text(aes(label = scales::percent(after_stat(count/sum(count)))), stat = 'count', color = 'black', vjust = -.5)+
  labs(title=  'Số người mắc bệnh mãn tính')

Có thể thấy, có sự chênh lệch lớn giữa những người không có bệnh mãn tính và nhóm người có bệnh mãn tính.

Số người không mắc bệnh mãn tính là 1025 người (chiếm 23%), còn lại là những người mắc nhiều bệnh mãn tính, có 3381 người (chiếm 77%).

Chương 3. Thống kê mô tả các cặp biến

3.1. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và thu nhập (income)

3.1.1. Bảng tần suất

addmargins(table(t$incomed,t$insurance))

##       
##        coBHCN khongBHCN  Sum
##   thap   2137       795 2932
##   cao    1284       190 1474
##   Sum    3421       985 4406

addmargins(prop.table(table(t$incomed,t$insurance)))

##       
##            coBHCN  khongBHCN        Sum
##   thap 0.48502043 0.18043577 0.66545620
##   cao  0.29142079 0.04312301 0.33454380
##   Sum  0.77644122 0.22355878 1.00000000

t %>% ggplot(aes(x = incomed, y = after_stat(count), fill = insurance)) +
  geom_bar(position = "dodge", color = "black") +
  geom_text(aes(label = after_stat(count)), stat = 'count', color = 'black', vjust = -0.25, position = position_dodge(width = 0.9)) +
  theme_classic() + 
  labs(x = "Thu nhập",
       y = "Số lượng",
       title = "Thu nhập và bảo hiểm") +
  scale_fill_manual(values = c("khongBHCN" = "dodgerblue", "coBHCN" = "blue"))

Bảng tần suất và biểu đồ về việc sở hữu bảo hiểm và thu nhập cho thấy:

Nhóm có thu nhập cao:

Có 190 người không có bảo hiểm cá nhân chiếm tỷ lệ là 4.31% tổng số người trong bộ dữ liệu

Có 1284 người có bảo hiểm cá nhân chiếm 29.14% tổng số người trong bộ dữ liệu

Nhóm có thu nhập thấp:

có 795 người không có bảo hiểm cá nhân chiếm 18.04% tổng số người trong bộ dữ liệu

có 2137 người có bảo hiểm cá nhân chiếm 48.50% tổng số người trong dữ liệu

Có thể thấy những người có thu nhập thấp thường có xu hướng mua bảo hiểm y tế nhiều hơn những người có thu nhập cao

3.1.2. Relative Risk

th8 <- table(t$incomed,t$insurance)
addmargins(th8)

##       
##        coBHCN khongBHCN  Sum
##   thap   2137       795 2932
##   cao    1284       190 1474
##   Sum    3421       985 4406

RelRisk(th8)

## [1] 0.8367063

Nhận xét:

Tỷ lệ sở hữu bảo hiểm cá nhân ở nhóm thu nhập thấp bằng 0.8367063 so với tỷ lệ sở hữu bảo hiểm của người thu nhập cao.

3.1.3. Tỷ lệ chênh - Oddratio

OddsRatio(th8)

## [1] 0.3977645

Kết quả Odd Ratio là 0.3977645 cho thấy rằng tỉ lệ người thu nhập thấp có khả năng sở hữu bảo hiểm cá nhân thấp hơn người thu nhập cao.

Vì sự chênh lệch giữa nhóm có BHCN và không có BHCN của “người thu nhập thấp” thấp hơn sự chênh lệch giữa nhóm có BHCN và không có BHCN của người thu nhập cao.

3.2. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và biến employed (có việc làm)

3.2.1. Bảng tần suất

addmargins(table(t$employed,t$insurance))

##               
##                coBHCN khongBHCN  Sum
##   covieclam       391        64  455
##   khongvieclam   3030       921 3951
##   Sum            3421       985 4406

round((table(t$employed,t$insurance)/sum(table(t$employed,t$insurance)))*100, 2)

##               
##                coBHCN khongBHCN
##   covieclam      8.87      1.45
##   khongvieclam  68.77     20.90

t %>% ggplot(aes(x = employed, y = after_stat(count), fill = insurance)) +
  geom_bar(position = "dodge", color = "black") +
  geom_text(aes(label = after_stat(count)), stat = 'count', color = 'black', vjust = -0.25, position = position_dodge(width = 0.9)) +
  theme_classic() +
  labs(x = "Việc làm",
       y = "Số lượng",
       title = "Việc làm và bảo hiểm") +
  scale_fill_manual(values = c("khongBHCN" = "darkorange", "coBHCN" = "red"))

Nhận xét:

Bảng tần suất và biểu đồ về việc sở hữu bảo hiểm và có việc làm cho thấy:

Nhóm không có việc làm là:

Có 921 người không có bảo hiểm cá nhân chiếm tỷ lệ là 20.9% tổng số người trong bộ dữ liệu

Có 3030 người có bảo hiểm cá nhân chiếm tỷ lệ là 68.77% tổng số người trong bộ dữ liệu.

Nhóm có việc làm

Có 64 người không có bảo hiểm cá nhân chiếm tỷ lệ là 1.45% tổng số người trong bộ dữ liệu

Có 331 người có bảo hiểm cá nhân chiếm 8.87% tổng số người trong bộ dữ liệu

Như vậy:

Tổng cộng có 985 người không có bảo hiểm y tế cá nhân, chiếm 22.35% tổng số người trong tập dữ liệu.

Tổng cộng có 3421 người có bảo hiểm y tế cá nhân, chiếm 77.65% tổng số người trong tập dữ liệu.

Tỷ lệ người không có việc làm có bảo hiểm cá nhân chiếm số lượng lớn hơn những người có việc làm. Cho thấy những người không có việc làm có xu hướng mua bảo hiểm nhiều hơn

3.2.2. Relative Risk

t3 <- table(t$employed,t$insurance)
t3

##               
##                coBHCN khongBHCN
##   covieclam       391        64
##   khongvieclam   3030       921

RelRisk(t3)

## [1] 1.120546

Nhận xét:

Tỷ lệ Relative Risk được tính là 1.120546 cho biết rằng khả năng của những người có việc làm và sở hữu bảo hiểm cá nhân cao hơn khoảng 12% so với những người có việc làm.

3.2.3. Tỷ lệ chênh - Oddratio

OddsRatio(t3)

## [1] 1.857008

Nhận xét:

Kết quả Odd Ratio của người có việc làm và không có việc làm khoảng 1.857008 điều này có nghĩa là odd của nhóm người sở hữu bảo hiểm, có việc làm cao hơn 1.857008 lần so với odd người không có việc làm.

Có sự khác biệt trong kết quả giữa Relative risk, odd ratios và bảng tần suất. Điều này có thể giải thích rằng do tỷ lệ sở hữu BHCN trong nhóm người có việc làm cao hơn tỷ lệ sở hữu BHCN của nhóm người không việc làm.

3.3. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và biến medicaid

3.3.1. Bảng tần suất

addmargins(table(t$medicaid,t$insurance))

##              
##               coBHCN khongBHCN  Sum
##   coBHMedi        61       341  402
##   khongBHMedi   3360       644 4004
##   Sum           3421       985 4406

round((table(t$medicaid,t$insurance)/sum(table(t$medicaid,t$insurance)))*100, 2)

##              
##               coBHCN khongBHCN
##   coBHMedi      1.38      7.74
##   khongBHMedi  76.26     14.62

t %>% ggplot(aes(x = medicaid, y = after_stat(count), fill = insurance)) +
  geom_bar(position = "dodge", color = "black") +
  geom_text(aes(label = after_stat(count)), stat = 'count', color = 'black', vjust = -0.25, position = position_dodge(width = 0.9)) +
  theme_classic() +
  labs(x = "Bảo hiểm medicaid",
       y = "Số lượng",
       title = "Bảo hiểm medicaid và bảo hiểm cá nhân") +
  scale_fill_manual(values = c("khongBHCN" = "tomato", "coBHCN" = "orange"))

Bảng tần suất và biểu đồ về việc sở hữu bảo hiểm và bảo hiểm medicaid cho thấy:

Nhóm không có bảo hiểm medicaid là:

Có 644 người không có bảo hiểm cá nhân chiếm tỷ lệ là 14.62% tổng số người trong bộ dữ liệu

Có 3360 người có bảo hiểm cá nhân chiếm tỷ lệ là 76.26% tổng số người trong bộ dữ liệu.

Nhóm có bảo hiểm medicaid

Có 341 người không có bảo hiểm cá nhân chiếm tỷ lệ là 7.74% tổng số người trong bộ dữ liệu

Có 61 người có bảo hiểm cá nhân chiếm 1.38% tổng số người trong bộ dữ liệu

Như vậy:

Tổng cộng có 985 người không có bảo hiểm y tế cá nhân, chiếm 22.35% tổng số người trong tập dữ liệu.

Tổng cộng có 3421 người có bảo hiểm y tế cá nhân, chiếm 77.65% tổng số người trong tập dữ liệu.

Điều này phản ánh rằng những người có bảo hiểm cá nhân ít khi có bảo hiểm Medicaid, trong khi ngược lại, những người không có bảo hiểm cá nhân có xu hướng có bảo hiểm Medicaid nhiều hơn.

3.3.2. Relative Risk

t4 <- table(t$medicaid,t$insurance)
t4

##              
##               coBHCN khongBHCN
##   coBHMedi        61       341
##   khongBHMedi   3360       644

RelRisk(t4)

## [1] 0.180825

Nhận xét:

Tỷ lệ Relative Risk được tính là 0.180825, cho biết rằng khả năng sở hữu bảo hiểm cá nhân của người có bảo hiểm medicaid bằng 0.180825 lần so với người không có bảo hiểm Medicaid.

3.3.3. Tỷ lệ chênh - Oddratio

OddsRatio(t4)

## [1] 0.03428641

Kết quả Odd Ratio của những người không có bảo hiểm medicaid so với người không có bảo hiểm medicaid khoảng 0.03428641 . Tỷ lệ sở hữu BHCN của người được tài trợ thấp hơn tỷ lệ sở hữu BHCN của người không được tài trợ.

3.4. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và biến hospital (số lần nhập viện)

Mục đích: Dùng biến hospital và insurance để nghiên cứu sự tương tác hoặc mối quan hệ giữa việc sử dụng dịch vụ y tế (đo bằng số lần nhập viện) và tình trạng bảo hiểm y tế của người dân vào năm 1988.

Phân tích: Có thể phân tích xem những người có bảo hiểm có xu hướng sử dụng dịch vụ y tế (nhập viện) như thế nào so với những người không có bảo hiểm.

3.4.1. Bảng tần suất

# Bảng tần suất của hospital và insurance:
addmargins(table(t$hospital1,t$insurance))

##        
##         coBHCN khongBHCN  Sum
##   it      3352       964 4316
##   nhieu     69        21   90
##   Sum     3421       985 4406

# Tỉ lệ phần trăm của bảng tần suất:
addmargins(round((table(t$hospital1,t$insurance)/sum(table(t$hospital,t$insurance)))*100, 2))

##        
##         coBHCN khongBHCN    Sum
##   it     76.08     21.88  97.96
##   nhieu   1.57      0.48   2.05
##   Sum    77.65     22.36 100.01

# Trực quan hóa dữ liệu 
t %>% ggplot(aes(x = hospital1, y = after_stat(count), fill = insurance)) +
  geom_bar(position = "dodge", color = "black") +
  geom_text(aes(label = after_stat(count)), stat = 'count', color = 'black', vjust = -0.25, position = position_dodge(width = 0.9)) +
  theme_classic() +
  labs(x = "tần suất nhập viện",
       y = "Số lượng",
       title = "tần suất nhập viện với bảo hiểm Medicaid") +
  scale_fill_manual(values = c("khongBHCN" = "lightblue", "coBHCN" = "lightgreen"))

Nhận xét:

Nhận xét về sự phân tích mối quan hệ giữa biến sở hữu bảo hiểm cá nhân(insurance) và số lần nhập viện (hospital) Phân bố số lần nhập viện theo tình trạng bảo hiểm:

Nhóm ít lần nhập viện (“it”):

Có 964 người không có bảo hiểm y tế tư nhân, chiếm 21.88% tổng số người trong tập dữ liệu.
Có 3352 người có bảo hiểm y tế tư nhân, chiếm 76.08% tổng số người trong tập dữ liệu.

Nhóm nhiều lần nhập viện (“nhieu”):

Có 21 người không có bảo hiểm y tế tư nhân, chiếm 0.48% tổng số người trong tập dữ liệu.
Có 69 người có bảo hiểm y tế tư nhân, chiếm 1.57% tổng số người trong tập dữ liệu.

Tổng cộng có 985 người không có bảo hiểm y tế tư nhân, chiếm 22.36% tổng số người trong tập dữ liệu.

Tổng cộng có 3421 người có bảo hiểm y tế tư nhân, chiếm 77.65% tổng số người trong tập dữ liệu.

\(\to\) Kết luận:

Bảo hiểm y tế tư nhân và số lần nhập viện:

Tỷ lệ người có bảo hiểm y tế tư nhân nhập viện nhiều lần (nhóm “nhieu”) là 69 người, trong khi tỷ lệ người không có bảo hiểm cá nhânnhập viện nhiều lần là 21 người.

Điều này cho thấy rằng những người có bảo hiểm y tế tư nhân có xu hướng sử dụng dịch vụ y tế nhiều hơn so với những người không có bảo hiểm.

Sự chênh lệch rõ ràng:

Số lượng người có bảo hiểm y tế tư nhân nhiều hơn hẳn so với số người không có bảo hiểm, chiếm gần 78% tổng số mẫu.

Chênh lệch này có thể phản ánh một phần sự phổ biến của bảo hiểm y tế tư nhân vào năm 1988, hoặc có thể liên quan đến các yếu tố kinh tế-xã hội khác chưa được phân tích ở đây.

3.4.2. Relative Risk

# Tính chỉ số liên quan đến rủi ro (Relative Risk và Risk Ratio)

th2 <- table(t$hospital1,t$insurance)
th2

##        
##         coBHCN khongBHCN
##   it      3352       964
##   nhieu     69        21

RelRisk(th2)

## [1] 1.013015

Tỷ lệ Relative Risk được tính là 1.013015, cho biết rằng khả năng sở hữu bảo hiểm cá nhân của người ít nhập viện cao hơn 1.013015 lần so với người nhập viện nhiều.

3.4.3. Tỷ lệ chênh - Oddratio

OddsRatio(th2)

## [1] 1.058272

Kết quả Odd Ratio 1.058272 (OR > 1 ) cho thấy rằng tỷ lệ sở hữu BHCN của người ít nhập viện cao hơn người nhập viện nhiều.

3.5. Phân tích biến insurance (sở hữu bảo hiểm tư nhân) và biến chronic (số bệnh mãn tính mà người được khảo sát mắc phải)

3.5.1. Bảng tần suất

addmargins(table(t$chronic1,t$insurance ))

##            
##             coBHCN khongBHCN  Sum
##   khongbenh    800       225 1025
##   cobenh      2621       760 3381
##   Sum         3421       985 4406

addmargins(round((table(t$chronic1,t$insurance)/sum(table(t$chronic1,t$insurance)))*100, 2))

##            
##             coBHCN khongBHCN    Sum
##   khongbenh  18.16      5.11  23.27
##   cobenh     59.49     17.25  76.74
##   Sum        77.65     22.36 100.01

# Trực quan hóa dữ liệu

t %>% ggplot(aes(x = chronic1, y = after_stat(count), fill = insurance)) +
  geom_bar(position = "dodge", color = "black") +
  geom_text(aes(label = after_stat(count)), stat = 'count', color = 'black', vjust = -0.25, position = position_dodge(width = 0.9)) +
  theme_classic() +
  labs(x = "Số bệnh mãn tính",
       y = "Số lượng",
       title = "Số bệnh mãn tính với Bảo hiểm cá nhân") +
  scale_fill_manual(values = c("khongBHCN"= "lightblue", "coBHCN" = "lightgreen"))

Nhận xét:

Nhóm không mắc bệnh mãn tính:

Có bảo hiểm: 800 người (18.16% trong tổng số người)

Không có bảo hiểm: 225 người (5.11% trong tổng số người)

Nhóm mắc bệnh mãn tính:

Có bảo hiểm: 2621 người ( 59.49% trong tổng số người)

Không có bảo hiểm: 760 người (17.25% trong tổng số người)

3.5.2. Relative Risk

th3 <- table(t$chronic1,t$insurance)
th3

##            
##             coBHCN khongBHCN
##   khongbenh    800       225
##   cobenh      2621       760

RelRisk(th3)

## [1] 1.006802

Nhận xét:

Nhận xét:

Tỷ lệ Relative Risk được tính là 1.006802. Gía trị này cho biết khả năng sở hữu bảo hiểm cá nhân của những người không có bệnh mãn tính bằng 1.006802 lần so với những người có bệnh mãn tính.

3.5.3. Tỷ lệ chênh - Oddratio

OddsRatio(th3)

## [1] 1.030989

Nhận xét:

Đối với nhóm không có bảo hiểm cá nhân, odds = 1.030989. Điều này nghĩa là đối với việc sở hữu bảo hiểm cá nhân, thì tỷ lệ chênh của những người có bệnh mãn tính và không bệnh mãn tính là 1.030989.

Chương 4. Thống kê suy diễn

4.1. Kiểm định tính độc lập của 2 biến

4.1.1. Kiểm tra tính độc lập cho biến insurance (sở hữu bảo hiểm cá nhân) và biến incomed (thu nhập)

Giả thuyết \(H_0\) : biến insurance (sở hữu bảo hiểm cá nhân) và biến incomed (thu nhập) độc lập. Giả thuyết \(H_1\): biến insurance (sở hữu bảo hiểm cá nhân) và biến incomed (thu nhập) không độc lập

chisq.test(th8)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  th8
## X-squared = 113.52, df = 1, p-value < 2.2e-16

Nhận xét:

Bác bỏ giả thuyết \(H_0\), chấp nhận giả thuyết \(H_1\), insurance(sở hữu bảo hiểm cá nhân) và thu nhập không độc lập.

4.1.2. Kiểm tra tính độc lập cho hai biến insurance (sở hữu bảo hiểm cá nhân) và biến employed (việc làm)

Giả thuyết \(H_0\) : biến insurance (sở hữu bảo hiểm cá nhân) và biến employed (việc làm) độc lập. Giả thuyết \(H_1\): biến insurance (sở hữu bảo hiểm cá nhân) và biến incomed (thu nhập) không độc lập

chisq.test(t3)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  t3
## X-squared = 19.56, df = 1, p-value = 9.751e-06

Nhận xét Bác bỏ giả thuyết \(H_0\), chấp nhận giả thuyết \(H_1\),biến insurance (sở hữu bảo hiểm cá nhân) và biếnemployed (việc làm) không độc lập

4.1.3. Kiểm tra tính độc lập cho hai biến insurance (sở hữu bảo hiểm cá nhân) và biến medicaid

chisq.test(t4)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  t4
## X-squared = 990.58, df = 1, p-value < 2.2e-16

Nhận xét Bác bỏ giả thuyết \(H_0\), chấp nhận giả thuyết \(H_1\), biến insurance (sở hữu bảo hiểm cá nhân) và biến medicaid không độc lập

4.1.4. Kiểm tra tính độc lập cho hai biến insurance (sở hữu bảo hiểm cá nhân) và biến hospital (số lần nhập viện )

Giả thuyết \(H_0\) : biến insurance (sở hữu bảo hiểm cá nhân) và biến hospital (số lần nhập viện ) độc lập. Giả thuyết \(H_1\): biến insurance (sở hữu bảo hiểm cá nhân) và biến hospital (số lần nhập viện ) không độc lập.

t5 <- table(t$hospital1, t$insurance)
chisq.test(t5)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  t5
## X-squared = 0.0094216, df = 1, p-value = 0.9227

Ta có p_value = 0.9227 > \(\alpha\). Bác bỏ giả thuyết \(H_1\), chấp nhận giả thuyết \(H_0\), biến insurance (sở hữu bảo hiểm cá nhân) và biến hospital (số lần nhập viện ) độc lập với nhau

4.1.5. Kiểm tra tính độc lập cho biến insurance (sở hữu bảo hiểm cá nhân) và biến chronic (số bệnh mãn tính mà người được khảo sát mắc phải)

Giả thuyết \(H_0\) : biến insurance (sở hữu bảo hiểm cá nhân) và biến chronic (số bệnh mãn tính mà người được khảo sát mắc phải) độc lập. Giả thuyết \(H_1\): biến insurance (sở hữu bảo hiểm cá nhân) và biến chronic (số bệnh mãn tính mà người được khảo sát mắc phải) không độc lập.

th5 <- table(t$chronic1,t$insurance)
chisq.test(th5)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  th5
## X-squared = 0.09746, df = 1, p-value = 0.7549

Ta có p_value = 0.7549 > \(\alpha\). Bác bỏ giả thuyết \(H_1\), chấp nhận giả thuyết \(H_0\), biến insurance (sở hữu bảo hiểm cá nhân) và biến chronic (số bệnh mãn tính mà người được khảo sát mắc phải) độc lập với nhau.

4.2. Khoảng ước lượng tỷ lệ

4.2.1. Ước lượng tỉ lệ người có và không có bảo hiểm cá nhân

p <- t[t$insurance == 'coBHCN',] 
prop.test(length(p$insurance), length(t$insurance))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(p$insurance) out of length(t$insurance), null probability 0.5
## X-squared = 1345.7, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.7637858 0.7886106
## sample estimates:
##         p 
## 0.7764412

t10 <- t[t$insurance == 'khongBHCN',] 
prop.test(length(t10$insurance), length(t$insurance))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(t10$insurance) out of length(t$insurance), null probability 0.5
## X-squared = 1345.7, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.2113894 0.2362142
## sample estimates:
##         p 
## 0.2235588

Với độ tin cậy 95%, ta có tỷ lệ người có bảo hiểm tư nhân so với tổng thể nằm trong khoảng từ 76.38% và 78.86%. Hay nói cách khác, tỉ lệ người không có bảo hiểm cá nhân nhân sẽ chiếm khoảng từ 21.11% đến 23.62%

4.2.2. Ước lượng tỉ lệ người có và không có việc làm

e <- t[t$employed == 'covieclam',] 
prop.test(length(e$employed == 'covieclam'), length(t$employed))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(e$employed == "covieclam") out of length(t$employed), null probability 0.5
## X-squared = 2772.4, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.09451665 0.11271981
## sample estimates:
##         p 
## 0.1032683

e1 <- t[t$insurance == 'khongvieclam',] 
prop.test(length(e1$employed == 'khongvieclam'), length(t$employed))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(e1$employed == "khongvieclam") out of length(t$employed), null probability 0.5
## X-squared = 4404, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.000000000 0.001086028
## sample estimates:
## p 
## 0

Với độ tin cậy 95%, ta có tỷ lệ người có việc làm so với tổng thể nằm trong khoảng từ 9.45% đến 11.27%. Hay nói cách khác, tỉ lệ người không có việc làm sẽ chiếm khoảng từ 0% đến 0.1086028%.

4.2.3. Ước lượng tỷ lệ người có và không có bảo hiểm medicaid

e3 <- t[t$medicaid == 'coBHMedi',] 
prop.test(length(e3$medicaid == 'coBHMedi'), length(t$medicaid))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(e3$medicaid == "coBHMedi") out of length(t$medicaid), null probability 0.5
## X-squared = 2943.1, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.08298051 0.10021957
## sample estimates:
##          p 
## 0.09123922

e4 <- t[t$medicaid == 'khongBHMedi',] 
prop.test(length(e4$medicaid == 'khongBHMedi'), length(t$medicaid))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(e4$medicaid == "khongBHMedi") out of length(t$medicaid), null probability 0.5
## X-squared = 2943.1, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.8997804 0.9170195
## sample estimates:
##         p 
## 0.9087608

Với độ tin cậy 95%, ta có tỷ lệ người có bảo hiểm medicaid với tổng thể nằm trong khoảng từ 8.3% và 10.02%. Hay nói cách khác, tỉ lệ người không có bảo hiểm medicaid sẽ chiếm khoảng từ 89.97% đến 91.70%

4.2.4. Ước lượng tỷ lệ người có thu nhập dưới 2.5 ngàn đô

i <- t[t$income == 2.5 | t$income < 2.5,]
prop.test(length(i$income == 30 | i$income <30), length(t$income))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(i$income == 30 | i$income < 30) out of length(t$income), null probability 0.5
## X-squared = 481.81, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.6512709 0.6793509
## sample estimates:
##         p 
## 0.6654562

Với độ tin cậy 95%, ta có tỷ lệ người có thu nhập dưới 2.5 ngàn đô chiếm từ 65.13% đến 67.94%.

4.2.5. Ước lượng tỷ lệ người không có bệnh mãn tính

i3 <- t[t$chronic == 1 | t$chronic < 1,]
prop.test(length(i3$chronic == 1 | i3$chronic <1), length( t$chronic))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(i3$chronic == 1 | i3$chronic < 1) out of length(t$chronic), null probability 0.5
## X-squared = 92.674, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.5578502 0.5872788
## sample estimates:
##         p 
## 0.5726282

Với độ tin cậy 95%, ta có tỷ lệ người có dưới một bệnh mãn tính là 55.79% đến 58.73%

4.2.6. Ứớc lượng người có dưới 2 lần nhập viện

i8 <- t[t$hospital == 2 | t$hospital < 2,]
prop.test(length(i8$hospital == 2 | i8$hospital <2), length( t$hospital))

## 
##  1-sample proportions test with continuity correction
## 
## data:  length(i8$hospital == 2 | i8$hospital < 2) out of length(t$hospital), null probability 0.5
## X-squared = 4051.4, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.9748351 0.9834535
## sample estimates:
##         p 
## 0.9795733

Với độ tin cậy 95%, ta có tỷ lệ người có dưới 2 lần nhập viện là 97.48% đến 98.35%

Chương 5. Hồi quy

5.1. Mô hình xác suất tuyến tính

5.1.1. Biến insurance theo income

t <- t %>% mutate(insurance1 = if_else(insurance == "khongBHCN", 0, 1))
levels(factor(t$insurance1))

## [1] "0" "1"

lm.insc <- glm(insurance1 ~ income, data = t)
summary(lm.insc)

## 
## Call:
## glm(formula = insurance1 ~ income, data = t)
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.719351   0.008193   87.81   <2e-16 ***
## income      0.022591   0.002120   10.66   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 0.1692928)
## 
##     Null deviance: 764.79  on 4405  degrees of freedom
## Residual deviance: 745.57  on 4404  degrees of freedom
## AIC: 4682.1
## 
## Number of Fisher Scoring iterations: 2

Mô hình xác suất tuyến tính của biến insurance và incomed có dạng:

\(\hat\pi = 0.719351 + 0.022591 \times income\)

Có thể thấy, cả hai hệ số trong mô hình đều có ý nghĩa thống kê rất cao (p-value < 0.01), cho thấy thu nhập ảnh hưởng đến việc mua bảo hiểm cá nhân hay không.

Hệ số chặn là 0.719351 cho biết nếu như không có sự ảnh hưởng của thu nhập thì xác suất sở hữu bảo hiểm cá nhân khá thấp (0.271146).
Hệ số của biểu hiện thu nhập cao là 0.022591. Cho biết thu nhập tăng 1 đơn vị thì khả năng sở hữu bảo hiểm cá nhân tăng 0.022591 lần.

5.1.2. Biến insurance theo employed

t$employed <- factor (t$employed)
lm.insem <- glm(insurance1 ~ employed, data = t)
summary(lm.insem)

## 
## Call:
## glm(formula = insurance1 ~ employed, data = t)
## 
## Coefficients:
##                      Estimate Std. Error t value Pr(>|t|)    
## (Intercept)           0.85934    0.01949  44.087  < 2e-16 ***
## employedkhongvieclam -0.09245    0.02058  -4.491 7.26e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 0.1728673)
## 
##     Null deviance: 764.79  on 4405  degrees of freedom
## Residual deviance: 761.31  on 4404  degrees of freedom
## AIC: 4774.1
## 
## Number of Fisher Scoring iterations: 2

Mô hình xác suất tuyến tính của biến insurance và employed có dạng:

\(\hat\pi = 0.85934 - 0.09245\times employedkhongvieclam\)

Có thể thấy, cả hai hệ số trong mô hình đều có ý nghĩa thống kê rất cao (p-value < 0.01), cho thấy tình trạng việc làm có ảnh hưởng đến việc mua bảo hiểm cá nhân hay không. Cụ thể:

Hệ số chặn là 0.85934, cho biết nếu như không có sự ảnh hưởng của tình trạng việc làm thì xác suất sở hữu bảo hiểm cá nhân là 0.85934
Hệ số của tình trạng không có việc làm là - 0.09245, cho biết xác suất để sở hữu bảo hiểm cá nhân của người không có việc làm thấp hơn 0.09245 lần so với người có việc làm.

5.1.3. Biến insurance theo medicaid

t$medicaid <- factor (t$medicaid)

lm.insm <- glm(insurance1 ~ medicaid, data = t)
summary(lm.insm)

## 
## Call:
## glm(formula = insurance1 ~ medicaid, data = t)
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          0.15174    0.01829   8.297   <2e-16 ***
## medicaidkhongBHMedi  0.68742    0.01918  35.831   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 0.1344603)
## 
##     Null deviance: 764.79  on 4405  degrees of freedom
## Residual deviance: 592.16  on 4404  degrees of freedom
## AIC: 3667.1
## 
## Number of Fisher Scoring iterations: 2

Mô hình xác suất tuyến tính của biến insurance và medicaid có dạng:

\(\hat\pi = 0.15174 + 0.68742 \times khongBHMedi\)

Có thể thấy, cả hai hệ số trong mô hình đều có ý nghĩa thống kê rất cao (p-value < \(\alpha\) = 0.001), cho thấy việc có sở hữu bảo hiểm Medicaid có ảnh hưởng đến việc mua bảo hiểm cá nhân. Cụ thể:

Hệ số chặn là 0.15174, cho biết nếu như không có sự ảnh hưởng của bảo hiểm Medicaid thì xác suất sở hữu bảo hiểm cá nhân là khá thấp, khoảng 0.15174
Hệ số của biểu hiện không có BH Medicaid là 0.68742, cho biết khi không sở hữu bảo hiểm của Medicaid thì xác suất để sở hữu bảo hiểm cá nhân là khá cao, khoảng 0.68742

5.1.4. Hồi quy đa biến

lm.in <- glm(insurance1 ~ income + employed + medicaid, data = t)
summary(lm.in)

## 
## Call:
## glm(formula = insurance1 ~ income + employed + medicaid, data = t)
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)           0.153295   0.025964   5.904 3.81e-09 ***
## income                0.012417   0.001926   6.448 1.25e-10 ***
## employedkhongvieclam -0.015943   0.018366  -0.868    0.385    
## medicaidkhongBHMedi   0.666913   0.019354  34.459  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 0.1331537)
## 
##     Null deviance: 764.79  on 4405  degrees of freedom
## Residual deviance: 586.14  on 4402  degrees of freedom
## AIC: 3626.1
## 
## Number of Fisher Scoring iterations: 2

\(\hat\pi = 0.153295 + 0.012417\times income + 0.666913 \times medicaidkhongBHMedi\)

Biến employed không có ý nghĩa do p-value = 0.315 lớn hơn mức ý nghĩa.

Hệ số chặn 0.153295 có nghĩa là nếu như không có tác động của các yếu tố thu nhập, bảo hiểm Medicaid thì xác suất sở hữu bảo hiểm cá nhân là 0.153295

Hệ số của biến thu nhập là 0.012417, cho biết khi thu nhập tăng 1 đơn vị thì xác suất sở hữu BHCN tăng 0.012417.

Hệ số của biểu hiện không có bảo hiểm Medicaid là 0.666913, cho biết người không có bảo hiểm Medicaid có xác suất mua bảo hiểm cá nhân cao hơn 0.666913 so với người có bảo hiểm Medicaid.

5.2. Mô hình hồi quy logit

5.2.1. Biến insurance và incomed

t$insurance <- factor(t$insurance, levels = c('khongBHCN','coBHCN'))
logit.insc <- glm(insurance ~ income, family = binomial(link = 'logit'), data = t)
summary(logit.insc)

## 
## Call:
## glm(formula = insurance ~ income, family = binomial(link = "logit"), 
##     data = t)
## 
## Coefficients:
##             Estimate Std. Error z value Pr(>|z|)    
## (Intercept)  0.63699    0.05828   10.93   <2e-16 ***
## income       0.29325    0.02550   11.50   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 4490.3  on 4404  degrees of freedom
## AIC: 4494.3
## 
## Number of Fisher Scoring iterations: 5

Mô hình hồi quy logit của insurance và incomed là:

\(ln(\frac{\hat\pi}{1-\hat\pi})= 0.63699 + 0.29325\times income\)

Hệ số chặn là 0.63699, có nghĩa là \(ln(\frac{\hat\pi}{1-\hat\pi})\) khi không có tác động của thu nhập thì xác suất sở hữu BHCN là 0.63699

Hệ số của biến thu nhập là 0.29325 . Hệ số này nghĩa là, đối với việc mua bảo hiểm cá nhân thì \(ln(\frac{\hat\pi}{1-\hat\pi})\) của thu nhập tăng thì khả năng sở hữu BHCN tăng.

5.2.2. Biến insurance và employed

logit.insem <- glm(insurance ~ employed, family = binomial(link = 'logit'), data = t)
summary(logit.insem)

## 
## Call:
## glm(formula = insurance ~ employed, family = binomial(link = "logit"), 
##     data = t)
## 
## Coefficients:
##                      Estimate Std. Error z value Pr(>|z|)    
## (Intercept)            1.8098     0.1348  13.422  < 2e-16 ***
## employedkhongvieclam  -0.6190     0.1400  -4.421 9.81e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 4660.4  on 4404  degrees of freedom
## AIC: 4664.4
## 
## Number of Fisher Scoring iterations: 4

Mô hình hồi quy logit của insurance và employ là:

\(ln(\frac{\hat\pi}{1-\hat\pi})= 1.8098 - 0.6190\times employedkhongvieclam\)

Hệ số chặn là 1.8098. Hệ số này có nghĩa là nếu không có sự ảnh hưởng của biểu hiện khác thì \(ln(\frac{\hat\pi}{1-\hat\pi})\) của người có việc làm, sở hữu bảo hiểm cá nhân là 1.8098

Hệ số của biểu hiện không có việc làm bằng - 0.6190. Hệ số này nghĩa là, đối với việc mua bảo hiểm cá nhân thì chênh lệch \(ln(\frac{\hat\pi}{1-\hat\pi})\) của người có không có việc làm và người có việc làm là -0.6190

5.2.3. Biến insurance và medicaid

logit.insm <- glm(insurance ~ medicaid , family = binomial(link = 'logit'), data = t)
summary(logit.insm)

## 
## Call:
## glm(formula = insurance ~ medicaid, family = binomial(link = "logit"), 
##     data = t)
## 
## Coefficients:
##                     Estimate Std. Error z value Pr(>|z|)    
## (Intercept)          -1.7210     0.1390  -12.38   <2e-16 ***
## medicaidkhongBHMedi   3.3730     0.1455   23.18   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 3874.3  on 4404  degrees of freedom
## AIC: 3878.3
## 
## Number of Fisher Scoring iterations: 3

Mô hình hồi quy logit của insurance và medicaid là:

\(ln(\frac{\hat\pi}{1-\hat\pi})= -1.7210 + 3.3730\times khongBHMedi\)

Hệ số chặn là - 1.7210 và có ý nghĩa thống kê. Hệ số này cho biết nếu không có sự ảnh hưởng của biểu hiện khác thì \(ln(\frac{\hat\pi}{1-\hat\pi})\) của người vừa có bảo hiểm Medi, vừa có bảo hiểm cá nhân là -1.7210

Hệ số của biểu hiện không có bảo hiểm Medicaid bằng 3.3730. Hệ số này nghĩa là, đối với việc mua bảo hiểm cá nhân thì chênh lệch \(ln(\frac{\hat\pi}{1-\hat\pi})\) của người không có bảo hiểm Medicaid so với người có bảo hiểm Medicaid là 3.3730

5.2.4. Hồi quy đa biến

logit.in <- glm(insurance ~ income + employed + medicaid , family = binomial(link = 'logit'), data = t)
summary(logit.in)

## 
## Call:
## glm(formula = insurance ~ income + employed + medicaid, family = binomial(link = "logit"), 
##     data = t)
## 
## Coefficients:
##                      Estimate Std. Error z value Pr(>|z|)    
## (Intercept)          -1.84778    0.20766  -8.898  < 2e-16 ***
## income                0.17004    0.02432   6.993  2.7e-12 ***
## employedkhongvieclam -0.08793    0.14853  -0.592    0.554    
## medicaidkhongBHMedi   3.18680    0.14779  21.564  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 3807.2  on 4402  degrees of freedom
## AIC: 3815.2
## 
## Number of Fisher Scoring iterations: 5

\(ln(\frac{\hat\pi}{1-\hat\pi})= -1.84778 + 0.17004 \times income + 3.18680\times medicaidkhongBHMedi\)

Hệ số của biến thu nhập là 0.17004 cho thấy khi thu nhập càng cao thì xác suất sở hữu bảo hiểm cá nhân càng cao.

Hệ số của biểu hiện không có bảo hiểm Medicaid là 3.18680 cho thấy xác suất sở hữu bảo hiểm cá nhân của người không có bảo hiểm Medicaid chênh lệch \(ln(\frac{\hat\pi}{1-\hat\pi})\) = 3.18680 so với người có bảo hiểm Meidcaid.

5.3. Mô hình probit

5.3.1. Biến insurance và incomed

probit.insc <- glm(insurance ~ income, family = binomial(link = 'probit'), data = t)

## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred

summary(probit.insc)

## 
## Call:
## glm(formula = insurance ~ income, family = binomial(link = "probit"), 
##     data = t)
## 
## Coefficients:
##             Estimate Std. Error z value Pr(>|z|)    
## (Intercept)  0.45580    0.03317   13.74   <2e-16 ***
## income       0.13983    0.01280   10.93   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 4508.5  on 4404  degrees of freedom
## AIC: 4512.5
## 
## Number of Fisher Scoring iterations: 6

Mô hình hồi quy như sau:

\(\hat\pi = \Phi(0.45580 + 0.13983\times income)\)

Hệ số chặn là 0.45580 và có ý nghĩa thống kê. Trong trường hợp thu nhập bằng 0, thì xác suất sở hữu bảo hiểm cá nhân là \(\Phi(0.45580)\)

Hệ số của biến income là 0.13983 và có ý nghĩa thống kê. khi thu nhập tăng thì xác suất sở hữu bao hiểm cá nhân tăng \(\Phi(0.13983)\).

5.3.2. Biến insurance và employed

probit.insem <- glm(insurance ~ employed, family = binomial(link = 'probit'), data = t)
summary(probit.insem)

## 
## Call:
## glm(formula = insurance ~ employed, family = binomial(link = "probit"), 
##     data = t)
## 
## Coefficients:
##                      Estimate Std. Error z value Pr(>|z|)    
## (Intercept)           1.07736    0.07300  14.759  < 2e-16 ***
## employedkhongvieclam -0.34870    0.07624  -4.574 4.78e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 4660.4  on 4404  degrees of freedom
## AIC: 4664.4
## 
## Number of Fisher Scoring iterations: 4

Mô hình hồi quy như sau: \(\hat\pi = \Phi(1.07736 - 0.34870 \times employedkhongvieclam)\)

Hệ số chặn là 1.07736, khi các yếu tố khác không đổi thì xác suất của một người có việc làm và có mua bảo hiểm cá nhân là \(\Phi(1.07736)\)

Đối với người có không có việc làm thì xác suất sở hữu bảo hiểm cá nhân là \(\Phi(-0.34870)\)

5.3.3. Biến insurance và medicaid

probit.insm <- glm(insurance ~ medicaid, family = binomial(link = 'probit'), data = t)
summary(probit.insm)

## 
## Call:
## glm(formula = insurance ~ medicaid, family = binomial(link = "probit"), 
##     data = t)
## 
## Coefficients:
##                     Estimate Std. Error z value Pr(>|z|)    
## (Intercept)         -1.02899    0.07616  -13.51   <2e-16 ***
## medicaidkhongBHMedi  2.02001    0.07978   25.32   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 3874.3  on 4404  degrees of freedom
## AIC: 3878.3
## 
## Number of Fisher Scoring iterations: 3

Mô hình hồi quy như sau:

\(\hat\pi = \Phi(-1.02899 + 2.02001 \times khongBHMedi)\)

Hệ số chặn là - 1.02899 và có ý nghĩa thống kê (p-value < \(\alpha\)). Hệ số này có nghĩa là với người có bảo hiểm Medicaid, xác suất để người đó mua thêm bảo hiểm cá nhân là \(\Phi(-1.02899)\)

Hệ số của biểu hiện không có bảo hiểm Medicaid là 2.02001 và có ý nghĩa thống kê. Hệ số cho biết nếu không có bảo hiểm của Medicaid thì xác suất để người đó sở hữu bảo hiểm cá nhân là \(\Phi(2.02001)\)

5.3.4. Hồi quy đa biến

probit.in <- glm(insurance ~ income + employed + medicaid, family = binomial(link = 'probit'), data = t)
summary(probit.in)

## 
## Call:
## glm(formula = insurance ~ income + employed + medicaid, family = binomial(link = "probit"), 
##     data = t)
## 
## Coefficients:
##                      Estimate Std. Error z value Pr(>|z|)    
## (Intercept)          -1.07661    0.11238  -9.580  < 2e-16 ***
## income                0.08466    0.01236   6.851 7.35e-12 ***
## employedkhongvieclam -0.05205    0.08037  -0.648    0.517    
## medicaidkhongBHMedi   1.91434    0.08104  23.623  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 4682.5  on 4405  degrees of freedom
## Residual deviance: 3810.3  on 4402  degrees of freedom
## AIC: 3818.3
## 
## Number of Fisher Scoring iterations: 5

\(\hat\pi = \Phi(-1.07661 + 0.08466 \times income - 1.91434 \times medicaidkhongBHMedi)\)

Hệ số chặn là \(\Phi(-1.07661)\) và có ý nghĩa thống kê. Điều này có nghĩa là nếu không có yếu tố tác động thì xác suất sở hữu bảo hiểm cá nhân là \(\Phi(-1.07661)\)

Xác suất sở hữu bảo hiểm cá nhân của người có thu nhập càng cao là \(\Phi(0.08466)\).

Xác suất sở hữu bảo hiểm cá nhân của người không có bảo hiểm Medicaid là \(\Phi(-1.91434)\).

5.4. Đánh giá mô hình

AIC(lm.insc)

## [1] 4682.078

AIC(logit.insc)

## [1] 4494.283

AIC(probit.insc)

## [1] 4512.45

Đối với mô hình hồi quy biến insurance theo biến incomed thì mô hình logit là tốt nhất vì chỉ số AIC của mô hình này nhỏ nhất.

AIC(lm.insem)

## [1] 4774.138

AIC(logit.insem)

## [1] 4664.403

AIC(probit.insem)

## [1] 4664.403

Trong việc giải thích mối quan hệ giữa tình trạng việc làm và sở hữu bảo hiểm cá nhân thì mô hình logit và probit có hiệu suất tốt hơn. Logit và Probit có hiệu suất tương đương nhau trong việc giải thích mối quan hệ giữa tình trạng việc làm và sở hữu bảo hiểm cá nhân vì chỉ số AIC bằng nhau.

AIC(lm.insm)

## [1] 3667.109

AIC(logit.insm)

## [1] 3878.276

AIC(probit.insm)

## [1] 3878.276

Trong việc giải thích mối quan hệ giữa việc sở hữu bảo hiểm Medicaid và sở hữu bảo hiểm cá nhân thì mô hình xác suất tuyến tính có hiệu suất tốt hơn.

AIC(lm.in)

## [1] 3626.083

AIC(logit.in)

## [1] 3815.234

AIC(probit.in)

## [1] 3818.286

Mô hình xác suất tuyến tính giải thích mối quan hệ giữa biến phụ thuộc insurance và các biến độc lập incomed, employed, medicaid tốt hơn mô hình logit, probit.

Chương 6. Kết luận

Sau quá trình, phân tính và chạy các tính toán liên quan cho biến phụ thuộc ‘insurance’ – bảo hiểm cá nhân, theo các biến độc lập: employed, medicaid, income, hospital và chronic. Ta thấy được rằng, trong 5 biến độc lập nghiên cứu thì chỉ có 3 biến employed – có đi làm hay không, medicaid – được tổ chức BH Medicaid hay không, income – thu nhập mang lại ý nghĩa thống kê, và cho thấy rằng chúng thật sự ảnh hưởng đến việc sở hữu bảo hiểm cá nhân.

Theo quan sát, cho thấy cá nhân có sở hữu bảo hiểm cá nhân nhiều hơn cá nhân không sở hữu bảo hiểm cá nhân. Trong đó, số người đang có việc làm chiếm tỷ lệ khá thấp so với cá nhân không có việc làm. Điều này khá đúng vì người được quan sát ở đây là người trên 66 tuổi. Họ được chính phủ cung cấp một loại bảo hiểm khác là Mediciad.

Đo lường thu nhập của cá nhân được khảo sát rơi vào khoảng 2.5271 ngàn đô là chủ yếu. Kết quả phân tích hồi quy, cũng cho thấy rằng việc thu nhập tỷ lệ thuận với việc họ có mua bảo hiểm cá nhân. Khi đó người có việc làm, mang lại nguồn thu nhập thấp có xu hướng sở hữu bảo hiểm cá nhân cao hơn là những người có việc làm nhưng lại có thu nhập cao. Mà những cá nhân đã sở hữu bảo hiểm Medicaid có xác suất sở hữu bảo hiểm cá nhân thấp hơn những cá nhân không sở hữu lại bảo hiểm do chính phủ cấp (Medicaid). Kết quả hồi quy, cho ta cái nhìn khẳng định rằng việc sở hũu bảo hiểm cá nhân có mối quan hệ cùng chiều với biểu hiện cá nhân không được tổ chức cung cấp bảo hiểm Mediciad.

Tuy nhiên, mối quan hệ của các biến này khi hồi quy đa biến lại cho thấy biến có việc làm hay không lại không có ý nghĩa thông kê. Chỉ có yếu tố thu nhập và bảo hiểm Medicaid (biểu hiện không có bảo hiểm Medicaid) tác động cùng chiều. Mặc khác, mô hình hồi quy logit, probit cho kết quả cho mối quan hệ giữa việc sỡ hữu bảo hiểm theo biến thu nhập tốt hơn. Trong việc giải thích mối quan hệ giữa việc sở hữu bảo hiểm Medicaid và sở hữu bảo hiểm cá nhân thì mô hình xác suất tuyến tính có hiệu suất tốt hơn.

LS0tDQp0aXRsZTogIlBow6JuIHTDrWNoIGPDoWMgeeG6v3UgdOG7kSDhuqNuaCBoxrDhu59uZyDEkeG6v24gdmnhu4djIG11YSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIg0Kb3V0cHV0Og0KICBodG1sX2RvY3VtZW50Og0KICAgIHRvYzogVFJVRQ0KICAgIHRvY19mbG9hdDogVFJVRQ0KICAgIGRmX3ByaW50OiBwYWdlZA0KICAgIGNvZGVfZG93bmxvYWQ6IHRydWUNCiAgICBjb2RlX2ZvbGRpbmc6IGhpZGUNCiAgcGRmX2RvY3VtZW50Og0KICAgIGV4dHJhX2RlcGVuZGVuY2llczoNCiAgICAgIHZpZXRuYW06IHV0ZjgNCiAgICB0b2M6IHllcw0KICAgIG51bWJlcl9zZWN0aW9uczogeWVzDQogIHdvcmRfZG9jdW1lbnQ6DQogICAgdG9jOiB5ZXMNCiAgICBudW1iZXJfc2VjdGlvbnM6IHllcw0KZ2VvbWV0cnk6DQogICAgICAtIGlubmVyPTNjbQ0KICAgICAgLSBvdXRlcj00Y20NCiAgICAgIC0gdG9wPTNjbQ0KICAgICAgLSBib3R0b209NGNtDQogICAgICAtIGhlYWRzZXA9MjJwdA0KICAgICAgLSBoZWFkaGVpZ2h0PTExcHQNCiAgICAgIC0gZm9vdHNraXA9MzNwdA0KICAgICAgLSBpZ25vcmVoZWFkDQogICAgICAtIGlnbm9yZWZvb3QNCiAgICAgIC0gaGVpZ2h0cm91bmRlZA0KYXV0aG9yOiAiTmjDs20gNSAtIE5ndXnhu4VuIFRo4buLIEh14buzbmggTmjGsCwgVHJp4buHdSBUcuG6p24gVGhhbmggTGFuLCBQaOG6oW0gVGjhu4sgVGh1IFRo4bqjbyINCmRhdGU6ICIyMDI0LTA2LTE0Ig0KLS0tDQoNCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQ0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KGVjaG8gPSBUUlVFKQ0KYGBgDQoNCiMgQ2jGsMahbmcgMS4gVOG7lW5nIHF1YW4gduG7gSDEkeG7gSB0w6BpDQoNCiMjIDEuMS4gR2nhu5tpIHRoaeG7h3UgZOG7ryBsaeG7h3UNCg0KTk1FUzE5ODggbMOgIG3hu5l0IGLhu5kgZOG7ryBsaeG7h3UgY2jDqW8sIHRyw61jaCB04burIGvhur90IHF14bqjIGPhu6dhIGN14buZYyBraOG6o28gc8OhdCBE4buLY2ggduG7pSBjaMSDbSBzw7NjIHPhu6ljIGto4buPZSB04bqhaSBIb2EgS+G7syAoTk1FUykuDQpOTUVTIGzDoCBjdeG7mWMga2jhuqNvIHPDoXQgbmjhu69uZyBuZ8aw4budaSDEkcOjIGtow7RuZyDEkcaw4bujYyBraMOhbSwgY2jhu69hIGLhu4duaCB04bqhaSBjw6FjIGPGoSBz4bufIGNoxINtIHPDs2Mgc+G7qWMga2jhu49lIHRyb25nIG3hu5l0IHRo4budaSBnaWFuIGTDoGkgdHJvbmcgc3Xhu5F0IG7Eg20gMTk4NyB2w6AgMTk4OC4NCg0KQuG7mSBk4buvIGxp4buHdSBOTUVTMTk4OCBjaOG7iSBs4bqleSBr4bq/dCBxdeG6oyBraOG6o28gc8OhdCBj4bunYSBuaOG7r25nIG5nxrDhu51pIHRyw6puIDY2IHR14buVaSB2w6AgdsOgbyBuxINtIDE5ODguDQoNClThuq1wIGThu68gbGnhu4d1IE5NRVMxOTg4IGJhbyBn4buTbSA0NDA2IHF1YW4gc8OhdCB24bubaSAxOSBiaeG6v246DQoNCnZpc2l0IDogc+G7kSBs4bqnbiBraMOhbSBi4buHbmggduG7m2kgYsOhYyBzxKkgdOG6oWkgYuG7h25oIHZp4buHbi4NCg0KbnZpc2l0IDogc+G7kSBs4bqnbiBraMOhbSBi4buHbmggduG7m2kgbmjDom4gdmnDqm4geSB04bq/IG5oxrAgeSB0w6EsIMSRaeG7gXUgZMaw4buhbmcgKGtow7RuZyBwaOG6o2kgYsOhYyBzxKkpIHThuqFpIGLhu4duaCB2aeG7h24uDQoNCm92aXNpdHMgOiBz4buRIGzhuqduIGtow6FtIGLhu4duaCB24bubaSBiw6FjIHPEqSB04bqhaSBwaMOybmcga2jDoW0gbmdv4bqhaSB0csO6Lg0KDQpub3Zpc2l0cyA6IHPhu5EgbOG6p24gdGjEg20ga2jDoW0gduG7m2kgbmjDom4gdmnDqm4geSB04bq/IHThuqFpIHBow7JuZyBraMOhbSBuZ2/huqFpIHRyw7ouDQoNCmVtZXJnZW5jeSA6IHPhu5EgbOG6p24gY+G6pXAgY+G7qXUuDQoNCmhvc3BpdGFsIDogc+G7kSBs4bqnbiBuaOG6rXAgdmnhu4duLg0KDQpoZWFsdGg6IHRy4bqhbmcgdGjDoWkgc+G7qWMga2jhu49lLg0KInBvb3IiIGzDoCBraMO0bmcga2jhu49lLCAiYXZlcmFnZSIgbMOgIGLDrG5oIHRoxrDhu51uZywgImV4Y2VsbGVudCIgbMOgIHLhuqV0IGto4buPZS4NCg0KY2hyb25pYyA6IHPhu5EgYuG7h25oIG3Do24gdMOtbmggbcOgIG5nxrDhu51pIMSRxrDhu6NjIGto4bqjbyBzw6F0IG3huq9jIHBo4bqjaS4NCg0KYWRsIDogbmfGsOG7nWkgxJHGsOG7o2Mga2jhuqNvIHPDoXQgY8OzIGfhurdwIGtow7Mga2jEg24gdHJvbmcgdmnhu4djIHRo4buxYyBoaeG7h24gY8OhYyBob+G6oXQgxJHhu5luZyBow6BuZyBuZ8OgeSBuaMawIHbhu4cgc2luaCBjw6EgbmjDom4sIGRpIGNodXnhu4NuLCDEg24gdeG7kW5nLi4uIGhheSBraMO0bmc/DQoibGltaXRlZCIgbMOgIGPDsyBn4bq3cCBraMOzIGtoxINuIHRyb25nIHZp4buHYyB0aOG7sWMgaGnhu4duIGPDoWMgaG/huqF0IMSR4buZbmcgaMOgbmcgbmfDoHksICJub3JtYWwiIGzDoCBraMO0bmcgY8OzLg0KDQpyZWdpb24gOiBuw7NpIHbhu4EgdsO5bmcgbWnhu4FuIGPhu6dhIG5nxrDhu51pIHRo4buxYyBoaeG7h24ga2jhuqNvIHPDoXQuDQoibm9ydGhlYXN0IiBsw6AgbWnhu4FuIMSRw7RuZyBi4bqvYywgIm1pZHdlc3QiIGzDoCBtaeG7gW4gdHJ1bmcgdMOieSBIb2EgS+G7sywgIndlc3QiIGzDoCBtaeG7gW4gdMOieSwgIm9ydGhlciIgbMOgIHbDuW5nIG1p4buBbiBraMOhYy4NCg0KYWdlIDogxJHhu5kgdHXhu5VpICjEkeG7mSB0deG7lWkgY2hpYSAxMCkuDQoNCmFmYW0gOiBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBjw7MgcGjhuqNpIG5nxrDhu51pIGfhu5FjIE3hu7kgZ8O0YyBQaGkgaGF5IGtow7RuZz8NCiJ5ZXMiIGzDoCBwaOG6o2ksICJubyIgbMOgIGtow7RuZy4NCg0KZ2VuZGVyIDogZ2nhu5tpIHTDrW5oLCAibWFsZSIgbMOgIG5hbSwgImZlbWFsZSIgbMOgIG7hu68uDQoNCm1hcnJpZWQgOiDEkcOjIGvhur90IGjDtG4gY2jGsGE/DQoieWVzIiBsw6Aga+G6v3QgaMO0biBy4buTaSwgIm5vIiBsw6AgY2jGsGEga+G6v3QgaMO0bi4NCg0Kc2Nob29sIDogbMOgIHPhu5EgbsSDbSDEkWkgaOG7jWMuDQoNCmluY29tZSA6IHRodSBuaOG6rXAgZ2lhIMSRw6xuaCAoxJHGoW4gduG7iyA6IDEwbmfDoG4gxJHDtCkuDQoNCmVtcGxveWVkIDogbmfGsOG7nWkgxJHGsOG7o2Mga2jhuqNvIHPDoXQgY8OzIMSRYW5nIMSRaSBsw6BtIGtow7RuZz8NCiJ5ZXMiIGzDoCBjw7MsICJubyIgbMOgIGtow7RuZy4NCg0KaW5zdXJhbmNlIDogbmfGsOG7nWkgxJHGsOG7o2Mga2jhuqNvIHPDoXQgY8OzIHThu7EgbXVhIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gaGF5IGtow7RuZz8NCiJ5ZXMiIGzDoCBjw7MsICJubyIgbMOgIGtow7RuZy4NCkLhuqNvIGhp4buDbSBjw6EgbmjDom4gaGF5IGPDsm4gZ+G7jWkgbMOgIGLhuqNvIGhp4buDbSB0xrAgbmjDom4gbMOgIGPDoWNoIGfhu41pIGNodW5nIGNobyBjw6FjIGxv4bqhaSBz4bqjbiBwaOG6qW0gYuG6o28gaGnhu4NtIGtow6FjIG5oYXUgbmjGsCBi4bqjbyBoaeG7g20gbmjDom4gdGjhu40sIGLhuqNvIGhp4buDbSBz4bupYyBraOG7j2UsIGLhuqNvIGhp4buDbSB5IHThur8sLi4uIELhuqNvIGhp4buDbSBjw6EgbmjDom4gYuG6r3QgYnXhu5ljIG5nxrDhu51pIMSRxrDhu6NjIGLhuqNvIGhp4buDbSBj4bqnbiBwaOG6o2kgYuG7jyB0aeG7gW4gcmEgbXVhIGLhuqNvIGhp4buDbS4NCg0KbWVkaWNhaWQgOiBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBjw7MgxJHGsOG7o2MgYuG6o28gaGnhu4NtIGLhu59pIHThu5UgY2jhu6ljIGLhuqNvIGhp4buDbSBNZWRpY2FpZCBraMO0bmc/DQoieWVzIiBsw6AgY8OzLCAibm8iIGzDoCBraMO0bmcuDQpNZWRpY2FpZCBsw6AgY2jGsMahbmcgdHLDrG5oIGLhuqNvIGhp4buDbSBjaMSDbSBzw7NjIHPhu6ljIGto4buPZSwgbcOgIG5nxrDhu51pIMSRxrDhu6NjIGLhuqNvIGhp4buDbSBraMO0bmcgY+G6p24gcGjhuqNpIGNoaSB0aeG7gW4gxJHhu4MgbXVhIGLhuqNvIGhp4buDbS4NCk1lZGljYWlkIGLhuqNvIGhp4buDbSBjaG8gY8O0bmcgZMOibiBN4bu5LCB0aOG7j2EgbcOjbiBt4buZdCBob+G6t2Mgbmhp4buBdSDEkWnhu4F1IGtp4buHbiB24buBIHRodSBuaOG6rXAsIMSR4buZIHR14buVaSBjxaluZyBuaMawIGdpw6EgdHLhu4sgdMOgaSBz4bqjbiBtw6AgbmfGsOG7nWkgxJHDsyBz4bufIGjhu691Lg0KDQojIyAxLjIuIEdp4bubaSB0aGnhu4d1IMSR4buBIHTDoGkNCg0KVHJvbmcgYsOgaSBwaMOibiB0w61jaCBuw6B5LCBuaMOzbSB0w6FjIGdp4bqjIHPhur0gcGjDom4gdMOtY2ggY8OhYyB54bq/dSB04buRIOG6o25oIGjGsOG7n25nIMSR4bq/biB2aeG7h2MgbXVhIGLhuqNvIGhp4buDbSBjw6EgbmjDom4uDQoNCk5ow7NtIMSRw6MgY2jhu41uIDYgYmnhur9uIMSR4buDIG5naGnDqm4gY+G7qXUgYmFvIGfhu5NtOg0KDQotICAgQmnhur9uIMSR4buLbmggdMOtbmggOiBpbnN1cmFuY2UsIGVtcGxveWVkLCBtZWRpY2FpZA0KDQotICAgQmnhur9uIMSR4buLbmggbMaw4bujbmcgOiBpbmNvbWUsIGhvc3BpdGFsLCBjaHJvbmljDQoNCk5naGnDqm4gY+G7qXUgbsOgeSBz4bq9IHBow6JuIHTDrWNoIGPDoWMgeeG6v3UgdOG7kSBuaMawIHRodSBuaOG6rXAsIHTDrG0gdHLhuqFuZyB2aeG7h2MgbMOgbSwgc+G7kSBs4bqnbiBuaOG6rXAgdmnhu4duLCBi4buHbmggbcOjbiB0w61uaCwgxJHGsOG7o2MgTWVkaWNhaWQgYuG6o28gaGnhu4NtIOG6o25oIGjGsOG7n25nIG5oxrAgdGjhur8gbsOgbyDEkeG6v24gdmnhu4djIG3hu5l0IGPDoSBuaMOibiByYSBxdXnhur90IMSR4buLbmggbXVhIGLhuqNvIGhp4buDbSBjw6EgbmjDom4uDQoNCi0gICBCaeG6v24gcGjhu6UgdGh14buZYyA6IGluc3VyYW5jZQ0KDQotICAgQmnhur9uIMSR4buZYyBs4bqtcCA6IGluY29tZSwgaG9zcGl0YWwsIGNocm9uaWMsIGVtcGxveWVkLCBtZWRpY2FpZA0KDQpgYGB7cixtZXNzYWdlPUZBTFNFLHdhcm5pbmc9RkFMU0V9DQpsaWJyYXJ5KHRpZHl2ZXJzZSkNCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkoQUVSKQ0KbGlicmFyeShEVCkNCmxpYnJhcnkoZHBseXIpDQpsaWJyYXJ5KERlc2NUb29scykNCmxpYnJhcnkoZXBpdG9vbHMpDQpsaWJyYXJ5KHNjYWxlcykNCmxpYnJhcnkoVkdBTSkNCmRhdGEoIk5NRVMxOTg4IikNCnQgPC0gTk1FUzE5ODgNCmRhdGF0YWJsZSh0KQ0KYGBgDQoNCiMgQ2jGsMahbmcgMi4gVGjhu5FuZyBrw6ogbcO0IHThuqMgMSBiaeG6v24NCg0KIyMgMi4xLiBCaeG6v24gxJHhu4tuaCB0w61uaA0KDQojIyMgMi4xLjEuIMSQ4buVaSB0w6puIGJp4buDdSBoaeG7h24gY+G7p2EgY8OhYyBiaeG6v24gxJHhu4tuaCB0w61uaA0KDQrEkOG7gyB0csOhbmggbmjhuqdtIHPhu7EgbmjhuqdtIGzhuqtuIHRyb25nIHF1w6EgdHLDrG5oIHBow6JuIHTDrWNoLCB0YSDEkeG7lWkgdMOqbiBjw6FjIGJp4buDdSBoaeG7h24gInllcyIsICJubyIgY+G7p2EgMyBiaeG6v24gxJHhu4tuaCB0w61uaCBpbnN1cmFuY2UsIGVtcGxveWVkLCBtZWRpY2FpZA0KDQotICAgxJDhu5FpIHbhu5tpIGJp4bq/biBpbnN1cmFuY2UgOiAieWVzIiDEkeG7lWkgdGjDoG5oICJjb0JIQ04iLCBuZ2jEqWEgbMOgIE5nxrDhu51pIGto4bqjbyBzw6F0IGPDsyBtdWEgYuG6o28gaGnhu4NtIGPDoSBuaMOibi4NCiAgICAibm8iIMSR4buVaSB0aMOgbmggImtob25nQkhDTiIsIG5naMSpYSBsw6AgTmfGsOG7nWkga2jhuqNvIHPDoXQga2jDtG5nIG11YSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuLg0KDQotICAgxJDhu5FpIHbhu5tpIGJp4bq/biBlbXBsb3llZCA6ICJ5ZXMiIMSR4buVaSB0aMOgbmggImNvdmllY2xhbSIsIG5naMSpYSBsw6AgTmfGsOG7nWkga2jhuqNvIHPDoXQgxJFhbmcgY8OzIHZp4buHYyBsw6BtLg0KICAgICJubyIgxJHhu5VpIHRow6BuaCAia2hvbmd2aWVjbGFtIiwgbmdoxKlhIGzDoCBOZ8aw4budaSBraOG6o28gc8OhdCDEkWFuZyBraMO0bmcgY8OzIHZp4buHYyBsw6BtLg0KDQotICAgxJDhu5FpIHbhu5tpIGJp4bq/biBtZWRpY2FpZCA6ICJ5ZXMiIMSR4buVaSB0aMOgbmggImNvQkhNZWRpIiwgbmdoxKlhIGzDoCBOZ8aw4budaSBraOG6o28gc8OhdCDEkcaw4bujYyBi4bqjbyBoaeG7g20gYuG7n2kgdOG7lSBjaOG7qWMgTWVkaWNhaWQuDQogICAgIm5vIiDEkeG7lWkgdGjDoG5oICJraG9uZ0JITWVkaSIsIG5naMSpYSBsw6AgTmfGsOG7nWkga2jhuqNvIHPDoXQga2jDtG5nIMSRxrDhu6NjIGLhuqNvIGhp4buDbSBi4bufaSB04buVIGNo4bupYyBNZWRpY2FpZC4NCg0KYGBge3J9DQp0JGluc3VyYW5jZSA8LSBmYWN0b3IodCRpbnN1cmFuY2UsIGxldmVscyA9IGMoJ3llcycsICdubycpLCBsYWJlbHMgPSBjKCdjb0JIQ04nLCdraG9uZ0JIQ04nKSkNCg0KdCRlbXBsb3llZCA8LSBmYWN0b3IodCRlbXBsb3llZCwgbGV2ZWxzID0gYygneWVzJywgJ25vJyksIGxhYmVscyA9IGMoJ2NvdmllY2xhbScsICdraG9uZ3ZpZWNsYW0nICkpDQoNCnQkbWVkaWNhaWQgPC0gZmFjdG9yKHQkbWVkaWNhaWQsIGxldmVscyA9IGMoJ3llcycsJ25vJyksIGxhYmVscyA9IGMoJ2NvQkhNZWRpJywna2hvbmdCSE1lZGknICkpDQpgYGANCg0KIyMjIDIuMS4yLiBCaeG6v24gaW5zdXJhbmNlDQoNCmBgYHtyfQ0KdGFibGUodCRpbnN1cmFuY2UpDQoNCnByb3AudGFibGUodGFibGUodCRpbnN1cmFuY2UpKQ0KDQpnZ3Bsb3QodCwgYWVzKHggPSBpbnN1cmFuY2UsIHk9IGFmdGVyX3N0YXQoY291bnQpKSkrDQogIGdlb21fYmFyKGZpbGwgPSAncGluaycpKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsID0gc2NhbGVzOjpwZXJjZW50KGFmdGVyX3N0YXQoY291bnQvc3VtKGNvdW50KSkpKSwgc3RhdCA9ICdjb3VudCcsIGNvbG9yID0gJ2JsYWNrJywgdmp1c3QgPSAtLjEpKw0KICBsYWJzKHRpdGxlPSAgJ1Phu5EgbMaw4bujbmcgbmfGsOG7nWkgbXVhIGLhuqNvIGhp4buDbSBjw6EgbmjDom4nKQ0KYGBgDQoNCkPDsyB0aOG7gyB0aOG6pXkgY8OzIHPhu7EgY2jDqm5oIGzhu4djaCBs4bubbiBnaeG7r2Egc+G7kSBsxrDhu6NuZyBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIGPDoSBuaMOibiB2w6Agc+G7kSBsxrDhu6NuZyBuZ8aw4budaSBraMO0bmcgY8OzLg0KVHJvbmcgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdGhhbSBnaWEga2jhuqNvIHPDoXQgY8OzIDM0MjEgbmfGsOG7nWkgY8OzIG11YSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIChjaGnhur9tIDc3LjY0JSksIGPDsm4gbOG6oWkgbMOgIDk4NSBuZ8aw4budaSBraMO0bmcgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIChjaGnhur9tIDIyLjM1JSkuDQoNCiMjIyAyLjEuMy4gQmnhur9uIGVtcGxveWVkDQoNCmBgYHtyfQ0KDQphZGRtYXJnaW5zKHRhYmxlKHQkZW1wbG95ZWQpKQ0KDQpnZ3Bsb3QodCwgYWVzKHggPSBlbXBsb3llZCwgeT0gYWZ0ZXJfc3RhdChjb3VudCkpKSsNCiAgZ2VvbV9iYXIoZmlsbCA9ICdibGFjaycpKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsID0gc2NhbGVzOjpwZXJjZW50KGFmdGVyX3N0YXQoY291bnQvc3VtKGNvdW50KSkpKSwgc3RhdCA9ICdjb3VudCcsIGNvbG9yID0gJ2JsYWNrJywgdmp1c3QgPSAtLjUpKw0KICBsYWJzKHRpdGxlPSAgJ1TDrG5oIHRy4bqhbmcgdmnhu4djIGzDoG0nKQ0KYGBgDQoNCkThu7FhIHbDoG8gYmnhu4N1IMSR4buTLCB0YSB0aOG6pXkgxJHGsOG7o2MgdHJvbmcgc+G7kSBuaOG7r25nIG5nxrDhu51pIHRoYW0gZ2lhIGto4bqjbyBzw6F0IHRow6wgcGjhuqduIGzhu5tuIGzDoCBuZ8aw4budaSBraMO0bmcgY8OzIHZp4buHYy4NCkJp4buDdSBoaeG7h24ga2jDtG5nIGPDsyB2aeG7h2MgbMOgbSBjw7Mgc+G7sSBjaMOqbmggbOG7h2NoIHLDtSByw6BuZyBzbyB24bubaSBiaeG7g3UgaGnhu4duIGPDsyB2aeG7h2MgbMOgbS4NCkPhu6UgdGjhu4MsIHRyb25nIHThu5VuZyBz4buRIDQ0MDYgbmfGsOG7nWksIHRow6wgY8OzIDQ1NSBuZ8aw4budaSDEkWFuZyBjw7Mgdmnhu4djIChjaOG7iSBjaGnhur9tIDEwJSkgdsOgIDM5NTEgbmfGsOG7nWkgY8OybiBs4bqhaSBsw6Aga2jDtG5nIGPDsyB2aeG7h2MgbMOgbSAoY2hp4bq/bSB04bqtbiA5MCUpLg0KDQojIyMgMi4xLjQuIEJp4bq/biBtZWRpY2FpZA0KDQpgYGB7cn0NCmFkZG1hcmdpbnModGFibGUodCRtZWRpY2FpZCkpDQoNCmdncGxvdCh0LCBhZXMoeCA9IG1lZGljYWlkLCB5PSBhZnRlcl9zdGF0KGNvdW50KSkpKw0KICBnZW9tX2JhcihmaWxsID0gJ2xpZ2h0Ymx1ZScpKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsID0gc2NhbGVzOjpwZXJjZW50KGFmdGVyX3N0YXQoY291bnQvc3VtKGNvdW50KSkpKSwgc3RhdCA9ICdjb3VudCcsIGNvbG9yID0gJ2JsYWNrJywgdmp1c3QgPSAtLjUpKw0KICBsYWJzKHRpdGxlPSAgJ0PDsyDEkcaw4bujYyB04buVIGNo4bupYyBNZWRpY2FpZCBi4bqjbyBoaeG7g20ga2jDtG5nPycpDQpgYGANCg0KQmnhur9uIG1lZGljYWlkIGPFqW5nIGPDsyBz4buxIGNow6puaCBs4buHY2ggbOG7m24gZ2nhu69hIGhhaSBiaeG7g3UgaGnhu4duLg0KVOG6oWkgdGjhu51pIMSRaeG7g20ga2jhuqNvIHPDoXQsIGNo4buJIGPDsyA0MDIgbmfGsOG7nWkgbMOgIMSRxrDhu6NjIHThu5UgY2jhu6ljIE1lZGljYWlkIGLhuqNvIGhp4buDbSAoY2hp4bq/bSBwaOG6p24gbmjhu48gOSUpLCA0MDA0IG5nxrDhu51pIGPDsm4gbOG6oWkga2jDtG5nIMSRxrDhu6NjIHThu5UgY2jhu6ljIE1lZGljYWlkIHTDoGkgdHLhu6MgYuG6o28gaGnhu4NtIChjaGnhur9tIDkxJSkuDQoNCiMjIDIuMi4gQmnhur9uIMSR4buLbmggbMaw4bujbmcNCg0KIyMjIDIuMi4xLiBCaeG6v24gaW5jb21lDQoNCkNodXnhu4NuIMSR4buVaSB04burIGJp4bq/biDEkeG7i25oIGzGsOG7o25nIHNhbmcgYmnhur9uIMSR4buLbmggdMOtbmggYuG6sW5nIGPDoWNoIGNoaWEga2hv4bqjbmcgdHJvbmcgZOG7ryBsaeG7h3UuDQoNCsSQ4bqndSB0acOqbiwgdGEgbmjDrG4gdOG7lW5nIHF1YW4gY8OhYyBiaeG6v24gaW5jb21lLg0KDQpgYGB7cn0NCnN1bW1hcnkodCRpbmNvbWUpDQpgYGANCg0KLSAgIFRydW5nIGLDrG5oIHRodSBuaOG6rXAgY+G7p2EgbmjDum5nIG5nxrDhu51pIMSRxrDhu6NjIGto4bqjbyBzw6F0IG7hurFtIOG7nyAyLjUyNzEuDQoNCi0gICBN4bupYyB0aHUgbmjhuq1wIHRo4bqlcCBuaOG6pXQgbMOgIC0xLjAxMjUsIGNhbyBuaOG6pXQgbMOgIDU0LjgzNTEuDQoNCi0gICBUcm9uZyDEkcOzIOG7nyBt4bupYyB0aHUgbmjhuq1wIDAuOTEyMiB0aMOsIGPDsyAyNSUgZOG7ryBsaeG7h3UgdGjhuqVwIGjGoW4gbMaw4bujbmcgdGh1IG5o4bqtcCBuw6B5IOG7niBt4bupYyAxLjY5ODIgdGjDrCBjw7MgNTAlIGThu68gbGnhu4d1IGPDsyB0aHUgbmjhuq1wIHRo4bqlcCBoxqFuIHbDoCA1MCUgZOG7ryBsaeG7h3UgY2FvIGjGoW4gbeG7qWMgbsOgeS4NCiAgICBDw7JuIOG7nyBt4bupYyAzLjE3MjggdGjDrCBjw7MgNzUlIGThu68gbGnhu4d1IGPDsyBsxrDhu6NuZyB0aHUgbmjhuq1wIHRo4bqlcCBoxqFuIG3hu6ljIG7DoHkuDQoNClRp4bq/cCB0aGVvLCB0YSB04bqhbyBt4buZdCBiaeG6v24gbeG7m2kgdMOqbiBsw6AgaW5jb21lMS4NCkJp4bq/biBuw6B5IHPhur0gY2jhu6lhIGThu68gbGnhu4d1IHbhu4EgdGh1IG5o4bqtcCBj4bunYSBuZ8aw4budaSB0aGFtIGdpYSBraOG6o28gc8OhdCB2w6AgY2hvIGJp4bq/dCDhu58gbeG7qWMgdGh1IG5o4bqtcCBuw6BvIGzDoCB0aOG6pXAsIG3hu6ljIG7DoG8gbMOgIGNhby4NCg0KLSAgIEtob+G6o25nIHThu6sgLTEuMSDEkeG6v24gMi41IHPhur0gY8OzIHTDqm4gbMOgICd0aGFwJyAsIMSRxrDhu6NjIHhlbSBt4bupYyB0aHUgbmjhuq1wIHRo4bqlcC4NCg0KLSAgIEtob+G6o25nIHThu6sgMi41IMSR4bq/biA1NSB0w6puIGzDoCAnY2FvJywgxJHGsOG7o2MgeGVtIG5oxrAgbeG7qWMgdGh1IG5o4bqtcCBjYW8uDQoNCmBgYHtyfQ0KdCRpbmNvbWVkIDwtIGN1dCh0JGluY29tZSwgYnJlYWtzID0gYygtMS4xICwyLjUsNTUpLCBsYWJlbHMgPSBjKCJ0aGFwIiwiY2FvIikpDQoNCnRhYmxlKHQkaW5jb21lZCkNCg0KZ2dwbG90KHQsIGFlcyh4ID0gaW5jb21lZCwgeT0gYWZ0ZXJfc3RhdChjb3VudCkpKSsNCiAgZ2VvbV9iYXIoZmlsbCA9ICdwaW5rJykrDQogIGdlb21fdGV4dChhZXMobGFiZWwgPSBzY2FsZXM6OnBlcmNlbnQoYWZ0ZXJfc3RhdChjb3VudC9zdW0oY291bnQpKSkpLCBzdGF0ID0gJ2NvdW50JywgY29sb3IgPSAnYmxhY2snLCB2anVzdCA9IC0uNSkrDQogIGxhYnModGl0bGU9ICAnVGh1IG5o4bqtcCBj4bunYSBuaOG7r25nIG5nxrDhu51pIGto4bqjbyBzw6F0JykNCmBgYA0KDQpE4buxYSB0csOqbiBiaeG7g3UgxJHhu5MsIHRhIGPDsyB0aOG7gyB0aOG6pXkgbmjhu69uZyBuZ8aw4budaSBjw7MgdGh1IG5o4bqtcCB0aOG6pXAgY2hp4bq/bSBz4buRIGzGsOG7o25nIG5oaeG7gXUgaMahbiBy4bqldCBuaGnhu4F1IHNvIHbhu5tpIG5nxrDhu51pIGPDsyB0aHUgbmjhuq1wIGNhby4NCkPhu6UgdGjhu4MsIGPDsyAyOTQzIG5nxrDhu51pIHRyb25nIHThu5VuZyBz4buRIDQ0MDYgbmfGsOG7nWkgY8OzIHRodSB0aOG6rXAgdGjhuqVwIChjaGnhur9tIDY3JSksIGPDsm4gbOG6oWkgbMOgIDE0NjMgbmfGsOG7nWkgY8OzIHRodSBuaOG6rXAgY2FvIChjaGnhur9tIDMzJSkuDQoNCiMjIyAyLjIuMi4gQmnhur9uIGhvc3BpdGFsDQoNClTGsMahbmcgdOG7sSBjaG8gY8OhYyBiaeG6v24gaG9zcGl0YWwNCg0KYGBge3J9DQpzdW1tYXJ5ICh0JGhvc3BpdGFsKQ0KYGBgDQoNCi0gICBUcnVuZyBiw6xuaCBuaOG7r25nIG5nxrDhu51pIGto4bqjbyBzw6F0IGtow7RuZyBuaOG6rXAgdmnhu4duIGzhuqduIG7DoG8sIChnacOhIHRy4buLIHRydW5nIGLDrG5oIGzDoCAwLjI5NikNCg0KLSAgIFPhu5EgbOG6p24gdGjhuqVwIG5o4bqldCBsw6AgMCwgY2FvIG5o4bqldCBsw6AgOCBs4bqnbi4NCg0KLSAgIEPDsyB0aOG7gyB0aOG6pXkg4bufIHThu6kgcGjDom4gduG7iyB0aOG7qSBuaOG6pXQsIHRo4bupIGhhaSB2w6AgdGjhu6kgYmEsIHRow6wgc+G7kSBs4bqnbiBuaOG6rXAgdmnhu4duIMSR4buBdSBsw6AgMCBs4bqnbi4NCiAgICDEkGnhu4F1IG7DoHkgY2hvIHRo4bqleSBz4buRIGzhuqduIG5o4bqtcCB2aeG7h24gY+G7p2Egbmjhu69uZyBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBy4bqldCDDrXQgdGjhuq1tIGNow60gbMOgIGtow7RuZyBjw7MuDQoNClRp4bq/cCB0aGVvLCB0YSB04bqhbyBt4buZdCBiaeG6v24gbeG7m2kgdMOqbiBsw6AgaG9zcGl0YWwxLg0KQmnhur9uIG7DoHkgc+G6vSBjaOG7qWEgZOG7ryBsaeG7h3UgduG7gSBz4buRIGzhuqduIG5o4bqtcCB2aeG7h24gY+G7p2EgbmfGsOG7nWkgdGhhbSBnaWEga2jhuqNvIHPDoXQgdsOgIGNobyBiaeG6v3Qg4bufIG3huqV5IGzhuqduIG7DoG8gbMOgIHRo4bqlcCwgbeG6pXkgbOG6p24gbMOgIGNhby4NCg0KLSAgIEtob+G6o25nIHThu6sgLTAuNSDEkeG6v24gMi41IHPhur0gY8OzIHTDqm4gbMOgICdpdCcgLCDEkcaw4bujYyB4ZW0gbMOgIGtow7RuZyBuaOG6rXAgdmnhu4duLg0KDQotICAgS2hv4bqjbmcgdOG7qyAyLjUgxJHhur9uIDkgdMOqbiBsw6AgJ25oaWV1JywgxJHGsOG7o2MgeGVtIG5oxrAgbmjhuq1wIHZp4buHbiBuaGnhu4F1Lg0KDQpgYGB7cn0NCnQkaG9zcGl0YWwxIDwtIGN1dCh0JGhvc3BpdGFsLCBicmVha3MgPSBjKC0wLjUsMi41LDkpLCBsYWJlbHMgPSBjKCJpdCIsIm5oaWV1IikpDQoNCmFkZG1hcmdpbnModGFibGUodCRob3NwaXRhbDEpKQ0KDQpnZ3Bsb3QodCwgYWVzKHggPSBob3NwaXRhbDEsIHk9IGFmdGVyX3N0YXQoY291bnQpKSkrDQogIGdlb21fYmFyKGZpbGwgPSAnYmx1ZScpKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsID0gc2NhbGVzOjpwZXJjZW50KGFmdGVyX3N0YXQoY291bnQvc3VtKGNvdW50KSkpKSwgc3RhdCA9ICdjb3VudCcsIGNvbG9yID0gJ2JsYWNrJywgdmp1c3QgPSAtLjUpKw0KICBsYWJzKHRpdGxlPSAgJ1Phu5EgbOG6p24gbmjhuq1wIHZp4buHbiBj4bunYSBuaOG7r25nIG5nxrDhu51pIGto4bqjbyBzw6F0JykNCmBgYA0KDQpE4buxYSB0csOqbiBiaeG7g3UgxJHhu5MsIHRhIGPDsyB0aOG7gyB0aOG6pXkgY8OzIGtob+G6o25nIGNow6puaCBs4buHY2ggcuG6pXQgbOG7m24gZ2nhu69hIGdp4buvYSBuZ8aw4budaSBuaOG6rXAgdmnhu4duIMOtdCB2w6AgbmfGsOG7nWkgbmjhuq1wIHZp4buHbiBuaGnhu4F1Lg0KQ+G7pSB0aOG7gywgc+G7kSBuZ8aw4budaSBraMO0bmcgbmjhuq1wIHZp4buHbiBjaGnhur9tIHThuq1uIDM1NDEgbmfGsOG7nWkgKDgwJSksIHRyb25nIGtoaSBuaMOzbSBuZ8aw4budaSBjw7Mgbmjhuq1wIHZp4buHbiBjaGnhur9tIHPhu5EgbMaw4bujbmcgcuG6pXQgw610LCBjaOG7iSBjw7MgODY1IG5nxrDhu51pKDIwJSkuDQoNCiMjIyAyLjIuMy4gQmnhur9uIGNocm9uaWMNCg0KxJDhu5FpIHbhu5tpIGJp4bq/biBjaHJvbmljOg0KDQpgYGB7cn0NCnN1bW1hcnkgKHQkY2hyb25pYykNCmBgYA0KDQotICAgVHJ1bmcgYsOsbmggc+G7kSBi4buHbmggbcOjbiB0w61uaCBtw6Agbmjhu69uZyBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBt4bqvYyBwaOG6o2kgbMOgIDEuNTQyDQoNCi0gICBT4buRIGLhu4duaCB0aOG6pXAgbmjhuqV0IGzDoCAwLCBjYW8gbmjhuqV0IGzDoCA4Lg0KDQotICAgVHJvbmcgxJHDsyDhu58gbeG7qWMgYuG7iyAxIGLhu4duaCBtw6NuIHTDrW5oIHRow6wgY8OzIDI1JSBk4buvIGxp4buHdSB0aOG6pXAgaMahbiBt4bupYyBuw6B5Lg0KICAgIOG7niBt4bupYyAxIGLhu4duaCB0aMOsIGPDsyA1MCUgZOG7ryBsaeG7h3UgY8OzIHPhu5EgbMaw4bujbmcgYuG7h25oIHRo4bqlcCBoxqFuIHbDoCA1MCUgZOG7ryBsaeG7h3UgY2FvIGjGoW4gbeG7qWMgbsOgeS4NCiAgICBDw7JuIOG7nyBt4bupYyBi4buLIDIgYuG7h25oIG3Do24gdMOtbmggdGjDrCBjw7MgNzUlIGThu68gbGnhu4d1IHRo4bqlcCBoxqFuIG3hu6ljIG7DoHkuDQoNClRp4bq/cCB0aGVvLCB0YSB04bqhbyBt4buZdCBiaeG6v24gbeG7m2kgdMOqbiBsw6AgY2hyb25pYzEuDQpCaeG6v24gbsOgeSBz4bq9IGNo4bupYSBk4buvIGxp4buHdSB24buBIHPhu5EgbMaw4bujbmcgYuG7h25oIG3Do24gdMOtbmggbcOgIG5nxrDhu51pIHRoYW0gZ2lhIGto4bqjbyBzw6F0IG3huq9jIHBo4bqjaS4NCg0KLSAgIEtob+G6o25nIHThu6sgLTAuNSDEkeG6v24gMC41IHPhur0gY8OzIHTDqm4gbMOgICdraG9uZ2JlbmgnICwgxJHGsOG7o2MgaGnhu4N1IGzDoCBuZ8aw4budaSBraOG6o28gc8OhdCBraMO0bmcgbeG6r2MgYuG7h25oIG3Do24gdMOtbmguDQoNCi0gICBLaG/huqNuZyB04burIDAuNSDEkeG6v24gOSB0w6puIGzDoCAnY29iZW5oJywgxJHGsOG7o2MgaGnhu4N1IGzDoCBuZ8aw4budaSBraOG6o28gc8OhdCBt4bqvYyAxIGhv4bq3YyBuaGnhu4F1IGjGoW4gMSBi4buHbmggbcOjbiB0w61uaC4NCg0KYGBge3J9DQp0JGNocm9uaWMxIDwtIGN1dCh0JGNocm9uaWMsIGJyZWFrcyA9IGMoLTAuNSwwLjUsOSksIGxhYmVscyA9IGMoImtob25nYmVuaCIsImNvYmVuaCIpKQ0KDQphZGRtYXJnaW5zKHRhYmxlKHQkY2hyb25pYzEpKQ0KDQpnZ3Bsb3QodCwgYWVzKHggPSBjaHJvbmljMSwgeT0gYWZ0ZXJfc3RhdChjb3VudCkpKSsNCiAgZ2VvbV9iYXIoZmlsbCA9ICdwaW5rJykrDQogIGdlb21fdGV4dChhZXMobGFiZWwgPSBzY2FsZXM6OnBlcmNlbnQoYWZ0ZXJfc3RhdChjb3VudC9zdW0oY291bnQpKSkpLCBzdGF0ID0gJ2NvdW50JywgY29sb3IgPSAnYmxhY2snLCB2anVzdCA9IC0uNSkrDQogIGxhYnModGl0bGU9ICAnU+G7kSBuZ8aw4budaSBt4bqvYyBi4buHbmggbcOjbiB0w61uaCcpDQoNCmBgYA0KDQpDw7MgdGjhu4MgdGjhuqV5LCBjw7Mgc+G7sSBjaMOqbmggbOG7h2NoIGzhu5tuIGdp4buvYSBuaOG7r25nIG5nxrDhu51pIGtow7RuZyBjw7MgYuG7h25oIG3Do24gdMOtbmggdsOgIG5ow7NtIG5nxrDhu51pIGPDsyBi4buHbmggbcOjbiB0w61uaC4NCg0KU+G7kSBuZ8aw4budaSBraMO0bmcgbeG6r2MgYuG7h25oIG3Do24gdMOtbmggbMOgIDEwMjUgbmfGsOG7nWkgKGNoaeG6v20gMjMlKSwgY8OybiBs4bqhaSBsw6Agbmjhu69uZyBuZ8aw4budaSBt4bqvYyBuaGnhu4F1IGLhu4duaCBtw6NuIHTDrW5oLCBjw7MgMzM4MSBuZ8aw4budaSAoY2hp4bq/bSA3NyUpLg0KDQojIENoxrDGoW5nIDMuIFRo4buRbmcga8OqIG3DtCB04bqjIGPDoWMgY+G6t3AgYmnhur9uDQoNCiMjIDMuMS4gUGjDom4gdMOtY2ggYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgdGh1IG5o4bqtcCAoaW5jb21lKQ0KDQojIyMgMy4xLjEuIELhuqNuZyB04bqnbiBzdeG6pXQNCg0KYGBge3J9DQoNCmFkZG1hcmdpbnModGFibGUodCRpbmNvbWVkLHQkaW5zdXJhbmNlKSkNCg0KYWRkbWFyZ2lucyhwcm9wLnRhYmxlKHRhYmxlKHQkaW5jb21lZCx0JGluc3VyYW5jZSkpKQ0KDQp0ICU+JSBnZ3Bsb3QoYWVzKHggPSBpbmNvbWVkLCB5ID0gYWZ0ZXJfc3RhdChjb3VudCksIGZpbGwgPSBpbnN1cmFuY2UpKSArDQogIGdlb21fYmFyKHBvc2l0aW9uID0gImRvZGdlIiwgY29sb3IgPSAiYmxhY2siKSArDQogIGdlb21fdGV4dChhZXMobGFiZWwgPSBhZnRlcl9zdGF0KGNvdW50KSksIHN0YXQgPSAnY291bnQnLCBjb2xvciA9ICdibGFjaycsIHZqdXN0ID0gLTAuMjUsIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjkpKSArDQogIHRoZW1lX2NsYXNzaWMoKSArIA0KICBsYWJzKHggPSAiVGh1IG5o4bqtcCIsDQogICAgICAgeSA9ICJT4buRIGzGsOG7o25nIiwNCiAgICAgICB0aXRsZSA9ICJUaHUgbmjhuq1wIHbDoCBi4bqjbyBoaeG7g20iKSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoImtob25nQkhDTiIgPSAiZG9kZ2VyYmx1ZSIsICJjb0JIQ04iID0gImJsdWUiKSkgDQoNCmBgYA0KDQoNCkLhuqNuZyB04bqnbiBzdeG6pXQgdsOgIGJp4buDdSDEkeG7kyB24buBIHZp4buHYyBz4bufIGjhu691IGLhuqNvIGhp4buDbSB2w6AgdGh1IG5o4bqtcCBjaG8gdGjhuqV5Og0KDQotIE5ow7NtIGPDsyB0aHUgbmjhuq1wIGNhbzoNCg0KQ8OzIDE5MCBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY2hp4bq/bSB04bu3IGzhu4cgbMOgIDQuMzElIHThu5VuZyBz4buRIG5nxrDhu51pIHRyb25nIGLhu5kgZOG7ryBsaeG7h3UNCg0KQ8OzIDEyODQgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY2hp4bq/bSAyOS4xNCUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgYuG7mSBk4buvIGxp4buHdQ0KDQotIE5ow7NtIGPDsyB0aHUgbmjhuq1wIHRo4bqlcDoNCg0KY8OzIDc5NSBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY2hp4bq/bSAxOC4wNCUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgYuG7mSBk4buvIGxp4buHdQ0KDQpjw7MgMjEzNyBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBjaGnhur9tIDQ4LjUwJSB04buVbmcgc+G7kSBuZ8aw4budaSB0cm9uZyBk4buvIGxp4buHdQ0KDQpDw7MgdGjhu4MgdGjhuqV5IG5o4buvbmcgbmfGsOG7nWkgY8OzIHRodSBuaOG6rXAgdGjhuqVwIHRoxrDhu51uZyBjw7MgeHUgaMaw4bubbmcgbXVhIGLhuqNvIGhp4buDbSB5IHThur8gbmhp4buBdSBoxqFuIG5o4buvbmcgbmfGsOG7nWkgY8OzIHRodSBuaOG6rXAgY2FvIA0KDQoNCiMjIyAzLjEuMi4gUmVsYXRpdmUgUmlzaw0KDQpgYGB7cn0NCnRoOCA8LSB0YWJsZSh0JGluY29tZWQsdCRpbnN1cmFuY2UpDQphZGRtYXJnaW5zKHRoOCkNClJlbFJpc2sodGg4KQ0KDQpgYGANCk5o4bqtbiB4w6l0Og0KDQpU4bu3IGzhu4cgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIOG7nyBuaMOzbSB0aHUgbmjhuq1wIHRo4bqlcCBi4bqxbmcgMC44MzY3MDYzIHNvIHbhu5tpIHThu7cgbOG7hyBz4bufIGjhu691IGLhuqNvIGhp4buDbSBj4bunYSBuZ8aw4budaSB0aHUgbmjhuq1wIGNhby4gDQoNCiMjIyAzLjEuMy4gVOG7tyBs4buHIGNow6puaCAtIE9kZHJhdGlvDQoNCmBgYHtyfQ0KT2Rkc1JhdGlvKHRoOCkNCmBgYA0KDQpL4bq/dCBxdeG6oyBPZGQgUmF0aW8gbMOgIDAuMzk3NzY0NSBjaG8gdGjhuqV5IHLhurFuZyB04buJIGzhu4cgbmfGsOG7nWkgdGh1IG5o4bqtcCB0aOG6pXAgY8OzIGto4bqjIG7Eg25nIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiB0aOG6pXAgaMahbiBuZ8aw4budaSB0aHUgbmjhuq1wIGNhby4NCg0KVsOsIHPhu7EgY2jDqm5oIGzhu4djaCBnaeG7r2EgbmjDs20gY8OzIEJIQ04gdsOgIGtow7RuZyBjw7MgQkhDTiBj4bunYSAibmfGsOG7nWkgdGh1IG5o4bqtcCB0aOG6pXAiIHRo4bqlcCBoxqFuIHPhu7EgY2jDqm5oIGzhu4djaCBnaeG7r2EgbmjDs20gY8OzIEJIQ04gdsOgIGtow7RuZyBjw7MgQkhDTiBj4bunYSBuZ8aw4budaSB0aHUgbmjhuq1wIGNhby4gDQoNCiMjIDMuMi4gUGjDom4gdMOtY2ggYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uIGVtcGxveWVkIChjw7Mgdmnhu4djIGzDoG0pDQoNCiMjIyAzLjIuMS4gQuG6o25nIHThuqduIHN14bqldA0KDQpgYGB7cn0NCmFkZG1hcmdpbnModGFibGUodCRlbXBsb3llZCx0JGluc3VyYW5jZSkpDQoNCnJvdW5kKCh0YWJsZSh0JGVtcGxveWVkLHQkaW5zdXJhbmNlKS9zdW0odGFibGUodCRlbXBsb3llZCx0JGluc3VyYW5jZSkpKSoxMDAsIDIpDQoNCnQgJT4lIGdncGxvdChhZXMoeCA9IGVtcGxveWVkLCB5ID0gYWZ0ZXJfc3RhdChjb3VudCksIGZpbGwgPSBpbnN1cmFuY2UpKSArDQogIGdlb21fYmFyKHBvc2l0aW9uID0gImRvZGdlIiwgY29sb3IgPSAiYmxhY2siKSArDQogIGdlb21fdGV4dChhZXMobGFiZWwgPSBhZnRlcl9zdGF0KGNvdW50KSksIHN0YXQgPSAnY291bnQnLCBjb2xvciA9ICdibGFjaycsIHZqdXN0ID0gLTAuMjUsIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjkpKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogIGxhYnMoeCA9ICJWaeG7h2MgbMOgbSIsDQogICAgICAgeSA9ICJT4buRIGzGsOG7o25nIiwNCiAgICAgICB0aXRsZSA9ICJWaeG7h2MgbMOgbSB2w6AgYuG6o28gaGnhu4NtIikgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKCJraG9uZ0JIQ04iID0gImRhcmtvcmFuZ2UiLCAiY29CSENOIiA9ICJyZWQiKSkgDQpgYGANCg0KTmjhuq1uIHjDqXQ6DQoNCkLhuqNuZyB04bqnbiBzdeG6pXQgdsOgIGJp4buDdSDEkeG7kyB24buBIHZp4buHYyBz4bufIGjhu691IGLhuqNvIGhp4buDbSB2w6AgY8OzIHZp4buHYyBsw6BtIGNobyB0aOG6pXk6DQoNCi0gTmjDs20ga2jDtG5nIGPDsyB2aeG7h2MgbMOgbSBsw6A6DQoNCkPDsyA5MjEgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGNoaeG6v20gdOG7tyBs4buHIGzDoCAyMC45JSB04buVbmcgc+G7kSBuZ8aw4budaSB0cm9uZyBi4buZIGThu68gbGnhu4d1DQoNCkPDsyAzMDMwIG5nxrDhu51pIGPDsyBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGNoaeG6v20gdOG7tyBs4buHIGzDoCA2OC43NyUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgYuG7mSBk4buvIGxp4buHdS4NCg0KLSBOaMOzbSBjw7Mgdmnhu4djIGzDoG0NCg0KQ8OzIDY0IG5nxrDhu51pIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBjaGnhur9tIHThu7cgbOG7hyBsw6AgMS40NSUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgYuG7mSBk4buvIGxp4buHdQ0KDQpDw7MgMzMxIG5nxrDhu51pIGPDsyBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGNoaeG6v20gOC44NyUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgYuG7mSBk4buvIGxp4buHdQ0KDQotIE5oxrAgduG6rXk6DQoNClThu5VuZyBj4buZbmcgY8OzIDk4NSBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSB5IHThur8gY8OhIG5ow6JuLCBjaGnhur9tIDIyLjM1JSB04buVbmcgc+G7kSBuZ8aw4budaSB0cm9uZyB04bqtcCBk4buvIGxp4buHdS4NCg0KVOG7lW5nIGPhu5luZyBjw7MgMzQyMSBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIHkgdOG6vyBjw6EgbmjDom4sIGNoaeG6v20gNzcuNjUlIHThu5VuZyBz4buRIG5nxrDhu51pIHRyb25nIHThuq1wIGThu68gbGnhu4d1Lg0KDQpU4bu3IGzhu4cgbmfGsOG7nWkga2jDtG5nIGPDsyB2aeG7h2MgbMOgbSBjw7MgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBjaGnhur9tIHPhu5EgbMaw4bujbmcgbOG7m24gaMahbiBuaOG7r25nIG5nxrDhu51pIGPDsyB2aeG7h2MgbMOgbS4gQ2hvIHRo4bqleSBuaOG7r25nIG5nxrDhu51pIGtow7RuZyBjw7Mgdmnhu4djIGzDoG0gY8OzIHh1IGjGsOG7m25nIG11YSBi4bqjbyBoaeG7g20gbmhp4buBdSBoxqFuDQoNCiMjIyAzLjIuMi4gUmVsYXRpdmUgUmlzaw0KDQpgYGB7cn0NCnQzIDwtIHRhYmxlKHQkZW1wbG95ZWQsdCRpbnN1cmFuY2UpDQp0Mw0KUmVsUmlzayh0MykNCmBgYA0KDQpOaOG6rW4geMOpdDoNCg0KVOG7tyBs4buHIFJlbGF0aXZlIFJpc2sgxJHGsOG7o2MgdMOtbmggbMOgIDEuMTIwNTQ2IGNobyBiaeG6v3QgcuG6sW5nIGto4bqjIG7Eg25nIGPhu6dhIG5o4buvbmcgbmfGsOG7nWkgY8OzIHZp4buHYyBsw6BtIHbDoCBz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY2FvIGjGoW4ga2hv4bqjbmcgMTIlIHNvIHbhu5tpIG5o4buvbmcgbmfGsOG7nWkgY8OzIHZp4buHYyBsw6BtLg0KDQojIyMgMy4yLjMuIFThu7cgbOG7hyBjaMOqbmggLSBPZGRyYXRpbw0KDQpgYGB7cn0NCg0KT2Rkc1JhdGlvKHQzKQ0KYGBgDQoNCk5o4bqtbiB4w6l0Og0KDQpL4bq/dCBxdeG6oyBPZGQgUmF0aW8gY+G7p2EgbmfGsOG7nWkgY8OzIHZp4buHYyBsw6BtIHbDoCBraMO0bmcgY8OzIHZp4buHYyBsw6BtIGtob+G6o25nIDEuODU3MDA4IMSRaeG7gXUgbsOgeSBjw7MgbmdoxKlhIGzDoCBvZGQgY+G7p2EgbmjDs20gbmfGsOG7nWkgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20sIGPDsyB2aeG7h2MgbMOgbSBjYW8gaMahbiAxLjg1NzAwOCBs4bqnbiBzbyB24bubaSBvZGQgbmfGsOG7nWkga2jDtG5nIGPDsyB2aeG7h2MgbMOgbS4NCg0KQ8OzIHPhu7Ega2jDoWMgYmnhu4d0IHRyb25nIGvhur90IHF14bqjIGdp4buvYSBSZWxhdGl2ZSByaXNrLCBvZGQgcmF0aW9zIHbDoCBi4bqjbmcgdOG6p24gc3XhuqV0LiDEkGnhu4F1IG7DoHkgY8OzIHRo4buDIGdp4bqjaSB0aMOtY2ggcuG6sW5nIGRvIHThu7cgbOG7hyBz4bufIGjhu691IEJIQ04gdHJvbmcgbmjDs20gbmfGsOG7nWkgY8OzIHZp4buHYyBsw6BtIGNhbyBoxqFuIHThu7cgbOG7hyBz4bufIGjhu691IEJIQ04gY+G7p2EgbmjDs20gbmfGsOG7nWkga2jDtG5nIHZp4buHYyBsw6BtLg0KDQojIyAzLjMuIFBow6JuIHTDrWNoIGJp4bq/biBpbnN1cmFuY2UgKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibikgdsOgIGJp4bq/biBtZWRpY2FpZA0KDQojIyMgMy4zLjEuIELhuqNuZyB04bqnbiBzdeG6pXQNCg0KYGBge3J9DQphZGRtYXJnaW5zKHRhYmxlKHQkbWVkaWNhaWQsdCRpbnN1cmFuY2UpKQ0KDQpyb3VuZCgodGFibGUodCRtZWRpY2FpZCx0JGluc3VyYW5jZSkvc3VtKHRhYmxlKHQkbWVkaWNhaWQsdCRpbnN1cmFuY2UpKSkqMTAwLCAyKQ0KDQp0ICU+JSBnZ3Bsb3QoYWVzKHggPSBtZWRpY2FpZCwgeSA9IGFmdGVyX3N0YXQoY291bnQpLCBmaWxsID0gaW5zdXJhbmNlKSkgKw0KICBnZW9tX2Jhcihwb3NpdGlvbiA9ICJkb2RnZSIsIGNvbG9yID0gImJsYWNrIikgKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsID0gYWZ0ZXJfc3RhdChjb3VudCkpLCBzdGF0ID0gJ2NvdW50JywgY29sb3IgPSAnYmxhY2snLCB2anVzdCA9IC0wLjI1LCBwb3NpdGlvbiA9IHBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC45KSkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICBsYWJzKHggPSAiQuG6o28gaGnhu4NtIG1lZGljYWlkIiwNCiAgICAgICB5ID0gIlPhu5EgbMaw4bujbmciLA0KICAgICAgIHRpdGxlID0gIkLhuqNvIGhp4buDbSBtZWRpY2FpZCB2w6AgYuG6o28gaGnhu4NtIGPDoSBuaMOibiIpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYygia2hvbmdCSENOIiA9ICJ0b21hdG8iLCAiY29CSENOIiA9ICJvcmFuZ2UiKSkgDQogICAgICAgICAgIA0KYGBgDQoNCkLhuqNuZyB04bqnbiBzdeG6pXQgdsOgIGJp4buDdSDEkeG7kyB24buBIHZp4buHYyBz4bufIGjhu691IGLhuqNvIGhp4buDbSB2w6AgYuG6o28gaGnhu4NtIG1lZGljYWlkIGNobyB0aOG6pXk6DQoNCi0gTmjDs20ga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gbWVkaWNhaWQgbMOgOg0KDQpDw7MgNjQ0IG5nxrDhu51pIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBjaGnhur9tIHThu7cgbOG7hyBsw6AgMTQuNjIlIHThu5VuZyBz4buRIG5nxrDhu51pIHRyb25nIGLhu5kgZOG7ryBsaeG7h3UNCg0KQ8OzIDMzNjAgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY2hp4bq/bSB04bu3IGzhu4cgbMOgIDc2LjI2JSB04buVbmcgc+G7kSBuZ8aw4budaSB0cm9uZyBi4buZIGThu68gbGnhu4d1Lg0KDQotIE5ow7NtIGPDsyBi4bqjbyBoaeG7g20gbWVkaWNhaWQNCg0KQ8OzIDM0MSBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY2hp4bq/bSB04bu3IGzhu4cgbMOgIDcuNzQlIHThu5VuZyBz4buRIG5nxrDhu51pIHRyb25nIGLhu5kgZOG7ryBsaeG7h3UNCg0KQ8OzIDYxIG5nxrDhu51pIGPDsyBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGNoaeG6v20gMS4zOCUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgYuG7mSBk4buvIGxp4buHdQ0KDQotIE5oxrAgduG6rXk6DQoNClThu5VuZyBj4buZbmcgY8OzIDk4NSBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSB5IHThur8gY8OhIG5ow6JuLCBjaGnhur9tIDIyLjM1JSB04buVbmcgc+G7kSBuZ8aw4budaSB0cm9uZyB04bqtcCBk4buvIGxp4buHdS4NCg0KVOG7lW5nIGPhu5luZyBjw7MgMzQyMSBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIHkgdOG6vyBjw6EgbmjDom4sIGNoaeG6v20gNzcuNjUlIHThu5VuZyBz4buRIG5nxrDhu51pIHRyb25nIHThuq1wIGThu68gbGnhu4d1Lg0KDQrEkGnhu4F1IG7DoHkgcGjhuqNuIMOhbmggcuG6sW5nIG5o4buvbmcgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gw610IGtoaSBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkLCB0cm9uZyBraGkgbmfGsOG7o2MgbOG6oWksIG5o4buvbmcgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGPDsyB4dSBoxrDhu5tuZyBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkIG5oaeG7gXUgaMahbi4NCg0KIyMjIDMuMy4yLiBSZWxhdGl2ZSBSaXNrDQoNCmBgYHtyfQ0KdDQgPC0gdGFibGUodCRtZWRpY2FpZCx0JGluc3VyYW5jZSkNCnQ0DQpSZWxSaXNrKHQ0KQ0KYGBgDQoNCk5o4bqtbiB4w6l0Og0KDQpU4bu3IGzhu4cgUmVsYXRpdmUgUmlzayDEkcaw4bujYyB0w61uaCBsw6AgMC4xODA4MjUsIGNobyBiaeG6v3QgcuG6sW5nIGto4bqjIG7Eg25nIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBj4bunYSBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIG1lZGljYWlkIGLhurFuZyAwLjE4MDgyNSBs4bqnbiBzbyB24bubaSBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBNZWRpY2FpZC4NCg0KIyMjIDMuMy4zLiBU4bu3IGzhu4cgY2jDqm5oIC0gT2RkcmF0aW8NCg0KYGBge3J9DQpPZGRzUmF0aW8odDQpDQoNCmBgYA0KS+G6v3QgcXXhuqMgT2RkIFJhdGlvIGPhu6dhIG5o4buvbmcgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gbWVkaWNhaWQgc28gduG7m2kgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gbWVkaWNhaWQga2hv4bqjbmcgMC4wMzQyODY0MSAuIFThu7cgbOG7hyBz4bufIGjhu691IEJIQ04gY+G7p2EgbmfGsOG7nWkgxJHGsOG7o2MgdMOgaSB0cuG7oyB0aOG6pXAgaMahbiB04bu3IGzhu4cgc+G7nyBo4buvdSBCSENOIGPhu6dhIG5nxrDhu51pIGtow7RuZyDEkcaw4bujYyB0w6BpIHRy4bujLg0KDQoNCiMjIDMuNC4gUGjDom4gdMOtY2ggYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uIGhvc3BpdGFsIChz4buRIGzhuqduIG5o4bqtcCB2aeG7h24pDQoNCk3hu6VjIMSRw61jaDogRMO5bmcgYmnhur9uIGhvc3BpdGFsIHbDoCBpbnN1cmFuY2UgxJHhu4MgbmdoacOqbiBj4bupdSBz4buxIHTGsMahbmcgdMOhYyBob+G6t2MgbeG7kWkgcXVhbiBo4buHIGdp4buvYSB2aeG7h2Mgc+G7rSBk4bulbmcgZOG7i2NoIHbhu6UgeSB04bq/ICjEkW8gYuG6sW5nIHPhu5EgbOG6p24gbmjhuq1wIHZp4buHbikgdsOgIHTDrG5oIHRy4bqhbmcgYuG6o28gaGnhu4NtIHkgdOG6vyBj4bunYSBuZ8aw4budaSBkw6JuIHbDoG8gbsSDbSAxOTg4Lg0KDQpQaMOibiB0w61jaDogQ8OzIHRo4buDIHBow6JuIHTDrWNoIHhlbSBuaOG7r25nIG5nxrDhu51pIGPDsyBi4bqjbyBoaeG7g20gY8OzIHh1IGjGsOG7m25nIHPhu60gZOG7pW5nIGThu4tjaCB24bulIHkgdOG6vyAobmjhuq1wIHZp4buHbikgbmjGsCB0aOG6vyBuw6BvIHNvIHbhu5tpIG5o4buvbmcgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20uDQoNCiMjIyAzLjQuMS4gQuG6o25nIHThuqduIHN14bqldA0KDQpgYGB7cn0NCiMgQuG6o25nIHThuqduIHN14bqldCBj4bunYSBob3NwaXRhbCB2w6AgaW5zdXJhbmNlOg0KYWRkbWFyZ2lucyh0YWJsZSh0JGhvc3BpdGFsMSx0JGluc3VyYW5jZSkpDQoNCiMgVOG7iSBs4buHIHBo4bqnbiB0csSDbSBj4bunYSBi4bqjbmcgdOG6p24gc3XhuqV0Og0KYWRkbWFyZ2lucyhyb3VuZCgodGFibGUodCRob3NwaXRhbDEsdCRpbnN1cmFuY2UpL3N1bSh0YWJsZSh0JGhvc3BpdGFsLHQkaW5zdXJhbmNlKSkpKjEwMCwgMikpDQoNCiMgVHLhu7FjIHF1YW4gaMOzYSBk4buvIGxp4buHdSANCnQgJT4lIGdncGxvdChhZXMoeCA9IGhvc3BpdGFsMSwgeSA9IGFmdGVyX3N0YXQoY291bnQpLCBmaWxsID0gaW5zdXJhbmNlKSkgKw0KICBnZW9tX2Jhcihwb3NpdGlvbiA9ICJkb2RnZSIsIGNvbG9yID0gImJsYWNrIikgKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsID0gYWZ0ZXJfc3RhdChjb3VudCkpLCBzdGF0ID0gJ2NvdW50JywgY29sb3IgPSAnYmxhY2snLCB2anVzdCA9IC0wLjI1LCBwb3NpdGlvbiA9IHBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC45KSkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICBsYWJzKHggPSAidOG6p24gc3XhuqV0IG5o4bqtcCB2aeG7h24iLA0KICAgICAgIHkgPSAiU+G7kSBsxrDhu6NuZyIsDQogICAgICAgdGl0bGUgPSAidOG6p24gc3XhuqV0IG5o4bqtcCB2aeG7h24gduG7m2kgYuG6o28gaGnhu4NtIE1lZGljYWlkIikgKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKCJraG9uZ0JIQ04iID0gImxpZ2h0Ymx1ZSIsICJjb0JIQ04iID0gImxpZ2h0Z3JlZW4iKSkgDQpgYGANCg0KDQoqKipOaOG6rW4geMOpdDoqKioNCg0KTmjhuq1uIHjDqXQgduG7gSBz4buxIHBow6JuIHTDrWNoIG3hu5FpIHF1YW4gaOG7hyBnaeG7r2EgYmnhur9uIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibihpbnN1cmFuY2UpIHbDoCBz4buRIGzhuqduIG5o4bqtcCB2aeG7h24gKGhvc3BpdGFsKQ0KKipQaMOibiBi4buRIHPhu5EgbOG6p24gbmjhuq1wIHZp4buHbiB0aGVvIHTDrG5oIHRy4bqhbmcgYuG6o28gaGnhu4NtOioqDQoNCk5ow7NtIMOtdCBs4bqnbiBuaOG6rXAgdmnhu4duICgiaXQiKToNCg0KLSBDw7MgOTY0IG5nxrDhu51pIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIHkgdOG6vyB0xrAgbmjDom4sIGNoaeG6v20gMjEuODglIHThu5VuZyBz4buRIG5nxrDhu51pIHRyb25nIHThuq1wIGThu68gbGnhu4d1Lg0KDQotIEPDsyAzMzUyIG5nxrDhu51pIGPDsyBi4bqjbyBoaeG7g20geSB04bq/IHTGsCBuaMOibiwgY2hp4bq/bSA3Ni4wOCUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgdOG6rXAgZOG7ryBsaeG7h3UuDQoNCk5ow7NtIG5oaeG7gXUgbOG6p24gbmjhuq1wIHZp4buHbiAoIm5oaWV1Iik6DQoNCi0gQ8OzIDIxIG5nxrDhu51pIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIHkgdOG6vyB0xrAgbmjDom4sIGNoaeG6v20gMC40OCUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgdOG6rXAgZOG7ryBsaeG7h3UuDQoNCi0gQ8OzIDY5IG5nxrDhu51pIGPDsyBi4bqjbyBoaeG7g20geSB04bq/IHTGsCBuaMOibiwgY2hp4bq/bSAxLjU3JSB04buVbmcgc+G7kSBuZ8aw4budaSB0cm9uZyB04bqtcCBk4buvIGxp4buHdS4NCg0KDQpU4buVbmcgY+G7mW5nIGPDsyA5ODUgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20geSB04bq/IHTGsCBuaMOibiwgY2hp4bq/bSAyMi4zNiUgdOG7lW5nIHPhu5EgbmfGsOG7nWkgdHJvbmcgdOG6rXAgZOG7ryBsaeG7h3UuDQoNClThu5VuZyBj4buZbmcgY8OzIDM0MjEgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSB5IHThur8gdMawIG5ow6JuLCBjaGnhur9tIDc3LjY1JSB04buVbmcgc+G7kSBuZ8aw4budaSB0cm9uZyB04bqtcCBk4buvIGxp4buHdS4NCg0KJFx0byQgKipL4bq/dCBsdeG6rW46KioNCg0KQuG6o28gaGnhu4NtIHkgdOG6vyB0xrAgbmjDom4gdsOgIHPhu5EgbOG6p24gbmjhuq1wIHZp4buHbjoNCg0KVOG7tyBs4buHIG5nxrDhu51pIGPDsyBi4bqjbyBoaeG7g20geSB04bq/IHTGsCBuaMOibiBuaOG6rXAgdmnhu4duIG5oaeG7gXUgbOG6p24gKG5ow7NtICJuaGlldSIpIGzDoCA2OSBuZ8aw4budaSwgdHJvbmcga2hpIHThu7cgbOG7hyBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBjw6EgbmjDom5uaOG6rXAgdmnhu4duIG5oaeG7gXUgbOG6p24gbMOgIDIxIG5nxrDhu51pLg0KDQrEkGnhu4F1IG7DoHkgY2hvIHRo4bqleSBy4bqxbmcgbmjhu69uZyBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIHkgdOG6vyB0xrAgbmjDom4gY8OzIHh1IGjGsOG7m25nIHPhu60gZOG7pW5nIGThu4tjaCB24bulIHkgdOG6vyBuaGnhu4F1IGjGoW4gc28gduG7m2kgbmjhu69uZyBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbS4NCg0KU+G7sSBjaMOqbmggbOG7h2NoIHLDtSByw6BuZzoNCg0KU+G7kSBsxrDhu6NuZyBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIHkgdOG6vyB0xrAgbmjDom4gbmhp4buBdSBoxqFuIGjhurNuIHNvIHbhu5tpIHPhu5EgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20sIGNoaeG6v20gZ+G6p24gNzglIHThu5VuZyBz4buRIG3huqt1Lg0KDQpDaMOqbmggbOG7h2NoIG7DoHkgY8OzIHRo4buDIHBo4bqjbiDDoW5oIG3hu5l0IHBo4bqnbiBz4buxIHBo4buVIGJp4bq/biBj4bunYSBi4bqjbyBoaeG7g20geSB04bq/IHTGsCBuaMOibiB2w6BvIG7Eg20gMTk4OCwgaG/hurdjIGPDsyB0aOG7gyBsacOqbiBxdWFuIMSR4bq/biBjw6FjIHnhur91IHThu5Ega2luaCB04bq/LXjDoyBo4buZaSBraMOhYyBjaMawYSDEkcaw4bujYyBwaMOibiB0w61jaCDhu58gxJHDonkuDQoNCg0KIyMjIDMuNC4yLiBSZWxhdGl2ZSBSaXNrDQoNCmBgYHtyfQ0KIyBUw61uaCBjaOG7iSBz4buRIGxpw6puIHF1YW4gxJHhur9uIHLhu6dpIHJvIChSZWxhdGl2ZSBSaXNrIHbDoCBSaXNrIFJhdGlvKQ0KDQp0aDIgPC0gdGFibGUodCRob3NwaXRhbDEsdCRpbnN1cmFuY2UpDQp0aDINClJlbFJpc2sodGgyKQ0KYGBgDQpU4bu3IGzhu4cgUmVsYXRpdmUgUmlzayDEkcaw4bujYyB0w61uaCBsw6AgMS4wMTMwMTUsIGNobyBiaeG6v3QgcuG6sW5nIGto4bqjIG7Eg25nIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBj4bunYSBuZ8aw4budaSDDrXQgbmjhuq1wIHZp4buHbiBjYW8gaMahbiAxLjAxMzAxNSBs4bqnbiBzbyB24bubaSBuZ8aw4budaSBuaOG6rXAgdmnhu4duIG5oaeG7gXUuDQoNCiMjIyAzLjQuMy4gVOG7tyBs4buHIGNow6puaCAtIE9kZHJhdGlvDQoNCmBgYHtyfQ0KT2Rkc1JhdGlvKHRoMikNCg0KYGBgDQpL4bq/dCBxdeG6oyBPZGQgUmF0aW8gMS4wNTgyNzIgKE9SID4gMSApIGNobyB0aOG6pXkgcuG6sW5nIHThu7cgbOG7hyBz4bufIGjhu691IEJIQ04gY+G7p2EgbmfGsOG7nWkgw610IG5o4bqtcCB2aeG7h24gY2FvIGjGoW4gbmfGsOG7nWkgbmjhuq1wIHZp4buHbiBuaGnhu4F1Lg0KDQojIyAzLjUuIFBow6JuIHTDrWNoIGJp4bq/biBpbnN1cmFuY2UgKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIHTGsCBuaMOibikgdsOgIGJp4bq/biBjaHJvbmljIChz4buRIGLhu4duaCBtw6NuIHTDrW5oIG3DoCBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBt4bqvYyBwaOG6o2kpDQoNCiMjIyAzLjUuMS4gQuG6o25nIHThuqduIHN14bqldA0KDQpgYGB7cn0NCg0KYWRkbWFyZ2lucyh0YWJsZSh0JGNocm9uaWMxLHQkaW5zdXJhbmNlICkpDQoNCmFkZG1hcmdpbnMocm91bmQoKHRhYmxlKHQkY2hyb25pYzEsdCRpbnN1cmFuY2UpL3N1bSh0YWJsZSh0JGNocm9uaWMxLHQkaW5zdXJhbmNlKSkpKjEwMCwgMikpDQoNCiMgVHLhu7FjIHF1YW4gaMOzYSBk4buvIGxp4buHdQ0KDQp0ICU+JSBnZ3Bsb3QoYWVzKHggPSBjaHJvbmljMSwgeSA9IGFmdGVyX3N0YXQoY291bnQpLCBmaWxsID0gaW5zdXJhbmNlKSkgKw0KICBnZW9tX2Jhcihwb3NpdGlvbiA9ICJkb2RnZSIsIGNvbG9yID0gImJsYWNrIikgKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsID0gYWZ0ZXJfc3RhdChjb3VudCkpLCBzdGF0ID0gJ2NvdW50JywgY29sb3IgPSAnYmxhY2snLCB2anVzdCA9IC0wLjI1LCBwb3NpdGlvbiA9IHBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC45KSkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICBsYWJzKHggPSAiU+G7kSBi4buHbmggbcOjbiB0w61uaCIsDQogICAgICAgeSA9ICJT4buRIGzGsOG7o25nIiwNCiAgICAgICB0aXRsZSA9ICJT4buRIGLhu4duaCBtw6NuIHTDrW5oIHbhu5tpIELhuqNvIGhp4buDbSBjw6EgbmjDom4iKSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoImtob25nQkhDTiI9ICJsaWdodGJsdWUiLCAiY29CSENOIiA9ICJsaWdodGdyZWVuIikpIA0KYGBgDQoNCg0KKioqTmjhuq1uIHjDqXQ6KioqDQoNCi0gTmjDs20ga2jDtG5nIG3huq9jIGLhu4duaCBtw6NuIHTDrW5oOg0KDQpDw7MgYuG6o28gaGnhu4NtOiA4MDAgbmfGsOG7nWkgKDE4LjE2JSB0cm9uZyB04buVbmcgc+G7kSBuZ8aw4budaSkNCg0KS2jDtG5nIGPDsyBi4bqjbyBoaeG7g206IDIyNSBuZ8aw4budaSAoNS4xMSUgdHJvbmcgdOG7lW5nIHPhu5EgbmfGsOG7nWkpDQoNCg0KLSBOaMOzbSBt4bqvYyBi4buHbmggbcOjbiB0w61uaDoNCg0KQ8OzIGLhuqNvIGhp4buDbTogMjYyMSBuZ8aw4budaSAoIDU5LjQ5JSB0cm9uZyB04buVbmcgc+G7kSBuZ8aw4budaSkNCg0KS2jDtG5nIGPDsyBi4bqjbyBoaeG7g206IDc2MCBuZ8aw4budaSAoMTcuMjUlIHRyb25nIHThu5VuZyBz4buRIG5nxrDhu51pKQ0KDQoNCg0KIyMjIDMuNS4yLiBSZWxhdGl2ZSBSaXNrDQoNCmBgYHtyfQ0KdGgzIDwtIHRhYmxlKHQkY2hyb25pYzEsdCRpbnN1cmFuY2UpDQp0aDMNClJlbFJpc2sodGgzKQ0KYGBgDQoNCioqKk5o4bqtbiB4w6l0OioqKg0KDQogTmjhuq1uIHjDqXQ6DQoNClThu7cgbOG7hyBSZWxhdGl2ZSBSaXNrIMSRxrDhu6NjIHTDrW5oIGzDoCAxLjAwNjgwMi4gR8OtYSB0cuG7iyBuw6B5IGNobyBiaeG6v3Qga2jhuqMgbsSDbmcgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGPhu6dhIG5o4buvbmcgbmfGsOG7nWkga2jDtG5nIGPDsyBi4buHbmggbcOjbiB0w61uaCBi4bqxbmcgMS4wMDY4MDIgbOG6p24gc28gduG7m2kgbmjhu69uZyBuZ8aw4budaSBjw7MgYuG7h25oIG3Do24gdMOtbmguDQoNCg0KIyMjIDMuNS4zLiBU4bu3IGzhu4cgY2jDqm5oIC0gT2RkcmF0aW8NCg0KYGBge3J9DQpPZGRzUmF0aW8odGgzKQ0KYGBgDQoNCioqKk5o4bqtbiB4w6l0OioqKg0KDQrEkOG7kWkgduG7m2kgbmjDs20ga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gY8OhIG5ow6JuLCBvZGRzID0gMS4wMzA5ODkuIMSQaeG7gXUgbsOgeSBuZ2jEqWEgbMOgIMSR4buRaSB24bubaSB2aeG7h2Mgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuLCB0aMOsIHThu7cgbOG7hyBjaMOqbmggY+G7p2Egbmjhu69uZyBuZ8aw4budaSBjw7MgYuG7h25oIG3Do24gdMOtbmggdsOgIGtow7RuZyBi4buHbmggbcOjbiB0w61uaCBsw6AgMS4wMzA5ODkuDQoNCiMgQ2jGsMahbmcgNC4gVGjhu5FuZyBrw6ogc3V5IGRp4buFbg0KDQojIyA0LjEuIEtp4buDbSDEkeG7i25oIHTDrW5oIMSR4buZYyBs4bqtcCBj4bunYSAyIGJp4bq/bg0KDQojIyMgNC4xLjEuIEtp4buDbSB0cmEgdMOtbmggxJHhu5ljIGzhuq1wIGNobyAgYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uIGluY29tZWQgKHRodSBuaOG6rXApDQoNCkdp4bqjIHRodXnhur90ICRIXzAkIDogYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uIGluY29tZWQgKHRodSBuaOG6rXApIMSR4buZYyBs4bqtcC4NCkdp4bqjIHRodXnhur90ICRIXzEkOiBiaeG6v24gaW5zdXJhbmNlIChz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4pIHbDoCBiaeG6v24gaW5jb21lZCAodGh1IG5o4bqtcCkga2jDtG5nIMSR4buZYyBs4bqtcA0KDQpgYGB7cn0NCmNoaXNxLnRlc3QodGg4KQ0KYGBgDQpOaOG6rW4geMOpdDogDQoNCkLDoWMgYuG7jyBnaeG6oyB0aHV54bq/dCAkSF8wJCwgY2jhuqVwIG5o4bqtbiBnaeG6oyB0aHV54bq/dCAkSF8xJCwgaW5zdXJhbmNlKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibikgdsOgIHRodSBuaOG6rXAga2jDtG5nIMSR4buZYyBs4bqtcC4NCg0KDQojIyMgNC4xLjIuIEtp4buDbSB0cmEgdMOtbmggxJHhu5ljIGzhuq1wIGNobyBoYWkgYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uIGVtcGxveWVkICh2aeG7h2MgbMOgbSkNCg0KR2nhuqMgdGh1eeG6v3QgJEhfMCQgOiBiaeG6v24gaW5zdXJhbmNlIChz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4pIHbDoCBiaeG6v24gZW1wbG95ZWQgKHZp4buHYyBsw6BtKSDEkeG7mWMgbOG6rXAuDQpHaeG6oyB0aHV54bq/dCAkSF8xJDogYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uIGluY29tZWQgKHRodSBuaOG6rXApIGtow7RuZyDEkeG7mWMgbOG6rXANCg0KYGBge3J9DQpjaGlzcS50ZXN0KHQzKQ0KYGBgDQoNCk5o4bqtbiB4w6l0IA0KQsOhYyBi4buPIGdp4bqjIHRodXnhur90ICRIXzAkLCBjaOG6pXAgbmjhuq1uIGdp4bqjIHRodXnhur90ICRIXzEkLGJp4bq/biBpbnN1cmFuY2UgKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibikgdsOgIGJp4bq/bmVtcGxveWVkICh2aeG7h2MgbMOgbSkga2jDtG5nIMSR4buZYyBs4bqtcA0KDQoNCiMjIyA0LjEuMy4gS2nhu4NtIHRyYSB0w61uaCDEkeG7mWMgbOG6rXAgY2hvIGhhaSBiaeG6v24gaW5zdXJhbmNlIChz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4pIHbDoCBiaeG6v24gbWVkaWNhaWQgDQoNCmBgYHtyfQ0KY2hpc3EudGVzdCh0NCkNCmBgYA0KDQoNCk5o4bqtbiB4w6l0IA0KQsOhYyBi4buPIGdp4bqjIHRodXnhur90ICRIXzAkLCBjaOG6pXAgbmjhuq1uIGdp4bqjIHRodXnhur90ICRIXzEkLCBiaeG6v24gaW5zdXJhbmNlIChz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4pIHbDoCBiaeG6v24gbWVkaWNhaWQgIGtow7RuZyDEkeG7mWMgbOG6rXANCg0KIyMjIDQuMS40LiBLaeG7g20gdHJhIHTDrW5oIMSR4buZYyBs4bqtcCBjaG8gaGFpIGJp4bq/biBpbnN1cmFuY2UgKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibikgdsOgIGJp4bq/biBob3NwaXRhbCAoc+G7kSBs4bqnbiBuaOG6rXAgdmnhu4duICkNCg0KR2nhuqMgdGh1eeG6v3QgJEhfMCQgOiBiaeG6v24gaW5zdXJhbmNlIChz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4pIHbDoCBiaeG6v24gIGhvc3BpdGFsIChz4buRIGzhuqduIG5o4bqtcCB2aeG7h24gKSDEkeG7mWMgbOG6rXAuDQpHaeG6oyB0aHV54bq/dCAkSF8xJDogYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uICBob3NwaXRhbCAoc+G7kSBs4bqnbiBuaOG6rXAgdmnhu4duICkga2jDtG5nIMSR4buZYyBs4bqtcC4NCg0KYGBge3J9DQp0NSA8LSB0YWJsZSh0JGhvc3BpdGFsMSwgdCRpbnN1cmFuY2UpDQpjaGlzcS50ZXN0KHQ1KQ0KYGBgDQoNCg0KVGEgY8OzIHBfdmFsdWUgPSAwLjkyMjcgPiAkXGFscGhhJC4gQsOhYyBi4buPIGdp4bqjIHRodXnhur90ICRIXzEkLCBjaOG6pXAgbmjhuq1uIGdp4bqjIHRodXnhur90ICRIXzAkLCBiaeG6v24gaW5zdXJhbmNlIChz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4pIHbDoCBiaeG6v24gIGhvc3BpdGFsIChz4buRIGzhuqduIG5o4bqtcCB2aeG7h24gKSDEkeG7mWMgbOG6rXAgduG7m2kgbmhhdQ0KDQoNCiMjIyA0LjEuNS4gS2nhu4NtIHRyYSB0w61uaCDEkeG7mWMgbOG6rXAgY2hvIGJp4bq/biBpbnN1cmFuY2UgKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibikgdsOgIGJp4bq/biBjaHJvbmljIChz4buRIGLhu4duaCBtw6NuIHTDrW5oIG3DoCBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBt4bqvYyBwaOG6o2kpDQoNCkdp4bqjIHRodXnhur90ICRIXzAkIDogYmnhur9uIGluc3VyYW5jZSAoc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuKSB2w6AgYmnhur9uIGNocm9uaWMgKHPhu5EgYuG7h25oIG3Do24gdMOtbmggbcOgIG5nxrDhu51pIMSRxrDhu6NjIGto4bqjbyBzw6F0IG3huq9jIHBo4bqjaSkgxJHhu5ljIGzhuq1wLg0KR2nhuqMgdGh1eeG6v3QgJEhfMSQ6IGJp4bq/biBpbnN1cmFuY2UgKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibikgdsOgIGJp4bq/biBjaHJvbmljIChz4buRIGLhu4duaCBtw6NuIHTDrW5oIG3DoCBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBt4bqvYyBwaOG6o2kpIGtow7RuZyDEkeG7mWMgbOG6rXAuDQoNCmBgYHtyfQ0KDQoNCnRoNSA8LSB0YWJsZSh0JGNocm9uaWMxLHQkaW5zdXJhbmNlKQ0KY2hpc3EudGVzdCh0aDUpDQoNCmBgYA0KDQpUYSBjw7MgcF92YWx1ZSA9IDAuNzU0OSA+ICRcYWxwaGEkLiBCw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3QgJEhfMSQsIGNo4bqlcCBuaOG6rW4gZ2nhuqMgdGh1eeG6v3QgJEhfMCQsIGJp4bq/biBpbnN1cmFuY2UgKHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibikgdsOgIGJp4bq/biBjaHJvbmljIChz4buRIGLhu4duaCBtw6NuIHTDrW5oIG3DoCBuZ8aw4budaSDEkcaw4bujYyBraOG6o28gc8OhdCBt4bqvYyBwaOG6o2kpIMSR4buZYyBs4bqtcCB24bubaSBuaGF1Lg0KDQojIyA0LjIuIEtob+G6o25nIMaw4bubYyBsxrDhu6NuZyB04bu3IGzhu4cgDQoNCiMjIyA0LjIuMS4gxq/hu5tjIGzGsOG7o25nIHThu4kgbOG7hyBuZ8aw4budaSBjw7MgdsOgIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIGPDoSBuaMOibg0KDQpgYGB7cn0NCnAgPC0gdFt0JGluc3VyYW5jZSA9PSAnY29CSENOJyxdIA0KcHJvcC50ZXN0KGxlbmd0aChwJGluc3VyYW5jZSksIGxlbmd0aCh0JGluc3VyYW5jZSkpIA0KdDEwIDwtIHRbdCRpbnN1cmFuY2UgPT0gJ2tob25nQkhDTicsXSANCnByb3AudGVzdChsZW5ndGgodDEwJGluc3VyYW5jZSksIGxlbmd0aCh0JGluc3VyYW5jZSkpIA0KYGBgDQoNClbhu5tpIMSR4buZIHRpbiBj4bqteSA5NSUsIHRhIGPDsyB04bu3IGzhu4cgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSB0xrAgIG5ow6JuIHNvIHbhu5tpIHThu5VuZyB0aOG7gyBu4bqxbSB0cm9uZyBraG/huqNuZyB04burIDc2LjM4JSB2w6AgNzguODYlLiBIYXkgbsOzaSBjw6FjaCBraMOhYywgdOG7iSBs4buHIG5nxrDhu51pIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBuaMOibiBz4bq9IGNoaeG6v20ga2hv4bqjbmcgdOG7qyAyMS4xMSUgxJHhur9uIDIzLjYyJQ0KDQojIyMgNC4yLjIuIMav4bubYyBsxrDhu6NuZyB04buJIGzhu4cgbmfGsOG7nWkgY8OzIHbDoCBraMO0bmcgY8OzIHZp4buHYyBsw6BtDQoNCmBgYHtyfQ0KDQplIDwtIHRbdCRlbXBsb3llZCA9PSAnY292aWVjbGFtJyxdIA0KcHJvcC50ZXN0KGxlbmd0aChlJGVtcGxveWVkID09ICdjb3ZpZWNsYW0nKSwgbGVuZ3RoKHQkZW1wbG95ZWQpKSANCmUxIDwtIHRbdCRpbnN1cmFuY2UgPT0gJ2tob25ndmllY2xhbScsXSANCnByb3AudGVzdChsZW5ndGgoZTEkZW1wbG95ZWQgPT0gJ2tob25ndmllY2xhbScpLCBsZW5ndGgodCRlbXBsb3llZCkpIA0KYGBgDQoNClbhu5tpIMSR4buZIHRpbiBj4bqteSA5NSUsIHRhIGPDsyB04bu3IGzhu4cgbmfGsOG7nWkgY8OzIHZp4buHYyBsw6BtIHNvIHbhu5tpIHThu5VuZyB0aOG7gyBu4bqxbSB0cm9uZyBraG/huqNuZyB04burIDkuNDUlIMSR4bq/biAxMS4yNyUuIEhheSBuw7NpIGPDoWNoIGtow6FjLCB04buJIGzhu4cgbmfGsOG7nWkga2jDtG5nIGPDsyB2aeG7h2MgbMOgbSBz4bq9IGNoaeG6v20ga2hv4bqjbmcgdOG7qyAwJSDEkeG6v24gMC4xMDg2MDI4JS4NCg0KDQoNCiMjIyA0LjIuMy4gxq/hu5tjIGzGsOG7o25nIHThu7cgbOG7hyBuZ8aw4budaSBjw7MgdsOgIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIG1lZGljYWlkIA0KDQpgYGB7cn0NCg0KZTMgPC0gdFt0JG1lZGljYWlkID09ICdjb0JITWVkaScsXSANCnByb3AudGVzdChsZW5ndGgoZTMkbWVkaWNhaWQgPT0gJ2NvQkhNZWRpJyksIGxlbmd0aCh0JG1lZGljYWlkKSkNCg0KZTQgPC0gdFt0JG1lZGljYWlkID09ICdraG9uZ0JITWVkaScsXSANCnByb3AudGVzdChsZW5ndGgoZTQkbWVkaWNhaWQgPT0gJ2tob25nQkhNZWRpJyksIGxlbmd0aCh0JG1lZGljYWlkKSkgDQpgYGANCg0KVuG7m2kgxJHhu5kgdGluIGPhuq15IDk1JSwgdGEgY8OzIHThu7cgbOG7hyBuZ8aw4budaSBjw7MgYuG6o28gaGnhu4NtIG1lZGljYWlkIHbhu5tpIHThu5VuZyB0aOG7gyBu4bqxbSB0cm9uZyBraG/huqNuZyB04burIDguMyUgdsOgIDEwLjAyJS4gSGF5IG7Ds2kgY8OhY2gga2jDoWMsIHThu4kgbOG7hyBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBtZWRpY2FpZCBz4bq9IGNoaeG6v20ga2hv4bqjbmcgdOG7qyA4OS45NyUgxJHhur9uIDkxLjcwJQ0KDQojIyMgNC4yLjQuIMav4bubYyBsxrDhu6NuZyB04bu3IGzhu4cgbmfGsOG7nWkgY8OzIHRodSBuaOG6rXAgZMaw4bubaSAyLjUgbmfDoG4gxJHDtA0KYGBge3J9DQppIDwtIHRbdCRpbmNvbWUgPT0gMi41IHwgdCRpbmNvbWUgPCAyLjUsXQ0KcHJvcC50ZXN0KGxlbmd0aChpJGluY29tZSA9PSAzMCB8IGkkaW5jb21lIDwzMCksIGxlbmd0aCh0JGluY29tZSkpDQoNCmBgYA0KDQpW4bubaSDEkeG7mSB0aW4gY+G6rXkgOTUlLCB0YSBjw7MgdOG7tyBs4buHIG5nxrDhu51pIGPDsyB0aHUgbmjhuq1wIGTGsOG7m2kgMi41IG5nw6BuIMSRw7QgY2hp4bq/bSB04burIDY1LjEzJSDEkeG6v24gNjcuOTQlLg0KDQoNCiMjIyA0LjIuNS4gxq/hu5tjIGzGsOG7o25nIHThu7cgbOG7hyBuZ8aw4budaSBraMO0bmcgY8OzIGLhu4duaCBtw6NuIHTDrW5oDQoNCmBgYHtyfQ0KaTMgPC0gdFt0JGNocm9uaWMgPT0gMSB8IHQkY2hyb25pYyA8IDEsXQ0KcHJvcC50ZXN0KGxlbmd0aChpMyRjaHJvbmljID09IDEgfCBpMyRjaHJvbmljIDwxKSwgbGVuZ3RoKCB0JGNocm9uaWMpKQ0KDQpgYGANClbhu5tpIMSR4buZIHRpbiBj4bqteSA5NSUsIHRhIGPDsyB04bu3IGzhu4cgbmfGsOG7nWkgY8OzIGTGsOG7m2kgbeG7mXQgYuG7h25oIG3Do24gdMOtbmggbMOgIDU1Ljc5JSDEkeG6v24gNTguNzMlDQoNCg0KIyMjIDQuMi42LiDhu6jhu5tjIGzGsOG7o25nIG5nxrDhu51pIGPDsyBkxrDhu5tpIDIgbOG6p24gbmjhuq1wIHZp4buHbg0KDQpgYGB7cn0NCmk4IDwtIHRbdCRob3NwaXRhbCA9PSAyIHwgdCRob3NwaXRhbCA8IDIsXQ0KcHJvcC50ZXN0KGxlbmd0aChpOCRob3NwaXRhbCA9PSAyIHwgaTgkaG9zcGl0YWwgPDIpLCBsZW5ndGgoIHQkaG9zcGl0YWwpKQ0KDQpgYGANClbhu5tpIMSR4buZIHRpbiBj4bqteSA5NSUsIHRhIGPDsyB04bu3IGzhu4cgbmfGsOG7nWkgY8OzIGTGsOG7m2kgMiBs4bqnbiBuaOG6rXAgdmnhu4duIGzDoCA5Ny40OCUgxJHhur9uIDk4LjM1JQ0KDQoNCiMgQ2jGsMahbmcgNS4gSOG7k2kgcXV5IA0KDQojIyA1LjEuIE3DtCBow6xuaCB4w6FjIHN14bqldCB0dXnhur9uIHTDrW5oIA0KDQojIyMgNS4xLjEuIEJp4bq/biBpbnN1cmFuY2UgdGhlbyBpbmNvbWUNCg0KYGBge3J9DQp0IDwtIHQgJT4lIG11dGF0ZShpbnN1cmFuY2UxID0gaWZfZWxzZShpbnN1cmFuY2UgPT0gImtob25nQkhDTiIsIDAsIDEpKQ0KbGV2ZWxzKGZhY3Rvcih0JGluc3VyYW5jZTEpKSAgIA0KbG0uaW5zYyA8LSBnbG0oaW5zdXJhbmNlMSB+IGluY29tZSwgZGF0YSA9IHQpDQpzdW1tYXJ5KGxtLmluc2MpDQpgYGANCk3DtCBow6xuaCB4w6FjIHN14bqldCB0dXnhur9uIHTDrW5oIGPhu6dhIGJp4bq/biBpbnN1cmFuY2UgdsOgIGluY29tZWQgY8OzIGThuqFuZzoNCg0KJFxoYXRccGkgPSAwLjcxOTM1MSArIDAuMDIyNTkxIFx0aW1lcyBpbmNvbWUkDQoNCkPDsyB0aOG7gyB0aOG6pXksIGPhuqMgaGFpIGjhu4cgc+G7kSB0cm9uZyBtw7QgaMOsbmggxJHhu4F1IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogcuG6pXQgY2FvIChwLXZhbHVlIDwgMC4wMSksIGNobyB0aOG6pXkgdGh1IG5o4bqtcCDhuqNuaCBoxrDhu59uZyDEkeG6v24gdmnhu4djIG11YSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGhheSBraMO0bmcuDQoNCi0gSOG7hyBz4buRIGNo4bq3biBsw6AgMC43MTkzNTEgY2hvIGJp4bq/dCBu4bq/dSBuaMawIGtow7RuZyBjw7Mgc+G7sSDhuqNuaCBoxrDhu59uZyBj4bunYSB0aHUgbmjhuq1wIHRow6wgeMOhYyBzdeG6pXQgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGtow6EgdGjhuqVwICgwLjI3MTE0NikuDQoNCi0gSOG7hyBz4buRIGPhu6dhIGJp4buDdSBoaeG7h24gdGh1IG5o4bqtcCBjYW8gbMOgIDAuMDIyNTkxLiBDaG8gYmnhur90IHRodSBuaOG6rXAgdMSDbmcgMSDEkcahbiB24buLIHRow6wga2jhuqMgbsSDbmcgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIHTEg25nIDAuMDIyNTkxIGzhuqduLg0KDQojIyMgNS4xLjIuIEJp4bq/biBpbnN1cmFuY2UgdGhlbyBlbXBsb3llZA0KDQpgYGB7cn0NCnQkZW1wbG95ZWQgPC0gZmFjdG9yICh0JGVtcGxveWVkKQ0KbG0uaW5zZW0gPC0gZ2xtKGluc3VyYW5jZTEgfiBlbXBsb3llZCwgZGF0YSA9IHQpDQpzdW1tYXJ5KGxtLmluc2VtKQ0KYGBgDQoNCk3DtCBow6xuaCB4w6FjIHN14bqldCB0dXnhur9uIHTDrW5oIGPhu6dhIGJp4bq/biBpbnN1cmFuY2UgdsOgIGVtcGxveWVkIGPDsyBk4bqhbmc6DQoNCiRcaGF0XHBpID0gMC44NTkzNCAtIDAuMDkyNDVcdGltZXMgZW1wbG95ZWRraG9uZ3ZpZWNsYW0kDQoNCkPDsyB0aOG7gyB0aOG6pXksIGPhuqMgaGFpIGjhu4cgc+G7kSB0cm9uZyBtw7QgaMOsbmggxJHhu4F1IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogcuG6pXQgY2FvIChwLXZhbHVlIDwgMC4wMSksIGNobyB0aOG6pXkgdMOsbmggdHLhuqFuZyB2aeG7h2MgbMOgbSBjw7Mg4bqjbmggaMaw4bufbmcgxJHhur9uIHZp4buHYyBtdWEgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBoYXkga2jDtG5nLiBD4bulIHRo4buDOg0KDQotIEjhu4cgc+G7kSBjaOG6t24gbMOgIDAuODU5MzQsIGNobyBiaeG6v3QgbuG6v3UgbmjGsCBraMO0bmcgY8OzIHPhu7Eg4bqjbmggaMaw4bufbmcgY+G7p2EgdMOsbmggdHLhuqFuZyB2aeG7h2MgbMOgbSB0aMOsIHjDoWMgc3XhuqV0IHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBsw6AgMC44NTkzNA0KDQotIEjhu4cgc+G7kSBj4bunYSB0w6xuaCB0cuG6oW5nIGtow7RuZyBjw7Mgdmnhu4djIGzDoG0gbMOgIC0gMC4wOTI0NSwgY2hvIGJp4bq/dCB4w6FjIHN14bqldCDEkeG7gyBz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY+G7p2EgbmfGsOG7nWkga2jDtG5nIGPDsyB2aeG7h2MgbMOgbSB0aOG6pXAgaMahbiAwLjA5MjQ1IGzhuqduIHNvIHbhu5tpIG5nxrDhu51pIGPDsyB2aeG7h2MgbMOgbS4NCg0KIyMjIDUuMS4zLiBCaeG6v24gaW5zdXJhbmNlIHRoZW8gbWVkaWNhaWQNCg0KYGBge3J9DQp0JG1lZGljYWlkIDwtIGZhY3RvciAodCRtZWRpY2FpZCkNCg0KbG0uaW5zbSA8LSBnbG0oaW5zdXJhbmNlMSB+IG1lZGljYWlkLCBkYXRhID0gdCkNCnN1bW1hcnkobG0uaW5zbSkNCmBgYA0KDQpNw7QgaMOsbmggeMOhYyBzdeG6pXQgdHV54bq/biB0w61uaCBj4bunYSBiaeG6v24gaW5zdXJhbmNlIHbDoCBtZWRpY2FpZCBjw7MgZOG6oW5nOg0KDQokXGhhdFxwaSA9ICAwLjE1MTc0ICsgMC42ODc0MiBcdGltZXMga2hvbmdCSE1lZGkkDQoNCkPDsyB0aOG7gyB0aOG6pXksIGPhuqMgaGFpIGjhu4cgc+G7kSB0cm9uZyBtw7QgaMOsbmggxJHhu4F1IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogcuG6pXQgY2FvIChwLXZhbHVlIDwgJFxhbHBoYSQgPSAwLjAwMSksIGNobyB0aOG6pXkgdmnhu4djIGPDsyBz4bufIGjhu691IGLhuqNvIGhp4buDbSBNZWRpY2FpZCBjw7Mg4bqjbmggaMaw4bufbmcgxJHhur9uIHZp4buHYyBtdWEgYuG6o28gaGnhu4NtIGPDoSBuaMOibi4gQ+G7pSB0aOG7gzoNCg0KLSBI4buHIHPhu5EgY2jhurduIGzDoCAwLjE1MTc0LCBjaG8gYmnhur90IG7hur91IG5oxrAga2jDtG5nIGPDsyBz4buxIOG6o25oIGjGsOG7n25nIGPhu6dhIGLhuqNvIGhp4buDbSBNZWRpY2FpZCB0aMOsIHjDoWMgc3XhuqV0IHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBsw6Aga2jDoSB0aOG6pXAsIGtob+G6o25nIDAuMTUxNzQgDQoNCi0gSOG7hyBz4buRIGPhu6dhIGJp4buDdSBoaeG7h24ga2jDtG5nIGPDsyBCSCBNZWRpY2FpZCBsw6AgMC42ODc0MiwgY2hvIGJp4bq/dCBraGkga2jDtG5nIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPhu6dhIE1lZGljYWlkIHRow6wgeMOhYyBzdeG6pXQgxJHhu4Mgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGzDoCBraMOhIGNhbywga2hv4bqjbmcgMC42ODc0Mg0KDQojIyMgNS4xLjQuIEjhu5NpIHF1eSDEkWEgYmnhur9uDQoNCmBgYHtyfQ0KbG0uaW4gPC0gZ2xtKGluc3VyYW5jZTEgfiBpbmNvbWUgKyBlbXBsb3llZCArIG1lZGljYWlkLCBkYXRhID0gdCkNCnN1bW1hcnkobG0uaW4pDQpgYGANCiRcaGF0XHBpID0gIDAuMTUzMjk1ICsgMC4wMTI0MTdcdGltZXMgaW5jb21lICsgMC42NjY5MTMgXHRpbWVzIG1lZGljYWlka2hvbmdCSE1lZGkkDQoNCkJp4bq/biBlbXBsb3llZCBraMO0bmcgY8OzIMO9IG5naMSpYSBkbyBwLXZhbHVlID0gMC4zMTUgbOG7m24gaMahbiBt4bupYyDDvSBuZ2jEqWEuDQoNCkjhu4cgc+G7kSBjaOG6t24gMC4xNTMyOTUgY8OzIG5naMSpYSBsw6AgbuG6v3UgbmjGsCBraMO0bmcgY8OzIHTDoWMgxJHhu5luZyBj4bunYSBjw6FjIHnhur91IHThu5EgdGh1IG5o4bqtcCwgYuG6o28gaGnhu4NtIE1lZGljYWlkIHRow6wgeMOhYyBzdeG6pXQgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGzDoCAwLjE1MzI5NQ0KDQpI4buHIHPhu5EgY+G7p2EgYmnhur9uIHRodSBuaOG6rXAgbMOgIDAuMDEyNDE3LCBjaG8gYmnhur90IGtoaSB0aHUgbmjhuq1wIHTEg25nIDEgxJHGoW4gduG7iyB0aMOsIHjDoWMgc3XhuqV0IHPhu58gaOG7r3UgQkhDTiB0xINuZyAwLjAxMjQxNy4gDQoNCkjhu4cgc+G7kSBj4bunYSBiaeG7g3UgaGnhu4duIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkIGzDoCAwLjY2NjkxMywgY2hvIGJp4bq/dCBuZ8aw4budaSBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBNZWRpY2FpZCBjw7MgeMOhYyBzdeG6pXQgbXVhIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gY2FvIGjGoW4gMC42NjY5MTMgc28gduG7m2kgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSBNZWRpY2FpZC4NCg0KIyMgNS4yLiBNw7QgaMOsbmggaOG7k2kgcXV5IGxvZ2l0DQoNCiMjIyA1LjIuMS4gQmnhur9uIGluc3VyYW5jZSB2w6AgaW5jb21lZA0KDQpgYGB7cn0NCnQkaW5zdXJhbmNlIDwtIGZhY3Rvcih0JGluc3VyYW5jZSwgbGV2ZWxzID0gYygna2hvbmdCSENOJywnY29CSENOJykpDQpsb2dpdC5pbnNjIDwtIGdsbShpbnN1cmFuY2UgfiBpbmNvbWUsIGZhbWlseSA9IGJpbm9taWFsKGxpbmsgPSAnbG9naXQnKSwgZGF0YSA9IHQpDQpzdW1tYXJ5KGxvZ2l0Lmluc2MpDQpgYGANCg0KTcO0IGjDrG5oIGjhu5NpIHF1eSBsb2dpdCBj4bunYSBpbnN1cmFuY2UgdsOgIGluY29tZWQgbMOgOg0KDQokbG4oXGZyYWN7XGhhdFxwaX17MS1caGF0XHBpfSk9IDAuNjM2OTkgKyAwLjI5MzI1XHRpbWVzIGluY29tZSQNCg0KSOG7hyBz4buRIGNo4bq3biBsw6AgMC42MzY5OSwgY8OzIG5naMSpYSBsw6AgJGxuKFxmcmFje1xoYXRccGl9ezEtXGhhdFxwaX0pJCBraGkga2jDtG5nIGPDsyB0w6FjIMSR4buZbmcgY+G7p2EgdGh1IG5o4bqtcCB0aMOsIHjDoWMgc3XhuqV0IHPhu58gaOG7r3UgQkhDTiBsw6AgMC42MzY5OQ0KDQpI4buHIHPhu5EgY+G7p2EgYmnhur9uIHRodSBuaOG6rXAgbMOgIDAuMjkzMjUgLiBI4buHIHPhu5EgbsOgeSBuZ2jEqWEgbMOgLCDEkeG7kWkgduG7m2kgdmnhu4djIG11YSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIHRow6wgJGxuKFxmcmFje1xoYXRccGl9ezEtXGhhdFxwaX0pJCBj4bunYSB0aHUgbmjhuq1wIHTEg25nIHRow6wga2jhuqMgbsSDbmcgc+G7nyBo4buvdSBCSENOIHTEg25nLg0KDQojIyMgNS4yLjIuIEJp4bq/biBpbnN1cmFuY2UgdsOgIGVtcGxveWVkDQoNCmBgYHtyfQ0KbG9naXQuaW5zZW0gPC0gZ2xtKGluc3VyYW5jZSB+IGVtcGxveWVkLCBmYW1pbHkgPSBiaW5vbWlhbChsaW5rID0gJ2xvZ2l0JyksIGRhdGEgPSB0KQ0Kc3VtbWFyeShsb2dpdC5pbnNlbSkNCmBgYA0KDQpNw7QgaMOsbmggaOG7k2kgcXV5IGxvZ2l0IGPhu6dhIGluc3VyYW5jZSB2w6AgZW1wbG95IGzDoDoNCg0KJGxuKFxmcmFje1xoYXRccGl9ezEtXGhhdFxwaX0pPSAxLjgwOTggLSAwLjYxOTBcdGltZXMgZW1wbG95ZWRraG9uZ3ZpZWNsYW0kDQoNCg0KSOG7hyBz4buRIGNo4bq3biBsw6AgMS44MDk4LiBI4buHIHPhu5EgbsOgeSBjw7MgbmdoxKlhIGzDoCBu4bq/dSBraMO0bmcgY8OzIHPhu7Eg4bqjbmggaMaw4bufbmcgY+G7p2EgYmnhu4N1IGhp4buHbiBraMOhYyB0aMOsICRsbihcZnJhY3tcaGF0XHBpfXsxLVxoYXRccGl9KSQgY+G7p2EgbmfGsOG7nWkgY8OzIHZp4buHYyBsw6BtLCBz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4gbMOgIDEuODA5OA0KDQpI4buHIHPhu5EgY+G7p2EgYmnhu4N1IGhp4buHbiBraMO0bmcgY8OzIHZp4buHYyBsw6BtIGLhurFuZyAtIDAuNjE5MC4gSOG7hyBz4buRIG7DoHkgbmdoxKlhIGzDoCwgxJHhu5FpIHbhu5tpIHZp4buHYyBtdWEgYuG6o28gaGnhu4NtIGPDoSBuaMOibiB0aMOsIGNow6puaCBs4buHY2ggJGxuKFxmcmFje1xoYXRccGl9ezEtXGhhdFxwaX0pJCBj4bunYSBuZ8aw4budaSBjw7Mga2jDtG5nIGPDsyB2aeG7h2MgbMOgbSB2w6AgbmfGsOG7nWkgY8OzIHZp4buHYyBsw6BtIGzDoCAtMC42MTkwDQoNCg0KIyMjIDUuMi4zLiBCaeG6v24gaW5zdXJhbmNlIHbDoCBtZWRpY2FpZA0KDQpgYGB7cn0NCmxvZ2l0Lmluc20gPC0gZ2xtKGluc3VyYW5jZSB+IG1lZGljYWlkICwgZmFtaWx5ID0gYmlub21pYWwobGluayA9ICdsb2dpdCcpLCBkYXRhID0gdCkNCnN1bW1hcnkobG9naXQuaW5zbSkNCmBgYA0KDQpNw7QgaMOsbmggaOG7k2kgcXV5IGxvZ2l0IGPhu6dhIGluc3VyYW5jZSB2w6AgbWVkaWNhaWQgbMOgOg0KDQokbG4oXGZyYWN7XGhhdFxwaX17MS1caGF0XHBpfSk9IC0xLjcyMTAgKyAzLjM3MzBcdGltZXMga2hvbmdCSE1lZGkkDQoNCkjhu4cgc+G7kSBjaOG6t24gbMOgIC0gMS43MjEwIHbDoCBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqLiBI4buHIHPhu5EgbsOgeSBjaG8gYmnhur90IG7hur91IGtow7RuZyBjw7Mgc+G7sSDhuqNuaCBoxrDhu59uZyBj4bunYSBiaeG7g3UgaGnhu4duIGtow6FjIHRow6wgJGxuKFxmcmFje1xoYXRccGl9ezEtXGhhdFxwaX0pJCBj4bunYSBuZ8aw4budaSB24burYSBjw7MgYuG6o28gaGnhu4NtIE1lZGksIHbhu6thIGPDsyBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGzDoCAtMS43MjEwDQoNCkjhu4cgc+G7kSBj4bunYSBiaeG7g3UgaGnhu4duIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkIGLhurFuZyAzLjM3MzAuIEjhu4cgc+G7kSBuw6B5IG5naMSpYSBsw6AsIMSR4buRaSB24bubaSB2aeG7h2MgbXVhIGLhuqNvIGhp4buDbSBjw6EgbmjDom4gdGjDrCBjaMOqbmggbOG7h2NoICRsbihcZnJhY3tcaGF0XHBpfXsxLVxoYXRccGl9KSQgY+G7p2EgbmfGsOG7nWkga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gTWVkaWNhaWQgc28gduG7m2kgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSBNZWRpY2FpZCBsw6AgMy4zNzMwDQoNCiMjIyA1LjIuNC4gSOG7k2kgcXV5IMSRYSBiaeG6v24NCg0KYGBge3J9DQpsb2dpdC5pbiA8LSBnbG0oaW5zdXJhbmNlIH4gaW5jb21lICsgZW1wbG95ZWQgKyBtZWRpY2FpZCAsIGZhbWlseSA9IGJpbm9taWFsKGxpbmsgPSAnbG9naXQnKSwgZGF0YSA9IHQpDQpzdW1tYXJ5KGxvZ2l0LmluKQ0KYGBgDQoNCiRsbihcZnJhY3tcaGF0XHBpfXsxLVxoYXRccGl9KT0gLTEuODQ3NzggKyAwLjE3MDA0IFx0aW1lcyBpbmNvbWUgKyAzLjE4NjgwXHRpbWVzIG1lZGljYWlka2hvbmdCSE1lZGkkDQoNCkjhu4cgc+G7kSBj4bunYSBiaeG6v24gdGh1IG5o4bqtcCBsw6AgMC4xNzAwNCBjaG8gdGjhuqV5IGtoaSB0aHUgbmjhuq1wIGPDoG5nIGNhbyB0aMOsIHjDoWMgc3XhuqV0IHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBjw6BuZyBjYW8uDQoNCkjhu4cgc+G7kSBj4bunYSBiaeG7g3UgaGnhu4duIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkIGzDoCAzLjE4NjgwIGNobyB0aOG6pXkgeMOhYyBzdeG6pXQgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGPhu6dhIG5nxrDhu51pIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkIGNow6puaCBs4buHY2ggJGxuKFxmcmFje1xoYXRccGl9ezEtXGhhdFxwaX0pJCA9IDMuMTg2ODAgc28gduG7m2kgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSBNZWlkY2FpZC4NCg0KIyMgNS4zLiBNw7QgaMOsbmggcHJvYml0DQoNCiMjIyA1LjMuMS4gQmnhur9uIGluc3VyYW5jZSB2w6AgaW5jb21lZA0KDQpgYGB7cn0NCg0KcHJvYml0Lmluc2MgPC0gZ2xtKGluc3VyYW5jZSB+IGluY29tZSwgZmFtaWx5ID0gYmlub21pYWwobGluayA9ICdwcm9iaXQnKSwgZGF0YSA9IHQpDQpzdW1tYXJ5KHByb2JpdC5pbnNjKQ0KYGBgDQoNCk3DtCBow6xuaCBo4buTaSBxdXkgbmjGsCBzYXU6DQoNCiRcaGF0XHBpID0gXFBoaSgwLjQ1NTgwICsgMC4xMzk4M1x0aW1lcyBpbmNvbWUpJA0KDQpI4buHIHPhu5EgY2jhurduIGzDoCAwLjQ1NTgwIHbDoCBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqLiBUcm9uZyB0csaw4budbmcgaOG7o3AgdGh1IG5o4bqtcCBi4bqxbmcgMCwgdGjDrCB4w6FjIHN14bqldCBz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4gbMOgICRcUGhpKDAuNDU1ODApJA0KDQpI4buHIHPhu5EgY+G7p2EgYmnhur9uIGluY29tZSBsw6AgMC4xMzk4MyB2w6AgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqi4ga2hpIHRodSBuaOG6rXAgdMSDbmcgdGjDrCB4w6FjIHN14bqldCBz4bufIGjhu691IGJhbyBoaeG7g20gY8OhIG5ow6JuIHTEg25nICRcUGhpKDAuMTM5ODMpJC4NCg0KIyMjIDUuMy4yLiBCaeG6v24gaW5zdXJhbmNlIHbDoCBlbXBsb3llZA0KDQpgYGB7cn0NCg0KcHJvYml0Lmluc2VtIDwtIGdsbShpbnN1cmFuY2UgfiBlbXBsb3llZCwgZmFtaWx5ID0gYmlub21pYWwobGluayA9ICdwcm9iaXQnKSwgZGF0YSA9IHQpDQpzdW1tYXJ5KHByb2JpdC5pbnNlbSkNCmBgYA0KDQpNw7QgaMOsbmggaOG7k2kgcXV5IG5oxrAgc2F1Og0KJFxoYXRccGkgPSBcUGhpKDEuMDc3MzYgLSAwLjM0ODcwIFx0aW1lcyBlbXBsb3llZGtob25ndmllY2xhbSkkDQoNCkjhu4cgc+G7kSBjaOG6t24gbMOgIDEuMDc3MzYsIGtoaSBjw6FjIHnhur91IHThu5Ega2jDoWMga2jDtG5nIMSR4buVaSB0aMOsIHjDoWMgc3XhuqV0IGPhu6dhIG3hu5l0IG5nxrDhu51pIGPDsyB2aeG7h2MgbMOgbSB2w6AgY8OzIG11YSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGzDoCAkXFBoaSgxLjA3NzM2KSQNCg0KxJDhu5FpIHbhu5tpIG5nxrDhu51pIGPDsyBraMO0bmcgY8OzIHZp4buHYyBsw6BtIHRow6wgeMOhYyBzdeG6pXQgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGzDoCAkXFBoaSgtMC4zNDg3MCkkDQoNCiMjIyA1LjMuMy4gQmnhur9uIGluc3VyYW5jZSB2w6AgbWVkaWNhaWQNCg0KYGBge3J9DQpwcm9iaXQuaW5zbSA8LSBnbG0oaW5zdXJhbmNlIH4gbWVkaWNhaWQsIGZhbWlseSA9IGJpbm9taWFsKGxpbmsgPSAncHJvYml0JyksIGRhdGEgPSB0KQ0Kc3VtbWFyeShwcm9iaXQuaW5zbSkNCmBgYA0KDQpNw7QgaMOsbmggaOG7k2kgcXV5IG5oxrAgc2F1Og0KDQokXGhhdFxwaSA9IFxQaGkoLTEuMDI4OTkgKyAyLjAyMDAxIFx0aW1lcyBraG9uZ0JITWVkaSkkDQoNCkjhu4cgc+G7kSBjaOG6t24gbMOgIC0gMS4wMjg5OSB2w6AgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiAocC12YWx1ZSA8ICRcYWxwaGEkKS4gSOG7hyBz4buRIG7DoHkgY8OzIG5naMSpYSBsw6AgduG7m2kgbmfGsOG7nWkgY8OzIGLhuqNvIGhp4buDbSBNZWRpY2FpZCwgeMOhYyBzdeG6pXQgxJHhu4MgbmfGsOG7nWkgxJHDsyBtdWEgdGjDqm0gYuG6o28gaGnhu4NtIGPDoSBuaMOibiBsw6AgJFxQaGkoLTEuMDI4OTkpJA0KDQpI4buHIHPhu5EgY+G7p2EgYmnhu4N1IGhp4buHbiBraMO0bmcgY8OzIGLhuqNvIGhp4buDbSBNZWRpY2FpZCBsw6AgMi4wMjAwMSB2w6AgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqi4gSOG7hyBz4buRIGNobyBiaeG6v3QgbuG6v3Uga2jDtG5nIGPDsyBi4bqjbyBoaeG7g20gY+G7p2EgTWVkaWNhaWQgdGjDrCB4w6FjIHN14bqldCDEkeG7gyBuZ8aw4budaSDEkcOzIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBsw6AgJFxQaGkoMi4wMjAwMSkkDQoNCiMjIyA1LjMuNC4gSOG7k2kgcXV5IMSRYSBiaeG6v24NCg0KYGBge3J9DQpwcm9iaXQuaW4gPC0gZ2xtKGluc3VyYW5jZSB+IGluY29tZSArIGVtcGxveWVkICsgbWVkaWNhaWQsIGZhbWlseSA9IGJpbm9taWFsKGxpbmsgPSAncHJvYml0JyksIGRhdGEgPSB0KQ0Kc3VtbWFyeShwcm9iaXQuaW4pDQpgYGANCg0KJFxoYXRccGkgPSBcUGhpKC0xLjA3NjYxICsgMC4wODQ2NiBcdGltZXMgaW5jb21lIC0gMS45MTQzNCBcdGltZXMgbWVkaWNhaWRraG9uZ0JITWVkaSkkDQoNCkjhu4cgc+G7kSBjaOG6t24gbMOgICRcUGhpKC0xLjA3NjYxKSQgdsOgIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ouIMSQaeG7gXUgbsOgeSBjw7MgbmdoxKlhIGzDoCBu4bq/dSBraMO0bmcgY8OzIHnhur91IHThu5EgdMOhYyDEkeG7mW5nIHRow6wgeMOhYyBzdeG6pXQgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGzDoCAkXFBoaSgtMS4wNzY2MSkkDQoNCljDoWMgc3XhuqV0IHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBj4bunYSBuZ8aw4budaSBjw7MgdGh1IG5o4bqtcCBjw6BuZyBjYW8gbMOgICRcUGhpKDAuMDg0NjYpJC4NCg0KWMOhYyBzdeG6pXQgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGPhu6dhIG5nxrDhu51pIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkIGzDoCAkXFBoaSgtMS45MTQzNCkkLg0KDQojIyA1LjQuIMSQw6FuaCBnacOhIG3DtCBow6xuaA0KDQpgYGB7cn0NCkFJQyhsbS5pbnNjKQ0KQUlDKGxvZ2l0Lmluc2MpDQpBSUMocHJvYml0Lmluc2MpDQpgYGANCg0KxJDhu5FpIHbhu5tpIG3DtCBow6xuaCBo4buTaSBxdXkgYmnhur9uIGluc3VyYW5jZSB0aGVvIGJp4bq/biBpbmNvbWVkIHRow6wgbcO0IGjDrG5oIGxvZ2l0IGzDoCB04buRdCBuaOG6pXQgdsOsIGNo4buJIHPhu5EgQUlDIGPhu6dhIG3DtCBow6xuaCBuw6B5IG5o4buPIG5o4bqldC4NCg0KYGBge3J9DQpBSUMobG0uaW5zZW0pDQpBSUMobG9naXQuaW5zZW0pDQpBSUMocHJvYml0Lmluc2VtKQ0KYGBgDQoNClRyb25nIHZp4buHYyBnaeG6o2kgdGjDrWNoIG3hu5FpIHF1YW4gaOG7hyBnaeG7r2EgdMOsbmggdHLhuqFuZyB2aeG7h2MgbMOgbSB2w6Agc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIHRow6wgbcO0IGjDrG5oIGxvZ2l0IHbDoCBwcm9iaXQgY8OzIGhp4buHdSBzdeG6pXQgdOG7kXQgaMahbi4gTG9naXQgdsOgIFByb2JpdCBjw7MgaGnhu4d1IHN14bqldCB0xrDGoW5nIMSRxrDGoW5nIG5oYXUgdHJvbmcgdmnhu4djIGdp4bqjaSB0aMOtY2ggbeG7kWkgcXVhbiBo4buHIGdp4buvYSB0w6xuaCB0cuG6oW5nIHZp4buHYyBsw6BtIHbDoCBz4bufIGjhu691IGLhuqNvIGhp4buDbSBjw6EgbmjDom4gdsOsIGNo4buJIHPhu5EgQUlDIGLhurFuZyBuaGF1Lg0KDQpgYGB7cn0NCkFJQyhsbS5pbnNtKQ0KQUlDKGxvZ2l0Lmluc20pDQpBSUMocHJvYml0Lmluc20pDQpgYGANCg0KVHJvbmcgdmnhu4djIGdp4bqjaSB0aMOtY2ggbeG7kWkgcXVhbiBo4buHIGdp4buvYSB2aeG7h2Mgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gTWVkaWNhaWQgdsOgIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiB0aMOsIG3DtCBow6xuaCB4w6FjIHN14bqldCB0dXnhur9uIHTDrW5oIGPDsyBoaeG7h3Ugc3XhuqV0IHThu5F0IGjGoW4uDQoNCmBgYHtyfQ0KQUlDKGxtLmluKQ0KQUlDKGxvZ2l0LmluKQ0KQUlDKHByb2JpdC5pbikNCmBgYA0KTcO0IGjDrG5oIHjDoWMgc3XhuqV0IHR1eeG6v24gdMOtbmggZ2nhuqNpIHRow61jaCBt4buRaSBxdWFuIGjhu4cgZ2nhu69hIGJp4bq/biBwaOG7pSB0aHXhu5ljIGluc3VyYW5jZSB2w6AgY8OhYyBiaeG6v24gxJHhu5ljIGzhuq1wIGluY29tZWQsIGVtcGxveWVkLCBtZWRpY2FpZCB04buRdCBoxqFuIG3DtCBow6xuaCBsb2dpdCwgcHJvYml0Lg0KDQojIENoxrDGoW5nIDYuIEvhur90IGx14bqtbg0KDQpTYXUgcXXDoSB0csOsbmgsIHBow6JuIHTDrW5oIHbDoCBjaOG6oXkgY8OhYyB0w61uaCB0b8OhbiBsacOqbiBxdWFuIGNobyBiaeG6v24gcGjhu6UgdGh14buZYyDigJhpbnN1cmFuY2XigJkg4oCTIGLhuqNvIGhp4buDbSBjw6EgbmjDom4sIHRoZW8gY8OhYyBiaeG6v24gxJHhu5ljIGzhuq1wOiBlbXBsb3llZCwgbWVkaWNhaWQsIGluY29tZSwgaG9zcGl0YWwgdsOgIGNocm9uaWMuIFRhIHRo4bqleSDEkcaw4bujYyBy4bqxbmcsIHRyb25nIDUgYmnhur9uIMSR4buZYyBs4bqtcCBuZ2hpw6puIGPhu6l1IHRow6wgY2jhu4kgY8OzIDMgYmnhur9uIGVtcGxveWVkIOKAkyBjw7MgxJFpIGzDoG0gaGF5IGtow7RuZywgbWVkaWNhaWQg4oCTIMSRxrDhu6NjIHThu5UgY2jhu6ljIEJIIE1lZGljYWlkIGhheSBraMO0bmcsIGluY29tZSDigJMgdGh1IG5o4bqtcCBtYW5nIGzhuqFpIMO9IG5naMSpYSB0aOG7kW5nIGvDqiwgdsOgIGNobyB0aOG6pXkgcuG6sW5nIGNow7puZyB0aOG6rXQgc+G7sSDhuqNuaCBoxrDhu59uZyDEkeG6v24gdmnhu4djIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibi4gDQoNClRoZW8gcXVhbiBzw6F0LCBjaG8gdGjhuqV5IGPDoSBuaMOibiBjw7Mgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIG5oaeG7gXUgaMahbiBjw6EgbmjDom4ga2jDtG5nIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibi4gVHJvbmcgxJHDsywgc+G7kSBuZ8aw4budaSDEkWFuZyBjw7Mgdmnhu4djIGzDoG0gY2hp4bq/bSB04bu3IGzhu4cga2jDoSB0aOG6pXAgc28gduG7m2kgY8OhIG5ow6JuIGtow7RuZyBjw7Mgdmnhu4djIGzDoG0uIMSQaeG7gXUgbsOgeSBraMOhIMSRw7puZyB2w6wgbmfGsOG7nWkgxJHGsOG7o2MgcXVhbiBzw6F0IOG7nyDEkcOieSBsw6AgbmfGsOG7nWkgdHLDqm4gNjYgdHXhu5VpLiBI4buNIMSRxrDhu6NjIGNow61uaCBwaOG7pyBjdW5nIGPhuqVwIG3hu5l0IGxv4bqhaSBi4bqjbyBoaeG7g20ga2jDoWMgbMOgIE1lZGljaWFkLg0KDQrEkG8gbMaw4budbmcgdGh1IG5o4bqtcCBj4bunYSBjw6EgbmjDom4gxJHGsOG7o2Mga2jhuqNvIHPDoXQgcsahaSB2w6BvIGtob+G6o25nIDIuNTI3MSBuZ8OgbiDEkcO0IGzDoCBjaOG7pyB54bq/dS4gS+G6v3QgcXXhuqMgcGjDom4gdMOtY2ggaOG7k2kgcXV5LCBjxaluZyBjaG8gdGjhuqV5IHLhurFuZyB2aeG7h2MgdGh1IG5o4bqtcCB04bu3IGzhu4cgdGh14bqtbiB24bubaSB2aeG7h2MgaOG7jSBjw7MgbXVhIGLhuqNvIGhp4buDbSBjw6EgbmjDom4uICBLaGkgxJHDsyBuZ8aw4budaSBjw7Mgdmnhu4djIGzDoG0sIG1hbmcgbOG6oWkgbmd14buTbiB0aHUgbmjhuq1wIHRo4bqlcCBjw7MgeHUgaMaw4bubbmcgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIGNhbyBoxqFuIGzDoCBuaOG7r25nIG5nxrDhu51pIGPDsyB2aeG7h2MgbMOgbSBuaMawbmcgbOG6oWkgY8OzIHRodSBuaOG6rXAgY2FvLiBNw6Agbmjhu69uZyBjw6EgbmjDom4gxJHDoyBz4bufIGjhu691IGLhuqNvIGhp4buDbSBNZWRpY2FpZCBjw7MgeMOhYyBzdeG6pXQgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gY8OhIG5ow6JuIHRo4bqlcCBoxqFuIG5o4buvbmcgY8OhIG5ow6JuIGtow7RuZyBz4bufIGjhu691IGzhuqFpIGLhuqNvIGhp4buDbSBkbyBjaMOtbmggcGjhu6cgY+G6pXAgKE1lZGljYWlkKS4gS+G6v3QgcXXhuqMgaOG7k2kgcXV5LCBjaG8gdGEgY8OhaSBuaMOsbiBraOG6s25nIMSR4buLbmggcuG6sW5nIHZp4buHYyBz4bufIGjFqXUgYuG6o28gaGnhu4NtIGPDoSBuaMOibiBjw7MgbeG7kWkgcXVhbiBo4buHIGPDuW5nIGNoaeG7gXUgduG7m2kgYmnhu4N1IGhp4buHbiBjw6EgbmjDom4ga2jDtG5nIMSRxrDhu6NjIHThu5UgY2jhu6ljIGN1bmcgY+G6pXAgYuG6o28gaGnhu4NtIE1lZGljaWFkLg0KDQpUdXkgbmhpw6puLCBt4buRaSBxdWFuIGjhu4cgY+G7p2EgY8OhYyBiaeG6v24gbsOgeSBraGkgaOG7k2kgcXV5IMSRYSBiaeG6v24gbOG6oWkgY2hvIHRo4bqleSBiaeG6v24gY8OzIHZp4buHYyBsw6BtIGhheSBraMO0bmcgbOG6oWkga2jDtG5nIGPDsyDDvSBuZ2jEqWEgdGjDtG5nIGvDqi4gQ2jhu4kgY8OzIHnhur91IHThu5EgdGh1IG5o4bqtcCB2w6AgYuG6o28gaGnhu4NtIE1lZGljYWlkIChiaeG7g3UgaGnhu4duIGtow7RuZyBjw7MgYuG6o28gaGnhu4NtIE1lZGljYWlkKSB0w6FjIMSR4buZbmcgY8O5bmcgY2hp4buBdS4gTeG6t2Mga2jDoWMsIG3DtCBow6xuaCBo4buTaSBxdXkgbG9naXQsIHByb2JpdCBjaG8ga+G6v3QgcXXhuqMgY2hvIG3hu5FpIHF1YW4gaOG7hyBnaeG7r2Egdmnhu4djIHPhu6EgaOG7r3UgYuG6o28gaGnhu4NtIHRoZW8gYmnhur9uIHRodSBuaOG6rXAgdOG7kXQgaMahbi4gVHJvbmcgdmnhu4djIGdp4bqjaSB0aMOtY2ggbeG7kWkgcXVhbiBo4buHIGdp4buvYSB2aeG7h2Mgc+G7nyBo4buvdSBi4bqjbyBoaeG7g20gTWVkaWNhaWQgdsOgIHPhu58gaOG7r3UgYuG6o28gaGnhu4NtIGPDoSBuaMOibiB0aMOsIG3DtCBow6xuaCB4w6FjIHN14bqldCB0dXnhur9uIHTDrW5oIGPDsyBoaeG7h3Ugc3XhuqV0IHThu5F0IGjGoW4uDQo=

Phân tích các yếu tố ảnh hưởng đến việc mua bảo hiểm cá nhân

Nhóm 5 - Nguyễn Thị Huỳnh Như, Triệu Trần Thanh Lan, Phạm Thị Thu Thảo

2024-06-14

Chương 1. Tổng quan về đề tài

1.1. Giới thiệu dữ liệu

1.2. Giới thiệu đề tài

Chương 2. Thống kê mô tả 1 biến

2.1. Biến định tính

2.1.1. Đổi tên biểu hiện của các biến định tính

2.1.2. Biến insurance

2.1.3. Biến employed

2.1.4. Biến medicaid

2.2. Biến định lượng

2.2.1. Biến income

2.2.2. Biến hospital

2.2.3. Biến chronic

Chương 3. Thống kê mô tả các cặp biến

3.1. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và thu nhập (income)

3.1.1. Bảng tần suất

3.1.2. Relative Risk

3.1.3. Tỷ lệ chênh - Oddratio

3.2. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và biến employed (có việc làm)

3.2.1. Bảng tần suất

3.2.2. Relative Risk

3.2.3. Tỷ lệ chênh - Oddratio

3.3. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và biến medicaid

3.3.1. Bảng tần suất

3.3.2. Relative Risk

3.3.3. Tỷ lệ chênh - Oddratio

3.4. Phân tích biến insurance (sở hữu bảo hiểm cá nhân) và biến hospital (số lần nhập viện)

3.4.1. Bảng tần suất

3.4.2. Relative Risk

3.4.3. Tỷ lệ chênh - Oddratio

3.5. Phân tích biến insurance (sở hữu bảo hiểm tư nhân) và biến chronic (số bệnh mãn tính mà người được khảo sát mắc phải)

3.5.1. Bảng tần suất

3.5.2. Relative Risk

3.5.3. Tỷ lệ chênh - Oddratio

Chương 4. Thống kê suy diễn

4.1. Kiểm định tính độc lập của 2 biến

4.1.1. Kiểm tra tính độc lập cho biến insurance (sở hữu bảo hiểm cá nhân) và biến incomed (thu nhập)

4.1.2. Kiểm tra tính độc lập cho hai biến insurance (sở hữu bảo hiểm cá nhân) và biến employed (việc làm)

4.1.3. Kiểm tra tính độc lập cho hai biến insurance (sở hữu bảo hiểm cá nhân) và biến medicaid

4.1.4. Kiểm tra tính độc lập cho hai biến insurance (sở hữu bảo hiểm cá nhân) và biến hospital (số lần nhập viện )

4.1.5. Kiểm tra tính độc lập cho biến insurance (sở hữu bảo hiểm cá nhân) và biến chronic (số bệnh mãn tính mà người được khảo sát mắc phải)

4.2. Khoảng ước lượng tỷ lệ

4.2.1. Ước lượng tỉ lệ người có và không có bảo hiểm cá nhân

4.2.2. Ước lượng tỉ lệ người có và không có việc làm

4.2.3. Ước lượng tỷ lệ người có và không có bảo hiểm medicaid

4.2.4. Ước lượng tỷ lệ người có thu nhập dưới 2.5 ngàn đô

4.2.5. Ước lượng tỷ lệ người không có bệnh mãn tính

4.2.6. Ứớc lượng người có dưới 2 lần nhập viện

Chương 5. Hồi quy

5.1. Mô hình xác suất tuyến tính

5.1.1. Biến insurance theo income

5.1.2. Biến insurance theo employed

5.1.3. Biến insurance theo medicaid

5.1.4. Hồi quy đa biến

5.2. Mô hình hồi quy logit

5.2.1. Biến insurance và incomed

5.2.2. Biến insurance và employed

5.2.3. Biến insurance và medicaid

5.2.4. Hồi quy đa biến

5.3. Mô hình probit

5.3.1. Biến insurance và incomed

5.3.2. Biến insurance và employed

5.3.3. Biến insurance và medicaid

5.3.4. Hồi quy đa biến

5.4. Đánh giá mô hình

Chương 6. Kết luận