NHIỆM VỤ TUẦN 2

PHẦN 1: TÌM HIỂU VÀ CHUẨN BỊ DỮ LIỆU

1.1 ĐỌC DỮ LIỆU

Đọc file dữ liệu gốc

data <- read.csv(file = file.choose(), header = T)

Các biến định tính

dat <- data[, dldt]: lấy trong bộ dữ liệu data các cột trong dldt vừa tạo ra (các biến định tính), và lấy tất cả các hàng

# Chọn các biến định tính
dldt <- c("Gender", "MaritalStatus", "Homeowner", "AnnualIncome", "City", "StateorProvince", "Country", "ProductFamily", "ProductDepartment", "ProductCategory")

# Tạo bộ dữ liệu mới chỉ chứa định tính
dat <- data[, dldt]

1.2. NỘI DUNG DỮ LIỆU

Bộ dữ liệu này gồm 14,059 giao dịch mua hàng tại siêu thị, với các thông tin liên quan đến khách hàng và sản phẩm bao gồm mã định danh khách hàng và các thông tin liên quan, vị trí địa lý, chi tiết các sản phẩm được mua.
Dữ liệu này có thể được sử dụng để phân tích hành vi mua sắm của khách hàng, phân khúc thị trường, hoặc đánh giá hiệu quả kinh doanh theo từng sản phẩm, từng khu vực và từng nhóm đối tượng.

Danh sách các biến và mô tả

Cấu trúc bộ dữ liệu Supermarket Transactions

str(data)

## 'data.frame':    14059 obs. of  16 variables:
##  $ X                : int  1 2 3 4 5 6 7 8 9 10 ...
##  $ PurchaseDate     : chr  "2007-12-18" "2007-12-20" "2007-12-21" "2007-12-21" ...
##  $ CustomerID       : int  7223 7841 8374 9619 1900 6696 9673 354 1293 7938 ...
##  $ Gender           : chr  "F" "M" "F" "M" ...
##  $ MaritalStatus    : chr  "S" "M" "M" "M" ...
##  $ Homeowner        : chr  "Y" "Y" "N" "Y" ...
##  $ Children         : int  2 5 2 3 3 3 2 2 3 1 ...
##  $ AnnualIncome     : chr  "$30K - $50K" "$70K - $90K" "$50K - $70K" "$30K - $50K" ...
##  $ City             : chr  "Los Angeles" "Los Angeles" "Bremerton" "Portland" ...
##  $ StateorProvince  : chr  "CA" "CA" "WA" "OR" ...
##  $ Country          : chr  "USA" "USA" "USA" "USA" ...
##  $ ProductFamily    : chr  "Food" "Food" "Food" "Food" ...
##  $ ProductDepartment: chr  "Snack Foods" "Produce" "Snack Foods" "Snacks" ...
##  $ ProductCategory  : chr  "Snack Foods" "Vegetables" "Snack Foods" "Candy" ...
##  $ UnitsSold        : int  5 5 3 4 4 3 4 6 1 2 ...
##  $ Revenue          : num  27.38 14.9 5.52 4.44 14 ...

Gender: Giới tính (F - Nữ, M - Nam)
MaritalStatus: Tình trạng hôn nhân (S - Độc thân, M - Đã kết hôn)
Homeowner: Có sở hữu nhà hay không (Y - Có, N - Không)
AnnualIncome: Thu nhập hàng năm (được biểu thị dưới dạng các khoảng)
City, StateorProvince, Country: Thông tin địa lý
ProductFamily, ProductDepartment, ProductCategory: Phân loại sản phẩm

Số biến và số quan sát

dim(dat)

## [1] 14059    10

Bộ dữ liệu dat chứa tổng cộng 14.059 quan sát cùng với 10 biến đặc trưng tương ứng là biến định tính.
Mỗi quan sát sẽ tưởng ứng với mỗi giao dịch mua hàng tại siêu thị.

Kiểm tra số lượng NA trong cột định tính

sapply(dat, …): áp dụng hàm cho từng cột trong dữ liệu dat.

sum(is.na(x)): đếm số lượng NA trong từng cột.

# Kiểm tra số lượng NA trong mỗi cột định tính
na_counts <- sapply(dat, function(x) sum(is.na(x)))
na_counts

##            Gender     MaritalStatus         Homeowner      AnnualIncome 
##                 0                 0                 0                 0 
##              City   StateorProvince           Country     ProductFamily 
##                 0                 0                 0                 0 
## ProductDepartment   ProductCategory 
##                 0                 0

Nhận xét - Kết quả thu được cho thấy rằng tất cả các cột định tính đều không có giá trị thiếu (NA) vì số lượng NA trong các cột đều hiển thị là 0

Chuyển đổi các cột sang factor

# Chuyển đổi các cột sang factor nếu chưa phải là factor
dat <- as.data.frame(lapply(dat, function(x) {
  if (!is.factor(x)) {
    return(as.factor(x))
  } else {
    return(x)
  }
}))

Kết quả thu được toàn bộ các biến đã được chuyển đổi thành dạng factor.

PHẦN 2: PHÂN TÍCH BIẾN ĐỊNH TÍNH

PRODUCTCATEGORY

Thống kê tuần suất và trực quan hóa

tab_ProductCategory <- table(dat$ProductCategory)
tab_ProductCategory

## 
##         Baking Goods    Bathroom Products        Beer and Wine 
##                  484                  365                  356 
##                Bread      Breakfast Foods              Candles 
##                  425                  417                   45 
##                Candy     Canned Anchovies         Canned Clams 
##                  352                   44                   53 
##       Canned Oysters      Canned Sardines        Canned Shrimp 
##                   35                   40                   38 
##          Canned Soup          Canned Tuna Carbonated Beverages 
##                  404                   87                  154 
##    Cleaning Supplies        Cold Remedies                Dairy 
##                  189                   93                  903 
##        Decongestants               Drinks                 Eggs 
##                   85                  135                  198 
##           Electrical      Frozen Desserts       Frozen Entrees 
##                  355                  323                  118 
##                Fruit             Hardware        Hot Beverages 
##                  765                  129                  226 
##              Hygiene     Jams and Jellies     Kitchen Products 
##                  197                  588                  217 
##            Magazines                 Meat        Miscellaneous 
##                  202                  761                   42 
##  Packaged Vegetables       Pain Relievers       Paper Products 
##                   48                  192                  345 
##                Pizza     Plastic Products Pure Juice Beverages 
##                  194                  141                  165 
##              Seafood          Side Dishes          Snack Foods 
##                  102                  153                 1600 
##            Specialty        Starchy Foods           Vegetables 
##                  289                  277                 1728

table(dat$ProductCategory)/sum(nrow(dat))

## 
##         Baking Goods    Bathroom Products        Beer and Wine 
##          0.034426346          0.025962017          0.025321858 
##                Bread      Breakfast Foods              Candles 
##          0.030229746          0.029660716          0.003200797 
##                Candy     Canned Anchovies         Canned Clams 
##          0.025037343          0.003129668          0.003769827 
##       Canned Oysters      Canned Sardines        Canned Shrimp 
##          0.002489508          0.002845153          0.002702895 
##          Canned Soup          Canned Tuna Carbonated Beverages 
##          0.028736041          0.006188207          0.010953837 
##    Cleaning Supplies        Cold Remedies                Dairy 
##          0.013443346          0.006614980          0.064229319 
##        Decongestants               Drinks                 Eggs 
##          0.006045949          0.009602390          0.014083505 
##           Electrical      Frozen Desserts       Frozen Entrees 
##          0.025250729          0.022974607          0.008393200 
##                Fruit             Hardware        Hot Beverages 
##          0.054413543          0.009175617          0.016075112 
##              Hygiene     Jams and Jellies     Kitchen Products 
##          0.014012376          0.041823743          0.015434953 
##            Magazines                 Meat        Miscellaneous 
##          0.014368020          0.054129028          0.002987410 
##  Packaged Vegetables       Pain Relievers       Paper Products 
##          0.003414183          0.013656732          0.024539441 
##                Pizza     Plastic Products Pure Juice Beverages 
##          0.013798990          0.010029163          0.011736254 
##              Seafood          Side Dishes          Snack Foods 
##          0.007255139          0.010882709          0.113806103 
##            Specialty        Starchy Foods           Vegetables 
##          0.020556227          0.019702682          0.122910591

barplot(tab_ProductCategory, main = "Danh mục sản phẩm - ProductCategory", col = "#FF0000", las = 2)

Nhận xét

Các loại thực phẩm phổ biến được quan sát từ biểu đồ:
- Vegetables (Rau): với số lượng là 1728 khách hàng và chiếm hết (12.2910591%)
- Snack Foods: 1600 khách hàng và chiếm hết (11.3806103%)
- Dairy (Sữa): 903 khách hàng và chiếm hết (6.4229319%)
- Fruit, Meat, Frozen Foods…
Đây là những mặt hàng tiêu dùng nhanh, có vòng đời ngắn.
Việc khách hàng mua chủ yếu các mặt hàng thiết yếu cho thấy siêu thị hoạt động tốt trong việc phục vụ nhu cầu cơ bản, chứ không phải sản phẩm xa xỉ hay không thường xuyên.

ANNUALINCOME

Thống kê tần suất và trực quan hóa

tab_AnnualIncome <- table(dat$AnnualIncome)
tab_AnnualIncome

## 
##   $10K - $30K $110K - $130K $130K - $150K       $150K +   $30K - $50K 
##          3090           643           760           273          4601 
##   $50K - $70K   $70K - $90K  $90K - $110K 
##          2370          1709           613

table(dat$AnnualIncome)/sum(nrow(dat))

## 
##   $10K - $30K $110K - $130K $130K - $150K       $150K +   $30K - $50K 
##    0.21978804    0.04573583    0.05405790    0.01941817    0.32726367 
##   $50K - $70K   $70K - $90K  $90K - $110K 
##    0.16857529    0.12155914    0.04360196

barplot(tab_AnnualIncome, main = "Thu nhập/tháng - AnnualIncome", col = "lightblue", las = 2)

Nhận xét

Nhóm $30K - $50K chiếm đa số với 3090 khách hàng (chiếm 32.7263675%), theo sau là:
- Nhóm $10K - $30K (21.9788036%)
- Nhóm $50K - $70K (16.857529%)
Rất ít khách hàng có thu nhập trên $150K (1.9418166%)
Điều này cho thấy rằng chuỗi siêu thị chủ yếu phục vụ nhóm thu nhập trung bình và thấp – phù hợp với chiến lược giá rẻ, hàng hóa thiết yếu.
Chiến lược marketing, khuyến mãi, và chính sách giá nên tập trung vào nhóm này (ví dụ: combo tiết kiệm, chiết khấu theo số lượng).

COUNTRY

Thống kê tần suất và trực quan hóa

tab_Country <- table(dat$Country)
tab_Country

## 
## Canada Mexico    USA 
##    809   3688   9562

table(dat$Country)/sum(nrow(dat))

## 
##     Canada     Mexico        USA 
## 0.05754321 0.26232307 0.68013372

barplot(tab_Country, main = "Quốc gia - Country", col = "orchid", las = 2)

Nhận xét

Tỷ lệ khách hàng theo quốc gia phân bổ theo trình tự giảm dần:
- Cao nhất là USA: 68.0133722%
- Kế tiếp là Mexico: 26.2323067%
- Thấp nhất là Canada: 5.7543211%
US Hoa Kỳ rõ ràng là thị trường trọng điểm với tỷ lệ khách hàng giao dịch cao nhất trong 3 quốc gia. Điều này chứng tỏ Hoa Kỳ là một quốc gia phát triển mạnh với nhu cầu tiêu dùng chiến phần lớn.
Tuy nhiên, sự hiện diện tại Mexico (hơn 1/4 số giao dịch) cho thấy đây là thị trường thứ cấp tiềm năng, có thể mở rộng hoặc tăng nhận diện thương hiệu.

GENDER

Thống kê tần suất và trực quan hóa

Sử dụng hàm table() để đếm số lần xuất hiện (tần suất) của từng giá trị trong biến định tính.
prop.table(tab_gender) tạo bảng tỷ lệ phần trăm (theo dạng thập phân). Nhân với 100 để ra phần trăm.
round(…, 2) làm tròn 2 chữ số sau dấu phẩy.
barplot() tạo biểu đồ cột từ bảng tần suất.
main: tiêu đề biểu đồ.
col: màu sắc cột (ở đây là màu “steelblue”).
las = 2: xoay nhãn trục hoành (X) theo chiều dọc để dễ đọc.

tab_gender <- table(dat$Gender)
tab_gender

## 
##    F    M 
## 7170 6889

table(dat$Gender)/sum(nrow(dat))

## 
##         F         M 
## 0.5099936 0.4900064

barplot(tab_gender, main = "Tần suất - Gender", col = "lightyellow", las = 2)

Nhận xét

Dữ liệu cho thấy có 7170 khách hàng là nữ (50.9993598%) và 6889 là khách hàng nam (49.0006402%).
Tỷ lệ giới tính giữa nam và nữ trong dữ liệu khá cân bằng, không có sự khác biệt đáng kể nào.
Điều này cho thấy chuỗi siêu thị hướng tới phục vụ cả hai giới một cách đồng đều, thay vì tập trung vào một nhóm cụ thể.
Mặt khác, việc nữ giới chiếm tỷ lệ nhỉnh hơn đôi chút có thể xuất phát từ vai trò thường thấy của họ trong việc mua sắm nhu yếu phẩm và quản lý chi tiêu gia đình, dẫn đến khả năng cao họ là người thực hiện các giao dịch được ghi nhận trong hóa đơn.

MARITALSTATUS

Thống kê tuần suất và trực quan hóa

tab_MaritalStatus <- table(dat$MaritalStatus)
tab_MaritalStatus

## 
##    M    S 
## 6866 7193

table(dat$MaritalStatus)/sum(nrow(dat))

## 
##         M         S 
## 0.4883704 0.5116296

barplot(tab_MaritalStatus, main = "Tình trạng hôn nhân - MaritalStatus", col = "skyblue", las = 2)

Nhận xét

Số lượng khách hàng độc thân (Single – ký hiệu S) là 6866 người (chiếm 48.8370439%), còn người đã kết hôn (Married – M) là 1, 0 người (chiếm 51.1629561%)
Mặc dù tỷ lệ khá cân bằng, nhưng nhóm độc thân chiếm tỷ lệ cao hơn một chút.
Giải thích hợp lý cho điều này có thể là:
- Những người độc thân có xu hướng tự đi mua sắm nhiều hơn, không dựa vào bạn đời.
- Ngoài ra, người trẻ độc thân thường sống ở thành thị và sử dụng dịch vụ siêu thị thường xuyên hơn, trong khi người đã kết hôn có thể mua theo nhóm hoặc theo hộ gia đình lớn và ít giao dịch hơn trên mỗi cá nhân.

HOMEOWNER

Thống kê tần suất và trực quan hóa

tab_Homeowner <- table(dat$Homeowner)
tab_Homeowner

## 
##    N    Y 
## 5615 8444

table(dat$Homeowner)/sum(nrow(dat))

## 
##         N         Y 
## 0.3993883 0.6006117

barplot(tab_Homeowner, main = "Sở hữu nhà - Homeowner", col = "darkseagreen", las = 2)

Nhận xét

Từ kết quả thu được ta thấy có 5615 khách hàng (chiếm 39.9388292%) không sở hữu nhà, còn 8444 khách hàng (chiếm 60.0611708%) là chủ sở hữu nhà.
Đây là một chỉ báo quan trọng về sự ổn định tài chính của khách hàng.
Người sở hữu nhà thường có:
- Chi tiêu ổn định hơn.
- Nhu cầu mua sắm nhiều hơn cho gia đình (thực phẩm, đồ gia dụng).
Trong khi đó, người không sở hữu nhà (có thể thuê nhà hoặc sống tạm thời) có thể:
- Chi tiêu dè dặt hơn.
- Ít mua hàng số lượng lớn, hoặc ưu tiên hàng nhỏ gọn, dễ vận chuyển.

CITY

Thống kê tần suất và trực quan hóa

tab_City <- table(dat$City)
tab_City

## 
##      Acapulco    Bellingham Beverly Hills     Bremerton       Camacho 
##           383           143           811           834           452 
##   Guadalajara       Hidalgo   Los Angeles        Merida   Mexico City 
##            75           845           926           654           194 
##       Orizaba      Portland         Salem    San Andres     San Diego 
##           464           876          1386           621           866 
## San Francisco       Seattle       Spokane        Tacoma     Vancouver 
##           130           922           875          1257           633 
##      Victoria   Walla Walla        Yakima 
##           176           160           376

table(dat$City)/sum(nrow(dat))

## 
##      Acapulco    Bellingham Beverly Hills     Bremerton       Camacho 
##   0.027242336   0.010171420   0.057685468   0.059321431   0.032150224 
##   Guadalajara       Hidalgo   Los Angeles        Merida   Mexico City 
##   0.005334661   0.060103848   0.065865282   0.046518245   0.013798990 
##       Orizaba      Portland         Salem    San Andres     San Diego 
##   0.033003770   0.062308841   0.098584537   0.044170994   0.061597553 
## San Francisco       Seattle       Spokane        Tacoma     Vancouver 
##   0.009246746   0.065580767   0.062237712   0.089408920   0.045024539 
##      Victoria   Walla Walla        Yakima 
##   0.012518671   0.011380610   0.026744434

barplot(tab_City, main = "Thành phố - City", col = "pink", las = 2)

Nhận xét

Các thành phố có tỷ lệ giao dịch cao là:
- Salem là thành phố có tỉ lệ khách hàng mua sản phẩm cao nhất với 1386 khách hàng và chiếm hết (9.8584537%)
- Theo sau là thành phố Tacoma (1.2518671%)
- Kế tiếp là thành phố Los Angeles (6.5865282%) và ta thấy thành phố Seattle với 6.5580767% xấp xỉ bằng với thành phố Los Angeles. Qua đó ta nhìn nhận được rằng 2 thành phố này có lượng khách hàng mua sắm gần như nhau.
Có thể nhận thấy sự tập trung mạnh tại các thành phố bờ Tây Hoa Kỳ, đặc biệt ở các bang Washington, Oregon và California.
Điều này gợi ý rằng:
- Các chi nhánh ở các thành phố này có lưu lượng khách hàng lớn, hoặc
- Các thành phố này có mật độ dân cư cao, nhu cầu tiêu dùng thường xuyên.

STATEORPROVINCE

Thống kê tần suất và trực quan hóa

tab_StateorProvince <- table(dat$StateorProvince)
tab_StateorProvince

## 
##        BC        CA        DF  Guerrero   Jalisco        OR  Veracruz        WA 
##       809      2733       815       383        75      2262       464      4567 
##   Yucatan Zacatecas 
##       654      1297

table(dat$StateorProvince)/sum(nrow(dat))

## 
##          BC          CA          DF    Guerrero     Jalisco          OR 
## 0.057543211 0.194395049 0.057969984 0.027242336 0.005334661 0.160893378 
##    Veracruz          WA     Yucatan   Zacatecas 
## 0.033003770 0.324845295 0.046518245 0.092254072

barplot(tab_StateorProvince, main = "Bang/Tỉnh - StateorProvince", col = "goldenrod", las = 2)

Nhận xét

Phân bố:
- Washington(WA): 32.4845295%
- California(CA): 19.4395049%
- Oregon(OR): 16.0893378%
Ba bang này đều ở bờ Tây Hoa Kỳ, và là những bang có:
- Dân số đông, thị trường bán lẻ phát triển.
- Thói quen tiêu dùng hiện đại, sử dụng dịch vụ siêu thị nhiều hơn so với các khu vực nông thôn.
Việc mở rộng kinh doanh hoặc đầu tư thêm vào các bang này là chiến lược khả thi.

PRODUCTFAMILY

Thống kê tần suất và trực quan hóa

tab_ProductFamily <- table(dat$ProductFamily)
tab_ProductFamily

## 
##          Drink           Food Non-Consumable 
##           1250          10153           2656

table(dat$ProductFamily)/sum(nrow(dat))

## 
##          Drink           Food Non-Consumable 
##     0.08891102     0.72217085     0.18891813

barplot(tab_ProductFamily, main = "Nhóm sản phẩm - ProductFamily", col = "darkblue", las = 2)

Nhận xét

Thực phẩm (Food) chiếm đa số (72.2170851%), cao vượt trội so với:
- Non-Consumables chiếm 18.8918131%
- Nhóm đồ uống chiếm 8.8911018%
Điều này cho thấy khách hàng mua sắm chủ yếu là để phục vụ nhu cầu ăn uống hàng ngày.
Siêu thị nên tiếp tục đầu tư mạnh vào:
- Nguồn cung thực phẩm tươi,
- Các sản phẩm thiết yếu giá tốt,
- Chính sách ưu đãi số lượng lớn cho thực phẩm.

PRODUCTDEPARTMENT

Thống kê tần suất và trực quan hóa

tab_ProductDepartment <- table(dat$ProductDepartment)
tab_ProductDepartment

## 
## Alcoholic Beverages         Baked Goods        Baking Goods           Beverages 
##                 356                 425                1072                 680 
##     Breakfast Foods        Canned Foods     Canned Products            Carousel 
##                 188                 977                 109                  59 
##            Checkout               Dairy                Deli                Eggs 
##                  82                 903                 699                 198 
##        Frozen Foods  Health and Hygiene           Household                Meat 
##                1382                 893                1420                  89 
##         Periodicals             Produce             Seafood         Snack Foods 
##                 202                1994                 102                1600 
##              Snacks       Starchy Foods 
##                 352                 277

table(dat$ProductDepartment)/sum(nrow(dat))

## 
## Alcoholic Beverages         Baked Goods        Baking Goods           Beverages 
##         0.025321858         0.030229746         0.076250089         0.048367594 
##     Breakfast Foods        Canned Foods     Canned Products            Carousel 
##         0.013372217         0.069492852         0.007753041         0.004196600 
##            Checkout               Dairy                Deli                Eggs 
##         0.005832563         0.064229319         0.049719041         0.014083505 
##        Frozen Foods  Health and Hygiene           Household                Meat 
##         0.098300021         0.063518031         0.101002916         0.006330464 
##         Periodicals             Produce             Seafood         Snack Foods 
##         0.014368020         0.141830856         0.007255139         0.113806103 
##              Snacks       Starchy Foods 
##         0.025037343         0.019702682

barplot(tab_ProductDepartment, main = "Bộ phận sản phẩm - ProductDepartment", col = "darkred", las = 2)

Nhận xét

Các bộ phân chiếm đa số trong toàn bộ các bộ phận được thống kê:
- Produce (Rau quả): 14.1830856%
- Snack Foods: 11.3806103%
- Household: 10.1002916%
Khách hàng mua nhiều các mặt hàng:
- Tươi sống, dễ hỏng (produce) -> yêu cầu bảo quản và chuỗi cung ứng tốt.
- Đồ ăn vặt và hàng gia dụng -> phản ánh thói quen tiêu dùng nhanh, tiện lợi.
Đây là cơ hội để tăng cường marketing nhóm sản phẩm tiện lợi, combo bữa ăn, sản phẩm cho hộ gia đình.

PHẦN 3: ƯỚC LƯỢNG KHOẢNG VÀ KIỂM ĐỊNH GIẢ THUYẾT CHO TỶ LỆ

table(): Thống kê số lượng từng loại trong biến định tính
sum(…): Đếm bao nhiêu dòng thỏa điều kiện (giá trị quan tâm)
length(…): Tổng số mẫu
prop.test(…): Tính khoảng tin cậy và kiểm định giả thuyết cho tỷ lệ
p = … : Tỷ lệ giả định trong giả thuyết H0
alternative = …: Kiểm định hai phía (“two.sided”), một phía lớn hơn (“greater”), nhỏ hơn (“less”)

3.1 Biến Gender

# Xác định hạng mục quan tâm
table(dat$Gender)

## 
##    F    M 
## 7170 6889

# Số lượng và tổng số 
slfemale <- sum(dat$Gender == "F")
totalgender <- length(dat$Gender)
# Ước lượng khoảng tin cậy 95% cho tỷ lệ "Female"
prop.test(slfemale, totalgender, conf.level = 0.95)

## 
##  1-sample proportions test with continuity correction
## 
## data:  slfemale out of totalgender, null probability 0.5
## X-squared = 5.5765, df = 1, p-value = 0.0182
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.5016931 0.5182886
## sample estimates:
##         p 
## 0.5099936

# Kiểm định giả thuyết: tỷ lệ Female = 0.5
prop.test(slfemale, totalgender, p = 0.5, alternative = "two.sided", conf.level = 0.95)

## 
##  1-sample proportions test with continuity correction
## 
## data:  slfemale out of totalgender, null probability 0.5
## X-squared = 5.5765, df = 1, p-value = 0.0182
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.5016931 0.5182886
## sample estimates:
##         p 
## 0.5099936

Giả thuyết

H0: p = 0.5
H1: p ≠ 0.5

Đọc kết quả và kết luận

Tổng số cá thể trong mẫu là 14.059 người, trong đó có 7.170 người là nữ và 6.889 người là nam. Tỷ lệ nữ giới trong mẫu là 50.99936%.
Kết quả kiểm định tỷ lệ một mẫu với giả thuyết tỷ lệ nữ trong tổng thể là 50% (H₀: p = 0.5) cho thấy giá trị thống kê Chi-squared = 5.5765 với 1 bậc tự do, và p-value = 0.0182. Với mức ý nghĩa thông thường α = 0.05, ta bác bỏ giả thuyết không, cho thấy tỷ lệ nữ giới trong mẫu có sự khác biệt có ý nghĩa thống kê so với tỷ lệ 50%.
Khoảng tin cậy 95% cho tỷ lệ nữ trong quần thể nằm trong khoảng từ 50.16931% đến 51.82886%, cho thấy tỷ lệ này chỉ hơi cao hơn 50%.

Tóm lại

Tuy tỷ lệ nữ giới trong mẫu chỉ nhỉnh hơn nam giới không quá cao (50.99936%), nhưng sự chênh lệch này là có ý nghĩa thống kê. Do đó, có thể kết luận rằng nữ giới chiếm ưu thế nhẹ nhưng đáng kể về mặt thống kê trong mẫu khảo sát này.

Kết luận

Tỷ lệ nữ trong tổng mẫu là 50.999%, khác có ý nghĩa thống kê so với 50%.
Với p-value = 0.0182, bạn bác bỏ giả thuyết rằng tỷ lệ nữ là 50% ở mức ý nghĩa 5%.
Tuy nhiên, sự chênh lệch dù có ý nghĩa thống kê, nhưng rất nhỏ về mặt thực tế (~1%)

3.2 Biến Homeowner

# Xác định hạng mục quan tâm
table(dat$Homeowner)

## 
##    N    Y 
## 5615 8444

slshn <- sum(dat$Homeowner == "Y")
totalhomeowner <- length(dat$Homeowner)

# Ước lượng khoảng tin cậy
prop.test(slshn, totalhomeowner, conf.level = 0.95)

## 
##  1-sample proportions test with continuity correction
## 
## data:  slshn out of totalhomeowner, null probability 0.5
## X-squared = 568.86, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.5924537 0.6087145
## sample estimates:
##         p 
## 0.6006117

# Kiểm định giả thuyết H0: tỷ lệ chủ nhà = 0.6
prop.test(slshn, totalhomeowner, p = 0.6, alternative = "two.sided", conf.level = 0.95)

## 
##  1-sample proportions test with continuity correction
## 
## data:  slshn out of totalhomeowner, null probability 0.6
## X-squared = 0.019445, df = 1, p-value = 0.8891
## alternative hypothesis: true p is not equal to 0.6
## 95 percent confidence interval:
##  0.5924537 0.6087145
## sample estimates:
##         p 
## 0.6006117

Giả thuyết 1

H0: p = 0.5
H1: p ≠ 0.5

Đọc kết quả và kết luận

Kết quả phân tích tỷ lệ một mẫu cho thấy tỷ lệ đối tượng chọn lựa trong tổng thể mẫu là 60.06%. Khi tiến hành kiểm định giả thuyết với tỷ lệ chuẩn là 50%, ta thu được giá trị thống kê Chi-squared là 568.86 với một bậc tự do, và p-value < 2.2e-16. Với mức ý nghĩa thông thường (α = 0.05), ta bác bỏ giả thuyết không (H₀: p = 0.5). Điều này cho thấy tỷ lệ thực tế khác biệt một cách có ý nghĩa thống kê so với 50%.

Giả thuyết 2

H0: p = 0.6
H1: p ≠ 0.6

Đọc kết quả và kết luận

Ngược lại, khi kiểm định với tỷ lệ giả định là 60%, giá trị thống kê Chi-squared chỉ là 0.019445 với p-value = 0.8891. Do p-value lớn hơn đáng kể so với mức ý nghĩa 0.05, ta không bác bỏ giả thuyết không (H₀: p = 0.6). Điều này cho thấy không có bằng chứng thống kê cho thấy tỷ lệ thực tế khác biệt so với 60%.

Tóm lại:

Tỷ lệ lựa chọn trong mẫu khảo sát là khoảng 60.06% và khác biệt có ý nghĩa thống kê so với tỷ lệ giả định 50%, nhưng không khác biệt đáng kể so với tỷ lệ 60%. Đây là cơ sở để giả định rằng tỷ lệ lựa chọn trong quần thể có thể vào khoảng 60%.

3.3 Biến ProductFamily

# Xác định hạng mục đề cập
table(dat$ProductFamily)

## 
##          Drink           Food Non-Consumable 
##           1250          10153           2656

slfood <- sum(dat$ProductFamily == "Food")
totalproduct <- length(dat$ProductFamily)
# Ước lượng khoảng tin cậy
prop.test(slfood, totalproduct, conf.level = 0.95)

## 
##  1-sample proportions test with continuity correction
## 
## data:  slfood out of totalproduct, null probability 0.5
## X-squared = 2774.9, df = 1, p-value < 2.2e-16
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.7146709 0.7295489
## sample estimates:
##         p 
## 0.7221709

# Kiểm định giả thuyết H0: tỷ lệ >= 0.7 (alternative = "greater": kiểm định 1 phía (>))
prop.test(slfood, totalproduct, p = 0.7, alternative = "greater", conf.level = 0.95)

## 
##  1-sample proportions test with continuity correction
## 
## data:  slfood out of totalproduct, null probability 0.7
## X-squared = 32.802, df = 1, p-value = 5.101e-09
## alternative hypothesis: true p is greater than 0.7
## 95 percent confidence interval:
##  0.7158789 1.0000000
## sample estimates:
##         p 
## 0.7221709

Giả thuyết 1

H0: p = 0.5
H1: p ≠ 0.5

Đọc kết quả và kết luận

Kết quả kiểm định tỷ lệ một mẫu cho giả thuyết rằng tỷ lệ sản phẩm thuộc nhóm Food bằng 50% (H₀: p = 0.5) cho thấy:
- Giá trị thống kê Chi-squared = 2774.9, df = 1
- p-value < 2.2e-16
Với mức ý nghĩa là 0.05, ta bác bỏ giả thuyết H0, nghĩa là tỷ lệ sản phẩm nhóm Food không chiếm 50% trong tổng số “ProductFamily”.
Khoảng tin cậy 95% cho tỷ lệ này nằm trong khoảng từ 71.47% đến 72.95%.

Giả thuyết 2

H0: p = 0.7
H1: p > 0.7

Đọc kết quả và kết luận

Kết quả kiểm định tiếp theo kiểm tra xem liệu tỷ lệ sản phẩm nhóm Food có cao hơn 70% hay không, cho thấy:
- Giá trị Chi-squared = 32.802, df = 1
- p-value = 5.101e-09
Với mức ý nghĩa 0.05, ta tiếp tục bác bỏ giả thuyết không, và chấp nhận giả thuyết đối (H₁: p > 0.7). Điều này chứng tỏ tỷ lệ sản phẩm nhóm Food lớn hơn 70% với ý nghĩa thống kê cao.
Khoảng tin cậy 95% cho tỷ lệ này nằm trong khoảng từ 71.58% đến 100%.

PHẦN 4: PHÂN TÍCH MỐI QUAN HỆ GIỮA HAI BIẾN ĐỊNH TÍNH

Chọn lần lượt các cặp biến có mối tương quan thú vị làm ảnh hưởng đến hành vi mua hàng, giao dịch của khách hàng.

install.packages("tidyverse", repos = "https://cran.r-project.org")

## Installing package into 'C:/Users/Khanh Hoa/AppData/Local/R/win-library/4.4'
## (as 'lib' is unspecified)

## package 'tidyverse' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\Khanh Hoa\AppData\Local\Temp\Rtmpcb0VMf\downloaded_packages

install.packages("janitor", repos = "https://cran.r-project.org")

## Installing package into 'C:/Users/Khanh Hoa/AppData/Local/R/win-library/4.4'
## (as 'lib' is unspecified)

## package 'janitor' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\Khanh Hoa\AppData\Local\Temp\Rtmpcb0VMf\downloaded_packages

install.packages("gmodels", repos = "https://cran.r-project.org")

## Installing package into 'C:/Users/Khanh Hoa/AppData/Local/R/win-library/4.4'
## (as 'lib' is unspecified)

## package 'gmodels' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\Khanh Hoa\AppData\Local\Temp\Rtmpcb0VMf\downloaded_packages

install.packages("ggplot2", repos = "https://cran.r-project.org")

## Installing package into 'C:/Users/Khanh Hoa/AppData/Local/R/win-library/4.4'
## (as 'lib' is unspecified)

## package 'ggplot2' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\Khanh Hoa\AppData\Local\Temp\Rtmpcb0VMf\downloaded_packages

install.packages("scales", repos = "https://cran.r-project.org")

## Installing package into 'C:/Users/Khanh Hoa/AppData/Local/R/win-library/4.4'
## (as 'lib' is unspecified)

## package 'scales' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\Khanh Hoa\AppData\Local\Temp\Rtmpcb0VMf\downloaded_packages

library(tidyverse)

## Warning: package 'tidyverse' was built under R version 4.4.3

## Warning: package 'ggplot2' was built under R version 4.4.3

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.2     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2

## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(janitor)

## Warning: package 'janitor' was built under R version 4.4.3

## 
## Attaching package: 'janitor'
## 
## The following objects are masked from 'package:stats':
## 
##     chisq.test, fisher.test

library(gmodels)

## Warning: package 'gmodels' was built under R version 4.4.3

library(ggplot2)
library(scales)

## Warning: package 'scales' was built under R version 4.4.3

## 
## Attaching package: 'scales'
## 
## The following object is masked from 'package:purrr':
## 
##     discard
## 
## The following object is masked from 'package:readr':
## 
##     col_factor

4.1 Gender và ProductFamily

Bảng tần suất chéo

tab1 <- table(dat$Gender, dat$ProductFamily)
tab1

##    
##     Drink Food Non-Consumable
##   F   669 5149           1352
##   M   581 5004           1304

prop.table(tab1, 1) # theo hàng

##    
##          Drink       Food Non-Consumable
##   F 0.09330544 0.71813110     0.18856346
##   M 0.08433735 0.72637538     0.18928727

Dựa trên bảng tần suất chéo giữa Giới tính (Gender) và Nhóm sản phẩm (ProductFamily):
Nữ giới (F) có xu hướng chọn sản phẩm Food nhiều nhất, chiếm khoảng 71.8% trong tổng số sản phẩm mà nữ giới chọn. Tiếp theo là nhóm Non-Consumable (18.9%) và Drink (9.3%).
Nam giới (M) cũng có xu hướng tương tự khi lựa chọn nhiều nhất là Food (72.6%), kế đến là Non-Consumable (18.9%) và ít nhất là Drink (8.4%).
So sánh giới tính:
- Cả hai giới đều ưu tiên nhóm sản phẩm Food với tỷ lệ cao nhất, trên 70%.
- Tỷ lệ chọn Drink của nữ cao hơn nam một chút (9.3% so với 8.4%).
- Tỷ lệ chọn Non-Consumable gần như tương đương giữa hai giới (18.9%).

Kết luận:

Sở thích về nhóm sản phẩm không có sự khác biệt đáng kể giữa nam và nữ. Nhóm sản phẩm Food chiếm ưu thế trong lựa chọn của cả hai giới tính, trong khi nhóm Drink là ít được chọn nhất.

Trực quan hóa

ggplot(dat, aes(x = Gender, fill = ProductFamily)) +
  geom_bar(position = "fill") +
  labs(title = "Tỷ lệ ProductFamily theo  Gender", y = "Tỷ lệ", x = "Gender") +
  scale_y_continuous(labels = scales::percent) +
  theme_minimal()

Nhận xét

Dựa trên bảng tần suất chéo giữa Giới tính (Gender) và Nhóm sản phẩm (ProductFamily):
Nữ giới (F) có xu hướng chọn sản phẩm Food nhiều nhất, chiếm khoảng 71.8% trong tổng số sản phẩm mà nữ giới chọn. Tiếp theo là nhóm Non-Consumable (18.9%) và Drink (9.3%).
Nam giới (M) cũng có xu hướng tương tự khi lựa chọn nhiều nhất là Food (72.6%), kế đến là Non-Consumable (18.9%) và ít nhất là Drink (8.4%).
So sánh giới tính:
- Cả hai giới đều ưu tiên nhóm sản phẩm Food với tỷ lệ cao nhất, trên 70%.
- Tỷ lệ chọn Drink của nữ cao hơn nam một chút (9.3% so với 8.4%).
- Tỷ lệ chọn Non-Consumable gần như tương đương giữa hai giới (18.9%).

Kết luận:

Sở thích về nhóm sản phẩm không có sự khác biệt đáng kể giữa nam và nữ. Nhóm sản phẩm Food chiếm ưu thế trong lựa chọn của cả hai giới, trong khi nhóm Drink là ít được chọn nhất.

Kiểm định Chi-bình phương

chi_result1 <- chisq.test(tab1)
chi_result1

## 
##  Pearson's Chi-squared test
## 
## data:  tab1
## X-squared = 3.5185, df = 2, p-value = 0.1722

# Giá trị thống kê Chi-squared
chi_result1$statistic

## X-squared 
##   3.51849

# Bậc tự do 
chi_result1$parameter

## df 
##  2

# Giá trị p
chi_result1$p.value

## [1] 0.1721748

Giả thuyết

Giả thuyết H0: Giới tính (Gender) và Nhóm sản phẩm (ProductFamily) là độc lập, không có mối liên hệ.
Giả thuyết H1: Giới tính và Nhóm sản phẩm có liên hệ với nhau.

Kết quả kiểm định:

Giá trị thống kê Chi-squared: X-squared = 3.51849
Bậc tự do: df = 2
Giá trị p: p = 0.1721748

Kết luận:

Với mức ý nghĩa 0.05, p_value = 0.1722 >= 0.05 → không đủ bằng chứng để bác bỏ H0, tức là không có mối liên hệ giữa 2 biến Gender và ProductFamily.

4.2 MaritalStatus và Homeowner

Bảng tần suất chéo

tab2 <- table(dat$MaritalStatus, dat$Homeowner)
tab2

##    
##        N    Y
##   M 1719 5147
##   S 3896 3297

prop.table(tab2, 1)

##    
##             N         Y
##   M 0.2503641 0.7496359
##   S 0.5416377 0.4583623

Trực quan hóa

ggplot(dat, aes(x = MaritalStatus, fill = Homeowner)) +
  geom_bar(position = "fill") +
  labs(title = "Tỷ lệ Homeowner theo Marital Status", y = "Tỷ lệ", x = "Marital Status") +
  scale_y_continuous(labels = scales::percent) +
  scale_fill_manual(values = c("Y" = "#FFDAB9", "N" = "#D8BFD8"))

  theme_minimal()

## List of 136
##  $ line                            :List of 6
##   ..$ colour       : chr "black"
##   ..$ linewidth    : num 0.5
##   ..$ linetype     : num 1
##   ..$ lineend      : chr "butt"
##   ..$ arrow        : logi FALSE
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_line" "element"
##  $ rect                            :List of 5
##   ..$ fill         : chr "white"
##   ..$ colour       : chr "black"
##   ..$ linewidth    : num 0.5
##   ..$ linetype     : num 1
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_rect" "element"
##  $ text                            :List of 11
##   ..$ family       : chr ""
##   ..$ face         : chr "plain"
##   ..$ colour       : chr "black"
##   ..$ size         : num 11
##   ..$ hjust        : num 0.5
##   ..$ vjust        : num 0.5
##   ..$ angle        : num 0
##   ..$ lineheight   : num 0.9
##   ..$ margin       : 'margin' num [1:4] 0points 0points 0points 0points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : logi FALSE
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ title                           : NULL
##  $ aspect.ratio                    : NULL
##  $ axis.title                      : NULL
##  $ axis.title.x                    :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : NULL
##   ..$ vjust        : num 1
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 2.75points 0points 0points 0points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.title.x.top                :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : NULL
##   ..$ vjust        : num 0
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 0points 0points 2.75points 0points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.title.x.bottom             : NULL
##  $ axis.title.y                    :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : NULL
##   ..$ vjust        : num 1
##   ..$ angle        : num 90
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 0points 2.75points 0points 0points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.title.y.left               : NULL
##  $ axis.title.y.right              :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : NULL
##   ..$ vjust        : num 1
##   ..$ angle        : num -90
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 0points 0points 0points 2.75points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.text                       :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : chr "grey30"
##   ..$ size         : 'rel' num 0.8
##   ..$ hjust        : NULL
##   ..$ vjust        : NULL
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : NULL
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.text.x                     :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : NULL
##   ..$ vjust        : num 1
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 2.2points 0points 0points 0points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.text.x.top                 :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : NULL
##   ..$ vjust        : num 0
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 0points 0points 2.2points 0points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.text.x.bottom              : NULL
##  $ axis.text.y                     :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : num 1
##   ..$ vjust        : NULL
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 0points 2.2points 0points 0points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.text.y.left                : NULL
##  $ axis.text.y.right               :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : num 0
##   ..$ vjust        : NULL
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 0points 0points 0points 2.2points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.text.theta                 : NULL
##  $ axis.text.r                     :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : num 0.5
##   ..$ vjust        : NULL
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : 'margin' num [1:4] 0points 2.2points 0points 2.2points
##   .. ..- attr(*, "unit")= int 8
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ axis.ticks                      : list()
##   ..- attr(*, "class")= chr [1:2] "element_blank" "element"
##  $ axis.ticks.x                    : NULL
##  $ axis.ticks.x.top                : NULL
##  $ axis.ticks.x.bottom             : NULL
##  $ axis.ticks.y                    : NULL
##  $ axis.ticks.y.left               : NULL
##  $ axis.ticks.y.right              : NULL
##  $ axis.ticks.theta                : NULL
##  $ axis.ticks.r                    : NULL
##  $ axis.minor.ticks.x.top          : NULL
##  $ axis.minor.ticks.x.bottom       : NULL
##  $ axis.minor.ticks.y.left         : NULL
##  $ axis.minor.ticks.y.right        : NULL
##  $ axis.minor.ticks.theta          : NULL
##  $ axis.minor.ticks.r              : NULL
##  $ axis.ticks.length               : 'simpleUnit' num 2.75points
##   ..- attr(*, "unit")= int 8
##  $ axis.ticks.length.x             : NULL
##  $ axis.ticks.length.x.top         : NULL
##  $ axis.ticks.length.x.bottom      : NULL
##  $ axis.ticks.length.y             : NULL
##  $ axis.ticks.length.y.left        : NULL
##  $ axis.ticks.length.y.right       : NULL
##  $ axis.ticks.length.theta         : NULL
##  $ axis.ticks.length.r             : NULL
##  $ axis.minor.ticks.length         : 'rel' num 0.75
##  $ axis.minor.ticks.length.x       : NULL
##  $ axis.minor.ticks.length.x.top   : NULL
##  $ axis.minor.ticks.length.x.bottom: NULL
##  $ axis.minor.ticks.length.y       : NULL
##  $ axis.minor.ticks.length.y.left  : NULL
##  $ axis.minor.ticks.length.y.right : NULL
##  $ axis.minor.ticks.length.theta   : NULL
##  $ axis.minor.ticks.length.r       : NULL
##  $ axis.line                       : list()
##   ..- attr(*, "class")= chr [1:2] "element_blank" "element"
##  $ axis.line.x                     : NULL
##  $ axis.line.x.top                 : NULL
##  $ axis.line.x.bottom              : NULL
##  $ axis.line.y                     : NULL
##  $ axis.line.y.left                : NULL
##  $ axis.line.y.right               : NULL
##  $ axis.line.theta                 : NULL
##  $ axis.line.r                     : NULL
##  $ legend.background               : list()
##   ..- attr(*, "class")= chr [1:2] "element_blank" "element"
##  $ legend.margin                   : 'margin' num [1:4] 5.5points 5.5points 5.5points 5.5points
##   ..- attr(*, "unit")= int 8
##  $ legend.spacing                  : 'simpleUnit' num 11points
##   ..- attr(*, "unit")= int 8
##  $ legend.spacing.x                : NULL
##  $ legend.spacing.y                : NULL
##  $ legend.key                      : list()
##   ..- attr(*, "class")= chr [1:2] "element_blank" "element"
##  $ legend.key.size                 : 'simpleUnit' num 1.2lines
##   ..- attr(*, "unit")= int 3
##  $ legend.key.height               : NULL
##  $ legend.key.width                : NULL
##  $ legend.key.spacing              : 'simpleUnit' num 5.5points
##   ..- attr(*, "unit")= int 8
##  $ legend.key.spacing.x            : NULL
##  $ legend.key.spacing.y            : NULL
##  $ legend.frame                    : NULL
##  $ legend.ticks                    : NULL
##  $ legend.ticks.length             : 'rel' num 0.2
##  $ legend.axis.line                : NULL
##  $ legend.text                     :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : 'rel' num 0.8
##   ..$ hjust        : NULL
##   ..$ vjust        : NULL
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : NULL
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ legend.text.position            : NULL
##  $ legend.title                    :List of 11
##   ..$ family       : NULL
##   ..$ face         : NULL
##   ..$ colour       : NULL
##   ..$ size         : NULL
##   ..$ hjust        : num 0
##   ..$ vjust        : NULL
##   ..$ angle        : NULL
##   ..$ lineheight   : NULL
##   ..$ margin       : NULL
##   ..$ debug        : NULL
##   ..$ inherit.blank: logi TRUE
##   ..- attr(*, "class")= chr [1:2] "element_text" "element"
##  $ legend.title.position           : NULL
##  $ legend.position                 : chr "right"
##  $ legend.position.inside          : NULL
##  $ legend.direction                : NULL
##  $ legend.byrow                    : NULL
##  $ legend.justification            : chr "center"
##  $ legend.justification.top        : NULL
##  $ legend.justification.bottom     : NULL
##  $ legend.justification.left       : NULL
##  $ legend.justification.right      : NULL
##  $ legend.justification.inside     : NULL
##  $ legend.location                 : NULL
##  $ legend.box                      : NULL
##  $ legend.box.just                 : NULL
##  $ legend.box.margin               : 'margin' num [1:4] 0cm 0cm 0cm 0cm
##   ..- attr(*, "unit")= int 1
##  $ legend.box.background           : list()
##   ..- attr(*, "class")= chr [1:2] "element_blank" "element"
##  $ legend.box.spacing              : 'simpleUnit' num 11points
##   ..- attr(*, "unit")= int 8
##   [list output truncated]
##  - attr(*, "class")= chr [1:2] "theme" "gg"
##  - attr(*, "complete")= logi TRUE
##  - attr(*, "validate")= logi TRUE

Nhận xét

Người đã kết hôn (M) có tỷ lệ sở hữu nhà cao hơn đáng kể so với người độc thân (S):
- Trong nhóm đã kết hôn (M): Gần 75% là chủ sở hữu nhà (Y).
- Trong nhóm độc thân (S): Chỉ khoảng 50% là chủ sở hữu nhà, phần còn lại là không sở hữu nhà (N).
Sự khác biệt rõ ràng trong cấu trúc tài sản theo tình trạng hôn nhân:
- Điều này phản ánh thực tế rằng những người kết hôn thường có động lực kinh tế và tài chính ổn định hơn, hoặc có xu hướng đầu tư chung vào tài sản như nhà ởở

Kiểm định Chi-bình phương

chi_result2 <- chisq.test(tab2)
chi_result2

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

# Giá trị thống kê Chi-squared
chi_result2$statistic

## X-squared 
##  1241.218

# Bậc tự do 
chi_result2$parameter

## df 
##  1

# Giá trị p
chi_result2$p.value

## [1] 6.724506e-272

Giả thuyết

Giả thuyết H0: Tình trạng hôn nhân (MaritalStatus) và sở hữu nhà (Homeowwner) là độc lập, không có mối liên hệ.
Giả thuyết H1: Tình trạng hôn nhân (MaritalStatus) và sở hữu nhà (Homeowwner) có liên hệ với nhau.

Kết quả kiểm định:

Giá trị thống kê Chi-squared: X-squared = 1241.22
Bậc tự do: df = 1
Giá trị p: p < 2.2e-16

Kết luận:

Với mức ý nghĩa 0.05, p_value < 0.05 → ta đủ điều kiện để bác bỏ H0, tức là có mối liên hệ giữa 2 biến MaritalStatus và Homeowner. Có thể hiểu rằng tình trạng hôn nhân có ảnh hưởng đến việc một người có sở hữu nhà hay không (hoặc ngược lại).

4.3 AnnualIncome và ProductCategory

Bảng tần suất chéo

tab3 <- table(dat$AnnualIncome, dat$ProductDepartment)
tab3

##                
##                 Alcoholic Beverages Baked Goods Baking Goods Beverages
##   $10K - $30K                    80         108          256       148
##   $110K - $130K                  14          23           49        35
##   $130K - $150K                  15          24           54        25
##   $150K +                         3          10           19        20
##   $30K - $50K                   121         134          336       228
##   $50K - $70K                    61          63          189        97
##   $70K - $90K                    39          50          130        91
##   $90K - $110K                   23          13           39        36
##                
##                 Breakfast Foods Canned Foods Canned Products Carousel Checkout
##   $10K - $30K                54          203              29        8       20
##   $110K - $130K               3           49               4        1        2
##   $130K - $150K              11           46               6        3        9
##   $150K +                     3           21               0        1        1
##   $30K - $50K                63          335              37       30       21
##   $50K - $70K                31          170              18        8       22
##   $70K - $90K                19          120              11        7        5
##   $90K - $110K                4           33               4        1        2
##                
##                 Dairy Deli Eggs Frozen Foods Health and Hygiene Household Meat
##   $10K - $30K     174  147   37          327                201       313   24
##   $110K - $130K    38   33    9           66                 33        63    5
##   $130K - $150K    49   45   15           84                 51        76    4
##   $150K +          17   13    4           26                 19        27    2
##   $30K - $50K     299  233   67          442                292       432   27
##   $50K - $70K     160  107   36          212                134       272   17
##   $70K - $90K     126   96   20          151                120       182    7
##   $90K - $110K     40   25   10           74                 43        55    3
##                
##                 Periodicals Produce Seafood Snack Foods Snacks Starchy Foods
##   $10K - $30K            49     420      17         329     76            70
##   $110K - $130K          17      77       3          85     16            18
##   $130K - $150K           9     109       8          83     19            15
##   $150K +                 1      33       5          35      8             5
##   $30K - $50K            65     650      44         533    127            85
##   $50K - $70K            36     344      11         274     60            48
##   $70K - $90K            18     264      12         184     29            28
##   $90K - $110K            7      97       2          77     17             8

prop.table(tab3, 1)

##                
##                 Alcoholic Beverages Baked Goods Baking Goods   Beverages
##   $10K - $30K           0.025889968 0.034951456  0.082847896 0.047896440
##   $110K - $130K         0.021772939 0.035769829  0.076205288 0.054432348
##   $130K - $150K         0.019736842 0.031578947  0.071052632 0.032894737
##   $150K +               0.010989011 0.036630037  0.069597070 0.073260073
##   $30K - $50K           0.026298631 0.029124103  0.073027603 0.049554445
##   $50K - $70K           0.025738397 0.026582278  0.079746835 0.040928270
##   $70K - $90K           0.022820363 0.029256875  0.076067876 0.053247513
##   $90K - $110K          0.037520392 0.021207178  0.063621533 0.058727569
##                
##                 Breakfast Foods Canned Foods Canned Products    Carousel
##   $10K - $30K       0.017475728  0.065695793     0.009385113 0.002588997
##   $110K - $130K     0.004665630  0.076205288     0.006220840 0.001555210
##   $130K - $150K     0.014473684  0.060526316     0.007894737 0.003947368
##   $150K +           0.010989011  0.076923077     0.000000000 0.003663004
##   $30K - $50K       0.013692676  0.072810259     0.008041730 0.006520322
##   $50K - $70K       0.013080169  0.071729958     0.007594937 0.003375527
##   $70K - $90K       0.011117613  0.070216501     0.006436513 0.004095963
##   $90K - $110K      0.006525285  0.053833605     0.006525285 0.001631321
##                
##                    Checkout       Dairy        Deli        Eggs Frozen Foods
##   $10K - $30K   0.006472492 0.056310680 0.047572816 0.011974110  0.105825243
##   $110K - $130K 0.003110420 0.059097978 0.051321928 0.013996890  0.102643857
##   $130K - $150K 0.011842105 0.064473684 0.059210526 0.019736842  0.110526316
##   $150K +       0.003663004 0.062271062 0.047619048 0.014652015  0.095238095
##   $30K - $50K   0.004564225 0.064985873 0.050641165 0.014562052  0.096066073
##   $50K - $70K   0.009282700 0.067510549 0.045147679 0.015189873  0.089451477
##   $70K - $90K   0.002925688 0.073727326 0.056173201 0.011702750  0.088355764
##   $90K - $110K  0.003262643 0.065252855 0.040783034 0.016313214  0.120717781
##                
##                 Health and Hygiene   Household        Meat Periodicals
##   $10K - $30K          0.065048544 0.101294498 0.007766990 0.015857605
##   $110K - $130K        0.051321928 0.097978227 0.007776050 0.026438569
##   $130K - $150K        0.067105263 0.100000000 0.005263158 0.011842105
##   $150K +              0.069597070 0.098901099 0.007326007 0.003663004
##   $30K - $50K          0.063464464 0.093892632 0.005868290 0.014127364
##   $50K - $70K          0.056540084 0.114767932 0.007172996 0.015189873
##   $70K - $90K          0.070216501 0.106495026 0.004095963 0.010532475
##   $90K - $110K         0.070146819 0.089722675 0.004893964 0.011419250
##                
##                     Produce     Seafood Snack Foods      Snacks Starchy Foods
##   $10K - $30K   0.135922330 0.005501618 0.106472492 0.024595469   0.022653722
##   $110K - $130K 0.119751166 0.004665630 0.132192846 0.024883359   0.027993779
##   $130K - $150K 0.143421053 0.010526316 0.109210526 0.025000000   0.019736842
##   $150K +       0.120879121 0.018315018 0.128205128 0.029304029   0.018315018
##   $30K - $50K   0.141273636 0.009563138 0.115844382 0.027602695   0.018474245
##   $50K - $70K   0.145147679 0.004641350 0.115611814 0.025316456   0.020253165
##   $70K - $90K   0.154476302 0.007021650 0.107665301 0.016968988   0.016383850
##   $90K - $110K  0.158238173 0.003262643 0.125611746 0.027732463   0.013050571

Trực quan hóa

ggplot(dat, aes(x = AnnualIncome, fill = ProductDepartment)) +
  geom_bar(position = "fill") +
  labs(title = "Tỷ lệ ProductDepartment theo Annual Income", y = "Tỷ lệ", x = "Annual Income") +
  scale_y_continuous(labels = scales::percent) +
  theme_minimal()

Nhận xét

Kiểm định chi - bình phương

chi_result3 <- chisq.test(tab3)

## Warning in stats::chisq.test(x, y, ...): Chi-squared approximation may be
## incorrect

chi_result3

## 
##  Pearson's Chi-squared test
## 
## data:  tab3
## X-squared = 170.16, df = 147, p-value = 0.09276

# Giá trị thống kê Chi-squared
chi_result3$statistic

## X-squared 
##  170.1572

# Bậc tự do 
chi_result3$parameter

##  df 
## 147

# Giá trị p
chi_result3$p.value

## [1] 0.0927647

Giả thuyết

Giả thuyết H0: Thu nhập hàng tháng (AnnualIncome) và ProductDepartment là độc lập, không có mối liên hệ.
Giả thuyết H1: Thu nhập hàng tháng (AnnualIncome) và ProductDepartment có liên hệ với nhau.

Kết quả kiểm định:

Giá trị thống kê Chi-squared: X-squared = 170.16
Bậc tự do: df = 147
Giá trị p: p = 0.0927647

Kết luận:

Với mức ý nghĩa 0.05 p_value = 0.09276 > 0.05 → ta không đủ điều kiện để bác bỏ H0, tức là không có mối liên hệ giữa 2 biến AnnualIncome và ProductDepartment.

Homeowner và AnnualIncome

Bảng tần suất chéo

tab4 <- table(dat$Homeowner, dat$AnnualIncome)
tab4

##    
##     $10K - $30K $110K - $130K $130K - $150K $150K + $30K - $50K $50K - $70K
##   N        1359           119           136      48        2087        1063
##   Y        1731           524           624     225        2514        1307
##    
##     $70K - $90K $90K - $110K
##   N         686          117
##   Y        1023          496

prop.table(tab4, 1)

##    
##     $10K - $30K $110K - $130K $130K - $150K     $150K + $30K - $50K $50K - $70K
##   N 0.242030276   0.021193232   0.024220837 0.008548531 0.371682992 0.189314337
##   Y 0.204997631   0.062055898   0.073898626 0.026646139 0.297726196 0.154784462
##    
##     $70K - $90K $90K - $110K
##   N 0.122172752  0.020837044
##   Y 0.121151113  0.058739934

Trực quan hóa

ggplot(dat, aes(x = AnnualIncome, fill = Homeowner)) +
  geom_bar(position = "fill") +
  labs(title = "Tỷ lệ Homeowner theo Annual Income", y = "Tỷ lệ", x = "Annual Income") +
  scale_y_continuous(labels = scales::percent) +
  theme_minimal()

Nhận xét

Kiểm định chi - bình phương

chi_result4 <- chisq.test(tab4)
chi_result4

## 
##  Pearson's Chi-squared test
## 
## data:  tab4
## X-squared = 546.37, df = 7, p-value < 2.2e-16

# Giá trị thống kê Chi-squared
chi_result4$statistic

## X-squared 
##  546.3723

# Bậc tự do 
chi_result4$parameter

## df 
##  7

# Giá trị p
chi_result4$p.value

## [1] 8.517274e-114

Giả thuyết

Giả thuyết H0: Thu nhập hàng tháng (AnnualIncome) và tình trạng sở hữu nhà (Homeowner) là độc lập, không có mối liên hệ.
Giả thuyết H1: Thu nhập hàng tháng (AnnualIncome) và tình trạng sở hữu nhà (Homeowner) có liên hệ với nhau.

Kết quả kiểm định:

Giá trị thống kê Chi-squared: X-squared = 546.37
Bậc tự do: df = 7
Giá trị p: p <2.2e-16

Kết luận:

Với mức ý nghĩa 0.05 p_value < 0.05 → ta có đủ điều kiện để bác bỏ H0, tức là có mối liên hệ giữa 2 biến thu nhập hàng tháng và tình trạng sở hữu nhà.

4.5 City và Gender

Bảng tần suất chéo

tab5 <- table(dat$Gender, dat$Country)
tab5

##    
##     Canada Mexico  USA
##   F    373   2056 4741
##   M    436   1632 4821

prop.table(tab5, 1)

##    
##         Canada     Mexico        USA
##   F 0.05202232 0.28675035 0.66122734
##   M 0.06328930 0.23689940 0.69981129

Trực quan hóa

ggplot(dat, aes(x = Country, fill = Gender)) +
  geom_bar(position = "fill") +
  labs(title = "Tỷ lệ giới tính theo quốc gia", y = "Tỷ lệ", x = "Country") +
  scale_y_continuous(labels = scales::percent) +
  scale_fill_manual(values = c("M" = "yellow", "F" = "pink")) +
  theme_minimal()

Nhận xét

Kiểm định chi - bình phương

chi_result5 <- chisq.test(tab5)
chi_result5

## 
##  Pearson's Chi-squared test
## 
## data:  tab5
## X-squared = 48.725, df = 2, p-value = 2.628e-11

# Giá trị thống kê Chi-squared
chi_result5$statistic

## X-squared 
##  48.72464

# Bậc tự do 
chi_result5$parameter

## df 
##  2

# Giá trị p
chi_result5$p.value

## [1] 2.627719e-11

Giả thuyết

Giả thuyết H0: Giới tính (Gender) và quốc gia (Country) là 2 biến độc lập, không có mối liên hệ.
Giả thuyết H1: Giới tính (Gender) và quốc gia (Country) là 2 biến có mối liên hệ với nhau.

Kết quả kiểm định:

Giá trị thống kê Chi-squared: X-squared = 48.725
Bậc tự do: df = 2
Giá trị p: p = 2.628e-11

Kết luận:

Với mức ý nghĩa 0.05 p_value = 2.628e-11 < 0.05 → ta có đủ điều kiện để bác bỏ H0, tức là có mối liên hệ giữa 2 biến là giới tính và quốc gia.

PHẦN 5: RELATIVE RISK

5.1. Giới thiệu

Trong phân tích thống kê định lượng, bảng ngẫu nhiên (contingency table) là một công cụ cơ bản để khảo sát mối liên hệ giữa hai biến phân loại. Trong trường hợp đơn giản nhất, bảng ngẫu nhiên 2x2 biểu diễn sự phân bố của hai biến nhị phân, giúp xác định và đo lường mối liên hệ giữa chúng. Các chỉ số như hiệu tỷ lệ (risk difference), tỷ số nguy cơ (Relative Risk - RR) và tỷ số chênh (Odds Ratio - OR) đóng vai trò quan trọng trong việc lượng hóa mối quan hệ đó. Bài viết này sẽ trình bày chi tiết cấu trúc xác suất sinh ra bảng ngẫu nhiên, phương pháp so sánh hai tỷ lệ, cách xây dựng khoảng tin cậy cho Odds Ratio, và kết thúc bằng một ví dụ thực tiễn trong lĩnh vực kinh doanh.

5.2. Cấu trúc và xác suất của bảng ngẫu nhiên

Một bảng ngẫu nhiên 2x2 là bảng tần suất đếm số quan sát thuộc vào từng tổ hợp của hai biến nhị phân:

\[ \begin{array}{|c|c|c|c|} \hline & \text{Kết quả (+)} & \text{Kết quả (–)} & \text{Tổng} \\ \hline \text{Phơi nhiễm (Yes)} & a & b & a + b \\ \hline \text{Không phơi nhiễm (No)} & c & d & c + d \\ \hline \text{Tổng cộng} & a + c & b + d & n = a + b + c + d \\ \hline \end{array} \]

Để hiểu được sự hình thành của bảng này, cần xác định mô hình xác suất sinh ra dữ liệu — trong đó phổ biến nhất là phân phối Poisson và Multinomial.

5.2.1. Phân phối Poisson

Phân phối Poisson thường được sử dụng để mô hình hóa số lượng sự kiện xảy ra trong một khoảng thời gian, không gian hoặc đơn vị cụ thể. Giả sử các sự kiện xảy ra độc lập và với một tỷ lệ trung bình không đổi, mỗi ô trong bảng có thể xem là biến ngẫu nhiên Poisson:

\[ X_{ij} \sim \text{Poisson}(\lambda_{ij}) \]

Ưu điểm:

Phù hợp khi tổng số quan sát không cố định.
Áp dụng trong phân tích số sự kiện như: số lượt truy cập, số đơn hàng lỗi, v.v.

Hạn chế:

Không kiểm soát được tổng số mẫu (n).

5.2.2. Phân phối Multinomial

Trong trường hợp tổng số quan sát $n$ là cố định, và mỗi quan sát rơi vào một trong bốn ô với xác suất $p_1, p_2, p_3, p_4$, thì bảng ngẫu nhiên 2x2 có thể được mô hình hóa theo phân phối đa thức Multinomial như sau:

\[ (a, b, c, d) \sim \text{Multinomial}(n; p_1, p_2, p_3, p_4) \]

Phân phối này thường được sử dụng trong các tình huống mà dữ liệu được thu thập từ khảo sát hoặc nghiên cứu xã hội học với kích thước mẫu cố định. Khi đó, các xác suất $p_1, p_2, p_3, p_4$ biểu diễn xác suất mà một cá thể rơi vào từng ô trong bảng 2x2.

Ưu điểm của mô hình Multinomial:

Kiểm soát được tổng số mẫu.
Phù hợp cho dữ liệu khảo sát, thử nghiệm.

So sánh giữa phân phối Poisson và Multinomial

Tiêu chí	Poisson	Multinomial
Tổng số mẫu	Không cố định	Cố định
Ứng dụng chính	Số sự kiện	Tần suất khảo sát
Dữ liệu phù hợp	Giao dịch, lỗi, tai nạn	Trả lời khảo sát, phân nhóm

5.3. Các chỉ số đo mối liên hệ

5.3.1. Hiệu tỷ lệ (Risk Difference - RD)

\[ RD = \frac{a}{a + b} - \frac{c}{c + d} \]

Ý nghĩa:

Đo lường chênh lệch tuyệt đối giữa xác suất xảy ra kết quả ở nhóm phơi nhiễm và nhóm không phơi nhiễm.
Thường được sử dụng trong đánh giá tác động của chính sách, can thiệp hoặc chương trình thí điểm, khi sự khác biệt về xác suất là quan trọng hơn so với tỷ số.

Khoảng tin cậy 95% (CI) cho RD có thể được ước lượng bằng:

\[ CI_{RD} = RD \pm Z_{1 - \alpha/2} \cdot SE_{RD} \]

Trong đó: - $Z_{1 - \alpha/2}$ là giá trị tới hạn (thường ≈ 1.96 với 95% CI),

$SE_{RD} = \sqrt{\frac{a(a+b-a)}{(a+b)^3} + \frac{c(c+d-c)}{(c+d)^3}}$

5.3.2. Tỷ số nguy cơ (Relative Risk - RR)

\[ RR = \frac{a / (a + b)}{c / (c + d)} \]

Ý nghĩa:

$RR > 1$: nguy cơ xảy ra kết quả cao hơn ở nhóm phơi nhiễm.
$RR < 1$: nguy cơ xảy ra kết quả thấp hơn ở nhóm phơi nhiễm.
$RR = 1$: không có mối liên hệ giữa phơi nhiễm và kết quả.

Ta thường lấy log của RR để tính khoảng tin cậy:

\[ CI_{RR} = \exp \left[ \ln(RR) \pm Z_{1 - \alpha/2} \cdot SE_{\ln(RR)} \right] \]

Trong đó:

\[ SE_{\ln(RR)} = \sqrt{ \frac{1}{a} - \frac{1}{a + b} + \frac{1}{c} - \frac{1}{c + d} } \]

Ví dụ trong kinh doanh:

Tỷ lệ khách mua hàng khi có khuyến mãi: $80/100 = 0.8$
Khi không có khuyến mãi: $40/100 = 0.4$
$RR = 0.8 / 0.4 = 2$

⇒ Khuyến mãi làm tăng gấp đôi khả năng khách mua hàng.

Hạn chế:

⚠️ Không sử dụng được trong nghiên cứu bệnh chứng (case-control) vì ta không biết nguy cơ tuyệt đối.

⚠️ Dễ gây hiểu nhầm nếu không đi kèm nguy cơ tuyệt đối. Ví dụ: RR = 2 có vẻ cao, nhưng nếu nguy cơ ban đầu là 1%, thì tăng lên 2% vẫn là rất thấp.

⚠️ Không đối xứng: Nếu đổi nhóm tham chiếu thì giá trị RR thay đổi, khác với OR.

5.3.3. Tỷ số chênh (Odds Ratio - OR)

5.3.3.1. Odds

Odds (tỷ số cơ hội):
\[ \text{Odds} = \frac{p}{1 - p} \] Trong đó, $p$ là xác suất xảy ra của một biến nhị phân (như việc sở hữu nhà).

5.3.3.1. Odds Ratio - OR

\[ OR = \frac{a / b}{c / d} = \frac{ad}{bc} \]

Hay

\[ \text{OR} = \frac{\text{Odds ở nhóm 1}}{\text{Odds ở nhóm 2}} \] OR thường dùng để đo lường mối liên hệ giữa một yếu tố và một kết quả trong bảng 2x2.

OR = 1: Không có sự khác biệt về odds giữa hai nhóm.
OR > 1: Nhóm 1 có odds xảy ra kết quả cao hơn nhóm 2.
OR < 1: Nhóm 1 có odds xảy ra kết quả thấp hơn nhóm 2.

Ví dụ: OR = 2 nghĩa là odds xảy ra kết quả ở nhóm 1 cao gấp 2 lần nhóm 2.

Lưu ý: OR khác với RR, đặc biệt khi biến phụ thuộc không hiếm gặp. Khi kết quả là hiếm (p nhỏ), thì OR ≈ RR.

Ý nghĩa:

OR cho biết tỷ lệ odds (tức là “xác suất chia cho 1 trừ xác suất”) giữa hai nhóm.
Thường được sử dụng phổ biến trong mô hình logistic regression.
Trong các nghiên cứu bệnh hiếm (rare disease assumption), OR xấp xỉ với RR.

Ưu điểm của Odds Ratio:

Dễ tính và dễ diễn giải.
Ổn định về mặt toán học khi sử dụng trong các mô hình hồi quy.
Không phụ thuộc vào tỷ lệ hiện diện trong mẫu (đặc biệt trong các thiết kế bệnh – chứng: case-control study).

5.4. Khoảng ước lượng tin cậy cho ODDS RATIO

Để xác định xem OR có ý nghĩa thống kê hay không, cần xây dựng khoảng tin cậy.

Logarithm tự nhiên của Odds Ratio:

\[ \log(OR) = \log\left( \frac{a \cdot d}{b \cdot c} \right) \]

Sai số chuẩn (Standard Error):

\[ SE = \sqrt{ \frac{1}{a} + \frac{1}{b} + \frac{1}{c} + \frac{1}{d} } \]

Khoảng tin cậy 95% cho $\log(OR)$:

\[ \log(OR) \pm 1.96 \cdot SE \]

Lấy mũ để có khoảng tin cậy 95% cho OR:

\[ CI_{95\%} = \left[ e^{\log(OR) - 1.96 \cdot SE},\quad e^{\log(OR) + 1.96 \cdot SE} \right] \]

Nếu khoảng tin cậy không chứa 1, ta có thể kết luận rằng mối liên hệ có ý nghĩa thống kê ở mức ý nghĩa 5%.

5.5. Ứng dụng trong kinh doanh: HIỆU QUẢ CỦA KHUYẾN MÃI

Một công ty thương mại điện tử thử nghiệm chiến dịch khuyến mãi để kiểm tra xem liệu có ảnh hưởng đến hành vi mua hàng không. Họ lấy mẫu 200 khách hàng ngẫu nhiên và thu được bảng sau:

	Mua hàng	Không mua	Tổng
Có khuyến mãi	80	20	100
Không khuyến mãi	40	60	100

Tính Odds Ratio (OR):

\[ OR = \frac{80 \cdot 60}{20 \cdot 40} = \frac{4800}{800} = 6 \]

→ Odds mua hàng khi có khuyến mãi cao gấp 6 lần so với khi không có khuyến mãi.

Logarithm tự nhiên của OR:

\[ \log(OR) = \log(6) \approx 1.79 \]

Sai số chuẩn (SE):

\[ SE = \sqrt{ \frac{1}{80} + \frac{1}{20} + \frac{1}{40} + \frac{1}{60} } \approx 0.35 \]

Khoảng tin cậy 95% cho OR:

\[ CI = \left[ e^{1.79 - 1.96 \cdot 0.35},\quad e^{1.79 + 1.96 \cdot 0.35} \right] \approx \left[ e^{1.11},\quad e^{2.47} \right] = [3.03,\ 11.81] \]

Vì khoảng tin cậy không chứa 1, có thể kết luận rằng khuyến mãi có ảnh hưởng có ý nghĩa thống kê đến hành vi mua hàng.
Với $OR = 6$ và $CI = (3.03,\ 11.81)$, doanh nghiệp có thể tự tin triển khai chiến lược khuyến mãi quy mô lớn hơn.

5.6 GENDER VÀ MARITALSTATUS

Tạo bảng chéo

# Tạo bảng chéo giữa Gender và MaritalStatus
table_gender_marital <- table(dat$Gender, dat$MaritalStatus)
addmargins(table_gender_marital)  # Thêm hàng và cột tổng

##      
##           M     S   Sum
##   F    3602  3568  7170
##   M    3264  3625  6889
##   Sum  6866  7193 14059

Tính toán Relative Risk

# Trích xuất số liệu từ bảng chéo
A <- table_gender_marital["F","S"] # nữ độc thân
B <- table_gender_marital["F","M"] # nữ kết hôn
C <- table_gender_marital["M","S"] # nam độc thân
D <- table_gender_marital["M","M"] # nam kết hôn

# Tính RR
RR <- (A/(A+B))/(C/(C+D))
RR

## [1] 0.945701

# Trong đó: A/(A+B) là tỷ lệ phụ nữ độc thân trong tổng số nữ
# Và: C/(C+D) là tỷ lệ nam giới độc thân trong tổng số nam

5.6.3 Ý nghĩa

RR = 1: Không có sự khác biệt về tỷ lệ độc thân giữa hai giới tính.
RR > 1: Nam giới có tỷ lệ độc thân cao hơn phụ nữ.
RR < 1: Phụ nữ có tỷ lệ độc thân cao hơn nam giới.
Tỷ lệ phụ nữ độc thân bằng khoảng 94.57% tỷ lệ nam giới độc thân. Điều này cho thấy không có sự khác biệt lớn về tỷ lệ độc thân giữa hai giới tính trong dữ liệu.

install.packages("epitools", repos = "https://cran.rstudio.com/")

## Installing package into 'C:/Users/Khanh Hoa/AppData/Local/R/win-library/4.4'
## (as 'lib' is unspecified)

## package 'epitools' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\Khanh Hoa\AppData\Local\Temp\Rtmpcb0VMf\downloaded_packages

library(epitools)
riskratio(table_gender_marital, method="wald")

## $data
##        
##            M    S Total
##   F     3602 3568  7170
##   M     3264 3625  6889
##   Total 6866 7193 14059
## 
## $measure
##    risk ratio with 95% C.I.
##     estimate    lower    upper
##   F 1.000000       NA       NA
##   M 1.057417 1.023813 1.092123
## 
## $p.value
##    two-sided
##       midp.exact fisher.exact   chi.square
##   F           NA           NA           NA
##   M 0.0007045092 0.0007373895 0.0007038837
## 
## $correction
## [1] FALSE
## 
## attr(,"method")
## [1] "Unconditional MLE & normal approximation (Wald) CI"

# Chỉ định phương pháp tính khoảng tin cậy (confidence interval) cho Risk Ratio.
# "wald" là phương pháp sử dụng ước lượng Wald, giả định phân phối chuẩn, để tính khoảng tin cậy.

Nhận xét

Relative Risk (RR)
- F (Female): 1.000000
- M (Male): 1.057417 (CI: 1.023813 - 1.092123)
Nhận xét về RR:
- RR của nhóm F (Phụ nữ) là 1 vì nhóm này được lấy làm nhóm tham chiếu.
- RR của nhóm M (Nam giới) là 1.0574, điều này cho thấy nam giới có tỷ lệ kết hôn (M) cao hơn một chút so với phụ nữ, với tỷ lệ 1.0574 lần so với nhóm phụ nữ. Sự khác biệt này có thể là không lớn, nhưng nó chỉ ra rằng nam giới có một khả năng kết hôn cao hơn phụ nữ trong mẫu nghiên cứu này.
- Confidence Interval (CI) của RR cho nhóm M là [1.0238, 1.0921], nghĩa là khoảng tin cậy của RR không bao gồm giá trị 1. Điều này cho thấy sự khác biệt giữa hai nhóm F và M là có ý nghĩa thống kê và không phải ngẫu nhiên.
Kết luận về RR:
- RR cho thấy nam giới có xu hướng kết hôn cao hơn so với phụ nữ, với RR = 1.0574.
- Confidence Interval không bao gồm giá trị 1, điều này hỗ trợ cho giả thuyết rằng sự khác biệt giữa tỷ lệ kết hôn của nam và nữ có ý nghĩa thống kê.
p-value từ các kiểm định:
- Chi-square test: p-value = 0.0007038837
- Fisher’s exact test: p-value = 0.0007373895
- Midpoint exact test: p-value = 0.0007045092
Nhận xét về p-value:
- Các p-value từ các kiểm định đều nhỏ hơn mức ý nghĩa thông thường là 0.05, nghĩa là có đủ bằng chứng để bác bỏ giả thuyết null hypothesis (giả thuyết không có sự khác biệt giữa nam và nữ).
- Điều này chỉ ra rằng có sự khác biệt có ý nghĩa thống kê giữa tỷ lệ kết hôn ở phụ nữ và nam giới.
Tóm tắt và nhận xét tổng quát
- Relative Risk (RR): Tỷ lệ kết hôn của nam giới là 1.0574 lần so với phụ nữ, điều này cho thấy nam giới có khả năng kết hôn cao hơn phụ nữ trong mẫu nghiên cứu này.
- Confidence Interval (C.I.): Khoảng tin cậy [1.0238, 1.0921] không bao gồm giá trị 1, điều này có nghĩa là sự khác biệt giữa tỷ lệ kết hôn của nam và nữ là có ý nghĩa thống kê.
- p-value: Các p-value từ Chi-square, Fisher’s exact test và midp.exact đều nhỏ hơn 0.05, cho thấy có sự khác biệt đáng kể giữa tỷ lệ kết hôn của nam và nữ.
Kết luận
- Tỷ lệ kết hôn giữa nam và nữ có sự khác biệt có ý nghĩa thống kê. Cụ thể, nam giới có tỷ lệ kết hôn cao hơn phụ nữ.
- Relative Risk cho thấy nam giới có khả năng kết hôn cao hơn phụ nữ (với RR = 1.0574), và kết quả này có ý nghĩa thống kê rõ ràng, vì p-value đều rất nhỏ.
- Khoảng tin cậy 95% của RR không bao gồm giá trị 1, điều này càng củng cố kết luận về sự khác biệt giữa tỷ lệ kết hôn của hai giới.

5.7 GENDER VÀ HOMEOWNER

Câu hỏi nghiên cứu: Liệu có sự khác biệt về tỷ lệ sở hữu nhà giữa nam và nữ không?

5.7.1 Tạo bảng chéo

# Tạo bảng chéo giữa Gender và Homeowner
table_gender_homeowner <- table(dat$Gender, dat$Homeowner)
addmargins(table_gender_homeowner)  # Thêm hàng và cột tổng

##      
##           N     Y   Sum
##   F    2826  4344  7170
##   M    2789  4100  6889
##   Sum  5615  8444 14059

5.7.2 Tính toán RR

# Trích xuất số liệu từ bảng chéo
A <- table_gender_homeowner["F", "Y"]  # nữ sở hữu nhà
B <- table_gender_homeowner["F", "N"]   # nữ không sở hữu nhà
C <- table_gender_homeowner["M", "Y"]  # nam sở hữu nhà
D <- table_gender_homeowner["M", "N"]   # nam không sở hữu nhà

# Tính RR
RR <- (A / (A + B)) / (C / (C + D))
RR

## [1] 1.017989

# Trong đó:
# A / (A + B): Tỷ lệ nữ sở hữu nhà trong tổng số nữ.
# C / (C + D): Tỷ lệ nam sở hữu nhà trong tổng số nam.

riskratio(table_gender_homeowner, method="wald")

## $data
##        
##            N    Y Total
##   F     2826 4344  7170
##   M     2789 4100  6889
##   Total 5615 8444 14059
## 
## $measure
##    risk ratio with 95% C.I.
##      estimate     lower    upper
##   F 1.0000000        NA       NA
##   M 0.9823291 0.9561812 1.009192
## 
## $p.value
##    two-sided
##     midp.exact fisher.exact chi.square
##   F         NA           NA         NA
##   M   0.195158    0.1964833  0.1950884
## 
## $correction
## [1] FALSE
## 
## attr(,"method")
## [1] "Unconditional MLE & normal approximation (Wald) CI"

5.7.3 Nhận xét

Nhóm nữ (F) được chọn làm nhóm tham chiếu, nên RR = 1. Điều này có nghĩa chúng ta so sánh nguy cơ sở hữu nhà của nhóm nam (M) với nhóm nữ.
Nhóm nam (M) có RR = 0.9823 nghĩa là:
- Tỷ lệ sở hữu nhà ở nam bằng khoảng 98.23% so với nữ. Hay nói cách khác, nguy cơ (xác suất) sở hữu nhà của nam thấp hơn nữ khoảng 1.77%.
Khoảng tin cậy 95% của RR là (0.9562 – 1.0092):
- Khoảng này bao gồm số 1, nghĩa là giá trị RR thực sự có thể là 1 (không có khác biệt).
- Khoảng tin cậy cho ta biết mức độ chính xác của ước lượng RR; khoảng rộng và bao gồm 1 cho thấy sự không chắc chắn và thiếu bằng chứng về sự khác biệt thực sự.
Ta có giả thuyết:
- H₀: Tỷ lệ sở hữu nhà ở Nam và Nữ không khác nhau, tức là RR = 1
- H₁: Tỷ lệ sở hữu nhà ở Nam và Nữ khác nhau, RR ≠ 1
p-value từ các kiểm định:
- Chi-square test: p-value = 0.1950884
- Fisher’s exact test: p-value = 0.1964833
- Midpoint exact test: p-value = 0.195158
Nhận xét về p-value:
- Các giá trị p-value từ chi-square, Fisher’s exact, và midp.exact đều lớn hơn mức ý nghĩa thông thường là 0.05.
- Điều này có nghĩa là không có đủ bằng chứng để bác bỏ giả thuyết null hypothesis, tức là không có sự khác biệt đáng kể về tỉ lệ tình trạng hôn nhân có (Y) giữa nhóm nam và nữ. Sự khác biệt giữa hai nhóm có thể là do ngẫu nhiên.

Nhận xét tổng quát

Relative Risk cho thấy rằng tỷ lệ sự kiện giữa hai nhóm giới tính F (phụ nữ) và M (nam giới) rất gần nhau, với RR của nhóm M là 0.9823, cho thấy không có sự khác biệt rõ rệt giữa tỷ lệ tình trạng hôn nhân có (Y) giữa hai nhóm.
Confidence Interval cho nhóm nam giới (M) bao gồm giá trị 1, điều này hỗ trợ giả thuyết rằng không có sự khác biệt thống kê đáng kể giữa tỷ lệ sự kiện ở nam và nữ.
Các kiểm định thống kê (Chi-square, Fisher’s exact test) đều cho p-value lớn hơn 0.05, điều này có nghĩa là không có sự khác biệt ý nghĩa giữa hai nhóm giới tính về tỷ lệ tình trạng hôn nhân có (Y).
Tóm lại, các kết quả này cho thấy không có sự khác biệt rõ rệt và có ý nghĩa giữa tỷ lệ tình trạng hôn nhân có (Y) ở phụ nữ và nam giới.

Kết luận

Dựa trên các kết quả từ Relative Risk và các kiểm định thống kê, chúng ta có thể kết luận rằng không có sự khác biệt đáng kể giữa nam và nữ về tỷ lệ tình trạng hôn nhân có (Y). Dữ liệu cho thấy rằng tỷ lệ hôn nhân có (Y) ở nam và nữ là khá tương đồng và không có sự khác biệt thống kê đáng kể.

5.8 MaritalStatus và Homeowner

table_mari_home <- table(dat$MaritalStatus, dat$Homeowner)
addmargins(table_mari_home)

##      
##           N     Y   Sum
##   M    1719  5147  6866
##   S    3896  3297  7193
##   Sum  5615  8444 14059

riskratio(table_mari_home)

## $data
##        
##            N    Y Total
##   M     1719 5147  6866
##   S     3896 3297  7193
##   Total 5615 8444 14059
## 
## $measure
##    risk ratio with 95% C.I.
##      estimate     lower     upper
##   M 1.0000000        NA        NA
##   S 0.6114466 0.5942071 0.6291862
## 
## $p.value
##    two-sided
##     midp.exact  fisher.exact    chi.square
##   M         NA            NA            NA
##   S          0 1.822183e-277 3.663022e-272
## 
## $correction
## [1] FALSE
## 
## attr(,"method")
## [1] "Unconditional MLE & normal approximation (Wald) CI"

Nhóm đã kết hôn (M) được dùng làm nhóm tham chiếu với RR = 1.

Nhóm độc thân (S) có:

RR = 0.6114 (khoảng tin cậy 95%: 0.5942 – 0.6292).Vì khoảng này không chứa 1, nên kết quả có ý nghĩa thống kê.
Điều này có nghĩa là người độc thân có nguy cơ sở hữu nhà thấp hơn khoảng 38.86% so với người đã kết hôn.
Khoảng tin cậy 95% khá hẹp và không bao gồm giá trị 1, cho thấy sự khác biệt này là rõ ràng và đáng tin cậy.

Ta thực hiện kiểm định giả thuyết như sau:

H₀: Tỷ lệ sở hữu nhà của người đã kết hôn và độc thân là như nhau (RR = 1).
H₁: Tỷ lệ sở hữu nhà của hai nhóm là khác nhau (RR ≠ 1).

Kết quả phân tích chỉ ra RR = 0.6114 với khoảng tin cậy 95% (0.5942 – 0.6292), và các kiểm định cho p-value ≈ 0 < mức ý nghĩa 0.05. Do đó, ta bác bỏ H₀ và kết luận rằng tình trạng hôn nhân có mối liên hệ thống kê rõ rệt với khả năng sở hữu nhà. Người độc thân có khả năng sở hữu nhà thấp hơn đáng kể so với người đã kết hôn.

PHẦN 6: ODDS RATIO

6.1 Gender và Homeowner

# Cài và gọi gói
install.packages("epitools")

## Warning: package 'epitools' is in use and will not be installed

library(epitools)

# Tính odds ratio
or_mh <- oddsratio(table_mari_home, method = "wald", conf.level = 0.95)
or_mh$measure

##    odds ratio with 95% C.I.
##      estimate     lower     upper
##   M 1.0000000        NA        NA
##   S 0.2826322 0.2630929 0.3036227

Kết quả phân tích Odds Ratio giữa giới tính và khả năng sở hữu nhà như sau:
- Odds Ratio (OR) = 0.956
- Khoảng tin cậy 95%: từ 0.894 đến 1.023
Diễn giải:
- Giá trị OR = 0.956 cho thấy odds (cơ hội) sở hữu nhà của nữ thấp hơn nam khoảng 4.4%. Tuy nhiên, sự khác biệt là nhỏ.
- Do khoảng tin cậy 95% chứa giá trị 1, nên:
  - Không có đủ bằng chứng thống kê để kết luận rằng giới tính ảnh hưởng đến odds sở hữu nhà.
Kết luận: Mặc dù nữ có vẻ có odds sở hữu nhà thấp hơn nam một chút, nhưng sự khác biệt không có ý nghĩa thống kê. Do đó, giới tính không phải là yếu tố ảnh hưởng rõ rệt đến khả năng sở hữu nhà trong dữ liệu hiện tại.

6.2 MaritalStatus và Homeowner

# Tính OR và khoảng tin cậy
or_mh <- oddsratio(table_mari_home, method = "wald", conf.level = 0.95)
or_mh$measure

##    odds ratio with 95% C.I.
##      estimate     lower     upper
##   M 1.0000000        NA        NA
##   S 0.2826322 0.2630929 0.3036227

Odds Ratio (OR) = 0.283
Khoảng tin cậy 95%: từ 0.263 đến 0.304
Diễn giải:
- Giá trị Odds Ratio 0.283 có nghĩa là odds sở hữu nhà của nhóm người độc thân chỉ bằng khoảng 28.3% odds của nhóm người đã kết hôn. Nói cách khác, người đã kết hôn có khả năng sở hữu nhà cao hơn gần 3.5 lần so với người độc thân (vì 1/0.283≈3.53).
- Khoảng tin cậy 95% cho Odds Ratio nằm hoàn toàn dưới 1, từ 0.263 đến 0.304, điều này cho thấy sự khác biệt về odds sở hữu nhà giữa hai nhóm là có ý nghĩa thống kê và không phải do ngẫu nhiên.
Kết luận:

Dữ liệu cho thấy tình trạng hôn nhân ảnh hưởng rõ ràng đến khả năng sở hữu nhà. Người đã kết hôn có odds sở hữu nhà cao hơn đáng kể so với người độc thân, và khác biệt này là có ý nghĩa thống kê với mức độ tin cậy 95%. Đây là bằng chứng cho thấy việc đã kết hôn có liên quan chặt chẽ đến khả năng sở hữu nhà trong mẫu khảo sát này.

6.3 Gender và MaritalStatus

# Tính OR và khoảng tin cậy
or_gm <- oddsratio(table_gender_marital, method = "wald", conf.level = 0.95)
or_gm$measure

##    odds ratio with 95% C.I.
##     estimate    lower    upper
##   F 1.000000       NA       NA
##   M 1.121184 1.049386 1.197893

Odds Ratio (OR) = 1.121
Khoảng tin cậy 95%: từ 1.049 đến 1.198
Diễn giải:
- Giá trị Odds Ratio 1.121 có nghĩa là odds (cơ hội) đã kết hôn của nữ cao hơn odds đã kết hôn của nam khoảng 12.1%, hay nói cách khác, nữ có khả năng đã kết hôn cao hơn nam khoảng 1.12 lần.
- Khoảng tin cậy 95% từ 1.049 đến 1.198, không chứa giá trị 1, chứng tỏ sự khác biệt này là có ý nghĩa thống kê. Điều này nghĩa là với mức tin cậy 95%, chúng ta chắc chắn rằng tỷ lệ đã kết hôn của nữ và nam không bằng nhau.
Kết luận:

Dữ liệu cho thấy nữ giới có odds đã kết hôn cao hơn nam giới khoảng 12%, tương đương với việc nữ có khả năng đã kết hôn gấp 1.12 lần so với nam giới.

LS0tDQp0aXRsZTogIk5ISeG7hk0gVuG7pCBUVeG6pk4gMiINCmF1dGhvcjogIlbDtSBIb8OgbmcgVMOtbmgiDQpkYXRlOiAiMjAyNS0wNS0yNSINCm91dHB1dDoNCiAgaHRtbF9kb2N1bWVudDogDQogICAgY29kZV9kb3dubG9hZDogdHJ1ZQ0KICAgIGNvZGVfZm9sZGluZzogaGlkZQ0KICAgIHRvY19kZXB0aDogNA0KICAgIHRvY19mbG9hdDogdHJ1ZQ0KICAgIHRvYzogdHJ1ZQ0KICBwZGZfZG9jdW1lbnQ6DQogICAgbGF0ZXhfZW5naW5lOiB4ZWxhdGV4DQotLS0NCg0KYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9DQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpDQpgYGANCg0KIyMgKipQSOG6pk4gMTogVMOMTSBISeG7glUgVsOAIENIVeG6qE4gQuG7iiBE4buuIExJ4buGVSoqDQoNCiMjIyAqKjEuMSDEkOG7jEMgROG7riBMSeG7hlUqKg0KDQoNCioqxJDhu41jIGZpbGUgZOG7ryBsaeG7h3UgZ+G7kWMqKg0KDQpgYGB7cn0NCmRhdGEgPC0gcmVhZC5jc3YoZmlsZSA9IGZpbGUuY2hvb3NlKCksIGhlYWRlciA9IFQpDQpgYGANCg0KKipDw6FjIGJp4bq/biDEkeG7i25oIHTDrW5oKioNCg0KLSBkYXQgPC0gZGF0YVssIGRsZHRdOiBs4bqleSB0cm9uZyBi4buZIGThu68gbGnhu4d1IGRhdGEgY8OhYyBj4buZdCB0cm9uZyBkbGR0IHbhu6thIHThuqFvIHJhIChjw6FjIGJp4bq/biDEkeG7i25oIHTDrW5oKSwgdsOgIGzhuqV5IHThuqV0IGPhuqMgY8OhYyBow6BuZyANCg0KYGBge3J9DQojIENo4buNbiBjw6FjIGJp4bq/biDEkeG7i25oIHTDrW5oDQpkbGR0IDwtIGMoIkdlbmRlciIsICJNYXJpdGFsU3RhdHVzIiwgIkhvbWVvd25lciIsICJBbm51YWxJbmNvbWUiLCAiQ2l0eSIsICJTdGF0ZW9yUHJvdmluY2UiLCAiQ291bnRyeSIsICJQcm9kdWN0RmFtaWx5IiwgIlByb2R1Y3REZXBhcnRtZW50IiwgIlByb2R1Y3RDYXRlZ29yeSIpDQoNCiMgVOG6oW8gYuG7mSBk4buvIGxp4buHdSBt4bubaSBjaOG7iSBjaOG7qWEgxJHhu4tuaCB0w61uaA0KZGF0IDwtIGRhdGFbLCBkbGR0XQ0KYGBgDQoNCg0KIyMjICoqMS4yLiBO4buYSSBEVU5HIEThu64gTEnhu4ZVKioNCg0KLSBC4buZIGThu68gbGnhu4d1IG7DoHkgZ+G7k20gMTQsMDU5IGdpYW8gZOG7i2NoIG11YSBow6BuZyB04bqhaSBzacOqdSB0aOG7iywgduG7m2kgY8OhYyB0aMO0bmcgdGluIGxpw6puIHF1YW4gxJHhur9uIGtow6FjaCBow6BuZyB2w6Agc+G6o24gcGjhuqltIGJhbyBn4buTbSAqKm3DoyDEkeG7i25oIGRhbmgga2jDoWNoIGjDoG5nIHbDoCBjw6FjIHRow7RuZyB0aW4gbGnDqm4gcXVhbioqLCAqKnbhu4sgdHLDrSDEkeG7i2EgbMO9KiosICoqY2hpIHRp4bq/dCBjw6FjIHPhuqNuIHBo4bqpbSDEkcaw4bujYyBtdWEqKi4NCg0KLSBE4buvIGxp4buHdSBuw6B5IGPDsyB0aOG7gyDEkcaw4bujYyBz4butIGThu6VuZyDEkeG7gyAqKnBow6JuIHTDrWNoIGjDoG5oIHZpIG11YSBz4bqvbSBj4bunYSBraMOhY2ggaMOgbmcqKiwgKipwaMOibiBraMO6YyB0aOG7iyB0csaw4budbmcqKiwgaG/hurdjICoqxJHDoW5oIGdpw6EgaGnhu4d1IHF14bqjIGtpbmggZG9hbmgqKiB0aGVvIHThu6tuZyBz4bqjbiBwaOG6qW0sIHThu6tuZyBraHUgduG7sWMgdsOgIHThu6tuZyBuaMOzbSDEkeG7kWkgdMaw4bujbmcuIA0KDQoqKkRhbmggc8OhY2ggY8OhYyBiaeG6v24gdsOgIG3DtCB04bqjKioNCg0KKipD4bqldSB0csO6YyBi4buZIGThu68gbGnhu4d1IFN1cGVybWFya2V0IFRyYW5zYWN0aW9ucyoqDQoNCmBgYHtyfQ0Kc3RyKGRhdGEpDQpgYGANCg0KICAtIEdlbmRlcjogR2nhu5tpIHTDrW5oIChGIC0gTuG7rywgTSAtIE5hbSkNCiAgLSBNYXJpdGFsU3RhdHVzOiBUw6xuaCB0cuG6oW5nIGjDtG4gbmjDom4gKFMgLSDEkOG7mWMgdGjDom4sIE0gLSDEkMOjIGvhur90IGjDtG4pDQogIC0gSG9tZW93bmVyOiBDw7Mgc+G7nyBo4buvdSBuaMOgIGhheSBraMO0bmcgKFkgLSBDw7MsIE4gLSBLaMO0bmcpDQogIC0gQW5udWFsSW5jb21lOiBUaHUgbmjhuq1wIGjDoG5nIG7Eg20gKMSRxrDhu6NjIGJp4buDdSB0aOG7iyBkxrDhu5tpIGThuqFuZyBjw6FjIGtob+G6o25nKQ0KICAtIENpdHksIFN0YXRlb3JQcm92aW5jZSwgQ291bnRyeTogVGjDtG5nIHRpbiDEkeG7i2EgbMO9DQogIC0gUHJvZHVjdEZhbWlseSwgUHJvZHVjdERlcGFydG1lbnQsIFByb2R1Y3RDYXRlZ29yeTogUGjDom4gbG/huqFpIHPhuqNuIHBo4bqpbQ0KICANCioqU+G7kSBiaeG6v24gdsOgIHPhu5EgcXVhbiBzw6F0KioNCg0KYGBge3J9DQpkaW0oZGF0KQ0KYGBgDQoNCi0gQuG7mSBk4buvIGxp4buHdSAqKmRhdCoqIGNo4bupYSB04buVbmcgY+G7mW5nICoqMTQuMDU5IHF1YW4gc8OhdCoqIGPDuW5nIHbhu5tpICoqMTAgYmnhur9uKiogxJHhurdjIHRyxrBuZyB0xrDGoW5nIOG7qW5nIGzDoCBiaeG6v24gxJHhu4tuaCB0w61uaC4NCi0gTeG7l2kgcXVhbiBzw6F0IHPhur0gdMaw4bufbmcg4bupbmcgduG7m2kgbeG7l2kgZ2lhbyBk4buLY2ggbXVhIGjDoG5nIHThuqFpIHNpw6p1IHRo4buLLg0KDQoNCioqS2nhu4NtIHRyYSBz4buRIGzGsOG7o25nIE5BIHRyb25nIGPhu5l0IMSR4buLbmggdMOtbmgqKg0KDQoNCnNhcHBseShkYXQsIC4uLik6IMOhcCBk4bulbmcgaMOgbSBjaG8gdOG7q25nIGPhu5l0IHRyb25nIGThu68gbGnhu4d1ICoqZGF0KiouDQoNCnN1bShpcy5uYSh4KSk6IMSR4bq/bSBz4buRIGzGsOG7o25nIE5BIHRyb25nIHThu6tuZyBj4buZdC4NCg0KYGBge3J9DQojIEtp4buDbSB0cmEgc+G7kSBsxrDhu6NuZyBOQSB0cm9uZyBt4buXaSBj4buZdCDEkeG7i25oIHTDrW5oDQpuYV9jb3VudHMgPC0gc2FwcGx5KGRhdCwgZnVuY3Rpb24oeCkgc3VtKGlzLm5hKHgpKSkNCm5hX2NvdW50cw0KYGBgDQoNCioqTmjhuq1uIHjDqXQqKiANCi0gS+G6v3QgcXXhuqMgdGh1IMSRxrDhu6NjIGNobyB0aOG6pXkgcuG6sW5nIHThuqV0IGPhuqMgY8OhYyBj4buZdCDEkeG7i25oIHTDrW5oIMSR4buBdSBraMO0bmcgY8OzIGdpw6EgdHLhu4sgdGhp4bq/dSAoTkEpIHbDrCBz4buRIGzGsOG7o25nIE5BIHRyb25nIGPDoWMgY+G7mXQgxJHhu4F1IGhp4buDbiB0aOG7iyBsw6AgMA0KDQojIyMjIyAqKkNodXnhu4NuIMSR4buVaSBjw6FjIGPhu5l0IHNhbmcgZmFjdG9yKioNCg0KDQpgYGB7cn0NCiMgQ2h1eeG7g24gxJHhu5VpIGPDoWMgY+G7mXQgc2FuZyBmYWN0b3IgbuG6v3UgY2jGsGEgcGjhuqNpIGzDoCBmYWN0b3INCmRhdCA8LSBhcy5kYXRhLmZyYW1lKGxhcHBseShkYXQsIGZ1bmN0aW9uKHgpIHsNCiAgaWYgKCFpcy5mYWN0b3IoeCkpIHsNCiAgICByZXR1cm4oYXMuZmFjdG9yKHgpKQ0KICB9IGVsc2Ugew0KICAgIHJldHVybih4KQ0KICB9DQp9KSkNCmBgYA0KICBL4bq/dCBxdeG6oyB0aHUgxJHGsOG7o2MgdG/DoG4gYuG7mSBjw6FjIGJp4bq/biDEkcOjIMSRxrDhu6NjIGNodXnhu4NuIMSR4buVaSB0aMOgbmggZOG6oW5nIGZhY3Rvci4NCg0KIyMgKipQSOG6pk4gMjogUEjDgk4gVMONQ0ggQknhur5OIMSQ4buKTkggVMONTkgqKg0KIyMjICoqUFJPRFVDVENBVEVHT1JZKioNCg0KKipUaOG7kW5nIGvDqiB0deG6p24gc3XhuqV0IHbDoCB0cuG7sWMgcXVhbiBow7NhKioNCg0KYGBge3J9DQp0YWJfUHJvZHVjdENhdGVnb3J5IDwtIHRhYmxlKGRhdCRQcm9kdWN0Q2F0ZWdvcnkpDQp0YWJfUHJvZHVjdENhdGVnb3J5DQp0YWJsZShkYXQkUHJvZHVjdENhdGVnb3J5KS9zdW0obnJvdyhkYXQpKQ0KYmFycGxvdCh0YWJfUHJvZHVjdENhdGVnb3J5LCBtYWluID0gIkRhbmggbeG7pWMgc+G6o24gcGjhuqltIC0gUHJvZHVjdENhdGVnb3J5IiwgY29sID0gIiNGRjAwMDAiLCBsYXMgPSAyKQ0KYGBgDQoNCioqTmjhuq1uIHjDqXQqKg0KDQotIEPDoWMgbG/huqFpIHRo4buxYyBwaOG6qW0gcGjhu5UgYmnhur9uIMSRxrDhu6NjIHF1YW4gc8OhdCB04burIGJp4buDdSDEkeG7kzoNCg0KICAtIFZlZ2V0YWJsZXMgKFJhdSk6IHbhu5tpIHPhu5EgbMaw4bujbmcgbMOgIGByIHRhYmxlKGRhdCRQcm9kdWN0Q2F0ZWdvcnkpWzQ1XWAga2jDoWNoIGjDoG5nIHbDoCBjaGnhur9tIGjhur90IChgciB0YWJsZShkYXQkUHJvZHVjdENhdGVnb3J5KVs0NV0vc3VtKG5yb3coZGF0KSkqMTAwYFwlKQ0KDQogIC0gU25hY2sgRm9vZHM6IGByIHRhYmxlKGRhdCRQcm9kdWN0Q2F0ZWdvcnkpWzQyXWAga2jDoWNoIGjDoG5nIHbDoCBjaGnhur9tIGjhur90IChgciB0YWJsZShkYXQkUHJvZHVjdENhdGVnb3J5KVs0Ml0vc3VtKG5yb3coZGF0KSkqMTAwYFwlKQ0KDQogIC0gRGFpcnkgKFPhu69hKTogYHIgdGFibGUoZGF0JFByb2R1Y3RDYXRlZ29yeSlbMThdYCBraMOhY2ggaMOgbmcgdsOgIGNoaeG6v20gaOG6v3QgKGByIHRhYmxlKGRhdCRQcm9kdWN0Q2F0ZWdvcnkpWzE4XS9zdW0obnJvdyhkYXQpKSoxMDBgXCUpDQoNCiAgLSBGcnVpdCwgTWVhdCwgRnJvemVuIEZvb2RzLi4uDQoNCi0gxJDDonkgbMOgIG5o4buvbmcgbeG6t3QgaMOgbmcgdGnDqnUgZMO5bmcgbmhhbmgsIGPDsyB2w7JuZyDEkeG7nWkgbmfhuq9uLg0KDQotIFZp4buHYyBraMOhY2ggaMOgbmcgbXVhIGNo4bunIHnhur91IGPDoWMgbeG6t3QgaMOgbmcgdGhp4bq/dCB54bq/dSBjaG8gdGjhuqV5IHNpw6p1IHRo4buLIGhv4bqhdCDEkeG7mW5nIHThu5F0IHRyb25nIHZp4buHYyBwaOG7pWMgduG7pSBuaHUgY+G6p3UgY8ahIGLhuqNuLCBjaOG7qSBraMO0bmcgcGjhuqNpIHPhuqNuIHBo4bqpbSB4YSB44buJIGhheSBraMO0bmcgdGjGsOG7nW5nIHh1ecOqbi4NCg0KIyMjICoqQU5OVUFMSU5DT01FKioNCg0KKipUaOG7kW5nIGvDqiB04bqnbiBzdeG6pXQgdsOgIHRy4buxYyBxdWFuIGjDs2EqKg0KDQpgYGB7cn0NCnRhYl9Bbm51YWxJbmNvbWUgPC0gdGFibGUoZGF0JEFubnVhbEluY29tZSkNCnRhYl9Bbm51YWxJbmNvbWUNCnRhYmxlKGRhdCRBbm51YWxJbmNvbWUpL3N1bShucm93KGRhdCkpDQpiYXJwbG90KHRhYl9Bbm51YWxJbmNvbWUsIG1haW4gPSAiVGh1IG5o4bqtcC90aMOhbmcgLSBBbm51YWxJbmNvbWUiLCBjb2wgPSAibGlnaHRibHVlIiwgbGFzID0gMikNCmBgYA0KDQoqKk5o4bqtbiB4w6l0KioNCg0KLSBOaMOzbSAkMzBLIC0gJDUwSyBjaGnhur9tIMSRYSBz4buRIHbhu5tpIGByIHRhYmxlKGRhdCRBbm51YWxJbmNvbWUpWzFdYCBraMOhY2ggaMOgbmcgIChjaGnhur9tIGByIHRhYmxlKGRhdCRBbm51YWxJbmNvbWUpWzVdL3N1bShucm93KGRhdCkpKjEwMGBcJSksIHRoZW8gc2F1IGzDoDoNCg0KICAtIE5ow7NtICQxMEsgLSAkMzBLIChgciB0YWJsZShkYXQkQW5udWFsSW5jb21lKVsxXS9zdW0obnJvdyhkYXQpKSoxMDBgXCUpDQoNCiAgLSBOaMOzbSAkNTBLIC0gJDcwSyAoYHIgdGFibGUoZGF0JEFubnVhbEluY29tZSlbNl0vc3VtKG5yb3coZGF0KSkqMTAwYFwlKQ0KDQotIFLhuqV0IMOtdCBraMOhY2ggaMOgbmcgY8OzIHRodSBuaOG6rXAgdHLDqm4gJDE1MEsgKGByIHRhYmxlKGRhdCRBbm51YWxJbmNvbWUpWzRdL3N1bShucm93KGRhdCkpKjEwMGBcJSkNCg0KLSDEkGnhu4F1IG7DoHkgY2hvIHRo4bqleSBy4bqxbmcgY2h14buXaSBzacOqdSB0aOG7iyBjaOG7pyB54bq/dSBwaOG7pWMgduG7pSBuaMOzbSB0aHUgbmjhuq1wIHRydW5nIGLDrG5oIHbDoCB0aOG6pXAg4oCTIHBow7kgaOG7o3AgduG7m2kgY2hp4bq/biBsxrDhu6NjIGdpw6EgcuG6uywgaMOgbmcgaMOzYSB0aGnhur90IHnhur91Lg0KDQotIENoaeG6v24gbMaw4bujYyBtYXJrZXRpbmcsIGtodXnhur9uIG3Do2ksIHbDoCBjaMOtbmggc8OhY2ggZ2nDoSBuw6puIHThuq1wIHRydW5nIHbDoG8gbmjDs20gbsOgeSAodsOtIGThu6U6IGNvbWJvIHRp4bq/dCBraeG7h20sIGNoaeG6v3Qga2jhuqV1IHRoZW8gc+G7kSBsxrDhu6NuZykuDQoNCiMjIyAqKkNPVU5UUlkqKg0KDQoqKlRo4buRbmcga8OqIHThuqduIHN14bqldCB2w6AgdHLhu7FjIHF1YW4gaMOzYSoqDQoNCmBgYHtyfQ0KdGFiX0NvdW50cnkgPC0gdGFibGUoZGF0JENvdW50cnkpDQp0YWJfQ291bnRyeQ0KdGFibGUoZGF0JENvdW50cnkpL3N1bShucm93KGRhdCkpDQpiYXJwbG90KHRhYl9Db3VudHJ5LCBtYWluID0gIlF14buRYyBnaWEgLSBDb3VudHJ5IiwgY29sID0gIm9yY2hpZCIsIGxhcyA9IDIpDQpgYGANCg0KKipOaOG6rW4geMOpdCoqDQoNCi0gVOG7tyBs4buHIGtow6FjaCBow6BuZyB0aGVvIHF14buRYyBnaWEgcGjDom4gYuG7lSB0aGVvIHRyw6xuaCB04buxIGdp4bqjbSBk4bqnbjoNCg0KICAtIENhbyBuaOG6pXQgbMOgIFVTQTogYHIgdGFibGUoZGF0JENvdW50cnkpWzNdL3N1bShucm93KGRhdCkpKjEwMGBcJQ0KDQogIC0gS+G6vyB0aeG6v3AgbMOgIE1leGljbzogYHIgdGFibGUoZGF0JENvdW50cnkpWzJdL3N1bShucm93KGRhdCkpKjEwMGBcJQ0KDQogIC0gVGjhuqVwIG5o4bqldCBsw6AgQ2FuYWRhOiBgciB0YWJsZShkYXQkQ291bnRyeSlbMV0vc3VtKG5yb3coZGF0KSkqMTAwYFwlDQogIA0KLSBVUyBIb2EgS+G7syByw7UgcsOgbmcgbMOgIHRo4buLIHRyxrDhu51uZyB0cuG7jW5nIMSRaeG7g20gduG7m2kgdOG7tyBs4buHIGtow6FjaCBow6BuZyBnaWFvIGThu4tjaCAgY2FvIG5o4bqldCB0cm9uZyAzIHF14buRYyBnaWEuIMSQaeG7gXUgbsOgeSBjaOG7qW5nIHThu48gSG9hIEvhu7MgbMOgIG3hu5l0IHF14buRYyBnaWEgcGjDoXQgdHJp4buDbiBt4bqhbmggduG7m2kgbmh1IGPhuqd1IHRpw6p1IGTDuW5nIGNoaeG6v24gcGjhuqduIGzhu5tuLiANCg0KLSBUdXkgbmhpw6puLCBz4buxIGhp4buHbiBkaeG7h24gdOG6oWkgTWV4aWNvIChoxqFuIDEvNCBz4buRIGdpYW8gZOG7i2NoKSBjaG8gdGjhuqV5IMSRw6J5IGzDoCB0aOG7iyB0csaw4budbmcgdGjhu6kgY+G6pXAgdGnhu4FtIG7Eg25nLCBjw7MgdGjhu4MgbeG7nyBy4buZbmcgaG/hurdjIHTEg25nIG5o4bqtbiBkaeG7h24gdGjGsMahbmcNCmhp4buHdS4NCg0KIyMjICoqR0VOREVSKioNCg0KKipUaOG7kW5nIGvDqiB04bqnbiBzdeG6pXQgdsOgIHRy4buxYyBxdWFuIGjDs2EqKg0KDQotIFPhu60gZOG7pW5nIGjDoG0gdGFibGUoKSDEkeG7gyDEkeG6v20gc+G7kSBs4bqnbiB4deG6pXQgaGnhu4duICh04bqnbiBzdeG6pXQpIGPhu6dhIHThu6tuZyBnacOhIHRy4buLIHRyb25nIGJp4bq/biDEkeG7i25oIHTDrW5oLg0KDQotIHByb3AudGFibGUodGFiX2dlbmRlcikgdOG6oW8gYuG6o25nIHThu7cgbOG7hyBwaOG6p24gdHLEg20gKHRoZW8gZOG6oW5nIHRo4bqtcCBwaMOibikuIE5ow6JuIHbhu5tpIDEwMCDEkeG7gyByYSBwaOG6p24gdHLEg20uDQoNCi0gcm91bmQoLi4uLCAyKSBsw6BtIHRyw7JuIDIgY2jhu68gc+G7kSBzYXUgZOG6pXUgcGjhuql5Lg0KDQotIGJhcnBsb3QoKSB04bqhbyBiaeG7g3UgxJHhu5MgY+G7mXQgdOG7qyBi4bqjbmcgdOG6p24gc3XhuqV0Lg0KDQotIG1haW46IHRpw6p1IMSR4buBIGJp4buDdSDEkeG7ky4NCg0KLSBjb2w6IG3DoHUgc+G6r2MgY+G7mXQgKOG7nyDEkcOieSBsw6AgbcOgdSAic3RlZWxibHVlIikuDQoNCi0gbGFzID0gMjogeG9heSBuaMOjbiB0cuG7pWMgaG/DoG5oIChYKSB0aGVvIGNoaeG7gXUgZOG7jWMgxJHhu4MgZOG7hSDEkeG7jWMuDQoNCmBgYHtyfQ0KdGFiX2dlbmRlciA8LSB0YWJsZShkYXQkR2VuZGVyKQ0KdGFiX2dlbmRlcg0KdGFibGUoZGF0JEdlbmRlcikvc3VtKG5yb3coZGF0KSkNCmJhcnBsb3QodGFiX2dlbmRlciwgbWFpbiA9ICJU4bqnbiBzdeG6pXQgLSBHZW5kZXIiLCBjb2wgPSAibGlnaHR5ZWxsb3ciLCBsYXMgPSAyKQ0KYGBgDQoNCioqTmjhuq1uIHjDqXQqKg0KDQotIEThu68gbGnhu4d1IGNobyB0aOG6pXkgY8OzIGByIHRhYmxlKGRhdCRHZW5kZXIpWzFdYCBraMOhY2ggaMOgbmcgbMOgIG7hu68gKGByIHRhYmxlKGRhdCRHZW5kZXIpWzFdL3N1bShucm93KGRhdCkpKjEwMGBcJSkgdsOgIGByIHRhYmxlKGRhdCRHZW5kZXIpWzJdYCBsw6Aga2jDoWNoIGjDoG5nIG5hbSAoYHIgdGFibGUoZGF0JEdlbmRlcilbMl0vc3VtKG5yb3coZGF0KSkqMTAwYFwlKS4NCg0KLSBU4bu3IGzhu4cgZ2nhu5tpIHTDrW5oIGdp4buvYSBuYW0gdsOgIG7hu68gdHJvbmcgZOG7ryBsaeG7h3Uga2jDoSBjw6JuIGLhurFuZywga2jDtG5nIGPDsyBz4buxIGtow6FjIGJp4buHdCDEkcOhbmcga+G7gyBuw6BvLg0KDQotIMSQaeG7gXUgbsOgeSBjaG8gdGjhuqV5IGNodeG7l2kgc2nDqnUgdGjhu4sgaMaw4bubbmcgdOG7m2kgcGjhu6VjIHbhu6UgY+G6oyBoYWkgZ2nhu5tpIG3hu5l0IGPDoWNoIMSR4buTbmcgxJHhu4F1LCB0aGF5IHbDrCB04bqtcCB0cnVuZyB2w6BvIG3hu5l0IG5ow7NtIGPhu6UgdGjhu4MuDQoNCi0gTeG6t3Qga2jDoWMsIHZp4buHYyBu4buvIGdp4bubaSBjaGnhur9tIHThu7cgbOG7hyBuaOG7iW5oIGjGoW4gxJHDtGkgY2jDunQgY8OzIHRo4buDIHh14bqldCBwaMOhdCB04burIHZhaSB0csOyIHRoxrDhu51uZyB0aOG6pXkgY+G7p2EgaOG7jSB0cm9uZyB2aeG7h2MgbXVhIHPhuq9tIG5odSB54bq/dSBwaOG6qW0gdsOgIHF14bqjbiBsw70gY2hpIHRpw6p1IGdpYSDEkcOsbmgsIGThuqtuIMSR4bq/biBraOG6oyBuxINuZyBjYW8gaOG7jSBsw6AgbmfGsOG7nWkgdGjhu7FjIGhp4buHbiBjw6FjIGdpYW8gZOG7i2NoIMSRxrDhu6NjIGdoaSBuaOG6rW4gdHJvbmcgaMOzYSDEkcahbi4NCg0KIyMjICoqTUFSSVRBTFNUQVRVUyoqDQoNCioqVGjhu5FuZyBrw6ogdHXhuqduIHN14bqldCB2w6AgdHLhu7FjIHF1YW4gaMOzYSoqDQoNCmBgYHtyfQ0KdGFiX01hcml0YWxTdGF0dXMgPC0gdGFibGUoZGF0JE1hcml0YWxTdGF0dXMpDQp0YWJfTWFyaXRhbFN0YXR1cw0KdGFibGUoZGF0JE1hcml0YWxTdGF0dXMpL3N1bShucm93KGRhdCkpDQpiYXJwbG90KHRhYl9NYXJpdGFsU3RhdHVzLCBtYWluID0gIlTDrG5oIHRy4bqhbmcgaMO0biBuaMOibiAtIE1hcml0YWxTdGF0dXMiLCBjb2wgPSAic2t5Ymx1ZSIsIGxhcyA9IDIpIA0KYGBgDQoNCioqTmjhuq1uIHjDqXQqKg0KDQotIFPhu5EgbMaw4bujbmcga2jDoWNoIGjDoG5nIMSR4buZYyB0aMOibiAoU2luZ2xlIOKAkyBrw70gaGnhu4d1IFMpIGzDoCBgciB0YWJsZShkYXQkTWFyaXRhbFN0YXR1cylbMV1gIG5nxrDhu51pIChjaGnhur9tIGByIHRhYmxlKGRhdCRNYXJpdGFsU3RhdHVzKVsxXS9zdW0obnJvdyhkYXQpKSoxMDBgXCUpLCBjw7JuIG5nxrDhu51pIMSRw6Mga+G6v3QgaMO0biAoTWFycmllZCDigJMgTSkgbMOgIGByIHRhYmxlKGRhdCRNYXJpdGFsU3RhdHVzWzJdKWAgbmfGsOG7nWkgKGNoaeG6v20gYHIgdGFibGUoZGF0JE1hcml0YWxTdGF0dXMpWzJdL3N1bShucm93KGRhdCkpKjEwMGBcJSkNCi0gTeG6t2MgZMO5IHThu7cgbOG7hyBraMOhIGPDom4gYuG6sW5nLCBuaMawbmcgbmjDs20gxJHhu5ljIHRow6JuIGNoaeG6v20gdOG7tyBs4buHIGNhbyBoxqFuIG3hu5l0IGNow7p0Lg0KDQotIEdp4bqjaSB0aMOtY2ggaOG7o3AgbMO9IGNobyDEkWnhu4F1IG7DoHkgY8OzIHRo4buDIGzDoDoNCg0KICAtIE5o4buvbmcgbmfGsOG7nWkgxJHhu5ljIHRow6JuIGPDsyB4dSBoxrDhu5tuZyB04buxIMSRaSBtdWEgc+G6r20gbmhp4buBdSBoxqFuLCBraMO0bmcgZOG7sWEgdsOgbyBi4bqhbiDEkeG7nWkuDQoNCiAgLSBOZ2/DoGkgcmEsIG5nxrDhu51pIHRy4bq7IMSR4buZYyB0aMOibiB0aMaw4budbmcgc+G7kW5nIOG7nyB0aMOgbmggdGjhu4sgdsOgIHPhu60gZOG7pW5nIGThu4tjaCB24bulIHNpw6p1IHRo4buLIHRoxrDhu51uZyB4dXnDqm4gaMahbiwgdHJvbmcga2hpIG5nxrDhu51pIMSRw6Mga+G6v3QgaMO0biBjw7MgdGjhu4MgbXVhIHRoZW8gbmjDs20gaG/hurdjIHRoZW8gaOG7mSBnaWEgxJHDrG5oIGzhu5tuIHbDoCDDrXQgZ2lhbyBk4buLY2ggaMahbiB0csOqbiBt4buXaSBjw6EgbmjDom4uDQogIA0KIyMjICoqSE9NRU9XTkVSKioNCg0KKipUaOG7kW5nIGvDqiB04bqnbiBzdeG6pXQgdsOgIHRy4buxYyBxdWFuIGjDs2EqKg0KDQpgYGB7cn0NCnRhYl9Ib21lb3duZXIgPC0gdGFibGUoZGF0JEhvbWVvd25lcikNCnRhYl9Ib21lb3duZXINCnRhYmxlKGRhdCRIb21lb3duZXIpL3N1bShucm93KGRhdCkpDQpiYXJwbG90KHRhYl9Ib21lb3duZXIsIG1haW4gPSAiU+G7nyBo4buvdSBuaMOgIC0gSG9tZW93bmVyIiwgY29sID0gImRhcmtzZWFncmVlbiIsIGxhcyA9IDIpDQpgYGANCg0KKipOaOG6rW4geMOpdCoqDQoNCi0gVOG7qyBr4bq/dCBxdeG6oyB0aHUgxJHGsOG7o2MgdGEgdGjhuqV5IGPDsyBgciB0YWJsZShkYXQkSG9tZW93bmVyKVsxXWAga2jDoWNoIGjDoG5nIChjaGnhur9tIGByIHRhYmxlKGRhdCRIb21lb3duZXIpWzFdL3N1bShucm93KGRhdCkpKjEwMGBcJSkga2jDtG5nIHPhu58gaOG7r3UgbmjDoCwgY8OybiBgciB0YWJsZShkYXQkSG9tZW93bmVyKVsyXWAga2jDoWNoIGjDoG5nIChjaGnhur9tIGByIHRhYmxlKGRhdCRIb21lb3duZXIpWzJdL3N1bShucm93KGRhdCkpKjEwMGBcJSkgbMOgIGNo4bunIHPhu58gaOG7r3UgbmjDoC4NCg0KLSDEkMOieSBsw6AgbeG7mXQgY2jhu4kgYsOhbyBxdWFuIHRy4buNbmcgduG7gSBz4buxIOG7lW4gxJHhu4tuaCB0w6BpIGNow61uaCBj4bunYSBraMOhY2ggaMOgbmcuDQoNCi0gTmfGsOG7nWkgc+G7nyBo4buvdSBuaMOgIHRoxrDhu51uZyBjw7M6DQoNCiAgLSBDaGkgdGnDqnUg4buVbiDEkeG7i25oIGjGoW4uDQoNCiAgLSBOaHUgY+G6p3UgbXVhIHPhuq9tIG5oaeG7gXUgaMahbiBjaG8gZ2lhIMSRw6xuaCAodGjhu7FjIHBo4bqpbSwgxJHhu5MgZ2lhIGThu6VuZykuDQoNCi0gVHJvbmcga2hpIMSRw7MsIG5nxrDhu51pIGtow7RuZyBz4bufIGjhu691IG5ow6AgKGPDsyB0aOG7gyB0aHXDqiBuaMOgIGhv4bq3YyBz4buRbmcgdOG6oW0gdGjhu51pKSBjw7MgdGjhu4M6DQoNCiAgLSBDaGkgdGnDqnUgZMOoIGThurd0IGjGoW4uDQoNCiAgLSDDjXQgbXVhIGjDoG5nIHPhu5EgbMaw4bujbmcgbOG7m24sIGhv4bq3YyDGsHUgdGnDqm4gaMOgbmcgbmjhu48gZ+G7jW4sIGThu4UgduG6rW4gY2h1eeG7g24uDQogIA0KIyMjICoqQ0lUWSoqDQoNCioqVGjhu5FuZyBrw6ogdOG6p24gc3XhuqV0IHbDoCB0cuG7sWMgcXVhbiBow7NhKioNCg0KYGBge3J9DQp0YWJfQ2l0eSA8LSB0YWJsZShkYXQkQ2l0eSkNCnRhYl9DaXR5DQp0YWJsZShkYXQkQ2l0eSkvc3VtKG5yb3coZGF0KSkNCmJhcnBsb3QodGFiX0NpdHksIG1haW4gPSAiVGjDoG5oIHBo4buRIC0gQ2l0eSIsIGNvbCA9ICJwaW5rIiwgbGFzID0gMikNCmBgYA0KDQoqKk5o4bqtbiB4w6l0KioNCg0KLSBDw6FjIHRow6BuaCBwaOG7kSBjw7MgdOG7tyBs4buHIGdpYW8gZOG7i2NoIGNhbyBsw6A6DQoNCiAgLSBTYWxlbSBsw6AgdGjDoG5oIHBo4buRIGPDsyB04buJIGzhu4cga2jDoWNoIGjDoG5nIG11YSBz4bqjbiBwaOG6qW0gY2FvIG5o4bqldCB24bubaSBgciB0YWJsZShkYXQkQ2l0eSlbMTNdYCBraMOhY2ggaMOgbmcgdsOgIGNoaeG6v20gaOG6v3QgKGByIHRhYmxlKGRhdCRDaXR5KVsxM10vc3VtKG5yb3coZGF0KSkqMTAwYFwlKQ0KDQogIC0gVGhlbyBzYXUgbMOgIHRow6BuaCBwaOG7kSBUYWNvbWEgKGByIHRhYmxlKGRhdCRDaXR5KVsyMV0vc3VtKG5yb3coZGF0KSkqMTAwYFwlKQ0KDQogIC0gS+G6vyB0aeG6v3AgbMOgIHRow6BuaCBwaOG7kSBMb3MgQW5nZWxlcyAoYHIgdGFibGUoZGF0JENpdHkpWzhdL3N1bShucm93KGRhdCkpKjEwMGBcJSkgdsOgIHRhIHRo4bqleSB0aMOgbmggcGjhu5EgU2VhdHRsZSB24bubaSBgciB0YWJsZShkYXQkQ2l0eSlbMTddL3N1bShucm93KGRhdCkpKjEwMGBcJSB44bqlcCB44buJIGLhurFuZyB24bubaSB0aMOgbmggcGjhu5EgTG9zIEFuZ2VsZXMuIFF1YSDEkcOzIHRhIG5ow6xuIG5o4bqtbiDEkcaw4bujYyBy4bqxbmcgMiB0aMOgbmggcGjhu5EgbsOgeSBjw7MgbMaw4bujbmcga2jDoWNoIGjDoG5nIG11YSBz4bqvbSBn4bqnbiBuaMawIG5oYXUuIA0KDQotIEPDsyB0aOG7gyBuaOG6rW4gdGjhuqV5IHPhu7EgdOG6rXAgdHJ1bmcgbeG6oW5oIHThuqFpIGPDoWMgdGjDoG5oIHBo4buRIGLhu50gVMOieSBIb2EgS+G7sywgxJHhurdjIGJp4buHdCDhu58gY8OhYyBiYW5nIFdhc2hpbmd0b24sIE9yZWdvbiB2w6AgQ2FsaWZvcm5pYS4NCg0KLSDEkGnhu4F1IG7DoHkgZ+G7o2kgw70gcuG6sW5nOg0KDQogIC0gQ8OhYyBjaGkgbmjDoW5oIOG7nyBjw6FjIHRow6BuaCBwaOG7kSBuw6B5IGPDsyBsxrB1IGzGsOG7o25nIGtow6FjaCBow6BuZyBs4bubbiwgaG/hurdjDQoNCiAgLSBDw6FjIHRow6BuaCBwaOG7kSBuw6B5IGPDsyBt4bqtdCDEkeG7mSBkw6JuIGPGsCBjYW8sIG5odSBj4bqndSB0acOqdSBkw7luZyB0aMaw4budbmcgeHV5w6puLg0KICANCiMjIyAqKlNUQVRFT1JQUk9WSU5DRSoqDQoNCioqVGjhu5FuZyBrw6ogdOG6p24gc3XhuqV0IHbDoCB0cuG7sWMgcXVhbiBow7NhKioNCg0KYGBge3J9DQp0YWJfU3RhdGVvclByb3ZpbmNlIDwtIHRhYmxlKGRhdCRTdGF0ZW9yUHJvdmluY2UpDQp0YWJfU3RhdGVvclByb3ZpbmNlDQp0YWJsZShkYXQkU3RhdGVvclByb3ZpbmNlKS9zdW0obnJvdyhkYXQpKQ0KYmFycGxvdCh0YWJfU3RhdGVvclByb3ZpbmNlLCBtYWluID0gIkJhbmcvVOG7iW5oIC0gU3RhdGVvclByb3ZpbmNlIiwgY29sID0gImdvbGRlbnJvZCIsIGxhcyA9IDIpDQpgYGANCg0KKipOaOG6rW4geMOpdCoqDQoNCi0gUGjDom4gYuG7kToNCg0KICAtIFdhc2hpbmd0b24oV0EpOiBgciB0YWJsZShkYXQkU3RhdGVvclByb3ZpbmNlKVs4XS9zdW0obnJvdyhkYXQpKSoxMDBgXCUNCg0KICAtIENhbGlmb3JuaWEoQ0EpOiBgciB0YWJsZShkYXQkU3RhdGVvclByb3ZpbmNlKVsyXS9zdW0obnJvdyhkYXQpKSoxMDBgXCUNCg0KICAtIE9yZWdvbihPUik6IGByIHRhYmxlKGRhdCRTdGF0ZW9yUHJvdmluY2UpWzZdL3N1bShucm93KGRhdCkpKjEwMGBcJQ0KDQotIEJhIGJhbmcgbsOgeSDEkeG7gXUg4bufIGLhu50gVMOieSBIb2EgS+G7sywgdsOgIGzDoCBuaOG7r25nIGJhbmcgY8OzOg0KDQogIC0gRMOibiBz4buRIMSRw7RuZywgdGjhu4sgdHLGsOG7nW5nIGLDoW4gbOG6uyBwaMOhdCB0cmnhu4NuLg0KDQogIC0gVGjDs2kgcXVlbiB0acOqdSBkw7luZyBoaeG7h24gxJHhuqFpLCBz4butIGThu6VuZyBk4buLY2ggduG7pSBzacOqdSB0aOG7iyBuaGnhu4F1IGjGoW4gc28gduG7m2kgY8OhYyBraHUgduG7sWMgbsO0bmcgdGjDtG4uDQoNCi0gVmnhu4djIG3hu58gcuG7mW5nIGtpbmggZG9hbmggaG/hurdjIMSR4bqndSB0xrAgdGjDqm0gdsOgbyBjw6FjIGJhbmcgbsOgeSBsw6AgY2hp4bq/biBsxrDhu6NjIGto4bqjIHRoaS4NCg0KDQojIyMgKipQUk9EVUNURkFNSUxZKioNCg0KKipUaOG7kW5nIGvDqiB04bqnbiBzdeG6pXQgdsOgIHRy4buxYyBxdWFuIGjDs2EqKg0KDQpgYGB7cn0NCnRhYl9Qcm9kdWN0RmFtaWx5IDwtIHRhYmxlKGRhdCRQcm9kdWN0RmFtaWx5KQ0KdGFiX1Byb2R1Y3RGYW1pbHkNCnRhYmxlKGRhdCRQcm9kdWN0RmFtaWx5KS9zdW0obnJvdyhkYXQpKQ0KYmFycGxvdCh0YWJfUHJvZHVjdEZhbWlseSwgbWFpbiA9ICJOaMOzbSBz4bqjbiBwaOG6qW0gLSBQcm9kdWN0RmFtaWx5IiwgY29sID0gImRhcmtibHVlIiwgbGFzID0gMikNCmBgYA0KDQoqKk5o4bqtbiB4w6l0KioNCg0KLSBUaOG7sWMgcGjhuqltIChGb29kKSBjaGnhur9tIMSRYSBz4buRIChgciB0YWJsZShkYXQkUHJvZHVjdEZhbWlseSlbMl0vc3VtKG5yb3coZGF0KSkqMTAwYFwlKSwgY2FvIHbGsOG7o3QgdHLhu5lpIHNvIHbhu5tpOg0KDQogIC0gTm9uLUNvbnN1bWFibGVzIGNoaeG6v20gYHIgdGFibGUoZGF0JFByb2R1Y3RGYW1pbHkpWzNdL3N1bShucm93KGRhdCkpKjEwMGBcJQ0KDQogIC0gTmjDs20gxJHhu5MgdeG7kW5nIGNoaeG6v20gYHIgdGFibGUoZGF0JFByb2R1Y3RGYW1pbHkpWzFdL3N1bShucm93KGRhdCkpKjEwMGBcJQ0KDQotIMSQaeG7gXUgbsOgeSBjaG8gdGjhuqV5IGtow6FjaCBow6BuZyBtdWEgc+G6r20gY2jhu6cgeeG6v3UgbMOgIMSR4buDIHBo4bulYyB24bulIG5odSBj4bqndSDEg24gdeG7kW5nIGjDoG5nIG5nw6B5Lg0KDQotIFNpw6p1IHRo4buLIG7Dqm4gdGnhur9wIHThu6VjIMSR4bqndSB0xrAgbeG6oW5oIHbDoG86DQoNCiAgLSBOZ3Xhu5NuIGN1bmcgdGjhu7FjIHBo4bqpbSB0xrDGoWksDQoNCiAgLSBDw6FjIHPhuqNuIHBo4bqpbSB0aGnhur90IHnhur91IGdpw6EgdOG7kXQsDQoNCiAgLSBDaMOtbmggc8OhY2ggxrB1IMSRw6NpIHPhu5EgbMaw4bujbmcgbOG7m24gY2hvIHRo4buxYyBwaOG6qW0uDQoNCiMjIyAqKlBST0RVQ1RERVBBUlRNRU5UKioNCg0KKipUaOG7kW5nIGvDqiB04bqnbiBzdeG6pXQgdsOgIHRy4buxYyBxdWFuIGjDs2EqKg0KDQpgYGB7cn0NCnRhYl9Qcm9kdWN0RGVwYXJ0bWVudCA8LSB0YWJsZShkYXQkUHJvZHVjdERlcGFydG1lbnQpDQp0YWJfUHJvZHVjdERlcGFydG1lbnQNCnRhYmxlKGRhdCRQcm9kdWN0RGVwYXJ0bWVudCkvc3VtKG5yb3coZGF0KSkNCmJhcnBsb3QodGFiX1Byb2R1Y3REZXBhcnRtZW50LCBtYWluID0gIkLhu5kgcGjhuq1uIHPhuqNuIHBo4bqpbSAtIFByb2R1Y3REZXBhcnRtZW50IiwgY29sID0gImRhcmtyZWQiLCBsYXMgPSAyKQ0KYGBgDQoNCioqTmjhuq1uIHjDqXQgKioNCg0KLSBDw6FjIGLhu5kgcGjDom4gY2hp4bq/bSDEkWEgc+G7kSB0cm9uZyB0b8OgbiBi4buZIGPDoWMgYuG7mSBwaOG6rW4gxJHGsOG7o2MgdGjhu5FuZyBrw6o6DQoNCiAgLSBQcm9kdWNlIChSYXUgcXXhuqMpOiBgciB0YWJsZShkYXQkUHJvZHVjdERlcGFydG1lbnQpWzE4XS9zdW0obnJvdyhkYXQpKSoxMDBgXCUNCg0KICAtIFNuYWNrIEZvb2RzOiBgciB0YWJsZShkYXQkUHJvZHVjdERlcGFydG1lbnQpWzIwXS9zdW0obnJvdyhkYXQpKSoxMDBgXCUNCg0KICAtIEhvdXNlaG9sZDogYHIgdGFibGUoZGF0JFByb2R1Y3REZXBhcnRtZW50KVsxNV0vc3VtKG5yb3coZGF0KSkqMTAwYFwlDQoNCi0gS2jDoWNoIGjDoG5nIG11YSBuaGnhu4F1IGPDoWMgbeG6t3QgaMOgbmc6DQoNCiAgLSBUxrDGoWkgc+G7kW5nLCBk4buFIGjhu49uZyAocHJvZHVjZSkgLT4gecOqdSBj4bqndSBi4bqjbyBxdeG6o24gdsOgIGNodeG7l2kgY3VuZyDhu6luZyB04buRdC4NCg0KICAtIMSQ4buTIMSDbiB24bq3dCB2w6AgaMOgbmcgZ2lhIGThu6VuZyAtPiBwaOG6o24gw6FuaCB0aMOzaSBxdWVuIHRpw6p1IGTDuW5nIG5oYW5oLCB0aeG7h24gbOG7o2kuDQoNCi0gxJDDonkgbMOgIGPGoSBo4buZaSDEkeG7gyB0xINuZyBjxrDhu51uZyBtYXJrZXRpbmcgbmjDs20gc+G6o24gcGjhuqltIHRp4buHbiBs4bujaSwgY29tYm8gYuG7r2EgxINuLCBz4bqjbiBwaOG6qW0gY2hvIGjhu5kgZ2lhIMSRw6xuaC4NCg0KDQojIyAqKlBI4bqmTiAzOiDGr+G7mkMgTMav4buiTkcgS0hP4bqiTkcgVsOAIEtJ4buCTSDEkOG7ik5IIEdJ4bqiIFRIVVnhur5UIENITyBU4bu2IEzhu4YqKg0KDQotIHRhYmxlKCk6CVRo4buRbmcga8OqIHPhu5EgbMaw4bujbmcgdOG7q25nIGxv4bqhaSB0cm9uZyBiaeG6v24gxJHhu4tuaCB0w61uaA0KLSBzdW0oLi4uKToJxJDhur9tIGJhbyBuaGnDqnUgZMOybmcgdGjhu49hIMSRaeG7gXUga2nhu4duIChnacOhIHRy4buLIHF1YW4gdMOibSkNCi0gbGVuZ3RoKC4uLik6CVThu5VuZyBz4buRIG3huqt1DQotIHByb3AudGVzdCguLi4pOglUw61uaCBraG/huqNuZyB0aW4gY+G6rXkgdsOgIGtp4buDbSDEkeG7i25oIGdp4bqjIHRodXnhur90IGNobyB04bu3IGzhu4cNCi0gcCA9IC4uLgk6IFThu7cgbOG7hyBnaeG6oyDEkeG7i25oIHRyb25nIGdp4bqjIHRodXnhur90IEgwDQotIGFsdGVybmF0aXZlID0gLi4uOglLaeG7g20gxJHhu4tuaCBoYWkgcGjDrWEgKCJ0d28uc2lkZWQiKSwgbeG7mXQgcGjDrWEgbOG7m24gaMahbiAoImdyZWF0ZXIiKSwgbmjhu48gaMahbiAoImxlc3MiKQ0KDQojIyMgKiozLjEgQmnhur9uIEdlbmRlcioqDQoNCmBgYHtyfQ0KIyBYw6FjIMSR4buLbmggaOG6oW5nIG3hu6VjIHF1YW4gdMOibQ0KdGFibGUoZGF0JEdlbmRlcikNCiMgU+G7kSBsxrDhu6NuZyB2w6AgdOG7lW5nIHPhu5EgDQpzbGZlbWFsZSA8LSBzdW0oZGF0JEdlbmRlciA9PSAiRiIpDQp0b3RhbGdlbmRlciA8LSBsZW5ndGgoZGF0JEdlbmRlcikNCiMgxq/hu5tjIGzGsOG7o25nIGtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIHThu7cgbOG7hyAiRmVtYWxlIg0KcHJvcC50ZXN0KHNsZmVtYWxlLCB0b3RhbGdlbmRlciwgY29uZi5sZXZlbCA9IDAuOTUpDQojIEtp4buDbSDEkeG7i25oIGdp4bqjIHRodXnhur90OiB04bu3IGzhu4cgRmVtYWxlID0gMC41DQpwcm9wLnRlc3Qoc2xmZW1hbGUsIHRvdGFsZ2VuZGVyLCBwID0gMC41LCBhbHRlcm5hdGl2ZSA9ICJ0d28uc2lkZWQiLCBjb25mLmxldmVsID0gMC45NSkNCmBgYA0KDQoqKkdp4bqjIHRodXnhur90KioNCg0KLSBIMDogcCA9IDAuNQ0KDQotIEgxOiBwIOKJoCAwLjUNCg0KKirEkOG7jWMga+G6v3QgcXXhuqMgdsOgIGvhur90IGx14bqtbiAqKg0KDQotIFThu5VuZyBz4buRIGPDoSB0aOG7gyB0cm9uZyBt4bqrdSBsw6AgMTQuMDU5IG5nxrDhu51pLCB0cm9uZyDEkcOzIGPDsyA3LjE3MCBuZ8aw4budaSBsw6AgbuG7ryB2w6AgNi44ODkgbmfGsOG7nWkgbMOgIG5hbS4gVOG7tyBs4buHIG7hu68gZ2nhu5tpIHRyb25nIG3huqt1IGzDoCA1MC45OTkzNiUuDQoNCi0gS+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmggdOG7tyBs4buHIG3hu5l0IG3huqt1IHbhu5tpIGdp4bqjIHRodXnhur90IHThu7cgbOG7hyBu4buvIHRyb25nIHThu5VuZyB0aOG7gyBsw6AgNTAlIChI4oKAOiBwID0gMC41KSBjaG8gdGjhuqV5IGdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQgPSA1LjU3NjUgduG7m2kgMSBi4bqtYyB04buxIGRvLCB2w6AgcC12YWx1ZSA9IDAuMDE4Mi4gVuG7m2kgbeG7qWMgw70gbmdoxKlhIHRow7RuZyB0aMaw4budbmcgzrEgPSAwLjA1LCB0YSBiw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3Qga2jDtG5nLCBjaG8gdGjhuqV5IHThu7cgbOG7hyBu4buvIGdp4bubaSB0cm9uZyBt4bqrdSBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiBzbyB24bubaSB04bu3IGzhu4cgNTAlLg0KDQotIEtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIHThu7cgbOG7hyBu4buvIHRyb25nIHF14bqnbiB0aOG7gyBu4bqxbSB0cm9uZyBraG/huqNuZyB04burIDUwLjE2OTMxJSDEkeG6v24gNTEuODI4ODYlLCBjaG8gdGjhuqV5IHThu7cgbOG7hyBuw6B5IGNo4buJIGjGoWkgY2FvIGjGoW4gNTAlLg0KDQoqKlTDs20gbOG6oWkqKg0KDQotIFR1eSB04bu3IGzhu4cgbuG7ryBnaeG7m2kgdHJvbmcgbeG6q3UgY2jhu4kgbmjhu4luaCBoxqFuIG5hbSBnaeG7m2kga2jDtG5nIHF1w6EgY2FvICAoNTAuOTk5MzYlKSwgbmjGsG5nIHPhu7EgY2jDqm5oIGzhu4djaCBuw6B5IGzDoCBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqLiBEbyDEkcOzLCBjw7MgdGjhu4Mga+G6v3QgbHXhuq1uIHLhurFuZyBu4buvIGdp4bubaSBjaGnhur9tIMawdSB0aOG6vyBuaOG6uSBuaMawbmcgxJHDoW5nIGvhu4MgduG7gSBt4bq3dCB0aOG7kW5nIGvDqiB0cm9uZyBt4bqrdSBraOG6o28gc8OhdCBuw6B5Lg0KDQoqKkvhur90IGx14bqtbioqDQoNCi0gVOG7tyBs4buHIG7hu68gdHJvbmcgdOG7lW5nIG3huqt1IGzDoCA1MC45OTklLCBraMOhYyBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqIHNvIHbhu5tpIDUwJS4NCi0gVuG7m2kgcC12YWx1ZSA9IDAuMDE4MiwgYuG6oW4gYsOhYyBi4buPIGdp4bqjIHRodXnhur90IHLhurFuZyB04bu3IGzhu4cgbuG7ryBsw6AgNTAlIOG7nyBt4bupYyDDvSBuZ2jEqWEgNSUuDQotIFR1eSBuaGnDqm4sIHPhu7EgY2jDqm5oIGzhu4djaCBkw7kgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiwgbmjGsG5nIHLhuqV0IG5o4buPIHbhu4EgbeG6t3QgdGjhu7FjIHThur8gKH4xJSkNCg0KIyMjICoqMy4yIEJp4bq/biBIb21lb3duZXIqKg0KDQpgYGB7cn0NCiMgWMOhYyDEkeG7i25oIGjhuqFuZyBt4bulYyBxdWFuIHTDom0NCnRhYmxlKGRhdCRIb21lb3duZXIpDQpzbHNobiA8LSBzdW0oZGF0JEhvbWVvd25lciA9PSAiWSIpDQp0b3RhbGhvbWVvd25lciA8LSBsZW5ndGgoZGF0JEhvbWVvd25lcikNCg0KIyDGr+G7m2MgbMaw4bujbmcga2hv4bqjbmcgdGluIGPhuq15DQpwcm9wLnRlc3Qoc2xzaG4sIHRvdGFsaG9tZW93bmVyLCBjb25mLmxldmVsID0gMC45NSkNCg0KIyBLaeG7g20gxJHhu4tuaCBnaeG6oyB0aHV54bq/dCBIMDogdOG7tyBs4buHIGNo4bunIG5ow6AgPSAwLjYNCnByb3AudGVzdChzbHNobiwgdG90YWxob21lb3duZXIsIHAgPSAwLjYsIGFsdGVybmF0aXZlID0gInR3by5zaWRlZCIsIGNvbmYubGV2ZWwgPSAwLjk1KQ0KYGBgDQoNCioqR2nhuqMgdGh1eeG6v3QgMSoqDQoNCi0gSDA6IHAgPSAwLjUNCg0KLSBIMTogcCDiiaAgMC41DQoNCioqxJDhu41jIGvhur90IHF14bqjIHbDoCBr4bq/dCBsdeG6rW4qKg0KDQotIEvhur90IHF14bqjIHBow6JuIHTDrWNoIHThu7cgbOG7hyBt4buZdCBt4bqrdSBjaG8gdGjhuqV5IHThu7cgbOG7hyDEkeG7kWkgdMaw4bujbmcgY2jhu41uIGzhu7FhIHRyb25nIHThu5VuZyB0aOG7gyBt4bqrdSBsw6AgNjAuMDYlLiBLaGkgdGnhur9uIGjDoG5oIGtp4buDbSDEkeG7i25oIGdp4bqjIHRodXnhur90IHbhu5tpIHThu7cgbOG7hyBjaHXhuqluIGzDoCA1MCUsIHRhIHRodSDEkcaw4bujYyBnacOhIHRy4buLIHRo4buRbmcga8OqIENoaS1zcXVhcmVkIGzDoCA1NjguODYgduG7m2kgbeG7mXQgYuG6rWMgdOG7sSBkbywgdsOgIHAtdmFsdWUgPCAyLjJlLTE2LiBW4bubaSBt4bupYyDDvSBuZ2jEqWEgdGjDtG5nIHRoxrDhu51uZyAozrEgPSAwLjA1KSwgdGEgYsOhYyBi4buPIGdp4bqjIHRodXnhur90IGtow7RuZyAoSOKCgDogcCA9IDAuNSkuIMSQaeG7gXUgbsOgeSBjaG8gdGjhuqV5IHThu7cgbOG7hyB0aOG7sWMgdOG6vyBraMOhYyBiaeG7h3QgbeG7mXQgY8OhY2ggY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiBzbyB24bubaSA1MCUuDQoNCioqR2nhuqMgdGh1eeG6v3QgMioqDQoNCi0gSDA6IHAgPSAwLjYNCg0KLSBIMTogcCDiiaAgMC42DQoNCioqxJDhu41jIGvhur90IHF14bqjIHbDoCBr4bq/dCBsdeG6rW4qKg0KDQotIE5nxrDhu6NjIGzhuqFpLCBraGkga2nhu4NtIMSR4buLbmggduG7m2kgdOG7tyBs4buHIGdp4bqjIMSR4buLbmggbMOgIDYwJSwgZ2nDoSB0cuG7iyB0aOG7kW5nIGvDqiBDaGktc3F1YXJlZCBjaOG7iSBsw6AgMC4wMTk0NDUgduG7m2kgcC12YWx1ZSA9IDAuODg5MS4gRG8gcC12YWx1ZSBs4bubbiBoxqFuIMSRw6FuZyBr4buDIHNvIHbhu5tpIG3hu6ljIMO9IG5naMSpYSAwLjA1LCB0YSBraMO0bmcgYsOhYyBi4buPIGdp4bqjIHRodXnhur90IGtow7RuZyAoSOKCgDogcCA9IDAuNikuIMSQaeG7gXUgbsOgeSBjaG8gdGjhuqV5IGtow7RuZyBjw7MgYuG6sW5nIGNo4bupbmcgdGjhu5FuZyBrw6ogY2hvIHRo4bqleSB04bu3IGzhu4cgdGjhu7FjIHThur8ga2jDoWMgYmnhu4d0IHNvIHbhu5tpIDYwJS4NCg0KKipUw7NtIGzhuqFpOioqDQoNCi0gVOG7tyBs4buHIGzhu7FhIGNo4buNbiB0cm9uZyBt4bqrdSBraOG6o28gc8OhdCBsw6Aga2hv4bqjbmcgNjAuMDYlIHbDoCBraMOhYyBiaeG7h3QgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiBzbyB24bubaSB04bu3IGzhu4cgZ2nhuqMgxJHhu4tuaCA1MCUsIG5oxrBuZyBraMO0bmcga2jDoWMgYmnhu4d0IMSRw6FuZyBr4buDIHNvIHbhu5tpIHThu7cgbOG7hyA2MCUuIMSQw6J5IGzDoCBjxqEgc+G7nyDEkeG7gyBnaeG6oyDEkeG7i25oIHLhurFuZyB04bu3IGzhu4cgbOG7sWEgY2jhu41uIHRyb25nIHF14bqnbiB0aOG7gyBjw7MgdGjhu4MgdsOgbyBraG/huqNuZyA2MCUuDQoNCiMjIyAqKjMuMyBCaeG6v24gUHJvZHVjdEZhbWlseSoqDQoNCmBgYHtyfQ0KIyBYw6FjIMSR4buLbmggaOG6oW5nIG3hu6VjIMSR4buBIGPhuq1wDQp0YWJsZShkYXQkUHJvZHVjdEZhbWlseSkNCnNsZm9vZCA8LSBzdW0oZGF0JFByb2R1Y3RGYW1pbHkgPT0gIkZvb2QiKQ0KdG90YWxwcm9kdWN0IDwtIGxlbmd0aChkYXQkUHJvZHVjdEZhbWlseSkNCiMgxq/hu5tjIGzGsOG7o25nIGtob+G6o25nIHRpbiBj4bqteQ0KcHJvcC50ZXN0KHNsZm9vZCwgdG90YWxwcm9kdWN0LCBjb25mLmxldmVsID0gMC45NSkNCiMgS2nhu4NtIMSR4buLbmggZ2nhuqMgdGh1eeG6v3QgSDA6IHThu7cgbOG7hyA+PSAwLjcgKGFsdGVybmF0aXZlID0gImdyZWF0ZXIiOiBraeG7g20gxJHhu4tuaCAxIHBow61hICg+KSkNCnByb3AudGVzdChzbGZvb2QsIHRvdGFscHJvZHVjdCwgcCA9IDAuNywgYWx0ZXJuYXRpdmUgPSAiZ3JlYXRlciIsIGNvbmYubGV2ZWwgPSAwLjk1KQ0KDQoNCmBgYA0KDQoqKkdp4bqjIHRodXnhur90IDEqKg0KDQotIEgwOiBwID0gMC41DQoNCi0gSDE6IHAg4omgIDAuNQ0KDQoqKsSQ4buNYyBr4bq/dCBxdeG6oyB2w6Aga+G6v3QgbHXhuq1uKioNCg0KLSBL4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaCB04bu3IGzhu4cgbeG7mXQgbeG6q3UgY2hvIGdp4bqjIHRodXnhur90IHLhurFuZyB04bu3IGzhu4cgc+G6o24gcGjhuqltIHRodeG7mWMgbmjDs20gRm9vZCBi4bqxbmcgNTAlIChI4oKAOiBwID0gMC41KSBjaG8gdGjhuqV5Og0KDQogIC0gR2nDoSB0cuG7iyB0aOG7kW5nIGvDqiBDaGktc3F1YXJlZCA9IDI3NzQuOSwgZGYgPSAxDQoNCiAgLSBwLXZhbHVlIDwgMi4yZS0xNg0KDQotIFbhu5tpIG3hu6ljIMO9IG5naMSpYSBsw6AgMC4wNSwgdGEgYsOhYyBi4buPIGdp4bqjIHRodXnhur90IEgwLCBuZ2jEqWEgbMOgIHThu7cgbOG7hyBz4bqjbiBwaOG6qW0gbmjDs20gRm9vZCBraMO0bmcgY2hp4bq/bSA1MCUgdHJvbmcgdOG7lW5nIHPhu5EgIlByb2R1Y3RGYW1pbHkiLg0KDQotIEtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIHThu7cgbOG7hyBuw6B5IG7hurFtIHRyb25nIGtob+G6o25nIHThu6sgNzEuNDclIMSR4bq/biA3Mi45NSUuDQoNCioqR2nhuqMgdGh1eeG6v3QgMioqDQoNCi0gSDA6IHAgPSAwLjcNCg0KLSBIMTogcCA+IDAuNw0KDQoqKsSQ4buNYyBr4bq/dCBxdeG6oyB2w6Aga+G6v3QgbHXhuq1uKioNCg0KLSBL4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaCB0aeG6v3AgdGhlbyBraeG7g20gdHJhIHhlbSBsaeG7h3UgdOG7tyBs4buHIHPhuqNuIHBo4bqpbSBuaMOzbSBGb29kIGPDsyBjYW8gaMahbiA3MCUgaGF5IGtow7RuZywgY2hvIHRo4bqleToNCg0KICAtIEdpw6EgdHLhu4sgQ2hpLXNxdWFyZWQgPSAzMi44MDIsIGRmID0gMQ0KDQogIC0gcC12YWx1ZSA9IDUuMTAxZS0wOQ0KDQotIFbhu5tpIG3hu6ljIMO9IG5naMSpYSAwLjA1LCB0YSB0aeG6v3AgdOG7pWMgYsOhYyBi4buPIGdp4bqjIHRodXnhur90IGtow7RuZywgdsOgIGNo4bqlcCBuaOG6rW4gZ2nhuqMgdGh1eeG6v3QgxJHhu5FpIChI4oKBOiBwID4gMC43KS4gxJBp4buBdSBuw6B5IGNo4bupbmcgdOG7jyB04bu3IGzhu4cgc+G6o24gcGjhuqltIG5ow7NtIEZvb2QgbOG7m24gaMahbiA3MCUgduG7m2kgw70gbmdoxKlhIHRo4buRbmcga8OqIGNhby4NCg0KLSBLaG/huqNuZyB0aW4gY+G6rXkgOTUlIGNobyB04bu3IGzhu4cgbsOgeSBu4bqxbSB0cm9uZyBraG/huqNuZyB04burIDcxLjU4JSDEkeG6v24gMTAwJS4NCg0KDQojIyAqKlBI4bqmTiA0OiBQSMOCTiBUw41DSCBN4buQSSBRVUFOIEjhu4YgR0nhu65BIEhBSSBCSeG6vk4gxJDhu4pOSCBUw41OSCoqDQoNCi0gQ2jhu41uIGzhuqduIGzGsOG7o3QgY8OhYyBj4bq3cCBiaeG6v24gY8OzIG3hu5FpIHTGsMahbmcgcXVhbiB0aMO6IHbhu4sgbMOgbSDhuqNuaCBoxrDhu59uZyDEkeG6v24gaMOgbmggdmkgbXVhIGjDoG5nLCBnaWFvIGThu4tjaCBj4bunYSBraMOhY2ggaMOgbmcuDQoNCg0KYGBge3J9DQppbnN0YWxsLnBhY2thZ2VzKCJ0aWR5dmVyc2UiLCByZXBvcyA9ICJodHRwczovL2NyYW4uci1wcm9qZWN0Lm9yZyIpDQppbnN0YWxsLnBhY2thZ2VzKCJqYW5pdG9yIiwgcmVwb3MgPSAiaHR0cHM6Ly9jcmFuLnItcHJvamVjdC5vcmciKQ0KaW5zdGFsbC5wYWNrYWdlcygiZ21vZGVscyIsIHJlcG9zID0gImh0dHBzOi8vY3Jhbi5yLXByb2plY3Qub3JnIikNCmluc3RhbGwucGFja2FnZXMoImdncGxvdDIiLCByZXBvcyA9ICJodHRwczovL2NyYW4uci1wcm9qZWN0Lm9yZyIpDQppbnN0YWxsLnBhY2thZ2VzKCJzY2FsZXMiLCByZXBvcyA9ICJodHRwczovL2NyYW4uci1wcm9qZWN0Lm9yZyIpDQpsaWJyYXJ5KHRpZHl2ZXJzZSkNCmxpYnJhcnkoamFuaXRvcikNCmxpYnJhcnkoZ21vZGVscykNCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkoc2NhbGVzKQ0KYGBgDQoNCg0KIyMjICoqNC4xIEdlbmRlciB2w6AgUHJvZHVjdEZhbWlseSoqDQoNCioqQuG6o25nIHThuqduIHN14bqldCBjaMOpbyoqDQoNCmBgYHtyfQ0KdGFiMSA8LSB0YWJsZShkYXQkR2VuZGVyLCBkYXQkUHJvZHVjdEZhbWlseSkNCnRhYjEgDQpwcm9wLnRhYmxlKHRhYjEsIDEpICMgdGhlbyBow6BuZw0KYGBgDQoNCg0KLSBE4buxYSB0csOqbiBi4bqjbmcgdOG6p24gc3XhuqV0IGNow6lvIGdp4buvYSBHaeG7m2kgdMOtbmggKEdlbmRlcikgdsOgIE5ow7NtIHPhuqNuIHBo4bqpbSAoUHJvZHVjdEZhbWlseSk6DQoNCi0gTuG7ryBnaeG7m2kgKEYpIGPDsyB4dSBoxrDhu5tuZyBjaOG7jW4gc+G6o24gcGjhuqltIEZvb2Qgbmhp4buBdSBuaOG6pXQsIGNoaeG6v20ga2hv4bqjbmcgNzEuOCUgdHJvbmcgdOG7lW5nIHPhu5Egc+G6o24gcGjhuqltIG3DoCBu4buvIGdp4bubaSBjaOG7jW4uIFRp4bq/cCB0aGVvIGzDoCBuaMOzbSBOb24tQ29uc3VtYWJsZSAoMTguOSUpIHbDoCBEcmluayAoOS4zJSkuDQoNCi0gTmFtIGdp4bubaSAoTSkgY8WpbmcgY8OzIHh1IGjGsOG7m25nIHTGsMahbmcgdOG7sSBraGkgbOG7sWEgY2jhu41uIG5oaeG7gXUgbmjhuqV0IGzDoCBGb29kICg3Mi42JSksIGvhur8gxJHhur9uIGzDoCBOb24tQ29uc3VtYWJsZSAoMTguOSUpIHbDoCDDrXQgbmjhuqV0IGzDoCBEcmluayAoOC40JSkuDQoNCi0gU28gc8OhbmggZ2nhu5tpIHTDrW5oOg0KDQogIC0gQ+G6oyBoYWkgZ2nhu5tpIMSR4buBdSDGsHUgdGnDqm4gbmjDs20gc+G6o24gcGjhuqltIEZvb2QgduG7m2kgdOG7tyBs4buHIGNhbyBuaOG6pXQsIHRyw6puIDcwJS4NCg0KICAtIFThu7cgbOG7hyBjaOG7jW4gRHJpbmsgY+G7p2EgbuG7ryBjYW8gaMahbiBuYW0gbeG7mXQgY2jDunQgKDkuMyUgc28gduG7m2kgOC40JSkuDQoNCiAgLSBU4bu3IGzhu4cgY2jhu41uIE5vbi1Db25zdW1hYmxlIGfhuqduIG5oxrAgdMawxqFuZyDEkcawxqFuZyBnaeG7r2EgaGFpIGdp4bubaSAoMTguOSUpLg0KDQoqKkvhur90IGx14bqtbjoqKg0KDQotIFPhu58gdGjDrWNoIHbhu4EgbmjDs20gc+G6o24gcGjhuqltIGtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgxJHDoW5nIGvhu4MgZ2nhu69hIG5hbSB2w6AgbuG7ry4gTmjDs20gc+G6o24gcGjhuqltIEZvb2QgY2hp4bq/bSDGsHUgdGjhur8gdHJvbmcgbOG7sWEgY2jhu41uIGPhu6dhIGPhuqMgaGFpIGdp4bubaSB0w61uaCwgdHJvbmcga2hpIG5ow7NtIERyaW5rIGzDoCDDrXQgxJHGsOG7o2MgY2jhu41uIG5o4bqldC4NCg0KKipUcuG7sWMgcXVhbiBow7NhKioNCg0KYGBge3J9DQpnZ3Bsb3QoZGF0LCBhZXMoeCA9IEdlbmRlciwgZmlsbCA9IFByb2R1Y3RGYW1pbHkpKSArDQogIGdlb21fYmFyKHBvc2l0aW9uID0gImZpbGwiKSArDQogIGxhYnModGl0bGUgPSAiVOG7tyBs4buHIFByb2R1Y3RGYW1pbHkgdGhlbyAgR2VuZGVyIiwgeSA9ICJU4bu3IGzhu4ciLCB4ID0gIkdlbmRlciIpICsNCiAgc2NhbGVfeV9jb250aW51b3VzKGxhYmVscyA9IHNjYWxlczo6cGVyY2VudCkgKw0KICB0aGVtZV9taW5pbWFsKCkNCmBgYA0KDQoNCioqTmjhuq1uIHjDqXQqKg0KDQotIEThu7FhIHRyw6puIGLhuqNuZyB04bqnbiBzdeG6pXQgY2jDqW8gZ2nhu69hIEdp4bubaSB0w61uaCAoR2VuZGVyKSB2w6AgTmjDs20gc+G6o24gcGjhuqltIChQcm9kdWN0RmFtaWx5KToNCg0KLSBO4buvIGdp4bubaSAoRikgY8OzIHh1IGjGsOG7m25nIGNo4buNbiBz4bqjbiBwaOG6qW0gRm9vZCBuaGnhu4F1IG5o4bqldCwgY2hp4bq/bSBraG/huqNuZyA3MS44JSB0cm9uZyB04buVbmcgc+G7kSBz4bqjbiBwaOG6qW0gbcOgIG7hu68gZ2nhu5tpIGNo4buNbi4gVGnhur9wIHRoZW8gbMOgIG5ow7NtIE5vbi1Db25zdW1hYmxlICgxOC45JSkgdsOgIERyaW5rICg5LjMlKS4NCg0KLSBOYW0gZ2nhu5tpIChNKSBjxaluZyBjw7MgeHUgaMaw4bubbmcgdMawxqFuZyB04buxIGtoaSBs4buxYSBjaOG7jW4gbmhp4buBdSBuaOG6pXQgbMOgIEZvb2QgKDcyLjYlKSwga+G6vyDEkeG6v24gbMOgIE5vbi1Db25zdW1hYmxlICgxOC45JSkgdsOgIMOtdCBuaOG6pXQgbMOgIERyaW5rICg4LjQlKS4NCg0KLSBTbyBzw6FuaCBnaeG7m2kgdMOtbmg6DQoNCiAgLSBD4bqjIGhhaSBnaeG7m2kgxJHhu4F1IMawdSB0acOqbiBuaMOzbSBz4bqjbiBwaOG6qW0gRm9vZCB24bubaSB04bu3IGzhu4cgY2FvIG5o4bqldCwgdHLDqm4gNzAlLg0KDQogIC0gVOG7tyBs4buHIGNo4buNbiBEcmluayBj4bunYSBu4buvIGNhbyBoxqFuIG5hbSBt4buZdCBjaMO6dCAoOS4zJSBzbyB24bubaSA4LjQlKS4NCg0KICAtIFThu7cgbOG7hyBjaOG7jW4gTm9uLUNvbnN1bWFibGUgZ+G6p24gbmjGsCB0xrDGoW5nIMSRxrDGoW5nIGdp4buvYSBoYWkgZ2nhu5tpICgxOC45JSkuDQoNCioqS+G6v3QgbHXhuq1uOioqDQoNCi0gU+G7nyB0aMOtY2ggduG7gSBuaMOzbSBz4bqjbiBwaOG6qW0ga2jDtG5nIGPDsyBz4buxIGtow6FjIGJp4buHdCDEkcOhbmcga+G7gyBnaeG7r2EgbmFtIHbDoCBu4buvLiBOaMOzbSBz4bqjbiBwaOG6qW0gRm9vZCBjaGnhur9tIMawdSB0aOG6vyB0cm9uZyBs4buxYSBjaOG7jW4gY+G7p2EgY+G6oyBoYWkgZ2nhu5tpLCB0cm9uZyBraGkgbmjDs20gRHJpbmsgbMOgIMOtdCDEkcaw4bujYyBjaOG7jW4gbmjhuqV0Lg0KDQoqKktp4buDbSDEkeG7i25oIENoaS1iw6xuaCBwaMawxqFuZyoqDQoNCmBgYHtyfQ0KY2hpX3Jlc3VsdDEgPC0gY2hpc3EudGVzdCh0YWIxKQ0KY2hpX3Jlc3VsdDENCg0KIyBHacOhIHRy4buLIHRo4buRbmcga8OqIENoaS1zcXVhcmVkDQpjaGlfcmVzdWx0MSRzdGF0aXN0aWMNCiMgQuG6rWMgdOG7sSBkbyANCmNoaV9yZXN1bHQxJHBhcmFtZXRlcg0KIyBHacOhIHRy4buLIHANCmNoaV9yZXN1bHQxJHAudmFsdWUNCmBgYA0KDQoqKkdp4bqjIHRodXnhur90KioNCg0KLSBHaeG6oyB0aHV54bq/dCBIMDogR2nhu5tpIHTDrW5oIChHZW5kZXIpIHbDoCBOaMOzbSBz4bqjbiBwaOG6qW0gKFByb2R1Y3RGYW1pbHkpIGzDoCDEkeG7mWMgbOG6rXAsIGtow7RuZyBjw7MgbeG7kWkgbGnDqm4gaOG7hy4NCg0KLSBHaeG6oyB0aHV54bq/dCBIMTogR2nhu5tpIHTDrW5oIHbDoCBOaMOzbSBz4bqjbiBwaOG6qW0gY8OzIGxpw6puIGjhu4cgduG7m2kgbmhhdS4NCg0KKipL4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaDoqKg0KDQotIEdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQ6IFgtc3F1YXJlZCA9IDMuNTE4NDkNCg0KLSBC4bqtYyB04buxIGRvOiBkZiA9IDINCg0KLSBHacOhIHRy4buLIHA6IHAgPSAwLjE3MjE3NDgNCg0KKipL4bq/dCBsdeG6rW46KioNCg0KLSBW4bubaSBt4bupYyDDvSBuZ2jEqWEgMC4wNSwgcF92YWx1ZSA9IDAuMTcyMiA+PSAwLjA1IOKGkiBraMO0bmcgxJHhu6cgYuG6sW5nIGNo4bupbmcgxJHhu4MgYsOhYyBi4buPIEgwLCB04bupYyBsw6Aga2jDtG5nIGPDsyBt4buRaSBsacOqbiBo4buHIGdp4buvYSAyIGJp4bq/biBHZW5kZXIgdsOgIFByb2R1Y3RGYW1pbHkuDQoNCg0KIyMjICoqNC4yIE1hcml0YWxTdGF0dXMgdsOgIEhvbWVvd25lcioqDQoNCioqQuG6o25nIHThuqduIHN14bqldCBjaMOpbyoqDQoNCmBgYHtyfQ0KdGFiMiA8LSB0YWJsZShkYXQkTWFyaXRhbFN0YXR1cywgZGF0JEhvbWVvd25lcikNCnRhYjIgDQpwcm9wLnRhYmxlKHRhYjIsIDEpDQpgYGANCg0KKipUcuG7sWMgcXVhbiBow7NhKioNCg0KYGBge3J9DQpnZ3Bsb3QoZGF0LCBhZXMoeCA9IE1hcml0YWxTdGF0dXMsIGZpbGwgPSBIb21lb3duZXIpKSArDQogIGdlb21fYmFyKHBvc2l0aW9uID0gImZpbGwiKSArDQogIGxhYnModGl0bGUgPSAiVOG7tyBs4buHIEhvbWVvd25lciB0aGVvIE1hcml0YWwgU3RhdHVzIiwgeSA9ICJU4bu3IGzhu4ciLCB4ID0gIk1hcml0YWwgU3RhdHVzIikgKw0KICBzY2FsZV95X2NvbnRpbnVvdXMobGFiZWxzID0gc2NhbGVzOjpwZXJjZW50KSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoIlkiID0gIiNGRkRBQjkiLCAiTiIgPSAiI0Q4QkZEOCIpKQ0KICB0aGVtZV9taW5pbWFsKCkNCmBgYA0KDQoqKk5o4bqtbiB4w6l0KioNCg0KLSBOZ8aw4budaSDEkcOjIGvhur90IGjDtG4gKE0pIGPDsyB04bu3IGzhu4cgc+G7nyBo4buvdSBuaMOgIGNhbyBoxqFuIMSRw6FuZyBr4buDIHNvIHbhu5tpIG5nxrDhu51pIMSR4buZYyB0aMOibiAoUyk6DQogIC0gVHJvbmcgbmjDs20gxJHDoyBr4bq/dCBow7RuIChNKTogR+G6p24gNzUlIGzDoCBjaOG7pyBz4bufIGjhu691IG5ow6AgKFkpLg0KDQogIC0gVHJvbmcgbmjDs20gxJHhu5ljIHRow6JuIChTKTogQ2jhu4kga2hv4bqjbmcgNTAlIGzDoCBjaOG7pyBz4bufIGjhu691IG5ow6AsIHBo4bqnbiBjw7JuIGzhuqFpIGzDoCBraMO0bmcgc+G7nyBo4buvdSBuaMOgIChOKS4NCg0KLSBT4buxIGtow6FjIGJp4buHdCByw7UgcsOgbmcgdHJvbmcgY+G6pXUgdHLDumMgdMOgaSBz4bqjbiB0aGVvIHTDrG5oIHRy4bqhbmcgaMO0biBuaMOibjoNCiAgLSDEkGnhu4F1IG7DoHkgcGjhuqNuIMOhbmggdGjhu7FjIHThur8gcuG6sW5nIG5o4buvbmcgbmfGsOG7nWkga+G6v3QgaMO0biB0aMaw4budbmcgY8OzIMSR4buZbmcgbOG7sWMga2luaCB04bq/IHbDoCB0w6BpIGNow61uaCDhu5VuIMSR4buLbmggaMahbiwgaG/hurdjIGPDsyB4dSBoxrDhu5tuZyDEkeG6p3UgdMawIGNodW5nIHbDoG8gdMOgaSBz4bqjbiBuaMawIG5ow6Ag4buf4bufDQoNCioqS2nhu4NtIMSR4buLbmggQ2hpLWLDrG5oIHBoxrDGoW5nKioNCg0KYGBge3J9DQpjaGlfcmVzdWx0MiA8LSBjaGlzcS50ZXN0KHRhYjIpDQpjaGlfcmVzdWx0Mg0KDQojIEdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQNCmNoaV9yZXN1bHQyJHN0YXRpc3RpYw0KIyBC4bqtYyB04buxIGRvIA0KY2hpX3Jlc3VsdDIkcGFyYW1ldGVyDQojIEdpw6EgdHLhu4sgcA0KY2hpX3Jlc3VsdDIkcC52YWx1ZQ0KYGBgDQoNCioqR2nhuqMgdGh1eeG6v3QqKg0KDQotIEdp4bqjIHRodXnhur90IEgwOiBUw6xuaCB0cuG6oW5nIGjDtG4gbmjDom4gKE1hcml0YWxTdGF0dXMpIHbDoCBz4bufIGjhu691IG5ow6AgKEhvbWVvd3duZXIpIGzDoCDEkeG7mWMgbOG6rXAsIGtow7RuZyBjw7MgbeG7kWkgbGnDqm4gaOG7hy4NCg0KLSBHaeG6oyB0aHV54bq/dCBIMTogVMOsbmggdHLhuqFuZyBow7RuIG5ow6JuIChNYXJpdGFsU3RhdHVzKSB2w6Agc+G7nyBo4buvdSBuaMOgIChIb21lb3d3bmVyKSBjw7MgbGnDqm4gaOG7hyB24bubaSBuaGF1Lg0KDQoqKkvhur90IHF14bqjIGtp4buDbSDEkeG7i25oOioqDQoNCi0gR2nDoSB0cuG7iyB0aOG7kW5nIGvDqiBDaGktc3F1YXJlZDogWC1zcXVhcmVkID0gMTI0MS4yMg0KDQotIELhuq1jIHThu7EgZG86IGRmID0gMQ0KDQotIEdpw6EgdHLhu4sgcDogcCA8IDIuMmUtMTYNCg0KKipL4bq/dCBsdeG6rW46KioNCg0KLSBW4bubaSBt4bupYyDDvSBuZ2jEqWEgMC4wNSwgcF92YWx1ZSA8IDAuMDUg4oaSIHRhIMSR4bunIMSRaeG7gXUga2nhu4duIMSR4buDIGLDoWMgYuG7jyBIMCwgdOG7qWMgbMOgIGPDsyBt4buRaSBsacOqbiBo4buHIGdp4buvYSAyIGJp4bq/biBNYXJpdGFsU3RhdHVzIHbDoCBIb21lb3duZXIuIEPDsyB0aOG7gyBoaeG7g3UgcuG6sW5nIHTDrG5oIHRy4bqhbmcgaMO0biBuaMOibiBjw7Mg4bqjbmggaMaw4bufbmcgxJHhur9uIHZp4buHYyBt4buZdCBuZ8aw4budaSBjw7Mgc+G7nyBo4buvdSBuaMOgIGhheSBraMO0bmcgKGhv4bq3YyBuZ8aw4bujYyBs4bqhaSkuDQoNCiMjIyAqKjQuMyBBbm51YWxJbmNvbWUgdsOgIFByb2R1Y3RDYXRlZ29yeSoqDQoNCioqQuG6o25nIHThuqduIHN14bqldCBjaMOpbyoqDQoNCmBgYHtyfQ0KdGFiMyA8LSB0YWJsZShkYXQkQW5udWFsSW5jb21lLCBkYXQkUHJvZHVjdERlcGFydG1lbnQpDQp0YWIzIA0KcHJvcC50YWJsZSh0YWIzLCAxKQ0KYGBgDQoNCioqVHLhu7FjIHF1YW4gaMOzYSoqDQoNCmBgYHtyfQ0KZ2dwbG90KGRhdCwgYWVzKHggPSBBbm51YWxJbmNvbWUsIGZpbGwgPSBQcm9kdWN0RGVwYXJ0bWVudCkpICsNCiAgZ2VvbV9iYXIocG9zaXRpb24gPSAiZmlsbCIpICsNCiAgbGFicyh0aXRsZSA9ICJU4bu3IGzhu4cgUHJvZHVjdERlcGFydG1lbnQgdGhlbyBBbm51YWwgSW5jb21lIiwgeSA9ICJU4bu3IGzhu4ciLCB4ID0gIkFubnVhbCBJbmNvbWUiKSArDQogIHNjYWxlX3lfY29udGludW91cyhsYWJlbHMgPSBzY2FsZXM6OnBlcmNlbnQpICsNCiAgdGhlbWVfbWluaW1hbCgpDQpgYGANCg0KKipOaOG6rW4geMOpdCoqDQoNCioqS2nhu4NtIMSR4buLbmggY2hpIC0gYsOsbmggcGjGsMahbmcqKg0KDQpgYGB7cn0NCmNoaV9yZXN1bHQzIDwtIGNoaXNxLnRlc3QodGFiMykNCmNoaV9yZXN1bHQzDQoNCiMgR2nDoSB0cuG7iyB0aOG7kW5nIGvDqiBDaGktc3F1YXJlZA0KY2hpX3Jlc3VsdDMkc3RhdGlzdGljDQojIELhuq1jIHThu7EgZG8gDQpjaGlfcmVzdWx0MyRwYXJhbWV0ZXINCiMgR2nDoSB0cuG7iyBwDQpjaGlfcmVzdWx0MyRwLnZhbHVlDQpgYGANCg0KKipHaeG6oyB0aHV54bq/dCoqDQoNCi0gR2nhuqMgdGh1eeG6v3QgSDA6IFRodSBuaOG6rXAgaMOgbmcgdGjDoW5nIChBbm51YWxJbmNvbWUpIHbDoCBQcm9kdWN0RGVwYXJ0bWVudCBsw6AgxJHhu5ljIGzhuq1wLCBraMO0bmcgY8OzIG3hu5FpIGxpw6puIGjhu4cuDQoNCi0gR2nhuqMgdGh1eeG6v3QgSDE6IFRodSBuaOG6rXAgaMOgbmcgdGjDoW5nIChBbm51YWxJbmNvbWUpIHbDoCBQcm9kdWN0RGVwYXJ0bWVudCBjw7MgbGnDqm4gaOG7hyB24bubaSBuaGF1Lg0KDQoqKkvhur90IHF14bqjIGtp4buDbSDEkeG7i25oOioqDQoNCi0gR2nDoSB0cuG7iyB0aOG7kW5nIGvDqiBDaGktc3F1YXJlZDogWC1zcXVhcmVkID0gMTcwLjE2DQoNCi0gQuG6rWMgdOG7sSBkbzogZGYgPSAxNDcNCg0KLSBHacOhIHRy4buLIHA6IHAgPSAwLjA5Mjc2NDcNCg0KKipL4bq/dCBsdeG6rW46KioNCg0KLSBW4bubaSBt4bupYyDDvSBuZ2jEqWEgMC4wNSBwX3ZhbHVlID0gMC4wOTI3NiA+IDAuMDUg4oaSIHRhIGtow7RuZyDEkeG7pyDEkWnhu4F1IGtp4buHbiDEkeG7gyBiw6FjIGLhu48gSDAsIHThu6ljIGzDoCBraMO0bmcgY8OzIG3hu5FpIGxpw6puIGjhu4cgZ2nhu69hIDIgYmnhur9uIEFubnVhbEluY29tZSB2w6AgUHJvZHVjdERlcGFydG1lbnQuIA0KDQojIyMgKipIb21lb3duZXIgdsOgIEFubnVhbEluY29tZSoqDQoNCioqQuG6o25nIHThuqduIHN14bqldCBjaMOpbyoqDQpgYGB7cn0NCnRhYjQgPC0gdGFibGUoZGF0JEhvbWVvd25lciwgZGF0JEFubnVhbEluY29tZSkNCnRhYjQNCnByb3AudGFibGUodGFiNCwgMSkNCg0KYGBgDQoNCioqVHLhu7FjIHF1YW4gaMOzYSoqDQoNCmBgYHtyfQ0KZ2dwbG90KGRhdCwgYWVzKHggPSBBbm51YWxJbmNvbWUsIGZpbGwgPSBIb21lb3duZXIpKSArDQogIGdlb21fYmFyKHBvc2l0aW9uID0gImZpbGwiKSArDQogIGxhYnModGl0bGUgPSAiVOG7tyBs4buHIEhvbWVvd25lciB0aGVvIEFubnVhbCBJbmNvbWUiLCB5ID0gIlThu7cgbOG7hyIsIHggPSAiQW5udWFsIEluY29tZSIpICsNCiAgc2NhbGVfeV9jb250aW51b3VzKGxhYmVscyA9IHNjYWxlczo6cGVyY2VudCkgKw0KICB0aGVtZV9taW5pbWFsKCkNCmBgYA0KDQoqKk5o4bqtbiB4w6l0KioNCg0KKipLaeG7g20gxJHhu4tuaCBjaGkgLSBiw6xuaCBwaMawxqFuZyoqDQoNCmBgYHtyfQ0KY2hpX3Jlc3VsdDQgPC0gY2hpc3EudGVzdCh0YWI0KQ0KY2hpX3Jlc3VsdDQNCg0KIyBHacOhIHRy4buLIHRo4buRbmcga8OqIENoaS1zcXVhcmVkDQpjaGlfcmVzdWx0NCRzdGF0aXN0aWMNCiMgQuG6rWMgdOG7sSBkbyANCmNoaV9yZXN1bHQ0JHBhcmFtZXRlcg0KIyBHacOhIHRy4buLIHANCmNoaV9yZXN1bHQ0JHAudmFsdWUNCmBgYA0KDQoqKkdp4bqjIHRodXnhur90KioNCg0KLSBHaeG6oyB0aHV54bq/dCBIMDogVGh1IG5o4bqtcCBow6BuZyB0aMOhbmcgKEFubnVhbEluY29tZSkgdsOgIHTDrG5oIHRy4bqhbmcgc+G7nyBo4buvdSBuaMOgIChIb21lb3duZXIpIGzDoCDEkeG7mWMgbOG6rXAsIGtow7RuZyBjw7MgbeG7kWkgbGnDqm4gaOG7hy4NCg0KLSBHaeG6oyB0aHV54bq/dCBIMTogVGh1IG5o4bqtcCBow6BuZyB0aMOhbmcgKEFubnVhbEluY29tZSkgdsOgIHTDrG5oIHRy4bqhbmcgc+G7nyBo4buvdSBuaMOgIChIb21lb3duZXIpIGPDsyBsacOqbiBo4buHIHbhu5tpIG5oYXUuDQoNCioqS+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmg6KioNCg0KLSBHacOhIHRy4buLIHRo4buRbmcga8OqIENoaS1zcXVhcmVkOiBYLXNxdWFyZWQgPSA1NDYuMzcNCg0KLSBC4bqtYyB04buxIGRvOiBkZiA9IDcNCg0KLSBHacOhIHRy4buLIHA6IHAgPDIuMmUtMTYNCg0KKipL4bq/dCBsdeG6rW46KioNCg0KLSBW4bubaSBt4bupYyDDvSBuZ2jEqWEgMC4wNSBwX3ZhbHVlIDwgMC4wNSDihpIgdGEgY8OzIMSR4bunIMSRaeG7gXUga2nhu4duIMSR4buDIGLDoWMgYuG7jyBIMCwgdOG7qWMgbMOgICBjw7MgbeG7kWkgbGnDqm4gaOG7hyBnaeG7r2EgMiBiaeG6v24gdGh1IG5o4bqtcCBow6BuZyB0aMOhbmcgdsOgIHTDrG5oIHRy4bqhbmcgc+G7nyBo4buvdSBuaMOgLiANCg0KIyMjICoqNC41IENpdHkgdsOgIEdlbmRlcioqDQoNCioqQuG6o25nIHThuqduIHN14bqldCBjaMOpbyoqDQoNCmBgYHtyfQ0KdGFiNSA8LSB0YWJsZShkYXQkR2VuZGVyLCBkYXQkQ291bnRyeSkNCnRhYjUNCnByb3AudGFibGUodGFiNSwgMSkNCmBgYA0KDQoqKlRy4buxYyBxdWFuIGjDs2EqKg0KDQpgYGB7cn0NCmdncGxvdChkYXQsIGFlcyh4ID0gQ291bnRyeSwgZmlsbCA9IEdlbmRlcikpICsNCiAgZ2VvbV9iYXIocG9zaXRpb24gPSAiZmlsbCIpICsNCiAgbGFicyh0aXRsZSA9ICJU4bu3IGzhu4cgZ2nhu5tpIHTDrW5oIHRoZW8gcXXhu5FjIGdpYSIsIHkgPSAiVOG7tyBs4buHIiwgeCA9ICJDb3VudHJ5IikgKw0KICBzY2FsZV95X2NvbnRpbnVvdXMobGFiZWxzID0gc2NhbGVzOjpwZXJjZW50KSArDQogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoIk0iID0gInllbGxvdyIsICJGIiA9ICJwaW5rIikpICsNCiAgdGhlbWVfbWluaW1hbCgpDQpgYGANCg0KKipOaOG6rW4geMOpdCoqDQoNCg0KDQoqKktp4buDbSDEkeG7i25oIGNoaSAtIGLDrG5oIHBoxrDGoW5nKioNCg0KYGBge3J9DQpjaGlfcmVzdWx0NSA8LSBjaGlzcS50ZXN0KHRhYjUpDQpjaGlfcmVzdWx0NQ0KDQojIEdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQNCmNoaV9yZXN1bHQ1JHN0YXRpc3RpYw0KIyBC4bqtYyB04buxIGRvIA0KY2hpX3Jlc3VsdDUkcGFyYW1ldGVyDQojIEdpw6EgdHLhu4sgcA0KY2hpX3Jlc3VsdDUkcC52YWx1ZQ0KYGBgDQoNCioqR2nhuqMgdGh1eeG6v3QqKg0KDQotIEdp4bqjIHRodXnhur90IEgwOiBHaeG7m2kgdMOtbmggKEdlbmRlcikgdsOgIHF14buRYyBnaWEgKENvdW50cnkpIGzDoCAyIGJp4bq/biDEkeG7mWMgbOG6rXAsIGtow7RuZyBjw7MgbeG7kWkgbGnDqm4gaOG7hy4NCg0KLSBHaeG6oyB0aHV54bq/dCBIMTogR2nhu5tpIHTDrW5oIChHZW5kZXIpIHbDoCBxdeG7kWMgZ2lhIChDb3VudHJ5KSBsw6AgMiBiaeG6v24gY8OzIG3hu5FpIGxpw6puIGjhu4cgduG7m2kgbmhhdS4NCg0KKipL4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaDoqKg0KDQotIEdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQ6IFgtc3F1YXJlZCA9IDQ4LjcyNQ0KDQotIELhuq1jIHThu7EgZG86IGRmID0gMg0KDQotIEdpw6EgdHLhu4sgcDogcCA9IDIuNjI4ZS0xMQ0KDQoqKkvhur90IGx14bqtbjoqKg0KDQotIFbhu5tpIG3hu6ljIMO9IG5naMSpYSAwLjA1IHBfdmFsdWUgPSAyLjYyOGUtMTEgPCAwLjA1IOKGkiB0YSBjw7MgxJHhu6cgxJFp4buBdSBraeG7h24gxJHhu4MgYsOhYyBi4buPIEgwLCB04bupYyBsw6AgIGPDsyBt4buRaSBsacOqbiBo4buHIGdp4buvYSAyIGJp4bq/biBsw6AgZ2nhu5tpIHTDrW5oIHbDoCBxdeG7kWMgZ2lhLiANCg0KDQojIyAqKlBI4bqmTiA1OiBSRUxBVElWRSBSSVNLICoqDQoNCiMjIyAqKjUuMS4gR2nhu5tpIHRoaeG7h3UqKg0KDQpUcm9uZyBwaMOibiB0w61jaCB0aOG7kW5nIGvDqiDEkeG7i25oIGzGsOG7o25nLCBi4bqjbmcgbmfhuqt1IG5oacOqbiAoY29udGluZ2VuY3kgdGFibGUpIGzDoCBt4buZdCBjw7RuZyBj4bulIGPGoSBi4bqjbiDEkeG7gyBraOG6o28gc8OhdCBt4buRaSBsacOqbiBo4buHIGdp4buvYSBoYWkgYmnhur9uIHBow6JuIGxv4bqhaS4gVHJvbmcgdHLGsOG7nW5nIGjhu6NwIMSRxqFuIGdp4bqjbiBuaOG6pXQsIGLhuqNuZyBuZ+G6q3Ugbmhpw6puIDJ4MiBiaeG7g3UgZGnhu4VuIHPhu7EgcGjDom4gYuG7kSBj4bunYSBoYWkgYmnhur9uIG5o4buLIHBow6JuLCBnacO6cCB4w6FjIMSR4buLbmggdsOgIMSRbyBsxrDhu51uZyBt4buRaSBsacOqbiBo4buHIGdp4buvYSBjaMO6bmcuIEPDoWMgY2jhu4kgc+G7kSBuaMawIGhp4buHdSB04bu3IGzhu4cgKHJpc2sgZGlmZmVyZW5jZSksIHThu7cgc+G7kSBuZ3V5IGPGoSAoUmVsYXRpdmUgUmlzayAtIFJSKSB2w6AgdOG7tyBz4buRIGNow6puaCAoT2RkcyBSYXRpbyAtIE9SKSDEkcOzbmcgdmFpIHRyw7IgcXVhbiB0cuG7jW5nIHRyb25nIHZp4buHYyBsxrDhu6NuZyBow7NhIG3hu5FpIHF1YW4gaOG7hyDEkcOzLiBCw6BpIHZp4bq/dCBuw6B5IHPhur0gdHLDrG5oIGLDoHkgY2hpIHRp4bq/dCBj4bqldSB0csO6YyB4w6FjIHN14bqldCBzaW5oIHJhIGLhuqNuZyBuZ+G6q3Ugbmhpw6puLCBwaMawxqFuZyBwaMOhcCBzbyBzw6FuaCBoYWkgdOG7tyBs4buHLCBjw6FjaCB4w6J5IGThu7FuZyBraG/huqNuZyB0aW4gY+G6rXkgY2hvIE9kZHMgUmF0aW8sIHbDoCBr4bq/dCB0aMO6YyBi4bqxbmcgbeG7mXQgdsOtIGThu6UgdGjhu7FjIHRp4buFbiB0cm9uZyBsxKluaCB24buxYyBraW5oIGRvYW5oLg0KDQojIyMgKio1LjIuIEPhuqV1IHRyw7pjIHbDoCB4w6FjIHN14bqldCBj4bunYSBi4bqjbmcgbmfhuqt1IG5oacOqbioqDQoNCk3hu5l0IGLhuqNuZyBuZ+G6q3Ugbmhpw6puIDJ4MiBsw6AgYuG6o25nIHThuqduIHN14bqldCDEkeG6v20gc+G7kSBxdWFuIHPDoXQgdGh14buZYyB2w6BvIHThu6tuZyB04buVIGjhu6NwIGPhu6dhIGhhaSBiaeG6v24gbmjhu4sgcGjDom46DQoNClxbDQpcYmVnaW57YXJyYXl9e3xjfGN8Y3xjfH0NClxobGluZQ0KICYgXHRleHR7S+G6v3QgcXXhuqMgKCspfSAmIFx0ZXh0e0vhur90IHF14bqjICjigJMpfSAmIFx0ZXh0e1Thu5VuZ30gXFwNClxobGluZQ0KXHRleHR7UGjGoWkgbmhp4buFbSAoWWVzKX0gJiBhICYgYiAmIGEgKyBiIFxcDQpcaGxpbmUNClx0ZXh0e0tow7RuZyBwaMahaSBuaGnhu4VtIChObyl9ICYgYyAmIGQgJiBjICsgZCBcXA0KXGhsaW5lDQpcdGV4dHtU4buVbmcgY+G7mW5nfSAmIGEgKyBjICYgYiArIGQgJiBuID0gYSArIGIgKyBjICsgZCBcXA0KXGhsaW5lDQpcZW5ke2FycmF5fQ0KXF0NCg0KxJDhu4MgaGnhu4N1IMSRxrDhu6NjIHPhu7EgaMOsbmggdGjDoG5oIGPhu6dhIGLhuqNuZyBuw6B5LCBj4bqnbiB4w6FjIMSR4buLbmggbcO0IGjDrG5oIHjDoWMgc3XhuqV0IHNpbmggcmEgZOG7ryBsaeG7h3Ug4oCUIHRyb25nIMSRw7MgcGjhu5UgYmnhur9uIG5o4bqldCBsw6AgcGjDom4gcGjhu5FpICoqUG9pc3NvbioqIHbDoCAqKk11bHRpbm9taWFsKiouDQoNCiMjIyMgKio1LjIuMS4gUGjDom4gcGjhu5FpIFBvaXNzb24qKg0KDQpQaMOibiBwaOG7kWkgUG9pc3NvbiB0aMaw4budbmcgxJHGsOG7o2Mgc+G7rSBk4bulbmcgxJHhu4MgbcO0IGjDrG5oIGjDs2Egc+G7kSBsxrDhu6NuZyBz4buxIGtp4buHbiB44bqjeSByYSB0cm9uZyBt4buZdCBraG/huqNuZyB0aOG7nWkgZ2lhbiwga2jDtG5nIGdpYW4gaG/hurdjIMSRxqFuIHbhu4sgY+G7pSB0aOG7gy4gR2nhuqMgc+G7rSBjw6FjIHPhu7Ega2nhu4duIHjhuqN5IHJhIMSR4buZYyBs4bqtcCB2w6AgduG7m2kgbeG7mXQgdOG7tyBs4buHIHRydW5nIGLDrG5oIGtow7RuZyDEkeG7lWksIG3hu5dpIMO0IHRyb25nIGLhuqNuZyBjw7MgdGjhu4MgeGVtIGzDoCBiaeG6v24gbmfhuqt1IG5oacOqbiBQb2lzc29uOg0KDQpcWw0KWF97aWp9IFxzaW0gXHRleHR7UG9pc3Nvbn0oXGxhbWJkYV97aWp9KQ0KXF0NCg0KKirGr3UgxJFp4buDbToqKg0KDQotIFBow7kgaOG7o3Aga2hpIHThu5VuZyBz4buRIHF1YW4gc8OhdCBraMO0bmcgY+G7kSDEkeG7i25oLg0KDQotIMOBcCBk4bulbmcgdHJvbmcgcGjDom4gdMOtY2ggc+G7kSBz4buxIGtp4buHbiBuaMawOiBz4buRIGzGsOG7o3QgdHJ1eSBj4bqtcCwgc+G7kSDEkcahbiBow6BuZyBs4buXaSwgdi52Lg0KDQoqKkjhuqFuIGNo4bq/OioqDQoNCi0gS2jDtG5nIGtp4buDbSBzb8OhdCDEkcaw4bujYyB04buVbmcgc+G7kSBt4bqrdSAobikuDQoNCiMjIyMgKio1LjIuMi4gUGjDom4gcGjhu5FpIE11bHRpbm9taWFsKioNCg0KVHJvbmcgdHLGsOG7nW5nIGjhu6NwIHThu5VuZyBz4buRIHF1YW4gc8OhdCBcKCBuIFwpIGzDoCBj4buRIMSR4buLbmgsIHbDoCBt4buXaSBxdWFuIHPDoXQgcsahaSB2w6BvIG3hu5l0IHRyb25nIGLhu5FuIMO0IHbhu5tpIHjDoWMgc3XhuqV0IFwoIHBfMSwgcF8yLCBwXzMsIHBfNCBcKSwgdGjDrCBi4bqjbmcgbmfhuqt1IG5oacOqbiAyeDIgY8OzIHRo4buDIMSRxrDhu6NjIG3DtCBow6xuaCBow7NhIHRoZW8gcGjDom4gcGjhu5FpIMSRYSB0aOG7qWMgTXVsdGlub21pYWwgbmjGsCBzYXU6DQoNClxbDQooYSwgYiwgYywgZCkgXHNpbSBcdGV4dHtNdWx0aW5vbWlhbH0objsgcF8xLCBwXzIsIHBfMywgcF80KQ0KXF0NCg0KUGjDom4gcGjhu5FpIG7DoHkgdGjGsOG7nW5nIMSRxrDhu6NjIHPhu60gZOG7pW5nIHRyb25nIGPDoWMgdMOsbmggaHXhu5FuZyBtw6AgZOG7ryBsaeG7h3UgxJHGsOG7o2MgdGh1IHRo4bqtcCB04burIGto4bqjbyBzw6F0IGhv4bq3YyBuZ2hpw6puIGPhu6l1IHjDoyBo4buZaSBo4buNYyB24bubaSBrw61jaCB0aMaw4bubYyBt4bqrdSBj4buRIMSR4buLbmguIEtoaSDEkcOzLCBjw6FjIHjDoWMgc3XhuqV0IFwoIHBfMSwgcF8yLCBwXzMsIHBfNCBcKSBiaeG7g3UgZGnhu4VuIHjDoWMgc3XhuqV0IG3DoCBt4buZdCBjw6EgdGjhu4MgcsahaSB2w6BvIHThu6tuZyDDtCB0cm9uZyBi4bqjbmcgMngyLg0KDQoqKsavdSDEkWnhu4NtIGPhu6dhIG3DtCBow6xuaCBNdWx0aW5vbWlhbDoqKg0KDQotIEtp4buDbSBzb8OhdCDEkcaw4bujYyB04buVbmcgc+G7kSBt4bqrdS4NCg0KLSBQaMO5IGjhu6NwIGNobyBk4buvIGxp4buHdSBraOG6o28gc8OhdCwgdGjhu60gbmdoaeG7h20uDQoNCioqU28gc8OhbmggZ2nhu69hIHBow6JuIHBo4buRaSBQb2lzc29uIHbDoCBNdWx0aW5vbWlhbCoqDQoNCnwgVGnDqnUgY2jDrSAgICAgICAgIHwgUG9pc3NvbiAgICAgICAgICAgICAgICAgIHwgTXVsdGlub21pYWwgICAgICAgICAgICAgICAgICAgICAgICAgIHwNCnwtLS0tLS0tLS0tLS0tLS0tLS18LS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tfC0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLXwNCnwgVOG7lW5nIHPhu5EgbeG6q3UgICAgICB8IEtow7RuZyBj4buRIMSR4buLbmggICAgICAgICAgICAgfCBD4buRIMSR4buLbmggICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgfA0KfCDhu6huZyBk4bulbmcgY2jDrW5oICAgfCBT4buRIHPhu7Ega2nhu4duICAgICAgICAgICAgICAgIHwgVOG6p24gc3XhuqV0IGto4bqjbyBzw6F0ICAgICAgICAgICAgICAgICAgICAgfA0KfCBE4buvIGxp4buHdSBwaMO5IGjhu6NwICB8IEdpYW8gZOG7i2NoLCBs4buXaSwgdGFpIG7huqFuICAgfCBUcuG6oyBs4budaSBraOG6o28gc8OhdCwgcGjDom4gbmjDs20gICAgICAgICAgfA0KDQojIyMgKio1LjMuIEPDoWMgY2jhu4kgc+G7kSDEkW8gbeG7kWkgbGnDqm4gaOG7hyoqDQoNCiMjIyMgKio1LjMuMS4gIEhp4buHdSB04bu3IGzhu4cgKFJpc2sgRGlmZmVyZW5jZSAtIFJEKSoqDQoNClxbDQpSRCA9IFxmcmFje2F9e2EgKyBifSAtIFxmcmFje2N9e2MgKyBkfQ0KXF0NCg0KKirDnSBuZ2jEqWEqKjoNCg0KLSDEkG8gbMaw4budbmcgKipjaMOqbmggbOG7h2NoIHR1eeG7h3QgxJHhu5FpKiogZ2nhu69hIHjDoWMgc3XhuqV0IHjhuqN5IHJhIGvhur90IHF14bqjIOG7nyBuaMOzbSBwaMahaSBuaGnhu4VtIHbDoCBuaMOzbSBraMO0bmcgcGjGoWkgbmhp4buFbS4NCg0KLSBUaMaw4budbmcgxJHGsOG7o2Mgc+G7rSBk4bulbmcgdHJvbmcgxJHDoW5oIGdpw6EgKip0w6FjIMSR4buZbmcgY+G7p2EgY2jDrW5oIHPDoWNoKiosIGNhbiB0aGnhu4dwIGhv4bq3YyBjaMawxqFuZyB0csOsbmggdGjDrSDEkWnhu4NtLCBraGkgc+G7sSBraMOhYyBiaeG7h3QgduG7gSB4w6FjIHN14bqldCBsw6AgcXVhbiB0cuG7jW5nIGjGoW4gc28gduG7m2kgdOG7tyBz4buRLg0KDQpLaG/huqNuZyB0aW4gY+G6rXkgOTUlIChDSSkgY2hvIFJEIGPDsyB0aOG7gyDEkcaw4bujYyDGsOG7m2MgbMaw4bujbmcgYuG6sW5nOg0KDQpcWw0KQ0lfe1JEfSA9IFJEIFxwbSBaX3sxIC0gXGFscGhhLzJ9IFxjZG90IFNFX3tSRH0NClxdDQoNClRyb25nIMSRw7M6DQotIFwoIFpfezEgLSBcYWxwaGEvMn0gXCkgbMOgIGdpw6EgdHLhu4sgdOG7m2kgaOG6oW4gKHRoxrDhu51uZyDiiYggMS45NiB24bubaSA5NSUgQ0kpLA0KDQotIFwoIFNFX3tSRH0gPSBcc3FydHtcZnJhY3thKGErYi1hKX17KGErYileM30gKyBcZnJhY3tjKGMrZC1jKX17KGMrZCleM319IFwpDQoNCiMjIyMgKio1LjMuMi4gIFThu7cgc+G7kSBuZ3V5IGPGoSAoUmVsYXRpdmUgUmlzayAtIFJSKSoqDQoNClxbDQpSUiA9IFxmcmFje2EgLyAoYSArIGIpfXtjIC8gKGMgKyBkKX0NClxdDQoNCioqw50gbmdoxKlhKio6DQoNCi0gXCggUlIgPiAxIFwpOiBuZ3V5IGPGoSB44bqjeSByYSBr4bq/dCBxdeG6oyBjYW8gaMahbiDhu58gbmjDs20gcGjGoWkgbmhp4buFbS4NCg0KLSBcKCBSUiA8IDEgXCk6IG5ndXkgY8ahIHjhuqN5IHJhIGvhur90IHF14bqjIHRo4bqlcCBoxqFuIOG7nyBuaMOzbSBwaMahaSBuaGnhu4VtLg0KDQotIFwoIFJSID0gMSBcKToga2jDtG5nIGPDsyBt4buRaSBsacOqbiBo4buHIGdp4buvYSBwaMahaSBuaGnhu4VtIHbDoCBr4bq/dCBxdeG6oy4NCg0KVGEgdGjGsOG7nW5nIGzhuqV5ICoqbG9nKiogY+G7p2EgUlIgxJHhu4MgdMOtbmgga2hv4bqjbmcgdGluIGPhuq15Og0KDQpcWw0KQ0lfe1JSfSA9IFxleHAgXGxlZnRbIFxsbihSUikgXHBtIFpfezEgLSBcYWxwaGEvMn0gXGNkb3QgU0Vfe1xsbihSUil9IFxyaWdodF0NClxdDQoNClRyb25nIMSRw7M6DQoNClxbDQpTRV97XGxuKFJSKX0gPSBcc3FydHsgXGZyYWN7MX17YX0gLSBcZnJhY3sxfXthICsgYn0gKyBcZnJhY3sxfXtjfSAtIFxmcmFjezF9e2MgKyBkfSB9DQpcXQ0KDQoqKlbDrSBk4bulIHRyb25nIGtpbmggZG9hbmgqKjoNCg0KLSBU4bu3IGzhu4cga2jDoWNoIG11YSBow6BuZyBraGkgY8OzIGtodXnhur9uIG3Do2k6IFwoIDgwLzEwMCA9IDAuOCBcKQ0KDQotIEtoaSBraMO0bmcgY8OzIGtodXnhur9uIG3Do2k6IFwoIDQwLzEwMCA9IDAuNCBcKQ0KDQotIFwoIFJSID0gMC44IC8gMC40ID0gMiBcKQ0KDQrih5IgS2h1eeG6v24gbcOjaSBsw6BtIHTEg25nIGfhuqVwICoqxJHDtGkqKiBraOG6oyBuxINuZyBraMOhY2ggbXVhIGjDoG5nLg0KDQpI4bqhbiBjaOG6vzoNCg0K4pqg77iPIEtow7RuZyBz4butIGThu6VuZyDEkcaw4bujYyB0cm9uZyBuZ2hpw6puIGPhu6l1IGLhu4duaCBjaOG7qW5nIChjYXNlLWNvbnRyb2wpIHbDrCB0YSBraMO0bmcgYmnhur90IG5ndXkgY8ahIHR1eeG7h3QgxJHhu5FpLg0KDQrimqDvuI8gROG7hSBnw6J5IGhp4buDdSBuaOG6p20gbuG6v3Uga2jDtG5nIMSRaSBrw6htIG5ndXkgY8ahIHR1eeG7h3QgxJHhu5FpLiBWw60gZOG7pTogUlIgPSAyIGPDsyB24bq7IGNhbywgbmjGsG5nIG7hur91IG5ndXkgY8ahIGJhbiDEkeG6p3UgbMOgIDElLCB0aMOsIHTEg25nIGzDqm4gMiUgduG6q24gbMOgIHLhuqV0IHRo4bqlcC4NCg0K4pqg77iPIEtow7RuZyDEkeG7kWkgeOG7qW5nOiBO4bq/dSDEkeG7lWkgbmjDs20gdGhhbSBjaGnhur91IHRow6wgZ2nDoSB0cuG7iyBSUiB0aGF5IMSR4buVaSwga2jDoWMgduG7m2kgT1IuDQoNCiMjIyMgKio1LjMuMy4gIFThu7cgc+G7kSBjaMOqbmggKE9kZHMgUmF0aW8gLSBPUikqKg0KDQojIyMjIyAqKjUuMy4zLjEuIE9kZHMqKg0KDQotICoqT2RkcyAodOG7tyBz4buRIGPGoSBo4buZaSkqKjogIA0KICBcWw0KICBcdGV4dHtPZGRzfSA9IFxmcmFje3B9ezEgLSBwfQ0KICBcXQ0KICBUcm9uZyDEkcOzLCBcKHBcKSBsw6AgeMOhYyBzdeG6pXQgeOG6o3kgcmEgY+G7p2EgbeG7mXQgYmnhur9uIG5o4buLIHBow6JuIChuaMawIHZp4buHYyBz4bufIGjhu691IG5ow6ApLg0KICANCiMjIyMjICoqNS4zLjMuMS4gT2RkcyBSYXRpbyAtIE9SKioNCg0KXFsNCk9SID0gXGZyYWN7YSAvIGJ9e2MgLyBkfSA9IFxmcmFje2FkfXtiY30NClxdDQoNCkhheSAgDQoNClxbDQogIFx0ZXh0e09SfSA9IFxmcmFje1x0ZXh0e09kZHMg4bufIG5ow7NtIDF9fXtcdGV4dHtPZGRzIOG7nyBuaMOzbSAyfX0NClxdDQogIE9SIHRoxrDhu51uZyBkw7luZyDEkeG7gyDEkW8gbMaw4budbmcgbeG7kWkgbGnDqm4gaOG7hyBnaeG7r2EgbeG7mXQgeeG6v3UgdOG7kSB2w6AgbeG7mXQga+G6v3QgcXXhuqMgdHJvbmcgYuG6o25nIDJ4Mi4NCiAgDQotICoqT1IgPSAxKio6IEtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgduG7gSBvZGRzIGdp4buvYSBoYWkgbmjDs20uICANCg0KLSAqKk9SID4gMSoqOiBOaMOzbSAxIGPDsyBvZGRzIHjhuqN5IHJhIGvhur90IHF14bqjIGNhbyBoxqFuIG5ow7NtIDIuICANCg0KLSAqKk9SIDwgMSoqOiBOaMOzbSAxIGPDsyBvZGRzIHjhuqN5IHJhIGvhur90IHF14bqjIHRo4bqlcCBoxqFuIG5ow7NtIDIuDQoNCj4gVsOtIGThu6U6IE9SID0gMiBuZ2jEqWEgbMOgICoqb2RkcyB44bqjeSByYSBr4bq/dCBxdeG6oyDhu58gbmjDs20gMSBjYW8gZ+G6pXAgMiBs4bqnbiBuaMOzbSAyKiouDQoNCi0gKipMxrB1IMO9Kio6IE9SIGtow6FjIHbhu5tpIFJSLCDEkeG6t2MgYmnhu4d0IGtoaSBiaeG6v24gcGjhu6UgdGh14buZYyBraMO0bmcgaGnhur9tIGfhurdwLiBLaGkga+G6v3QgcXXhuqMgbMOgIGhp4bq/bSAocCBuaOG7jyksIHRow6wgT1Ig4omIIFJSLg0KDQoqKsOdIG5naMSpYSoqOg0KDQotIE9SIGNobyBiaeG6v3QgdOG7tyBs4buHICoqb2RkcyoqICh04bupYyBsw6AgInjDoWMgc3XhuqV0IGNoaWEgY2hvIDEgdHLhu6sgeMOhYyBzdeG6pXQiKSBnaeG7r2EgaGFpIG5ow7NtLg0KDQotIFRoxrDhu51uZyDEkcaw4bujYyBz4butIGThu6VuZyBwaOG7lSBiaeG6v24gdHJvbmcgbcO0IGjDrG5oICoqbG9naXN0aWMgcmVncmVzc2lvbioqLg0KDQotIFRyb25nIGPDoWMgbmdoacOqbiBj4bupdSBi4buHbmggaGnhur9tIChyYXJlIGRpc2Vhc2UgYXNzdW1wdGlvbiksIE9SIHjhuqVwIHjhu4kgduG7m2kgUlIuDQoNCioqxq91IMSRaeG7g20gY+G7p2EgT2RkcyBSYXRpbyoqOg0KDQotIEThu4UgdMOtbmggdsOgIGThu4UgZGnhu4VuIGdp4bqjaS4NCg0KLSDhu5RuIMSR4buLbmggduG7gSBt4bq3dCB0b8OhbiBo4buNYyBraGkgc+G7rSBk4bulbmcgdHJvbmcgY8OhYyBtw7QgaMOsbmggaOG7k2kgcXV5Lg0KDQotIEtow7RuZyBwaOG7pSB0aHXhu5ljIHbDoG8gdOG7tyBs4buHIGhp4buHbiBkaeG7h24gdHJvbmcgbeG6q3UgKMSR4bq3YyBiaeG7h3QgdHJvbmcgY8OhYyB0aGnhur90IGvhur8gKipi4buHbmgg4oCTIGNo4bupbmcqKjogKmNhc2UtY29udHJvbCBzdHVkeSopLg0KDQojIyMgKio1LjQuIEtob+G6o25nIMaw4bubYyBsxrDhu6NuZyB0aW4gY+G6rXkgY2hvIE9ERFMgUkFUSU8qKg0KDQrEkOG7gyB4w6FjIMSR4buLbmggeGVtIE9SIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogaGF5IGtow7RuZywgY+G6p24geMOieSBk4buxbmcga2hv4bqjbmcgdGluIGPhuq15Lg0KDQoqKkxvZ2FyaXRobSB04buxIG5oacOqbiBj4bunYSBPZGRzIFJhdGlvOioqDQoNClxbDQpcbG9nKE9SKSA9IFxsb2dcbGVmdCggXGZyYWN7YSBcY2RvdCBkfXtiIFxjZG90IGN9IFxyaWdodCkNClxdDQoNCioqU2FpIHPhu5EgY2h14bqpbiAoU3RhbmRhcmQgRXJyb3IpOioqDQoNClxbDQpTRSA9IFxzcXJ0eyBcZnJhY3sxfXthfSArIFxmcmFjezF9e2J9ICsgXGZyYWN7MX17Y30gKyBcZnJhY3sxfXtkfSB9DQpcXQ0KDQoqKktob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIFwoXGxvZyhPUilcKToqKg0KDQpcWw0KXGxvZyhPUikgXHBtIDEuOTYgXGNkb3QgU0UNClxdDQoNCioqTOG6pXkgbcWpIMSR4buDIGPDsyBraG/huqNuZyB0aW4gY+G6rXkgOTUlIGNobyBPUjoqKg0KDQpcWw0KQ0lfezk1XCV9ID0gXGxlZnRbIGVee1xsb2coT1IpIC0gMS45NiBcY2RvdCBTRX0sXHF1YWQgZV57XGxvZyhPUikgKyAxLjk2IFxjZG90IFNFfSBccmlnaHRdDQpcXQ0KDQotIE7hur91IGtob+G6o25nIHRpbiBj4bqteSAqKmtow7RuZyBjaOG7qWEgMSoqLCB0YSBjw7MgdGjhu4Mga+G6v3QgbHXhuq1uIHLhurFuZyAqKm3hu5FpIGxpw6puIGjhu4cgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqioqIOG7nyBt4bupYyDDvSBuZ2jEqWEgNSUuDQoNCiMjIyAqKjUuNS4g4buobmcgZOG7pW5nIHRyb25nIGtpbmggZG9hbmg6IEhJ4buGVSBRVeG6oiBD4bumQSBLSFVZ4bq+TiBNw4NJKioNCg0KTeG7mXQgY8O0bmcgdHkgdGjGsMahbmcgbeG6oWkgxJFp4buHbiB04butIHRo4butIG5naGnhu4dtIGNoaeG6v24gZOG7i2NoIGtodXnhur9uIG3Do2kgxJHhu4Mga2nhu4NtIHRyYSB4ZW0gbGnhu4d1IGPDsyDhuqNuaCBoxrDhu59uZyDEkeG6v24gaMOgbmggdmkgbXVhIGjDoG5nIGtow7RuZy4gSOG7jSBs4bqleSBt4bqrdSAyMDAga2jDoWNoIGjDoG5nIG5n4bqrdSBuaGnDqm4gdsOgIHRodSDEkcaw4bujYyBi4bqjbmcgc2F1Og0KDQp8ICAgICAgICAgICAgICAgICAgfCBNdWEgaMOgbmcgfCBLaMO0bmcgbXVhIHwgVOG7lW5nIHwNCnwtLS0tLS0tLS0tLS0tLS0tLS18LS0tLS0tLS0tLXwtLS0tLS0tLS0tLS18LS0tLS0tLXwNCnwgQ8OzIGtodXnhur9uIG3Do2kgICAgfCA4MCAgICAgICB8IDIwICAgICAgICAgfCAxMDAgICB8DQp8IEtow7RuZyBraHV54bq/biBtw6NpIHwgNDAgICAgICAgfCA2MCAgICAgICAgIHwgMTAwICAgfA0KDQotIFTDrW5oIE9kZHMgUmF0aW8gKE9SKToNCg0KXFsNCk9SID0gXGZyYWN7ODAgXGNkb3QgNjB9ezIwIFxjZG90IDQwfSA9IFxmcmFjezQ4MDB9ezgwMH0gPSA2DQpcXQ0KDQrihpIgT2RkcyBtdWEgaMOgbmcga2hpIGPDsyBraHV54bq/biBtw6NpICoqY2FvIGfhuqVwIDYgbOG6p24qKiBzbyB24bubaSBraGkga2jDtG5nIGPDsyBraHV54bq/biBtw6NpLg0KDQotIExvZ2FyaXRobSB04buxIG5oacOqbiBj4bunYSBPUjoNCg0KXFsNClxsb2coT1IpID0gXGxvZyg2KSBcYXBwcm94IDEuNzkNClxdDQoNCi0gU2FpIHPhu5EgY2h14bqpbiAoU0UpOg0KDQpcWw0KU0UgPSBcc3FydHsgXGZyYWN7MX17ODB9ICsgXGZyYWN7MX17MjB9ICsgXGZyYWN7MX17NDB9ICsgXGZyYWN7MX17NjB9IH0gXGFwcHJveCAwLjM1DQpcXQ0KDQotIEtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIE9SOg0KDQpcWw0KQ0kgPSBcbGVmdFsgZV57MS43OSAtIDEuOTYgXGNkb3QgMC4zNX0sXHF1YWQgZV57MS43OSArIDEuOTYgXGNkb3QgMC4zNX0gXHJpZ2h0XSANClxhcHByb3ggXGxlZnRbIGVeezEuMTF9LFxxdWFkIGVeezIuNDd9IFxyaWdodF0gDQo9IFszLjAzLFwgMTEuODFdDQpcXQ0KDQotIFbDrCBraG/huqNuZyB0aW4gY+G6rXkgKipraMO0bmcgY2jhu6lhIDEqKiwgY8OzIHRo4buDIGvhur90IGx14bqtbiBy4bqxbmcgKipraHV54bq/biBtw6NpIGPDsyDhuqNuaCBoxrDhu59uZyBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqKiogxJHhur9uIGjDoG5oIHZpIG11YSBow6BuZy4NCg0KLSBW4bubaSBcKCBPUiA9IDYgXCkgdsOgIFwoIENJID0gKDMuMDMsXCAxMS44MSkgXCksIGRvYW5oIG5naGnhu4dwIGPDsyB0aOG7gyB04buxIHRpbiB0cmnhu4NuIGtoYWkgKipjaGnhur9uIGzGsOG7o2Mga2h1eeG6v24gbcOjaSBxdXkgbcO0IGzhu5tuIGjGoW4qKi4NCg0KIyMjICoqNS42IEdFTkRFUiBWw4AgTUFSSVRBTFNUQVRVUyoqDQoNCiMjIyMgKipU4bqhbyBi4bqjbmcgY2jDqW8qKg0KDQpgYGB7cn0NCiMgVOG6oW8gYuG6o25nIGNow6lvIGdp4buvYSBHZW5kZXIgdsOgIE1hcml0YWxTdGF0dXMNCnRhYmxlX2dlbmRlcl9tYXJpdGFsIDwtIHRhYmxlKGRhdCRHZW5kZXIsIGRhdCRNYXJpdGFsU3RhdHVzKQ0KYWRkbWFyZ2lucyh0YWJsZV9nZW5kZXJfbWFyaXRhbCkgICMgVGjDqm0gaMOgbmcgdsOgIGPhu5l0IHThu5VuZw0KYGBgDQoNCiMjIyMgKipUw61uaCB0b8OhbiBSZWxhdGl2ZSBSaXNrKioNCg0KYGBge3J9DQojIFRyw61jaCB4deG6pXQgc+G7kSBsaeG7h3UgdOG7qyBi4bqjbmcgY2jDqW8NCkEgPC0gdGFibGVfZ2VuZGVyX21hcml0YWxbIkYiLCJTIl0gIyBu4buvIMSR4buZYyB0aMOibg0KQiA8LSB0YWJsZV9nZW5kZXJfbWFyaXRhbFsiRiIsIk0iXSAjIG7hu68ga+G6v3QgaMO0bg0KQyA8LSB0YWJsZV9nZW5kZXJfbWFyaXRhbFsiTSIsIlMiXSAjIG5hbSDEkeG7mWMgdGjDom4NCkQgPC0gdGFibGVfZ2VuZGVyX21hcml0YWxbIk0iLCJNIl0gIyBuYW0ga+G6v3QgaMO0bg0KDQojIFTDrW5oIFJSDQpSUiA8LSAoQS8oQStCKSkvKEMvKEMrRCkpDQpSUg0KIyBUcm9uZyDEkcOzOiBBLyhBK0IpIGzDoCB04bu3IGzhu4cgcGjhu6UgbuG7ryDEkeG7mWMgdGjDom4gdHJvbmcgdOG7lW5nIHPhu5EgbuG7rw0KIyBWw6A6IEMvKEMrRCkgbMOgIHThu7cgbOG7hyBuYW0gZ2nhu5tpIMSR4buZYyB0aMOibiB0cm9uZyB04buVbmcgc+G7kSBuYW0NCmBgYA0KDQojIyMjICoqNS42LjMgw50gbmdoxKlhKioNCg0KLSBSUiA9IDE6IEtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgduG7gSB04bu3IGzhu4cgxJHhu5ljIHRow6JuIGdp4buvYSBoYWkgZ2nhu5tpIHTDrW5oLg0KLSBSUiA+IDE6IE5hbSBnaeG7m2kgY8OzIHThu7cgbOG7hyDEkeG7mWMgdGjDom4gY2FvIGjGoW4gcGjhu6UgbuG7ry4NCi0gUlIgPCAxOiBQaOG7pSBu4buvIGPDsyB04bu3IGzhu4cgxJHhu5ljIHRow6JuIGNhbyBoxqFuIG5hbSBnaeG7m2kuDQoNCi0gVOG7tyBs4buHIHBo4bulIG7hu68gxJHhu5ljIHRow6JuIGLhurFuZyBraG/huqNuZyA5NC41NyUgdOG7tyBs4buHIG5hbSBnaeG7m2kgxJHhu5ljIHRow6JuLiDEkGnhu4F1IG7DoHkgY2hvIHRo4bqleSBraMO0bmcgY8OzIHPhu7Ega2jDoWMgYmnhu4d0IGzhu5tuIHbhu4EgdOG7tyBs4buHIMSR4buZYyB0aMOibiBnaeG7r2EgaGFpIGdp4bubaSB0w61uaCB0cm9uZyBk4buvIGxp4buHdS4gDQoNCmBgYHtyfQ0KaW5zdGFsbC5wYWNrYWdlcygiZXBpdG9vbHMiLCByZXBvcyA9ICJodHRwczovL2NyYW4ucnN0dWRpby5jb20vIikNCmxpYnJhcnkoZXBpdG9vbHMpDQpyaXNrcmF0aW8odGFibGVfZ2VuZGVyX21hcml0YWwsIG1ldGhvZD0id2FsZCIpDQojIENo4buJIMSR4buLbmggcGjGsMahbmcgcGjDoXAgdMOtbmgga2hv4bqjbmcgdGluIGPhuq15IChjb25maWRlbmNlIGludGVydmFsKSBjaG8gUmlzayBSYXRpby4NCiMgIndhbGQiIGzDoCBwaMawxqFuZyBwaMOhcCBz4butIGThu6VuZyDGsOG7m2MgbMaw4bujbmcgV2FsZCwgZ2nhuqMgxJHhu4tuaCBwaMOibiBwaOG7kWkgY2h14bqpbiwgxJHhu4MgdMOtbmgga2hv4bqjbmcgdGluIGPhuq15Lg0KYGBgDQoqKk5o4bqtbiB4w6l0KioNCg0KLSBSZWxhdGl2ZSBSaXNrIChSUikNCg0KICAtIEYgKEZlbWFsZSk6IDEuMDAwMDAwDQogIC0gTSAoTWFsZSk6IDEuMDU3NDE3IChDSTogMS4wMjM4MTMgLSAxLjA5MjEyMykNCiAgDQotIE5o4bqtbiB4w6l0IHbhu4EgUlI6DQoNCiAgLSBSUiBj4bunYSBuaMOzbSBGIChQaOG7pSBu4buvKSBsw6AgMSB2w6wgbmjDs20gbsOgeSDEkcaw4bujYyBs4bqleSBsw6BtIG5ow7NtIHRoYW0gY2hp4bq/dS4NCiAgLSBSUiBj4bunYSBuaMOzbSBNIChOYW0gZ2nhu5tpKSBsw6AgMS4wNTc0LCDEkWnhu4F1IG7DoHkgY2hvIHRo4bqleSBuYW0gZ2nhu5tpIGPDsyB04bu3IGzhu4cga+G6v3QgaMO0biAoTSkgY2FvIGjGoW4gbeG7mXQgY2jDunQgc28gduG7m2kgcGjhu6UgbuG7rywgduG7m2kgdOG7tyBs4buHIDEuMDU3NCBs4bqnbiBzbyB24bubaSBuaMOzbSBwaOG7pSBu4buvLiBT4buxIGtow6FjIGJp4buHdCBuw6B5IGPDsyB0aOG7gyBsw6Aga2jDtG5nIGzhu5tuLCBuaMawbmcgbsOzIGNo4buJIHJhIHLhurFuZyBuYW0gZ2nhu5tpIGPDsyBt4buZdCBraOG6oyBuxINuZyBr4bq/dCBow7RuIGNhbyBoxqFuIHBo4bulIG7hu68gdHJvbmcgbeG6q3UgbmdoacOqbiBj4bupdSBuw6B5Lg0KICAtIENvbmZpZGVuY2UgSW50ZXJ2YWwgKENJKSBj4bunYSBSUiBjaG8gbmjDs20gTSBsw6AgWzEuMDIzOCwgMS4wOTIxXSwgbmdoxKlhIGzDoCBraG/huqNuZyB0aW4gY+G6rXkgY+G7p2EgUlIga2jDtG5nIGJhbyBn4buTbSBnacOhIHRy4buLIDEuIMSQaeG7gXUgbsOgeSBjaG8gdGjhuqV5IHPhu7Ega2jDoWMgYmnhu4d0IGdp4buvYSBoYWkgbmjDs20gRiB2w6AgTSBsw6AgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiB2w6Aga2jDtG5nIHBo4bqjaSBuZ+G6q3Ugbmhpw6puLg0KLSBL4bq/dCBsdeG6rW4gduG7gSBSUjoNCg0KICAtIFJSIGNobyB0aOG6pXkgbmFtIGdp4bubaSBjw7MgeHUgaMaw4bubbmcga+G6v3QgaMO0biBjYW8gaMahbiBzbyB24bubaSBwaOG7pSBu4buvLCB24bubaSBSUiA9IDEuMDU3NC4NCiAgLSBDb25maWRlbmNlIEludGVydmFsIGtow7RuZyBiYW8gZ+G7k20gZ2nDoSB0cuG7iyAxLCDEkWnhu4F1IG7DoHkgaOG7lyB0cuG7oyBjaG8gZ2nhuqMgdGh1eeG6v3QgcuG6sW5nIHPhu7Ega2jDoWMgYmnhu4d0IGdp4buvYSB04bu3IGzhu4cga+G6v3QgaMO0biBj4bunYSBuYW0gdsOgIG7hu68gY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqi4NCiAgDQotIHAtdmFsdWUgdOG7qyBjw6FjIGtp4buDbSDEkeG7i25oOg0KDQogIC0gQ2hpLXNxdWFyZSB0ZXN0OiBwLXZhbHVlID0gMC4wMDA3MDM4ODM3DQogIC0gRmlzaGVyJ3MgZXhhY3QgdGVzdDogcC12YWx1ZSA9IDAuMDAwNzM3Mzg5NQ0KICAtIE1pZHBvaW50IGV4YWN0IHRlc3Q6IHAtdmFsdWUgPSAwLjAwMDcwNDUwOTINCiAgDQotIE5o4bqtbiB4w6l0IHbhu4EgcC12YWx1ZToNCiAgLSBDw6FjIHAtdmFsdWUgdOG7qyBjw6FjIGtp4buDbSDEkeG7i25oIMSR4buBdSBuaOG7jyBoxqFuIG3hu6ljIMO9IG5naMSpYSB0aMO0bmcgdGjGsOG7nW5nIGzDoCAwLjA1LCBuZ2jEqWEgbMOgIGPDsyDEkeG7pyBi4bqxbmcgY2jhu6luZyDEkeG7gyBiw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3QgbnVsbCBoeXBvdGhlc2lzIChnaeG6oyB0aHV54bq/dCBraMO0bmcgY8OzIHPhu7Ega2jDoWMgYmnhu4d0IGdp4buvYSBuYW0gdsOgIG7hu68pLg0KICAtIMSQaeG7gXUgbsOgeSBjaOG7iSByYSBy4bqxbmcgY8OzIHPhu7Ega2jDoWMgYmnhu4d0IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogZ2nhu69hIHThu7cgbOG7hyBr4bq/dCBow7RuIOG7nyBwaOG7pSBu4buvIHbDoCBuYW0gZ2nhu5tpLg0KICANCi0gVMOzbSB04bqvdCB2w6Agbmjhuq1uIHjDqXQgdOG7lW5nIHF1w6F0DQogIC0gUmVsYXRpdmUgUmlzayAoUlIpOiBU4bu3IGzhu4cga+G6v3QgaMO0biBj4bunYSBuYW0gZ2nhu5tpIGzDoCAxLjA1NzQgbOG6p24gc28gduG7m2kgcGjhu6UgbuG7rywgxJFp4buBdSBuw6B5IGNobyB0aOG6pXkgbmFtIGdp4bubaSBjw7Mga2jhuqMgbsSDbmcga+G6v3QgaMO0biBjYW8gaMahbiBwaOG7pSBu4buvIHRyb25nIG3huqt1IG5naGnDqm4gY+G7qXUgbsOgeS4NCiAgLSBDb25maWRlbmNlIEludGVydmFsIChDLkkuKTogS2hv4bqjbmcgdGluIGPhuq15IFsxLjAyMzgsIDEuMDkyMV0ga2jDtG5nIGJhbyBn4buTbSBnacOhIHRy4buLIDEsIMSRaeG7gXUgbsOgeSBjw7MgbmdoxKlhIGzDoCBz4buxIGtow6FjIGJp4buHdCBnaeG7r2EgdOG7tyBs4buHIGvhur90IGjDtG4gY+G7p2EgbmFtIHbDoCBu4buvIGzDoCBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqLg0KICAtIHAtdmFsdWU6IEPDoWMgcC12YWx1ZSB04burIENoaS1zcXVhcmUsIEZpc2hlcidzIGV4YWN0IHRlc3QgdsOgIG1pZHAuZXhhY3QgxJHhu4F1IG5o4buPIGjGoW4gMC4wNSwgY2hvIHRo4bqleSBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgxJHDoW5nIGvhu4MgZ2nhu69hIHThu7cgbOG7hyBr4bq/dCBow7RuIGPhu6dhIG5hbSB2w6AgbuG7ry4NCiAgDQotIEvhur90IGx14bqtbg0KICAtIFThu7cgbOG7hyBr4bq/dCBow7RuIGdp4buvYSBuYW0gdsOgIG7hu68gY8OzIHPhu7Ega2jDoWMgYmnhu4d0IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ouIEPhu6UgdGjhu4MsIG5hbSBnaeG7m2kgY8OzIHThu7cgbOG7hyBr4bq/dCBow7RuIGNhbyBoxqFuIHBo4bulIG7hu68uDQogIC0gUmVsYXRpdmUgUmlzayBjaG8gdGjhuqV5IG5hbSBnaeG7m2kgY8OzIGto4bqjIG7Eg25nIGvhur90IGjDtG4gY2FvIGjGoW4gcGjhu6UgbuG7ryAoduG7m2kgUlIgPSAxLjA1NzQpLCB2w6Aga+G6v3QgcXXhuqMgbsOgeSBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqIHLDtSByw6BuZywgdsOsIHAtdmFsdWUgxJHhu4F1IHLhuqV0IG5o4buPLg0KICAtIEtob+G6o25nIHRpbiBj4bqteSA5NSUgY+G7p2EgUlIga2jDtG5nIGJhbyBn4buTbSBnacOhIHRy4buLIDEsIMSRaeG7gXUgbsOgeSBjw6BuZyBj4bunbmcgY+G7kSBr4bq/dCBsdeG6rW4gduG7gSBz4buxIGtow6FjIGJp4buHdCBnaeG7r2EgdOG7tyBs4buHIGvhur90IGjDtG4gY+G7p2EgaGFpIGdp4bubaS4NCiAgDQojIyMgKio1LjcgR0VOREVSIFbDgCBIT01FT1dORVIqKg0KDQpDw6J1IGjhu49pIG5naGnDqm4gY+G7qXU6IExp4buHdSBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgduG7gSB04bu3IGzhu4cgc+G7nyBo4buvdSBuaMOgIGdp4buvYSBuYW0gdsOgIG7hu68ga2jDtG5nPw0KDQojIyMjICoqNS43LjEgVOG6oW8gYuG6o25nIGNow6lvKioNCg0KYGBge3J9DQojIFThuqFvIGLhuqNuZyBjaMOpbyBnaeG7r2EgR2VuZGVyIHbDoCBIb21lb3duZXINCnRhYmxlX2dlbmRlcl9ob21lb3duZXIgPC0gdGFibGUoZGF0JEdlbmRlciwgZGF0JEhvbWVvd25lcikNCmFkZG1hcmdpbnModGFibGVfZ2VuZGVyX2hvbWVvd25lcikgICMgVGjDqm0gaMOgbmcgdsOgIGPhu5l0IHThu5VuZw0KYGBgDQoNCiMjIyMgKio1LjcuMiBUw61uaCB0b8OhbiBSUioqDQpgYGB7cn0NCiMgVHLDrWNoIHh14bqldCBz4buRIGxp4buHdSB04burIGLhuqNuZyBjaMOpbw0KQSA8LSB0YWJsZV9nZW5kZXJfaG9tZW93bmVyWyJGIiwgIlkiXSAgIyBu4buvIHPhu58gaOG7r3UgbmjDoA0KQiA8LSB0YWJsZV9nZW5kZXJfaG9tZW93bmVyWyJGIiwgIk4iXSAgICMgbuG7ryBraMO0bmcgc+G7nyBo4buvdSBuaMOgDQpDIDwtIHRhYmxlX2dlbmRlcl9ob21lb3duZXJbIk0iLCAiWSJdICAjIG5hbSBz4bufIGjhu691IG5ow6ANCkQgPC0gdGFibGVfZ2VuZGVyX2hvbWVvd25lclsiTSIsICJOIl0gICAjIG5hbSBraMO0bmcgc+G7nyBo4buvdSBuaMOgDQoNCiMgVMOtbmggUlINClJSIDwtIChBIC8gKEEgKyBCKSkgLyAoQyAvIChDICsgRCkpDQpSUg0KDQojIFRyb25nIMSRw7M6DQojIEEgLyAoQSArIEIpOiBU4bu3IGzhu4cgbuG7ryBz4bufIGjhu691IG5ow6AgdHJvbmcgdOG7lW5nIHPhu5EgbuG7ry4NCiMgQyAvIChDICsgRCk6IFThu7cgbOG7hyBuYW0gc+G7nyBo4buvdSBuaMOgIHRyb25nIHThu5VuZyBz4buRIG5hbS4NCg0KDQpgYGANCg0KYGBge3J9DQpyaXNrcmF0aW8odGFibGVfZ2VuZGVyX2hvbWVvd25lciwgbWV0aG9kPSJ3YWxkIikNCmBgYA0KDQoNCiMjIyMgKio1LjcuMyBOaOG6rW4geMOpdCAqKg0KDQotIE5ow7NtIG7hu68gKEYpIMSRxrDhu6NjIGNo4buNbiBsw6BtIG5ow7NtIHRoYW0gY2hp4bq/dSwgbsOqbiBSUiA9IDEuIMSQaeG7gXUgbsOgeSBjw7MgbmdoxKlhIGNow7puZyB0YSBzbyBzw6FuaCBuZ3V5IGPGoSBz4bufIGjhu691IG5ow6AgY+G7p2EgbmjDs20gbmFtIChNKSB24bubaSBuaMOzbSBu4buvLg0KDQotIE5ow7NtIG5hbSAoTSkgY8OzIFJSID0gMC45ODIzIG5naMSpYSBsw6A6DQoNCiAgLSBU4bu3IGzhu4cgc+G7nyBo4buvdSBuaMOgIOG7nyBuYW0gYuG6sW5nIGtob+G6o25nIDk4LjIzJSBzbyB24bubaSBu4buvLg0KSGF5IG7Ds2kgY8OhY2gga2jDoWMsIG5ndXkgY8ahICh4w6FjIHN14bqldCkgc+G7nyBo4buvdSBuaMOgIGPhu6dhIG5hbSB0aOG6pXAgaMahbiBu4buvIGtob+G6o25nIDEuNzclLg0KDQotIEtob+G6o25nIHRpbiBj4bqteSA5NSUgY+G7p2EgUlIgbMOgICgwLjk1NjIg4oCTIDEuMDA5Mik6DQoNCiAgLSBLaG/huqNuZyBuw6B5IGJhbyBn4buTbSBz4buRIDEsIG5naMSpYSBsw6AgZ2nDoSB0cuG7iyBSUiB0aOG7sWMgc+G7sSBjw7MgdGjhu4MgbMOgIDEgKGtow7RuZyBjw7Mga2jDoWMgYmnhu4d0KS4NCg0KICAtIEtob+G6o25nIHRpbiBj4bqteSBjaG8gdGEgYmnhur90IG3hu6ljIMSR4buZIGNow61uaCB4w6FjIGPhu6dhIMaw4bubYyBsxrDhu6NuZyBSUjsga2hv4bqjbmcgcuG7mW5nIHbDoCBiYW8gZ+G7k20gMSBjaG8gdGjhuqV5IHPhu7Ega2jDtG5nIGNo4bqvYyBjaOG6r24gdsOgIHRoaeG6v3UgYuG6sW5nIGNo4bupbmcgduG7gSBz4buxIGtow6FjIGJp4buHdCB0aOG7sWMgc+G7sS4NCg0KLSBUYSBjw7MgZ2nhuqMgdGh1eeG6v3Q6DQoNCiAgLSBI4oKAOiBU4bu3IGzhu4cgc+G7nyBo4buvdSBuaMOgIOG7nyBOYW0gdsOgIE7hu68ga2jDtG5nIGtow6FjIG5oYXUsIHThu6ljIGzDoCBSUiA9IDENCg0KICAtIEjigoE6IFThu7cgbOG7hyBz4bufIGjhu691IG5ow6Ag4bufIE5hbSB2w6AgTuG7ryBraMOhYyBuaGF1LCBSUiDiiaAgMQ0KDQotIHAtdmFsdWUgdOG7qyBjw6FjIGtp4buDbSDEkeG7i25oOg0KDQogIC0gQ2hpLXNxdWFyZSB0ZXN0OiBwLXZhbHVlID0gMC4xOTUwODg0DQogIC0gRmlzaGVyJ3MgZXhhY3QgdGVzdDogcC12YWx1ZSA9IDAuMTk2NDgzMw0KICAtIE1pZHBvaW50IGV4YWN0IHRlc3Q6IHAtdmFsdWUgPSAwLjE5NTE1OA0KICANCi0gTmjhuq1uIHjDqXQgduG7gSBwLXZhbHVlOg0KDQogIC0gQ8OhYyBnacOhIHRy4buLIHAtdmFsdWUgdOG7qyBjaGktc3F1YXJlLCBGaXNoZXIncyBleGFjdCwgdsOgIG1pZHAuZXhhY3QgxJHhu4F1IGzhu5tuIGjGoW4gbeG7qWMgw70gbmdoxKlhIHRow7RuZyB0aMaw4budbmcgbMOgIDAuMDUuDQogIC0gxJBp4buBdSBuw6B5IGPDsyBuZ2jEqWEgbMOgIGtow7RuZyBjw7MgxJHhu6cgYuG6sW5nIGNo4bupbmcgxJHhu4MgYsOhYyBi4buPIGdp4bqjIHRodXnhur90IG51bGwgaHlwb3RoZXNpcywgdOG7qWMgbMOgIGtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgxJHDoW5nIGvhu4MgduG7gSB04buJIGzhu4cgdMOsbmggdHLhuqFuZyBow7RuIG5ow6JuIGPDsyAoWSkgZ2nhu69hIG5ow7NtIG5hbSB2w6AgbuG7ry4gU+G7sSBraMOhYyBiaeG7h3QgZ2nhu69hIGhhaSBuaMOzbSBjw7MgdGjhu4MgbMOgIGRvIG5n4bqrdSBuaGnDqm4uDQoNCioqTmjhuq1uIHjDqXQgdOG7lW5nIHF1w6F0KioNCg0KLSBSZWxhdGl2ZSBSaXNrIGNobyB0aOG6pXkgcuG6sW5nIHThu7cgbOG7hyBz4buxIGtp4buHbiBnaeG7r2EgaGFpIG5ow7NtIGdp4bubaSB0w61uaCBGIChwaOG7pSBu4buvKSB2w6AgTSAobmFtIGdp4bubaSkgcuG6pXQgZ+G6p24gbmhhdSwgduG7m2kgUlIgY+G7p2EgbmjDs20gTSBsw6AgMC45ODIzLCBjaG8gdGjhuqV5IGtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgcsO1IHLhu4d0IGdp4buvYSB04bu3IGzhu4cgdMOsbmggdHLhuqFuZyBow7RuIG5ow6JuIGPDsyAoWSkgZ2nhu69hIGhhaSBuaMOzbS4NCg0KLSBDb25maWRlbmNlIEludGVydmFsIGNobyBuaMOzbSBuYW0gZ2nhu5tpIChNKSBiYW8gZ+G7k20gZ2nDoSB0cuG7iyAxLCDEkWnhu4F1IG7DoHkgaOG7lyB0cuG7oyBnaeG6oyB0aHV54bq/dCBy4bqxbmcga2jDtG5nIGPDsyBz4buxIGtow6FjIGJp4buHdCB0aOG7kW5nIGvDqiDEkcOhbmcga+G7gyBnaeG7r2EgdOG7tyBs4buHIHPhu7Ega2nhu4duIOG7nyBuYW0gdsOgIG7hu68uDQoNCi0gQ8OhYyBraeG7g20gxJHhu4tuaCB0aOG7kW5nIGvDqiAoQ2hpLXNxdWFyZSwgRmlzaGVyJ3MgZXhhY3QgdGVzdCkgxJHhu4F1IGNobyBwLXZhbHVlIGzhu5tuIGjGoW4gMC4wNSwgxJFp4buBdSBuw6B5IGPDsyBuZ2jEqWEgbMOgIGtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3Qgw70gbmdoxKlhIGdp4buvYSBoYWkgbmjDs20gZ2nhu5tpIHTDrW5oIHbhu4EgdOG7tyBs4buHIHTDrG5oIHRy4bqhbmcgaMO0biBuaMOibiBjw7MgKFkpLg0KDQotIFTDs20gbOG6oWksIGPDoWMga+G6v3QgcXXhuqMgbsOgeSBjaG8gdGjhuqV5IGtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgcsO1IHLhu4d0IHbDoCBjw7Mgw70gbmdoxKlhIGdp4buvYSB04bu3IGzhu4cgdMOsbmggdHLhuqFuZyBow7RuIG5ow6JuIGPDsyAoWSkg4bufIHBo4bulIG7hu68gdsOgIG5hbSBnaeG7m2kuDQoNCioqS+G6v3QgbHXhuq1uKioNCg0KROG7sWEgdHLDqm4gY8OhYyBr4bq/dCBxdeG6oyB04burIFJlbGF0aXZlIFJpc2sgdsOgIGPDoWMga2nhu4NtIMSR4buLbmggdGjhu5FuZyBrw6osIGNow7puZyB0YSBjw7MgdGjhu4Mga+G6v3QgbHXhuq1uIHLhurFuZyBraMO0bmcgY8OzIHPhu7Ega2jDoWMgYmnhu4d0IMSRw6FuZyBr4buDIGdp4buvYSBuYW0gdsOgIG7hu68gduG7gSB04bu3IGzhu4cgdMOsbmggdHLhuqFuZyBow7RuIG5ow6JuIGPDsyAoWSkuIEThu68gbGnhu4d1IGNobyB0aOG6pXkgcuG6sW5nIHThu7cgbOG7hyBow7RuIG5ow6JuIGPDsyAoWSkg4bufIG5hbSB2w6AgbuG7ryBsw6Aga2jDoSB0xrDGoW5nIMSR4buTbmcgdsOgIGtow7RuZyBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgdGjhu5FuZyBrw6ogxJHDoW5nIGvhu4MuDQoNCiMjIyAqKjUuOCBNYXJpdGFsU3RhdHVzIHbDoCBIb21lb3duZXIqKg0KDQpgYGB7cn0NCnRhYmxlX21hcmlfaG9tZSA8LSB0YWJsZShkYXQkTWFyaXRhbFN0YXR1cywgZGF0JEhvbWVvd25lcikNCmFkZG1hcmdpbnModGFibGVfbWFyaV9ob21lKQ0KYGBgDQoNCmBgYHtyfQ0Kcmlza3JhdGlvKHRhYmxlX21hcmlfaG9tZSkNCmBgYA0KDQpOaMOzbSDEkcOjIGvhur90IGjDtG4gKE0pIMSRxrDhu6NjIGTDuW5nIGzDoG0gbmjDs20gdGhhbSBjaGnhur91IHbhu5tpIFJSID0gMS4NCg0KTmjDs20gxJHhu5ljIHRow6JuIChTKSBjw7M6DQoNCi0gUlIgPSAwLjYxMTQgKGtob+G6o25nIHRpbiBj4bqteSA5NSU6IDAuNTk0MiDigJMgMC42MjkyKS5Ww6wga2hv4bqjbmcgbsOgeSBraMO0bmcgY2jhu6lhIDEsIG7Dqm4ga+G6v3QgcXXhuqMgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqi4gDQoNCi0gxJBp4buBdSBuw6B5IGPDsyBuZ2jEqWEgbMOgIG5nxrDhu51pIMSR4buZYyB0aMOibiBjw7Mgbmd1eSBjxqEgc+G7nyBo4buvdSBuaMOgIHRo4bqlcCBoxqFuIGtob+G6o25nIDM4Ljg2JSBzbyB24bubaSBuZ8aw4budaSDEkcOjIGvhur90IGjDtG4uDQoNCi0gS2hv4bqjbmcgdGluIGPhuq15IDk1JSBraMOhIGjhurlwIHbDoCBraMO0bmcgYmFvIGfhu5NtIGdpw6EgdHLhu4sgMSwgY2hvIHRo4bqleSBz4buxIGtow6FjIGJp4buHdCBuw6B5IGzDoCByw7UgcsOgbmcgdsOgIMSRw6FuZyB0aW4gY+G6rXkuDQoNClRhIHRo4buxYyBoaeG7h24ga2nhu4NtIMSR4buLbmggZ2nhuqMgdGh1eeG6v3QgbmjGsCBzYXU6DQoNCi0gKipI4oKAOioqIFThu7cgbOG7hyBz4bufIGjhu691IG5ow6AgY+G7p2EgbmfGsOG7nWkgxJHDoyBr4bq/dCBow7RuIHbDoCDEkeG7mWMgdGjDom4gbMOgIG5oxrAgbmhhdSAoUlIgPSAxKS4NCg0KLSAqKkjigoE6KiogVOG7tyBs4buHIHPhu58gaOG7r3UgbmjDoCBj4bunYSBoYWkgbmjDs20gbMOgIGtow6FjIG5oYXUgKFJSIOKJoCAxKS4NCg0KS+G6v3QgcXXhuqMgcGjDom4gdMOtY2ggY2jhu4kgcmEgUlIgPSAwLjYxMTQgduG7m2kga2hv4bqjbmcgdGluIGPhuq15IDk1JSAoMC41OTQyIOKAkyAwLjYyOTIpLCB2w6AgY8OhYyBraeG7g20gxJHhu4tuaCBjaG8gcC12YWx1ZSDiiYggMCA8IG3hu6ljIMO9IG5naMSpYSAwLjA1LiBEbyDEkcOzLCB0YSBiw6FjIGLhu48gSOKCgCB2w6Aga+G6v3QgbHXhuq1uIHLhurFuZyB0w6xuaCB0cuG6oW5nIGjDtG4gbmjDom4gY8OzIG3hu5FpIGxpw6puIGjhu4cgdGjhu5FuZyBrw6ogcsO1IHLhu4d0IHbhu5tpIGto4bqjIG7Eg25nIHPhu58gaOG7r3UgbmjDoC4gTmfGsOG7nWkgxJHhu5ljIHRow6JuIGPDsyBraOG6oyBuxINuZyBz4bufIGjhu691IG5ow6AgdGjhuqVwIGjGoW4gxJHDoW5nIGvhu4Mgc28gduG7m2kgbmfGsOG7nWkgxJHDoyBr4bq/dCBow7RuLg0KDQojIyAqKlBI4bqmTiA2OiBPRERTIFJBVElPKioNCg0KIyMjICoqNi4xIEdlbmRlciB2w6AgSG9tZW93bmVyKioNCg0KYGBge3J9DQojIEPDoGkgdsOgIGfhu41pIGfDs2kNCmluc3RhbGwucGFja2FnZXMoImVwaXRvb2xzIikNCmxpYnJhcnkoZXBpdG9vbHMpDQpgYGANCg0KDQpgYGB7cn0NCiMgVMOtbmggb2RkcyByYXRpbw0Kb3JfbWggPC0gb2Rkc3JhdGlvKHRhYmxlX21hcmlfaG9tZSwgbWV0aG9kID0gIndhbGQiLCBjb25mLmxldmVsID0gMC45NSkNCm9yX21oJG1lYXN1cmUNCmBgYA0KDQoNCi0gS+G6v3QgcXXhuqMgcGjDom4gdMOtY2ggT2RkcyBSYXRpbyBnaeG7r2EgZ2nhu5tpIHTDrW5oIHbDoCBraOG6oyBuxINuZyBz4bufIGjhu691IG5ow6AgbmjGsCBzYXU6DQoNCiAgLSBPZGRzIFJhdGlvIChPUikgPSAwLjk1Ng0KDQogIC0gS2hv4bqjbmcgdGluIGPhuq15IDk1JTogdOG7qyAwLjg5NCDEkeG6v24gMS4wMjMNCg0KLSBEaeG7hW4gZ2nhuqNpOg0KDQogIC0gR2nDoSB0cuG7iyBPUiA9IDAuOTU2IGNobyB0aOG6pXkgb2RkcyAoY8ahIGjhu5lpKSBz4bufIGjhu691IG5ow6AgY+G7p2EgbuG7ryB0aOG6pXAgaMahbiBuYW0ga2hv4bqjbmcgNC40JS4gVHV5IG5oacOqbiwgc+G7sSBraMOhYyBiaeG7h3QgbMOgIG5o4buPLg0KDQogIC0gRG8ga2hv4bqjbmcgdGluIGPhuq15IDk1JSBjaOG7qWEgZ2nDoSB0cuG7iyAxLCBuw6puOg0KDQogICAgLSBLaMO0bmcgY8OzIMSR4bunIGLhurFuZyBjaOG7qW5nIHRo4buRbmcga8OqIMSR4buDIGvhur90IGx14bqtbiBy4bqxbmcgZ2nhu5tpIHTDrW5oIOG6o25oIGjGsOG7n25nIMSR4bq/biBvZGRzIHPhu58gaOG7r3UgbmjDoC4NCg0KLSAqKkvhur90IGx14bqtbjoqKiBN4bq3YyBkw7kgbuG7ryBjw7MgduG6uyBjw7Mgb2RkcyBz4bufIGjhu691IG5ow6AgdGjhuqVwIGjGoW4gbmFtIG3hu5l0IGNow7p0LCBuaMawbmcgc+G7sSBraMOhYyBiaeG7h3Qga2jDtG5nIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ouIERvIMSRw7MsIGdp4bubaSB0w61uaCBraMO0bmcgcGjhuqNpIGzDoCB54bq/dSB04buRIOG6o25oIGjGsOG7n25nIHLDtSBy4buHdCDEkeG6v24ga2jhuqMgbsSDbmcgc+G7nyBo4buvdSBuaMOgIHRyb25nIGThu68gbGnhu4d1IGhp4buHbiB04bqhaS4NCg0KIyMjICoqNi4yIE1hcml0YWxTdGF0dXMgdsOgIEhvbWVvd25lcioqDQoNCmBgYHtyfQ0KIyBUw61uaCBPUiB2w6Aga2hv4bqjbmcgdGluIGPhuq15DQpvcl9taCA8LSBvZGRzcmF0aW8odGFibGVfbWFyaV9ob21lLCBtZXRob2QgPSAid2FsZCIsIGNvbmYubGV2ZWwgPSAwLjk1KQ0Kb3JfbWgkbWVhc3VyZQ0KYGBgDQoNCi0gT2RkcyBSYXRpbyAoT1IpID0gMC4yODMNCg0KLSBLaG/huqNuZyB0aW4gY+G6rXkgOTUlOiB04burIDAuMjYzIMSR4bq/biAwLjMwNA0KDQotIERp4buFbiBnaeG6o2k6DQoNCiAgLSBHacOhIHRy4buLIE9kZHMgUmF0aW8gMC4yODMgY8OzIG5naMSpYSBsw6Agb2RkcyBz4bufIGjhu691IG5ow6AgY+G7p2EgbmjDs20gbmfGsOG7nWkgxJHhu5ljIHRow6JuIGNo4buJIGLhurFuZyBraG/huqNuZyAyOC4zJSBvZGRzIGPhu6dhIG5ow7NtIG5nxrDhu51pIMSRw6Mga+G6v3QgaMO0bi4gTsOzaSBjw6FjaCBraMOhYywgbmfGsOG7nWkgxJHDoyBr4bq/dCBow7RuIGPDsyBraOG6oyBuxINuZyBz4bufIGjhu691IG5ow6AgY2FvIGjGoW4gZ+G6p24gMy41IGzhuqduIHNvIHbhu5tpIG5nxrDhu51pIMSR4buZYyB0aMOibiAodsOsIDEvMC4yODPiiYgzLjUzKS4NCg0KICAtIEtob+G6o25nIHRpbiBj4bqteSA5NSUgY2hvIE9kZHMgUmF0aW8gbuG6sW0gaG/DoG4gdG/DoG4gZMaw4bubaSAxLCB04burIDAuMjYzIMSR4bq/biAwLjMwNCwgxJFp4buBdSBuw6B5IGNobyB0aOG6pXkgc+G7sSBraMOhYyBiaeG7h3QgduG7gSBvZGRzIHPhu58gaOG7r3UgbmjDoCBnaeG7r2EgaGFpIG5ow7NtIGzDoCBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqIHbDoCBraMO0bmcgcGjhuqNpIGRvIG5n4bqrdSBuaGnDqm4uDQoNCi0gS+G6v3QgbHXhuq1uOg0KDQpE4buvIGxp4buHdSBjaG8gdGjhuqV5IHTDrG5oIHRy4bqhbmcgaMO0biBuaMOibiDhuqNuaCBoxrDhu59uZyByw7UgcsOgbmcgxJHhur9uIGto4bqjIG7Eg25nIHPhu58gaOG7r3UgbmjDoC4gTmfGsOG7nWkgxJHDoyBr4bq/dCBow7RuIGPDsyBvZGRzIHPhu58gaOG7r3UgbmjDoCBjYW8gaMahbiDEkcOhbmcga+G7gyBzbyB24bubaSBuZ8aw4budaSDEkeG7mWMgdGjDom4sIHbDoCBraMOhYyBiaeG7h3QgbsOgeSBsw6AgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiB24bubaSBt4bupYyDEkeG7mSB0aW4gY+G6rXkgOTUlLiDEkMOieSBsw6AgYuG6sW5nIGNo4bupbmcgY2hvIHRo4bqleSB2aeG7h2MgxJHDoyBr4bq/dCBow7RuIGPDsyBsacOqbiBxdWFuIGNo4bq3dCBjaOG6vSDEkeG6v24ga2jhuqMgbsSDbmcgc+G7nyBo4buvdSBuaMOgIHRyb25nIG3huqt1IGto4bqjbyBzw6F0IG7DoHkuDQoNCiMjIyAqKjYuMyBHZW5kZXIgdsOgIE1hcml0YWxTdGF0dXMqKg0KDQpgYGB7cn0NCiMgVMOtbmggT1IgdsOgIGtob+G6o25nIHRpbiBj4bqteQ0Kb3JfZ20gPC0gb2Rkc3JhdGlvKHRhYmxlX2dlbmRlcl9tYXJpdGFsLCBtZXRob2QgPSAid2FsZCIsIGNvbmYubGV2ZWwgPSAwLjk1KQ0Kb3JfZ20kbWVhc3VyZQ0KYGBgDQoNCi0gT2RkcyBSYXRpbyAoT1IpID0gMS4xMjENCg0KLSBLaG/huqNuZyB0aW4gY+G6rXkgOTUlOiB04burIDEuMDQ5IMSR4bq/biAxLjE5OA0KDQotIERp4buFbiBnaeG6o2k6DQoNCiAgLSBHacOhIHRy4buLIE9kZHMgUmF0aW8gMS4xMjEgY8OzIG5naMSpYSBsw6Agb2RkcyAoY8ahIGjhu5lpKSDEkcOjIGvhur90IGjDtG4gY+G7p2EgbuG7ryBjYW8gaMahbiBvZGRzIMSRw6Mga+G6v3QgaMO0biBj4bunYSBuYW0ga2hv4bqjbmcgMTIuMSUsIGhheSBuw7NpIGPDoWNoIGtow6FjLCBu4buvIGPDsyBraOG6oyBuxINuZyDEkcOjIGvhur90IGjDtG4gY2FvIGjGoW4gbmFtIGtob+G6o25nIDEuMTIgbOG6p24uDQoNCiAgLSBLaG/huqNuZyB0aW4gY+G6rXkgOTUlIHThu6sgMS4wNDkgxJHhur9uIDEuMTk4LCBraMO0bmcgY2jhu6lhIGdpw6EgdHLhu4sgMSwgY2jhu6luZyB04buPIHPhu7Ega2jDoWMgYmnhu4d0IG7DoHkgbMOgIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ouIMSQaeG7gXUgbsOgeSBuZ2jEqWEgbMOgIHbhu5tpIG3hu6ljIHRpbiBj4bqteSA5NSUsIGNow7puZyB0YSBjaOG6r2MgY2jhuq9uIHLhurFuZyB04bu3IGzhu4cgxJHDoyBr4bq/dCBow7RuIGPhu6dhIG7hu68gdsOgIG5hbSBraMO0bmcgYuG6sW5nIG5oYXUuDQoNCi0gS+G6v3QgbHXhuq1uOg0KDQpE4buvIGxp4buHdSBjaG8gdGjhuqV5IG7hu68gZ2nhu5tpIGPDsyBvZGRzIMSRw6Mga+G6v3QgaMO0biBjYW8gaMahbiBuYW0gZ2nhu5tpIGtob+G6o25nIDEyJSwgdMawxqFuZyDEkcawxqFuZyB24bubaSB2aeG7h2MgbuG7ryBjw7Mga2jhuqMgbsSDbmcgxJHDoyBr4bq/dCBow7RuIGfhuqVwIDEuMTIgbOG6p24gc28gduG7m2kgbmFtIGdp4bubaS4=

NHIỆM VỤ TUẦN 2

Võ Hoàng Tính

2025-05-25

PHẦN 1: TÌM HIỂU VÀ CHUẨN BỊ DỮ LIỆU

1.1 ĐỌC DỮ LIỆU

1.2. NỘI DUNG DỮ LIỆU

Chuyển đổi các cột sang factor

PHẦN 2: PHÂN TÍCH BIẾN ĐỊNH TÍNH

PRODUCTCATEGORY

ANNUALINCOME

COUNTRY

GENDER

MARITALSTATUS

HOMEOWNER

CITY

STATEORPROVINCE

PRODUCTFAMILY

PRODUCTDEPARTMENT

PHẦN 3: ƯỚC LƯỢNG KHOẢNG VÀ KIỂM ĐỊNH GIẢ THUYẾT CHO TỶ LỆ

3.1 Biến Gender

3.2 Biến Homeowner

3.3 Biến ProductFamily

PHẦN 4: PHÂN TÍCH MỐI QUAN HỆ GIỮA HAI BIẾN ĐỊNH TÍNH

4.1 Gender và ProductFamily

4.2 MaritalStatus và Homeowner

4.3 AnnualIncome và ProductCategory

Homeowner và AnnualIncome

4.5 City và Gender

PHẦN 5: RELATIVE RISK

5.1. Giới thiệu

5.2. Cấu trúc và xác suất của bảng ngẫu nhiên

5.2.1. Phân phối Poisson

5.2.2. Phân phối Multinomial

5.3. Các chỉ số đo mối liên hệ

5.3.1. Hiệu tỷ lệ (Risk Difference - RD)

5.3.2. Tỷ số nguy cơ (Relative Risk - RR)

5.3.3. Tỷ số chênh (Odds Ratio - OR)

5.3.3.1. Odds

5.3.3.1. Odds Ratio - OR

5.4. Khoảng ước lượng tin cậy cho ODDS RATIO

5.5. Ứng dụng trong kinh doanh: HIỆU QUẢ CỦA KHUYẾN MÃI

5.6 GENDER VÀ MARITALSTATUS

Tạo bảng chéo

Tính toán Relative Risk

5.6.3 Ý nghĩa

5.7 GENDER VÀ HOMEOWNER

5.7.1 Tạo bảng chéo

5.7.2 Tính toán RR

5.7.3 Nhận xét

5.8 MaritalStatus và Homeowner

PHẦN 6: ODDS RATIO

6.1 Gender và Homeowner

6.2 MaritalStatus và Homeowner

6.3 Gender và MaritalStatus