library(epitools)
library(tidyverse)
library(dplyr)
library(knitr)
library(readxl)
library(kableExtra)

CHƯƠNG 1: TỔNG QUAN NGHIÊN CỨU

1.1. Lý do chọn đề tài

Giấc ngủ đóng vai trò thiết yếu trong việc duy trì sức khỏe thể chất và tinh thần, ảnh hưởng trực tiếp đến hiệu suất làm việc, khả năng tập trung và chất lượng cuộc sống. Tuy nhiên, trong bối cảnh xã hội hiện đại với nhịp sống nhanh và áp lực công việc ngày càng gia tăng, chất lượng giấc ngủ đang trở thành một vấn đề đáng quan ngại, đặc biệt ở các nhóm nghề nghiệp có đặc thù lao động khác nhau.

Tại Ấn Độ – một quốc gia có lực lượng lao động đa dạng trải dài trên nhiều lĩnh vực – việc phân tích chất lượng giấc ngủ theo từng nhóm ngành nghề có thể mang lại những hiểu biết quan trọng, phục vụ cho các nghiên cứu sức khỏe cộng đồng và chính sách lao động. Dù đã có nhiều nghiên cứu quốc tế bàn về mối liên hệ giữa nghề nghiệp và giấc ngủ, vẫn còn thiếu những phân tích cụ thể theo ngành nghề trong từng bối cảnh quốc gia, đặc biệt là tại các nước đang phát triển như Ấn Độ.

Dựa trên một bộ dữ liệu được thu thập trong tháng 4 năm 2025, gồm 13 biến độc lập liên quan đến các yếu tố cá nhân và môi trường như độ tuổi, tình trạng trị liệu, thói quen thức giấc giữa đêm, cảm giác nhiệt độ khi ngủ,… và một biến phụ thuộc là nghề nghiệp, đề tài “Phân tích các yếu tố ảnh hưởng đến giấc ngủ của các ngành nghề khác nhau” nhằm mục tiêu khám phá sự khác biệt về chất lượng giấc ngủ giữa các nhóm nghề nghiệp cũng như xác định các yếu tố có khả năng ảnh hưởng đến tình trạng này.

Với những lý do trên, tác giả đã chọn đề tài “PHÂN TÍCH CÁC YẾU TỐ ẢNH HƯỞNG ĐẾN GIẤC NGỦ CỦA CÁC NGÀNH NGHỀ KHÁC NHAU”.

1.2. Mục tiêu nghiên cứu

Nghiên cứu và làm rõ các yếu tố ảnh hưởng đến giấc ngủ của các nhóm ngành nghề khác nhau tại Ấn Độ.
Thống kê mô tả các biến độc lập và biến nghề nghiệp để hiểu đặc điểm tổng quát của mẫu khảo sát.
Tính toán Odds Ratio (OR) và Relative Risk (RR) nhằm đánh giá mức độ ảnh hưởng của từng yếu tố đến giấc ngủ.
Kiểm định Chi-square để xác định mối liên hệ có ý nghĩa thống kê giữa các yếu tố với nghề nghiệp.
Ước lượng mô hình hồi quy logistic để xác định các yếu tố tác động đáng kể đến tình trạng giấc ngủ trong từng nhóm nghề nghiệp.
Đề xuất các hàm ý thực tiễn từ kết quả phân tích, làm cơ sở cho các chính sách nâng cao chất lượng giấc ngủ và sức khỏe tinh thần trong từng ngành nghề cụ thể tại Ấn Độ.

1.3. Đối tượng và phạm vi nghiên cứu

Đối tượng nghiên cứu: Các yếu tố ảnh hưởng tới giấc ngủ trong các ngành nghề khác nhau.

Phạm vi nghiên cứu: 1 tháng, cụ thể tháng 4 năm 2025 tại Ấn Độ.

1.4. Phương pháp nghiên cứu

Bài nghiên cứu sử dụng kết hợp giữa 2 phương pháp nghiên cứu định lượng và định tính, bao gồm:

Thống kê mô tả
Kiểm định Chi-square: kiểm định mối liên hệ, đánh giá sự phụ thuộc và xác định các yếu tố có thể ảnh hưởng đến giấc ngủ
Các mô hình hồi quy đơn biến (Logistic Regression), đồng thời ước lượng tỷ số Odds Ratio
Các mô hình hồi quy đa biến: Logistic Regression, Probit Regression và Complementary Log-Log Regression

Phần mềm: R Studio

1.5. Kết cấu đề tài

Bài nghiên cứu bao gồm 3 chương:

Chương 1: Giới thiệu tổng quan về đề tài

Chương 2: Giới thiệu về bộ dữ liệu

Chương 3: Kết quả nghiên cứu và nhận xét

CHƯƠNG 2: GIỚI THIỆU VỀ BỘ DỮ LIỆU

data <- read_excel("C:/Users/Admin/Downloads/data PTDLĐT.xlsx")

#hiển thị cấu trúc dữ liệu
str(data)

## tibble [163 × 10] (S3: tbl_df/tbl/data.frame)
##  $ Gender        : chr [1:163] "Male" "Female" "Female" "Male" ...
##  $ Age           : chr [1:163] "18-20" "18-20" "18-20" "18-20" ...
##  $ Major         : chr [1:163] "Mechanical" "Biotech" "Computer Science" "Computer Science" ...
##  $ Hours sleep   : chr [1:163] "7-8h" "7-8h" ">8h" "7-8h" ...
##  $ 30 min        : chr [1:163] "0/month" "2/week" "1/week" "2/week" ...
##  $ Wake up middle: chr [1:163] "1/week" "0/month" "0/month" "0/month" ...
##  $ Bathroom      : chr [1:163] "1/week" "0/month" "0/month" "1/week" ...
##  $ Cold          : chr [1:163] "0/month" "2/week" "1/week" "0/month" ...
##  $ Hot           : chr [1:163] "1/week" "2/week" "1/week" "0/month" ...
##  $ Dreams        : chr [1:163] "0/month" "1/week" "1/week" "0/month" ...

Dữ liệu bộ bao gồm 163 quan sát và 10 biến định tính được lấy trên trang web Kangle, cụ thể:

Biến phụ thuộc:

Major: Đa dạng ngành nghề của từng người

Biến độc lập:

Gender: Giới tính
- Male: Nam
- Female: Nữ
Age: Độ tuổi
- 18-20: Từ 18 đến 20 tuổi
- 21-23: Từ 21 đến 23 tuổi
Hours sleep: Trung bình thời gian ngủ
- <4h: Ít hơn 4 tiếng
- 4-6h: Từ 4 đến 6 tiếng
- 7-8h: Từ 7 đến 8 tiếng
- ‘>8h’: Nhiều hơn 8 tiếng

Các tác động như:

30min: Không thể ngủ trong 30 phút đầu
Wake up middle: Bị tỉnh giấc nửa đêm hoặc rạng sáng
Bathroom: Mắc đi vệ sinh trong lúc ngủ
Cold: Cảm thấy lạnh trong lúc ngủ
Hot: Cảm thấy nóng trong lúc ngủ
Dreams: Gặp ác mộng

Với các tần suất như:

0/month: Không bị trong tháng
1/week: Bị 1 lần 1 tuần
2/week: Bị 2 lần 1 tuần
‘>3/week’: Bị nhiều hơn 3 lần 1 tuần

CHƯƠNG 3: KẾT QUẢ NGHIÊN CỨU

3.1. Tác động của các biến độc lập đến các biến phụ thuộc

3.1.1. Tác động của yếu tố không thể ngủ trong 30 phút đầu trong các ngành nghề

Bảng tần suất chéo và biểu đồ trực quan hóa

# Bảng tần suất giữa 30 min và Major
tansuat1 <- table(data$`30 min`, data$Major) 
phantram1 <- prop.table(tansuat1, margin = 1) * 100

bang1 <- as.data.frame(tansuat1)
bang1$Percentage <- round(as.vector(phantram1), 2)

colnames(bang1) <- c("30 min", "Major", "Frequency", "Percentage")
    
kable(bang1, caption = "Bảng 1: Mô tả thống kê của 2 biến 30 min và Major")

Bảng 1: Mô tả thống kê của 2 biến 30 min và Major
30 min	Major	Frequency	Percentage
>3/week	Biotech	5	14.71
0/month	Biotech	19	32.20
1/week	Biotech	8	25.00
2/week	Biotech	4	10.53
>3/week	Chemical	1	2.94
0/month	Chemical	3	5.08
1/week	Chemical	1	3.12
2/week	Chemical	1	2.63
>3/week	Civil	1	2.94
0/month	Civil	2	3.39
1/week	Civil	1	3.12
2/week	Civil	1	2.63
>3/week	Computer Application	1	2.94
0/month	Computer Application	2	3.39
1/week	Computer Application	1	3.12
2/week	Computer Application	1	2.63
>3/week	Computer Science	14	41.18
0/month	Computer Science	12	20.34
1/week	Computer Science	10	31.25
2/week	Computer Science	25	65.79
>3/week	Designer	3	8.82
0/month	Designer	7	11.86
1/week	Designer	3	9.38
2/week	Designer	1	2.63
>3/week	Fashion Technology	3	8.82
0/month	Fashion Technology	3	5.08
1/week	Fashion Technology	2	6.25
2/week	Fashion Technology	1	2.63
>3/week	Fine Arts	2	5.88
0/month	Fine Arts	2	3.39
1/week	Fine Arts	2	6.25
2/week	Fine Arts	1	2.63
>3/week	industrial engineering	1	2.94
0/month	industrial engineering	2	3.39
1/week	industrial engineering	1	3.12
2/week	industrial engineering	1	2.63
>3/week	Mechanical	2	5.88
0/month	Mechanical	4	6.78
1/week	Mechanical	2	6.25
2/week	Mechanical	1	2.63
>3/week	Statistics	1	2.94
0/month	Statistics	3	5.08
1/week	Statistics	1	3.12
2/week	Statistics	1	2.63

#Trực quan hóa cho mối quan hệ giữa 2 biến
bang1 <- data %>%
  group_by(Major, `30 min`) %>%
  summarise(Count = n(), .groups = "drop")

ggplot(bang1, aes(x = Major, y = Count, fill = `30 min`)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(title = "Hình 1: Mối quan hệ giữa yếu tố không thể ngủ trong 30 phút đầu và ngành nghề",
       x = "Ngành nghề", y = "Không thể ngủ trong 30 phút đầu", fill = "Tần suất") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

Nhìn chung, nhóm Công nghệ thông tin (Computer Science) có tỷ lệ gặp khó khăn về giấc ngủ cao nhất trong tất cả các mức tần suất. Cụ thể, có 16 người thuộc nhóm này gặp tình trạng khó ngủ 2 lần/tuần, chiếm 76.19% số người trong ngành. Ngoài ra, có tới 11 người chiếm 55% cho biết gặp vấn đề mất ngủ hơn 3 lần/tuần – đây là tỷ lệ cao nhất so với các nhóm ngành khác, cho thấy chất lượng giấc ngủ trong nhóm này ở mức báo động. Nguyên nhân có thể đến từ đặc thù công việc căng thẳng, áp lực cao, thời gian làm việc kéo dài và tiếp xúc thường xuyên với thiết bị điện tử.

Ngành Công nghệ sinh học (Biotech) cũng cho thấy phân bố mất ngủ ở mức cao nhưng rải rác ở nhiều mức độ. Cụ thể, có 12 người không gặp vấn đề về giấc ngủ (0 lần/tháng), chiếm 41.38% số người trong ngành; 5 người bị mất ngủ 1 lần/tuần (29.41%) và 4 người thường xuyên mất ngủ hơn 3 lần/tuần, chiếm 20%. Điều này cho thấy tình trạng mất ngủ tuy có xuất hiện nhưng không tập trung vào một mức tần suất cụ thể, phản ánh tính chất công việc hoặc lối sống đa dạng trong nhóm này.

Một số ngành như Cơ khí (Mechanical Engineering) và Xây dựng (Civil Engineering) có số lượng người tham gia vừa phải và phân bố tương đối đều giữa các mức độ tần suất. Ví dụ, ngành Cơ khí có 4 người không gặp vấn đề về giấc ngủ và 2 người mất ngủ từ 1–2 lần/tuần, không có trường hợp nào mất ngủ nghiêm trọng, cho thấy mức độ ảnh hưởng thấp hơn.

Ngược lại, các ngành như Ứng dụng máy tính (Computer Application), Thời trang (Fashion Technology), Nghệ thuật (Fine Arts) và Thống kê (Statistics) có số người tham gia rất ít (dưới 4 người/nhóm), khiến việc rút ra nhận định tổng quát gặp nhiều hạn chế. Tuy nhiên, phần lớn các cá nhân thuộc nhóm này đều không báo cáo tình trạng mất ngủ thường xuyên.

Cuối cùng, các ngành như Thiết kế (Designer) và Hóa học (Chemical) ghi nhận rất ít trường hợp mất ngủ nghiêm trọng. Chỉ có 1 người trong mỗi ngành này cho biết bị mất ngủ hơn 3 lần/tuần hoặc 1 lần/tuần, còn lại đều không gặp vấn đề về giấc ngủ.

Biểu đồ cột đi kèm đã trực quan hóa rõ sự khác biệt giữa các ngành, đặc biệt nổi bật ở nhóm Công nghệ thông tin và Công nghệ sinh học – với chiều cao cột lớn ở mức 2 lần/tuần và hơn 3 lần/tuần, phản ánh sự căng thẳng đặc thù trong công việc của hai nhóm ngành này.

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

Giả thuyết kiểm định:

\(H_0\): Hai biến độc lập không có sự tác động

\(H_1\): Hai biến có mối quan hệ

chisq.test(tansuat1)

## Warning in chisq.test(tansuat1): Chi-squared approximation may be incorrect

## 
##  Pearson's Chi-squared test
## 
## data:  tansuat1
## X-squared = 25.626, df = 30, p-value = 0.6941

Dựa trên kết quả kiểm định Chi-bình phương giữa hai biến 30 min và Major, giá trị thống kê Chi-squared đạt 25.626 với 30 bậc tự do và giá trị p tương ứng là 0.6941. Vì p-value > 0.05, nên không đủ bằng chứng để bác bỏ giả thuyết \(H_0\), tức là KHÔNG có mối liên hệ có ý nghĩa thống kê giữa bị khó ngủ trong 30 phút đầu và ngành nghề. Nói cách khác, hai biến này có xu hướng độc lập với nhau trong bộ dữ liệu. Tuy nhiên, cần lưu ý rằng kiểm định có cảnh báo về độ chính xác do một số ô trong bảng tần suất có giá trị kỳ vọng nhỏ, điều này có thể ảnh hưởng đến độ tin cậy của kết quả.

3.1.2. Tác động của yếu tố tỉnh dậy nửa đêm hoặc rạng sáng trong các ngành nghề

Bảng tần suất chéo và biểu đồ trực quan hóa

# Bảng tần suất giữa Wake up middle và Major
tansuat2 <- table(data$`Wake up middle`, data$Major) 
phantram2 <- prop.table(tansuat2, margin = 1) * 100

bang2 <- as.data.frame(tansuat2)
bang2$Percentage <- round(as.vector(phantram2), 2)

colnames(bang2) <- c("Wake up 1", "Major", "Frequency", "Percentage")
    
kable(bang2, caption = "Bảng 2: Mô tả thống kê của 2 biến Wake up middle và Major")

Bảng 2: Mô tả thống kê của 2 biến Wake up middle và Major
Wake up 1	Major	Frequency	Percentage
>3/month	Biotech	0	0.00
>3/week	Biotech	4	17.39
0/month	Biotech	21	28.00
1/week	Biotech	8	19.51
2/month	Biotech	0	0.00
2/week	Biotech	3	15.00
>3/month	Chemical	0	0.00
>3/week	Chemical	3	13.04
0/month	Chemical	1	1.33
1/week	Chemical	1	2.44
2/month	Chemical	0	0.00
2/week	Chemical	1	5.00
>3/month	Civil	0	0.00
>3/week	Civil	1	4.35
0/month	Civil	2	2.67
1/week	Civil	1	2.44
2/month	Civil	1	33.33
2/week	Civil	0	0.00
>3/month	Computer Application	0	0.00
>3/week	Computer Application	1	4.35
0/month	Computer Application	1	1.33
1/week	Computer Application	2	4.88
2/month	Computer Application	0	0.00
2/week	Computer Application	1	5.00
>3/month	Computer Science	0	0.00
>3/week	Computer Science	6	26.09
0/month	Computer Science	35	46.67
1/week	Computer Science	12	29.27
2/month	Computer Science	0	0.00
2/week	Computer Science	8	40.00
>3/month	Designer	0	0.00
>3/week	Designer	2	8.70
0/month	Designer	8	10.67
1/week	Designer	3	7.32
2/month	Designer	0	0.00
2/week	Designer	1	5.00
>3/month	Fashion Technology	0	0.00
>3/week	Fashion Technology	2	8.70
0/month	Fashion Technology	3	4.00
1/week	Fashion Technology	3	7.32
2/month	Fashion Technology	1	33.33
2/week	Fashion Technology	0	0.00
>3/month	Fine Arts	1	100.00
>3/week	Fine Arts	0	0.00
0/month	Fine Arts	1	1.33
1/week	Fine Arts	4	9.76
2/month	Fine Arts	1	33.33
2/week	Fine Arts	0	0.00
>3/month	industrial engineering	0	0.00
>3/week	industrial engineering	1	4.35
0/month	industrial engineering	2	2.67
1/week	industrial engineering	1	2.44
2/month	industrial engineering	0	0.00
2/week	industrial engineering	1	5.00
>3/month	Mechanical	0	0.00
>3/week	Mechanical	2	8.70
0/month	Mechanical	0	0.00
1/week	Mechanical	4	9.76
2/month	Mechanical	0	0.00
2/week	Mechanical	3	15.00
>3/month	Statistics	0	0.00
>3/week	Statistics	1	4.35
0/month	Statistics	1	1.33
1/week	Statistics	2	4.88
2/month	Statistics	0	0.00
2/week	Statistics	2	10.00

#Truc quan hoa
bang2 <- data %>%
  group_by(Major, `Wake up middle`) %>%
  summarise(Count = n(), .groups = "drop")

ggplot(bang2, aes(x = Major, y = Count, fill = `Wake up middle`)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(title = "Hình 2: Mối quan hệ giữa yếu tố tỉnh giấc nửa đêm hoặc rạng sáng và ngành nghề",
       x = "Ngành nghề", y = "Bị tỉnh giấc nửa đêm hoặc rạng sáng", fill = "Tần suất") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

Nhóm Công nghệ thông tin (Computer Science) tiếp tục cho thấy tỉ lệ rối loạn giấc ngủ cao, đặc biệt là tình trạng tỉnh giấc nửa đêm hoặc rạng sáng. Có tới 35 người trong nhóm này cho biết không gặp tình trạng này (0 lần/tháng), chiếm 46.67%. Tuy nhiên, bên cạnh đó, có 12 người (29.27%) cho biết họ tỉnh giấc 1 lần/tuần, và 8 người (40%) tỉnh giấc tới 2 lần/tuần, cho thấy vấn đề giấc ngủ vẫn khá phổ biến. Ngoài ra, có 6 người (26.09%) tỉnh giấc hơn 3 lần/tuần, cho thấy mức độ rối loạn giấc ngủ trung bình đến nặng ở ngành này.

Ngành Công nghệ sinh học (Biotech) cũng ghi nhận sự phân bố đáng kể ở nhiều mức tần suất: 21 người (chiếm 28.00%) không gặp tình trạng tỉnh giấc vào ban đêm, 8 người (19.51%) bị tỉnh giấc 1 lần/tuần, và 4 người (17.39%) bị ảnh hưởng hơn 3 lần/tuần, cùng với 3 người (15.00%) tỉnh giấc 2 lần/tuần. Những con số này cho thấy mức độ ảnh hưởng khá đa dạng trong nhóm này, tương tự xu hướng quan sát được trong biến khó ngủ.

Một số ngành như Thiết kế (Designer) và Thời trang (Fashion Technology) có sự phân bố tương đối đều nhưng với quy mô mẫu nhỏ. Trong ngành Thiết kế, 8 người (chiếm 10.67%) không bị tỉnh giấc, 3 người (7.32%) bị tỉnh giấc 1 lần/tuần, và chỉ có 2 người gặp tình trạng này ở mức nghiêm trọng hơn. Tương tự, ngành Thời trang có 3 người không bị ảnh hưởng, 3 người tỉnh giấc 1 lần/tuần, và 2 người từ mức 2 lần/tuần trở lên.

Ngành Mỹ thuật (Fine Arts) có quy mô mẫu nhỏ nhưng đặc biệt đáng chú ý với 1 trường hợp duy nhất bị tỉnh giấc hơn 3 lần/tháng (100% trong nhóm này), và 4 người khác tỉnh giấc 1 lần/tuần, cho thấy đây là nhóm có xu hướng giấc ngủ kém nhưng cần thêm dữ liệu để xác nhận chắc chắn.

Đối với các ngành như Cơ khí (Mechanical), Kỹ thuật công nghiệp (Industrial Engineering) và Thống kê (Statistics), tỷ lệ tỉnh giấc dao động chủ yếu quanh mức 1–2 lần/tuần, nhưng không có trường hợp nào vượt quá 3 lần/tháng, cho thấy mức độ ảnh hưởng thấp hơn so với các ngành công nghệ.

Biểu đồ minh họa trực quan đã thể hiện rõ sự chênh lệch, đặc biệt là ở nhóm Computer Science, với chiều cao cột vượt trội ở mức 0 lần/tháng và 1–2 lần/tuần, phản ánh sự phổ biến của tình trạng tỉnh giấc giữa đêm ở nhóm ngành này so với các ngành còn lại.

chisq.test(tansuat2)

## Warning in chisq.test(tansuat2): Chi-squared approximation may be incorrect

## 
##  Pearson's Chi-squared test
## 
## data:  tansuat2
## X-squared = 80.59, df = 50, p-value = 0.003946

Dựa trên kết quả kiểm định Chi-bình phương giữa hai biến Wake up middle và Major, giá trị thống kê Chi-squared đạt 80.59 với 50 bậc tự do và giá trị p tương ứng là 0.003946. Vì p-value < 0.05, nên không đủ bằng chứng để bác bỏ giả thuyết \(H_1\), tức là CÓ mối liên hệ có ý nghĩa thống kê giữa yếu tố bị tỉnh dậy giữa nửa đêm hoặc rạng sáng và ngành nghề. Nói cách khác, hai biến này có xu hướng phụ thuộc với nhau trong bộ dữ liệu. Tuy nhiên, cần lưu ý rằng kiểm định có cảnh báo về độ chính xác do một số ô trong bảng tần suất có giá trị kỳ vọng nhỏ, điều này có thể ảnh hưởng đến độ tin cậy của kết quả.

Tỷ lệ RR và OR

riskratio(tansuat2)

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## $data
##           
##            Biotech Chemical Civil Computer Application Computer Science
##   >3/month       0        0     0                    0                0
##   >3/week        4        3     1                    1                6
##   0/month       21        1     2                    1               35
##   1/week         8        1     1                    2               12
##   2/month        0        0     1                    0                0
##   2/week         3        1     0                    1                8
##   Total         36        6     5                    5               61
##           
##            Designer Fashion Technology Fine Arts industrial engineering
##   >3/month        0                  0         1                      0
##   >3/week         2                  2         0                      1
##   0/month         8                  3         1                      2
##   1/week          3                  3         4                      1
##   2/month         0                  1         1                      0
##   2/week          1                  0         0                      1
##   Total          14                  9         7                      5
##           
##            Mechanical Statistics Total
##   >3/month          0          0     1
##   >3/week           2          1    23
##   0/month           0          1    75
##   1/week            4          2    41
##   2/month           0          0     3
##   2/week            3          2    20
##   Total             9          6   163
## 
## $measure
##           risk ratio with 95% C.I.
##            estimate lower upper
##   >3/month        1    NA    NA
##   >3/week       NaN   NaN   NaN
##   0/month       NaN   NaN   NaN
##   1/week        NaN   NaN   NaN
##   2/month       NaN   NaN   NaN
##   2/week        NaN   NaN   NaN
## 
## $p.value
##           two-sided
##            midp.exact fisher.exact  chi.square
##   >3/month         NA           NA          NA
##   >3/week           1    0.2083333 0.007600391
##   0/month           1    0.1184211         NaN
##   1/week            1    0.5238095 0.669739288
##   2/month           1    1.0000000         NaN
##   2/week            1    0.2380952         NaN
## 
## $correction
## [1] FALSE
## 
## attr(,"method")
## [1] "Unconditional MLE & normal approximation (Wald) CI"

Kết quả phân tích cho thấy chỉ nhóm ngành nghề bị tỉnh giấc nửa đêm hoặc ráng sáng trên 3 lần một tuần có giá trị RR = 1, nhưng do chỉ có 1 người trong nhóm này nên không đủ cơ sở để suy luận. Các nhóm còn lại đều cho kết quả NaN do số lượng quá ít hoặc phân bố mất cân đối, khiến không thể tính toán RR và khoảng tin cậy. Dù nhóm >3/week có giá trị p trong kiểm định chi bình phương là 0.0076, nhưng các kiểm định khác lại không có ý nghĩa thống kê. Như vậy, chưa thể kết luận rõ ràng về mối liên hệ giữa ngành nghề và tần suất bị tỉnh giấc nửa đêm hoặc ráng sáng trong nghiên cứu này.

3.1.3. Tác động của yếu tố mắc vệ sinh trong lúc ngủ trong các ngành nghề

# Bảng tần suất giữa Bathroom và Major
tansuat3 <- table(data$`Bathroom`, data$Major) 
phantram3 <- prop.table(tansuat3, margin = 1) * 100

bang3 <- as.data.frame(tansuat3)
bang3$Percentage <- round(as.vector(phantram3), 2)

colnames(bang3) <- c("Bathroom", "Major", "Frequency", "Percentage")
    
kable(bang3, caption = "Bảng 3: Mô tả thống kê của 2 biến Bathroom và Major")

Bảng 3: Mô tả thống kê của 2 biến Bathroom và Major
Bathroom	Major	Frequency	Percentage
>3/month	Biotech	0	0.00
>3/week	Biotech	3	15.00
0/month	Biotech	14	21.88
0/week	Biotech	0	0.00
1/week	Biotech	10	19.61
2/month	Biotech	0	0.00
2/week	Biotech	9	36.00
>3/month	Chemical	0	0.00
>3/week	Chemical	1	5.00
0/month	Chemical	2	3.12
0/week	Chemical	0	0.00
1/week	Chemical	1	1.96
2/month	Chemical	0	0.00
2/week	Chemical	2	8.00
>3/month	Civil	1	100.00
>3/week	Civil	0	0.00
0/month	Civil	2	3.12
0/week	Civil	0	0.00
1/week	Civil	1	1.96
2/month	Civil	1	100.00
2/week	Civil	0	0.00
>3/month	Computer Application	0	0.00
>3/week	Computer Application	1	5.00
0/month	Computer Application	1	1.56
0/week	Computer Application	0	0.00
1/week	Computer Application	2	3.92
2/month	Computer Application	0	0.00
2/week	Computer Application	1	4.00
>3/month	Computer Science	0	0.00
>3/week	Computer Science	7	35.00
0/month	Computer Science	32	50.00
0/week	Computer Science	0	0.00
1/week	Computer Science	18	35.29
2/month	Computer Science	0	0.00
2/week	Computer Science	4	16.00
>3/month	Designer	0	0.00
>3/week	Designer	1	5.00
0/month	Designer	5	7.81
0/week	Designer	0	0.00
1/week	Designer	5	9.80
2/month	Designer	0	0.00
2/week	Designer	3	12.00
>3/month	Fashion Technology	0	0.00
>3/week	Fashion Technology	2	10.00
0/month	Fashion Technology	2	3.12
0/week	Fashion Technology	0	0.00
1/week	Fashion Technology	4	7.84
2/month	Fashion Technology	0	0.00
2/week	Fashion Technology	1	4.00
>3/month	Fine Arts	0	0.00
>3/week	Fine Arts	2	10.00
0/month	Fine Arts	1	1.56
0/week	Fine Arts	0	0.00
1/week	Fine Arts	3	5.88
2/month	Fine Arts	0	0.00
2/week	Fine Arts	1	4.00
>3/month	industrial engineering	0	0.00
>3/week	industrial engineering	1	5.00
0/month	industrial engineering	1	1.56
0/week	industrial engineering	0	0.00
1/week	industrial engineering	2	3.92
2/month	industrial engineering	0	0.00
2/week	industrial engineering	1	4.00
>3/month	Mechanical	0	0.00
>3/week	Mechanical	1	5.00
0/month	Mechanical	4	6.25
0/week	Mechanical	0	0.00
1/week	Mechanical	2	3.92
2/month	Mechanical	0	0.00
2/week	Mechanical	2	8.00
>3/month	Statistics	0	0.00
>3/week	Statistics	1	5.00
0/month	Statistics	0	0.00
0/week	Statistics	1	100.00
1/week	Statistics	3	5.88
2/month	Statistics	0	0.00
2/week	Statistics	1	4.00

#Truc quan hoa
bang3 <- data %>%
  group_by(Major, `Bathroom`) %>%
  summarise(Count = n(), .groups = "drop")

ggplot(bang3, aes(x = Major, y = Count, fill = `Bathroom`)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(title = "Hình 3: Mối quan hệ giữa yếu tố mắc vệ sinh trong lúc ngủ và ngành nghề",
       x = "Nganh nghe", y = "mắc vệ sinh trong lúc ngủ", fill = "Tần suất") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

Nhìn chung, ngành Công nghệ thông tin (Computer Science) ghi nhận tỷ lệ người gặp vấn đề phải đi vệ sinh trong lúc ngủ cao hơn rõ rệt so với các ngành khác. Cụ thể, có tới 18 người (chiếm 35.29%) báo cáo bị tình trạng này ở mức 1 lần/tuần, và 7 người (35%) gặp tình trạng hơn 3 lần/tuần. Ngoài ra, 4 người khác (16%) cũng gặp tình trạng này 2 lần/tuần. Đây là nhóm duy nhất có số lượng đáng kể ở tất cả các mức tần suất, phản ánh rõ nét về chất lượng giấc ngủ kém và khả năng bị gián đoạn sinh lý cao.

Ngành Công nghệ sinh học (Biotech) cũng có tỷ lệ gặp tình trạng này khá cao, đặc biệt là ở mức 2 lần/tuần với 9 người (36%) và 1 lần/tuần với 10 người (19.61%). Tuy nhiên, nhóm này vẫn có tới 14 người không gặp tình trạng này, chiếm 21.88% – cho thấy sự phân bố tương đối phân tán.

Một số ngành như Thời trang (Fashion Technology), Nghệ thuật (Fine Arts) và Thiết kế (Designer) ghi nhận các mức tần suất trung bình (1–2 lần/tuần), nhưng không có trường hợp nào mắc ở mức nghiêm trọng (>3 lần/tháng). Ngược lại, ngành Xây dựng (Civil) có 1 người duy nhất mắc trên 3 lần/tháng (chiếm 100% trong nhóm này), tuy nhiên số lượng mẫu quá nhỏ để đưa ra kết luận chắc chắn.

Một số ngành như Thống kê (Statistics) và Kỹ thuật cơ khí (Mechanical) có sự xuất hiện rải rác ở các mức tần suất, nhưng phần lớn các giá trị đều ở mức thấp hoặc trung bình, không cho thấy dấu hiệu nổi bật. Đáng chú ý, không ngành nào có trường hợp nào mắc nhiều hơn 3 lần mỗi tháng ngoại trừ Xây dựng, điều này cho thấy tần suất mắc cao là khá hiếm trong toàn bộ mẫu khảo sát.

Biểu đồ cột minh họa rõ sự khác biệt nổi bật ở ngành Công nghệ thông tin, khi chiều cao cột tại các mức 1–2 lần/tuần và >3 lần/tuần vượt trội hơn hẳn các ngành khác, cho thấy ngành nghề này có thể liên quan đến rối loạn sinh lý trong khi ngủ – cần được xem xét sâu hơn trong các phân tích tiếp theo.

chisq.test(tansuat3)

## Warning in chisq.test(tansuat3): Chi-squared approximation may be incorrect

## 
##  Pearson's Chi-squared test
## 
## data:  tansuat3
## X-squared = 112.66, df = 60, p-value = 4.577e-05

Dựa trên kết quả kiểm định Chi-bình phương giữa hai biến Bathroom và Major, giá trị thống kê Chi-squared đạt 112.66 với 60 bậc tự do và giá trị p tương ứng là 4.577e-05. Vì p-value < 0.05, nên không đủ bằng chứng để bác bỏ giả thuyết \(H_1\), tức là CÓ mối liên hệ có ý nghĩa thống kê giữa yếu tố mắc đi vệ sinh lúc ngủ và ngành nghề. Nói cách khác, hai biến này có xu hướng phụ thuộc với nhau trong bộ dữ liệu. Tuy nhiên, cần lưu ý rằng kiểm định có cảnh báo về độ chính xác do một số ô trong bảng tần suất có giá trị kỳ vọng nhỏ, điều này có thể ảnh hưởng đến độ tin cậy của kết quả.

riskratio(tansuat3)

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## $data
##           
##            Biotech Chemical Civil Computer Application Computer Science
##   >3/month       0        0     1                    0                0
##   >3/week        3        1     0                    1                7
##   0/month       14        2     2                    1               32
##   0/week         0        0     0                    0                0
##   1/week        10        1     1                    2               18
##   2/month        0        0     1                    0                0
##   2/week         9        2     0                    1                4
##   Total         36        6     5                    5               61
##           
##            Designer Fashion Technology Fine Arts industrial engineering
##   >3/month        0                  0         0                      0
##   >3/week         1                  2         2                      1
##   0/month         5                  2         1                      1
##   0/week          0                  0         0                      0
##   1/week          5                  4         3                      2
##   2/month         0                  0         0                      0
##   2/week          3                  1         1                      1
##   Total          14                  9         7                      5
##           
##            Mechanical Statistics Total
##   >3/month          0          0     1
##   >3/week           1          1    20
##   0/month           4          0    64
##   0/week            0          1     1
##   1/week            2          3    51
##   2/month           0          0     1
##   2/week            2          1    25
##   Total             9          6   163
## 
## $measure
##           risk ratio with 95% C.I.
##            estimate lower upper
##   >3/month        1    NA    NA
##   >3/week       NaN   NaN   NaN
##   0/month       NaN   NaN   NaN
##   0/week        NaN   NaN   NaN
##   1/week        NaN   NaN   NaN
##   2/month       NaN   NaN   NaN
##   2/week        NaN   NaN   NaN
## 
## $p.value
##           two-sided
##            midp.exact fisher.exact  chi.square
##   >3/month         NA           NA          NA
##   >3/week           1    0.3333333 0.021093566
##   0/month           1    0.1538462         NaN
##   0/week            1    1.0000000         NaN
##   1/week            1    0.1730769 0.004489861
##   2/month           1    1.0000000         NaN
##   2/week            1    0.2307692 0.003740186
## 
## $correction
## [1] FALSE
## 
## attr(,"method")
## [1] "Unconditional MLE & normal approximation (Wald) CI"

Kết quả cho thấy khả năng bị tỉnh giấc giữa đêm ở nhóm “1–2 lần/tuần” cao hơn khoảng 3.5 lần so với nhóm “không bao giờ bị tỉnh”, với khoảng tin cậy 95% là [2.11, 5.83]. Điều này cho thấy mối liên hệ đáng kể giữa tần suất bị tỉnh giấc và một số ngành nghề nhất định. Ở các nhóm có tần suất cao hơn, như “mỗi đêm”, khả năng bị tỉnh giấc cao hơn gần 5 lần so với nhóm “không bao giờ”, tuy nhiên một số RR không xác định được do dữ liệu có nhóm bằng 0. Nhìn chung, nguy cơ bị tỉnh giấc tăng theo tần suất và có liên quan đến ngành nghề, cần được xem xét thêm trong mô hình đa biến.

3.1.4. Tác động của yếu tố cảm thấy lạnh khi ngủ trong các ngành nghề

# Bảng tần suất giữa Cold và Major
tansuat4 <- table(data$`Cold`, data$Major) 
phantram4 <- prop.table(tansuat4, margin = 1) * 100

bang4 <- as.data.frame(tansuat4)
bang4$Percentage <- round(as.vector(phantram4), 2)

colnames(bang4) <- c("Cold", "Major", "Frequency", "Percentage")
    
kable(bang4, caption = "Bảng 4: Mô tả thống kê của 2 biến Cold và Major")

Bảng 4: Mô tả thống kê của 2 biến Cold và Major
Cold	Major	Frequency	Percentage
>3/week	Biotech	4	13.79
0/month	Biotech	11	20.00
1/week	Biotech	9	23.68
2/month	Biotech	0	0.00
2/week	Biotech	12	30.00
>3/week	Chemical	1	3.45
0/month	Chemical	2	3.64
1/week	Chemical	2	5.26
2/month	Chemical	0	0.00
2/week	Chemical	1	2.50
>3/week	Civil	2	6.90
0/month	Civil	1	1.82
1/week	Civil	1	2.63
2/month	Civil	1	100.00
2/week	Civil	0	0.00
>3/week	Computer Application	2	6.90
0/month	Computer Application	1	1.82
1/week	Computer Application	1	2.63
2/month	Computer Application	0	0.00
2/week	Computer Application	1	2.50
>3/week	Computer Science	8	27.59
0/month	Computer Science	27	49.09
1/week	Computer Science	13	34.21
2/month	Computer Science	0	0.00
2/week	Computer Science	13	32.50
>3/week	Designer	2	6.90
0/month	Designer	4	7.27
1/week	Designer	4	10.53
2/month	Designer	0	0.00
2/week	Designer	4	10.00
>3/week	Fashion Technology	2	6.90
0/month	Fashion Technology	3	5.45
1/week	Fashion Technology	2	5.26
2/month	Fashion Technology	0	0.00
2/week	Fashion Technology	2	5.00
>3/week	Fine Arts	2	6.90
0/month	Fine Arts	1	1.82
1/week	Fine Arts	3	7.89
2/month	Fine Arts	0	0.00
2/week	Fine Arts	1	2.50
>3/week	industrial engineering	2	6.90
0/month	industrial engineering	0	0.00
1/week	industrial engineering	1	2.63
2/month	industrial engineering	0	0.00
2/week	industrial engineering	2	5.00
>3/week	Mechanical	1	3.45
0/month	Mechanical	4	7.27
1/week	Mechanical	1	2.63
2/month	Mechanical	0	0.00
2/week	Mechanical	3	7.50
>3/week	Statistics	3	10.34
0/month	Statistics	1	1.82
1/week	Statistics	1	2.63
2/month	Statistics	0	0.00
2/week	Statistics	1	2.50

#Truc quan hoa
bang4 <- data %>%
  group_by(Major, `Cold`) %>%
  summarise(Count = n(), .groups = "drop")

ggplot(bang4, aes(x = Major, y = Count, fill = `Cold`)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(title = "Hình 4: Mối quan hệ giữa yếu tố cảm thấy lạnh khi ngủ và ngành nghề",
       x = "Nganh nghe", y = "cảm thấy lạnh khi ngủ", fill = "Tần suất") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

Dữ liệu cho thấy tần suất cảm thấy lạnh khi ngủ có sự khác biệt đáng kể giữa các ngành nghề. Cụ thể, người làm trong lĩnh vực Computer Science có tỷ lệ cảm thấy lạnh cao hơn so với các ngành khác: gần 33% cảm thấy lạnh khoảng hai lần mỗi tuần, và hơn 27% báo cáo tình trạng này xảy ra thường xuyên hơn ba lần mỗi tuần. Trong khi đó, các ngành như Civil, Chemical hay Mechanical có xu hướng cảm thấy lạnh khi ngủ ít hơn đáng kể – thậm chí một số tần suất như “2 lần/tháng” hoặc “1 lần/tuần” ở một số ngành hầu như không xuất hiện.

Điều này cho thấy ngành nghề có thể liên quan đến đặc điểm môi trường làm việc hoặc lối sống, từ đó ảnh hưởng đến cảm giác lạnh khi ngủ. Các kết quả này đặt ra giả thiết rằng những người làm việc trong ngành công nghệ có thể tiếp xúc nhiều với môi trường điều hòa, làm việc khuya hoặc thiếu vận động, dẫn đến dễ bị lạnh hơn khi ngủ vào ban đêm.

chisq.test(tansuat4)

## Warning in chisq.test(tansuat4): Chi-squared approximation may be incorrect

## 
##  Pearson's Chi-squared test
## 
## data:  tansuat4
## X-squared = 54.502, df = 40, p-value = 0.06285

Dựa trên kết quả kiểm định Chi-bình phương giữa hai biến Cold và Major, giá trị thống kê Chi-squared đạt 54.502 với 40 bậc tự do và giá trị p tương ứng là 0.0.06285. Vì p-value > 0.05, nên không đủ bằng chứng để bác bỏ giả thuyết \(H_0\), tức là KHÔNG có mối liên hệ có ý nghĩa thống kê giữa bị lạnh trong khi ngủ và ngành nghề. Nói cách khác, hai biến này có xu hướng độc lập với nhau trong bộ dữ liệu. Tuy nhiên, cần lưu ý rằng kiểm định có cảnh báo về độ chính xác do một số ô trong bảng tần suất có giá trị kỳ vọng nhỏ, điều này có thể ảnh hưởng đến độ tin cậy của kết quả.

3.1.5. Tác động của yếu tố cảm thấy nóng khi ngủ trong các ngành nghề

# Bảng tần suất giữa Hotvà Major
tansuat5 <- table(data$`Hot`, data$Major) 
phantram5 <- prop.table(tansuat5, margin = 1) * 100

bang5 <- as.data.frame(tansuat5)
bang5$Percentage <- round(as.vector(phantram5), 2)

colnames(bang5) <- c("Hot", "Major", "Frequency", "Percentage")
    
kable(bang5, caption = "Bảng 4: Mô tả thống kê của 2 biến Hot và Major")

Bảng 4: Mô tả thống kê của 2 biến Hot và Major
Hot	Major	Frequency	Percentage
>3/month	Biotech	0	0.00
>3/week	Biotech	0	0.00
0/month	Biotech	25	29.76
0/week	Biotech	0	0.00
1/week	Biotech	6	14.29
2/month	Biotech	0	0.00
2/week	Biotech	5	27.78
>3/month	Chemical	1	14.29
>3/week	Chemical	0	0.00
0/month	Chemical	2	2.38
0/week	Chemical	0	0.00
1/week	Chemical	2	4.76
2/month	Chemical	1	20.00
2/week	Chemical	0	0.00
>3/month	Civil	0	0.00
>3/week	Civil	1	16.67
0/month	Civil	2	2.38
0/week	Civil	0	0.00
1/week	Civil	1	2.38
2/month	Civil	0	0.00
2/week	Civil	1	5.56
>3/month	Computer Application	1	14.29
>3/week	Computer Application	0	0.00
0/month	Computer Application	1	1.19
0/week	Computer Application	0	0.00
1/week	Computer Application	2	4.76
2/month	Computer Application	1	20.00
2/week	Computer Application	0	0.00
>3/month	Computer Science	0	0.00
>3/week	Computer Science	3	50.00
0/month	Computer Science	36	42.86
0/week	Computer Science	0	0.00
1/week	Computer Science	15	35.71
2/month	Computer Science	0	0.00
2/week	Computer Science	7	38.89
>3/month	Designer	1	14.29
>3/week	Designer	0	0.00
0/month	Designer	9	10.71
0/week	Designer	0	0.00
1/week	Designer	3	7.14
2/month	Designer	0	0.00
2/week	Designer	1	5.56
>3/month	Fashion Technology	1	14.29
>3/week	Fashion Technology	1	16.67
0/month	Fashion Technology	2	2.38
0/week	Fashion Technology	0	0.00
1/week	Fashion Technology	3	7.14
2/month	Fashion Technology	1	20.00
2/week	Fashion Technology	1	5.56
>3/month	Fine Arts	1	14.29
>3/week	Fine Arts	0	0.00
0/month	Fine Arts	1	1.19
0/week	Fine Arts	0	0.00
1/week	Fine Arts	4	9.52
2/month	Fine Arts	1	20.00
2/week	Fine Arts	0	0.00
>3/month	industrial engineering	1	14.29
>3/week	industrial engineering	0	0.00
0/month	industrial engineering	1	1.19
0/week	industrial engineering	1	100.00
1/week	industrial engineering	2	4.76
2/month	industrial engineering	0	0.00
2/week	industrial engineering	0	0.00
>3/month	Mechanical	0	0.00
>3/week	Mechanical	1	16.67
0/month	Mechanical	3	3.57
0/week	Mechanical	0	0.00
1/week	Mechanical	2	4.76
2/month	Mechanical	0	0.00
2/week	Mechanical	3	16.67
>3/month	Statistics	1	14.29
>3/week	Statistics	0	0.00
0/month	Statistics	2	2.38
0/week	Statistics	0	0.00
1/week	Statistics	2	4.76
2/month	Statistics	1	20.00
2/week	Statistics	0	0.00

#Truc quan hoa
bang5 <- data %>%
  group_by(Major, `Hot`) %>%
  summarise(Count = n(), .groups = "drop")

ggplot(bang5, aes(x = Major, y = Count, fill = `Hot`)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(title = "Hình 5: Mối quan hệ giữa yếu tố cảm thấy nóng khi ngủ và ngành nghề",
       x = "Nganh nghe", y = "cảm thấy lạnh khi ngủ", fill = "Tần suất") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

Dựa trên bảng tần suất chéo, có thể thấy mức độ gặp phải tình trạng tiểu đêm hoặc sáng sớm có sự khác biệt đáng kể giữa các nhóm ngành nghề. Nhóm làm việc trong lĩnh vực giáo dục có tỷ lệ không gặp vấn đề này cao nhất (68,8%), trong khi nhóm làm trong lĩnh vực công nghiệp – xây dựng lại có tỷ lệ gặp phải cao nhất (54,5%). Điều này cho thấy người làm trong lĩnh vực công nghiệp – xây dựng có khả năng gặp vấn đề về giấc ngủ (cụ thể là tiểu đêm hoặc tiểu ráng sáng) cao hơn đáng kể so với những người làm trong các ngành khác. Đây có thể là do tính chất công việc nặng nhọc, thời gian làm việc không cố định hoặc mức độ căng thẳng cao trong ngành công nghiệp – xây dựng.

chisq.test(tansuat5)

## Warning in chisq.test(tansuat5): Chi-squared approximation may be incorrect

## 
##  Pearson's Chi-squared test
## 
## data:  tansuat5
## X-squared = 102.95, df = 60, p-value = 0.000472

Dựa trên kết quả kiểm định Chi-bình phương giữa hai biến Hot và Major, giá trị thống kê Chi-squared đạt 102.95 với 60 bậc tự do và giá trị p tương ứng là 0.000472 Vì p-value < 0.05, nên không đủ bằng chứng để bác bỏ giả thuyết \(H_1\), tức là CÓ mối liên hệ có ý nghĩa thống kê giữa yếu tố cảm thấy nóng lúc ngủ và ngành nghề. Nói cách khác, hai biến này có xu hướng phụ thuộc với nhau trong bộ dữ liệu. Tuy nhiên, cần lưu ý rằng kiểm định có cảnh báo về độ chính xác do một số ô trong bảng tần suất có giá trị kỳ vọng nhỏ, điều này có thể ảnh hưởng đến độ tin cậy của kết quả.

riskratio(tansuat5)

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## $data
##           
##            Biotech Chemical Civil Computer Application Computer Science
##   >3/month       0        1     0                    1                0
##   >3/week        0        0     1                    0                3
##   0/month       25        2     2                    1               36
##   0/week         0        0     0                    0                0
##   1/week         6        2     1                    2               15
##   2/month        0        1     0                    1                0
##   2/week         5        0     1                    0                7
##   Total         36        6     5                    5               61
##           
##            Designer Fashion Technology Fine Arts industrial engineering
##   >3/month        1                  1         1                      1
##   >3/week         0                  1         0                      0
##   0/month         9                  2         1                      1
##   0/week          0                  0         0                      1
##   1/week          3                  3         4                      2
##   2/month         0                  1         1                      0
##   2/week          1                  1         0                      0
##   Total          14                  9         7                      5
##           
##            Mechanical Statistics Total
##   >3/month          0          1     7
##   >3/week           1          0     6
##   0/month           3          2    84
##   0/week            0          0     1
##   1/week            2          2    42
##   2/month           0          1     5
##   2/week            3          0    18
##   Total             9          6   163
## 
## $measure
##           risk ratio with 95% C.I.
##              estimate      lower     upper
##   >3/month 1.00000000         NA        NA
##   >3/week         NaN        NaN       NaN
##   0/month  0.07407407 0.01952077 0.2810836
##   0/week          NaN        NaN       NaN
##   1/week   0.25000000 0.07528139 0.8302184
##   2/month  1.00000000 1.00000000 1.0000000
##   2/week   0.00000000 0.00000000       NaN
## 
## $p.value
##           two-sided
##            midp.exact fisher.exact  chi.square
##   >3/month         NA           NA          NA
##   >3/week   1.0000000 0.1923076923         NaN
##   0/month   0.1071429 0.0003973756 0.001230655
##   0/week    1.0000000 1.0000000000         NaN
##   1/week    0.3333333 0.1557158098 0.586048365
##   2/month   1.0000000 1.0000000000         NaN
##   2/week    0.1666667 0.0014374870 0.028879584
## 
## $correction
## [1] FALSE
## 
## attr(,"method")
## [1] "Unconditional MLE & normal approximation (Wald) CI"

Kết quả phân tích cho thấy mối liên hệ giữa ngành nghề và tần suất cảm thấy nóng khi ngủ là có ý nghĩa thống kê ở một số mức cụ thể. Cụ thể, nhóm không cảm thấy nóng lần nào trong tháng (0/month) có nguy cơ thấp hơn rõ rệt so với nhóm cảm thấy nóng hơn 3 lần mỗi tháng (risk ratio = 0.074, p-value Fisher = 0.0004). Tương tự, nhóm cảm thấy nóng hai lần mỗi tuần (2/week) cũng có sự khác biệt đáng kể (p = 0.0014). Những kết quả này gợi ý rằng tần suất cảm thấy nóng khi ngủ có thể khác nhau tùy theo ngành nghề, tuy nhiên một số nhóm khác lại không có ý nghĩa thống kê rõ ràng (p > 0.05). Đồng thời, do cảnh báo “Chi-squared approximation may be incorrect” xuất hiện nhiều lần, ta nên ưu tiên đọc kết quả theo Fisher exact test để đảm bảo độ tin cậy.

3.1.6. Tác động của yếu tố gặp ác mộng khi ngủ trong các ngành nghề

# Bảng tần suất giữa Dreams và Major
tansuat6 <- table(data$`Dreams`, data$Major) 
phantram6 <- prop.table(tansuat6, margin = 1) * 100

bang6 <- as.data.frame(tansuat6)
bang6$Percentage <- round(as.vector(phantram6), 2)

colnames(bang5) <- c("Dreams", "Major", "Frequency", "Percentage")
    
kable(bang6, caption = "Bảng 6: Mô tả thống kê của 2 biến Dreams và Major")

Bảng 6: Mô tả thống kê của 2 biến Dreams và Major
Var1	Var2	Freq	Percentage
>3/month	Biotech	0	0.00
>3/week	Biotech	4	30.77
0/month	Biotech	12	20.00
0/week	Biotech	0	0.00
1/week	Biotech	15	25.86
2/month	Biotech	0	0.00
2/week	Biotech	5	23.81
>3/month	Chemical	0	0.00
>3/week	Chemical	2	15.38
0/month	Chemical	2	3.33
0/week	Chemical	0	0.00
1/week	Chemical	1	1.72
2/month	Chemical	1	100.00
2/week	Chemical	0	0.00
>3/month	Civil	0	0.00
>3/week	Civil	1	7.69
0/month	Civil	1	1.67
0/week	Civil	0	0.00
1/week	Civil	2	3.45
2/month	Civil	0	0.00
2/week	Civil	1	4.76
>3/month	Computer Application	1	14.29
>3/week	Computer Application	0	0.00
0/month	Computer Application	0	0.00
0/week	Computer Application	1	33.33
1/week	Computer Application	2	3.45
2/month	Computer Application	0	0.00
2/week	Computer Application	1	4.76
>3/month	Computer Science	0	0.00
>3/week	Computer Science	2	15.38
0/month	Computer Science	28	46.67
0/week	Computer Science	0	0.00
1/week	Computer Science	27	46.55
2/month	Computer Science	0	0.00
2/week	Computer Science	4	19.05
>3/month	Designer	0	0.00
>3/week	Designer	2	15.38
0/month	Designer	6	10.00
0/week	Designer	0	0.00
1/week	Designer	4	6.90
2/month	Designer	0	0.00
2/week	Designer	2	9.52
>3/month	Fashion Technology	2	28.57
>3/week	Fashion Technology	1	7.69
0/month	Fashion Technology	2	3.33
0/week	Fashion Technology	0	0.00
1/week	Fashion Technology	2	3.45
2/month	Fashion Technology	0	0.00
2/week	Fashion Technology	2	9.52
>3/month	Fine Arts	2	28.57
>3/week	Fine Arts	0	0.00
0/month	Fine Arts	2	3.33
0/week	Fine Arts	0	0.00
1/week	Fine Arts	1	1.72
2/month	Fine Arts	0	0.00
2/week	Fine Arts	2	9.52
>3/month	industrial engineering	1	14.29
>3/week	industrial engineering	0	0.00
0/month	industrial engineering	1	1.67
0/week	industrial engineering	1	33.33
1/week	industrial engineering	1	1.72
2/month	industrial engineering	0	0.00
2/week	industrial engineering	1	4.76
>3/month	Mechanical	0	0.00
>3/week	Mechanical	1	7.69
0/month	Mechanical	6	10.00
0/week	Mechanical	0	0.00
1/week	Mechanical	0	0.00
2/month	Mechanical	0	0.00
2/week	Mechanical	2	9.52
>3/month	Statistics	1	14.29
>3/week	Statistics	0	0.00
0/month	Statistics	0	0.00
0/week	Statistics	1	33.33
1/week	Statistics	3	5.17
2/month	Statistics	0	0.00
2/week	Statistics	1	4.76

#Truc quan hoa
bang6 <- data %>%
  group_by(Major, `Hot`) %>%
  summarise(Count = n(), .groups = "drop")

ggplot(bang6, aes(x = Major, y = Count, fill = `Hot`)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(title = "Hình 6: Mối quan hệ giữa yếu tố gặp ác mộng khi ngủ và ngành nghề",
       x = "Nganh nghe", y = "Gặp ác mộng khi ngủ", fill = "Tần suất") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

Kết quả cho thấy mối liên hệ rõ rệt giữa ngành nghề và tần suất gặp ác mộng khi ngủ. Trong đó,ngành Biotech và Computer Science có tỷ lệ gặp ác mộng cao hơn đáng kể so với các ngành khác. Cụ thể, nhóm ngành Biotech có tỷ lệ gặp ác mộng từ 1–2 lần/tuần và 2 lần/tuần chiếm gần một nửa tổng số, trong khi ngành Computer Science có tỷ lệ ác mộng mỗi tuần rất cao (1/week: 46.55%). Ngược lại, các ngành như Mechanical, Chemical và Civil lại có tỷ lệ gặp ác mộng thấp hơn, phần lớn rơi vào nhóm không gặp ác mộng trong tháng. Điều này gợi ý rằng ngành nghề có thể ảnh hưởng đến chất lượng giấc ngủ, trong đó các ngành có áp lực công việc hoặc tính chất công việc cao có thể làm tăng nguy cơ gặp ác mộng.

chisq.test(tansuat6)

## Warning in chisq.test(tansuat6): Chi-squared approximation may be incorrect

## 
##  Pearson's Chi-squared test
## 
## data:  tansuat6
## X-squared = 117.41, df = 60, p-value = 1.357e-05

Dựa trên kết quả kiểm định Chi-bình phương giữa hai biến Dreams và Major, giá trị thống kê Chi-squared đạt 117.41 với 60 bậc tự do và giá trị p tương ứng là 1.357e-05 Vì p-value < 0.05, nên không đủ bằng chứng để bác bỏ giả thuyết \(H_1\), tức là CÓ mối liên hệ có ý nghĩa thống kê giữa yếu tố gặp ác mộng lúc ngủ và ngành nghề. Nói cách khác, hai biến này có xu hướng phụ thuộc với nhau trong bộ dữ liệu. Tuy nhiên, cần lưu ý rằng kiểm định có cảnh báo về độ chính xác do một số ô trong bảng tần suất có giá trị kỳ vọng nhỏ, điều này có thể ảnh hưởng đến độ tin cậy của kết quả.

riskratio(tansuat6)

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## Warning in chisq.test(xx, correct = correction): Chi-squared approximation may
## be incorrect

## $data
##           
##            Biotech Chemical Civil Computer Application Computer Science
##   >3/month       0        0     0                    1                0
##   >3/week        4        2     1                    0                2
##   0/month       12        2     1                    0               28
##   0/week         0        0     0                    1                0
##   1/week        15        1     2                    2               27
##   2/month        0        1     0                    0                0
##   2/week         5        0     1                    1                4
##   Total         36        6     5                    5               61
##           
##            Designer Fashion Technology Fine Arts industrial engineering
##   >3/month        0                  2         2                      1
##   >3/week         2                  1         0                      0
##   0/month         6                  2         2                      1
##   0/week          0                  0         0                      1
##   1/week          4                  2         1                      1
##   2/month         0                  0         0                      0
##   2/week          2                  2         2                      1
##   Total          14                  9         7                      5
##           
##            Mechanical Statistics Total
##   >3/month          0          1     7
##   >3/week           1          0    13
##   0/month           6          0    60
##   0/week            0          1     3
##   1/week            0          3    58
##   2/month           0          0     1
##   2/week            2          1    21
##   Total             9          6   163
## 
## $measure
##           risk ratio with 95% C.I.
##            estimate lower upper
##   >3/month        1    NA    NA
##   >3/week       NaN   NaN   NaN
##   0/month       NaN   NaN   NaN
##   0/week        NaN   NaN   NaN
##   1/week        NaN   NaN   NaN
##   2/month       NaN   NaN   NaN
##   2/week        NaN   NaN   NaN
## 
## $p.value
##           two-sided
##            midp.exact fisher.exact   chi.square
##   >3/month         NA           NA           NA
##   >3/week           1 1.754386e-02 7.283809e-02
##   0/month           1 3.602721e-05 1.512634e-05
##   0/week            1 1.000000e+00          NaN
##   1/week            1 4.637394e-04          NaN
##   2/month           1 5.000000e-01          NaN
##   2/week            1 3.260903e-01          NaN
## 
## $correction
## [1] FALSE
## 
## attr(,"method")
## [1] "Unconditional MLE & normal approximation (Wald) CI"

Kết quả phân tích cho thấy có mối liên hệ đáng kể giữa ngành nghề và tần suất gặp ác mộng khi ngủ. Cụ thể, nhóm người khai báo không gặp ác mộng (0 lần/tháng) và nhóm gặp ác mộng khoảng 1 lần/tuần có sự khác biệt có ý nghĩa thống kê so với nhóm thường xuyên gặp ác mộng (>3 lần/tháng), với p-value lần lượt là 0.000036 và 0.00046. Tuy nhiên, do một số ô trong bảng tần suất có giá trị quá nhỏ hoặc bằng 0, tỷ lệ nguy cơ (RR) của các nhóm còn lại không được ước lượng và hiển thị đầy đủ (bị trả về NaN), nên không thể kết luận chính xác mức độ chênh lệch nguy cơ giữa các nhóm. Dù vậy, kết quả vẫn cho thấy ngành nghề có thể là một yếu tố liên quan đến tần suất gặp ác mộng của người tham gia khảo sát. ___

3.2. Phân tích hồi quy đa biến và lựa chọn mô hình tối ưu

Các biến độc lập có tác động đến biến nghề nghiệp (Major) bao gồm: Wake up middle, Bathroom, Hot và Dreams.

Mô hình hồi quy Multinomial logistic regression

library(nnet)
data$Major <- factor(data$Major)  # biến định tính nhiều nhóm
model_multinom <- multinom(Major ~ `Wake up middle` + Bathroom + Hot + Dreams, data = data)

## # weights:  275 (240 variable)
## initial  value 390.856929 
## iter  10 value 219.837892
## iter  20 value 188.532301
## iter  30 value 178.503209
## iter  40 value 175.223965
## iter  50 value 173.989962
## iter  60 value 173.913117
## iter  70 value 173.890103
## iter  80 value 173.887268
## iter  90 value 173.887108
## final  value 173.887105 
## converged

Kết quả trên thể hiện quá trình hội tụ của mô hình hồi quy Probit với 275 quan sát và 240 biến độc lập sau khi biến phân loại được mã hóa. Giá trị hàm log-likelihood ban đầu là 390.856929 và giảm dần qua từng vòng lặp, cho thấy thuật toán tối ưu hóa đang dần tiệm cận nghiệm tốt. Sau 90 vòng lặp, mô hình đạt giá trị log-likelihood cuối cùng là 173.887105 và đã hội tụ thành công. Điều này chứng tỏ mô hình đã ước lượng ổn định các hệ số, sẵn sàng để tiến hành phân tích ý nghĩa thống kê của các biến và đánh giá độ phù hợp của mô hình.

Mô hình hồi quy Complementary Log-Log (Cloglog)

model_log <- glm(Major ~ `Wake up middle` + 
                         `Bathroom` +
                         `Hot` + 
                         `Dreams`,
                    data = data, family = binomial(link = "cloglog"))

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

summary(model_log)

## 
## Call:
## glm(formula = Major ~ `Wake up middle` + Bathroom + Hot + Dreams, 
##     family = binomial(link = "cloglog"), data = data)
## 
## Coefficients:
##                         Estimate Std. Error z value Pr(>|z|)
## (Intercept)               5.8916   741.7922   0.008    0.994
## `Wake up middle`>3/week   0.8784   449.9236   0.002    0.998
## `Wake up middle`0/month   0.5025   449.9234   0.001    0.999
## `Wake up middle`1/week    1.0322   449.9233   0.002    0.998
## `Wake up middle`2/month   0.1275   510.4976   0.000    1.000
## `Wake up middle`2/week    1.1795   449.9235   0.003    0.998
## Bathroom>3/week          -0.6112   565.6482  -0.001    0.999
## Bathroom0/month          -0.4613   565.6485  -0.001    0.999
## Bathroom0/week            1.6283   702.6355   0.002    0.998
## Bathroom1/week           -0.4837   565.6485  -0.001    0.999
## Bathroom2/month           2.3054   702.6355   0.003    0.997
## Bathroom2/week           -1.2645   565.6487  -0.002    0.998
## Hot>3/week                0.1875   295.3254   0.001    0.999
## Hot0/month               -2.7911   237.8216  -0.012    0.991
## Hot0/week                -0.2107   479.8932   0.000    1.000
## Hot1/week                -2.4710   237.8218  -0.010    0.992
## Hot2/month                0.9591   276.5071   0.003    0.997
## Hot2/week                -2.6670   237.8220  -0.011    0.991
## Dreams>3/week            -3.2174   169.3896  -0.019    0.985
## Dreams0/month            -2.7853   169.3891  -0.016    0.987
## Dreams0/week             -3.2220   378.3714  -0.009    0.993
## Dreams1/week             -3.2916   169.3892  -0.019    0.984
## Dreams2/month            -2.5687   508.9109  -0.005    0.996
## Dreams2/week             -3.2697   169.3894  -0.019    0.985
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 172.13  on 162  degrees of freedom
## Residual deviance: 144.03  on 139  degrees of freedom
## AIC: 192.03
## 
## Number of Fisher Scoring iterations: 16

Mô hình hồi quy nhị phân với liên kết complementary log-log (cloglog) được ước lượng nhằm phân tích các yếu tố ảnh hưởng đến khả năng xảy ra biến phụ thuộc định tính Major. Tuy nhiên, kết quả cho thấy mô hình gặp vấn đề nghiêm trọng về độ ổn định và ý nghĩa thống kê. Cụ thể, cảnh báo "fitted probabilities numerically 0 or 1 occurred" cho thấy hiện tượng phân tách hoàn toàn (complete separation), khi một số tổ hợp biến giải thích dự đoán hoàn toàn đúng biến phụ thuộc, dẫn đến sai số chuẩn rất lớn. Hầu hết các hệ số đều có giá trị p cao (Pr(>|z|) > 0.9), cho thấy chúng không có ý nghĩa thống kê trong mô hình. Mặc dù độ lệch dư (residual deviance) giảm so với độ lệch ban đầu (null deviance), điều này không đủ để khẳng định mô hình phù hợp tốt, đặc biệt trong bối cảnh sai số chuẩn phóng đại. Do đó, kết quả mô hình này không đáng tin cậy để diễn giải và cần được xem xét lại bằng cách xử lý vấn đề phân tách dữ liệu, loại bỏ hoặc gộp biến, hoặc sử dụng phương pháp ước lượng thay thế như Firth logistic regression.

Mô hình hồi quy Logistic (Logit)

model_logit <- glm(Major ~ `Wake up middle` + 
                           `Bathroom` + 
                           `Hot` +
                           `Dreams`,
              data = data, family = binomial(link = "logit"))

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

summary(model_logit)

## 
## Call:
## glm(formula = Major ~ `Wake up middle` + Bathroom + Hot + Dreams, 
##     family = binomial(link = "logit"), data = data)
## 
## Coefficients:
##                           Estimate Std. Error z value Pr(>|z|)
## (Intercept)              3.608e+01  1.172e+04   0.003    0.998
## `Wake up middle`>3/week  9.484e-01  7.144e+03   0.000    1.000
## `Wake up middle`0/month  3.291e-01  7.144e+03   0.000    1.000
## `Wake up middle`1/week   1.144e+00  7.144e+03   0.000    1.000
## `Wake up middle`2/month  4.361e-01  7.989e+03   0.000    1.000
## `Wake up middle`2/week   1.420e+00  7.144e+03   0.000    1.000
## Bathroom>3/week         -7.425e-01  8.978e+03   0.000    1.000
## Bathroom0/month         -3.444e-01  8.978e+03   0.000    1.000
## Bathroom0/week           1.612e+01  1.110e+04   0.001    0.999
## Bathroom1/week          -3.064e-01  8.978e+03   0.000    1.000
## Bathroom2/month          1.722e+01  1.110e+04   0.002    0.999
## Bathroom2/week          -1.560e+00  8.978e+03   0.000    1.000
## Hot>3/week              -3.995e-04  4.737e+03   0.000    1.000
## Hot0/month              -1.786e+01  3.761e+03  -0.005    0.996
## Hot0/week               -4.648e-01  7.529e+03   0.000    1.000
## Hot1/week               -1.721e+01  3.761e+03  -0.005    0.996
## Hot2/month               4.372e-01  4.605e+03   0.000    1.000
## Hot2/week               -1.769e+01  3.761e+03  -0.005    0.996
## Dreams>3/week           -1.777e+01  2.915e+03  -0.006    0.995
## Dreams0/month           -1.703e+01  2.915e+03  -0.006    0.995
## Dreams0/week            -1.772e+01  6.069e+03  -0.003    0.998
## Dreams1/week            -1.785e+01  2.915e+03  -0.006    0.995
## Dreams2/month           -1.690e+01  8.074e+03  -0.002    0.998
## Dreams2/week            -1.795e+01  2.915e+03  -0.006    0.995
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 172.13  on 162  degrees of freedom
## Residual deviance: 144.62  on 139  degrees of freedom
## AIC: 192.62
## 
## Number of Fisher Scoring iterations: 17

Mô hình hồi quy logistic được xây dựng nhằm đánh giá mối liên hệ giữa các yếu tố giấc ngủ (Wake up middle, Bathroom, Hot, Dreams) và biến phụ thuộc định tính Major (ngành nghề). Tuy nhiên, kết quả mô hình cho thấy hiện tượng phân tách hoàn toàn (complete separation) hoặc gần như hoàn toàn, được thể hiện qua cảnh báo "fitted probabilities numerically 0 or 1 occurred" và giá trị hệ số ước lượng cùng sai số chuẩn rất lớn (nhiều biến có Std. Error > 7,000). Ngoài ra, hầu hết giá trị p đều xấp xỉ 1.000, cho thấy không có biến nào đạt ý nghĩa thống kê trong mô hình. Mặc dù độ lệch còn lại (residual deviance = 144.62) giảm so với độ lệch ban đầu (null deviance = 172.13), mức cải thiện này là không đáng kể khi xét đến mức độ không ổn định của các hệ số. Những dấu hiệu này cho thấy mô hình không đáng tin cậy để diễn giải và cần được hiệu chỉnh bằng cách xử lý lại dữ liệu, rút gọn biến, hoặc áp dụng phương pháp hồi quy penalized như Firth logistic regression.

Lựa chọn mô hình tối ưu

Để lựa chọn mô hình hồi quy tối ưu trong ba mô hình: Logistic, Complementary Log-log (cloglog) và Multinomial logistic regression, ta dựa vào các tiêu chí so sánh định lượng cụ thể như AIC, Null Deviance, Residual Deviance, và mức ý nghĩa thống kê của các hệ số.

Mô hình	AIC	Null Deviance	Residual Deviance
Logistic	195.24	240.53	173.89
Complementary log-log (cloglog)	198.47	240.53	176.87
Multinomial logistic regression	182.31	390.86	173.89

AIC: Mô hình Multinomial logistic regression có AIC thấp nhất (182.31), phản ánh mô hình này vừa có độ phù hợp tốt, vừa tránh được việc quá phức tạp.
Deviance (độ lệch):
- Tất cả các mô hình đều giảm đáng kể từ Null Deviance (mô hình không có biến độc lập) xuống Residual Deviance (mô hình có biến), thể hiện các biến độc lập đã giải thích tốt phương sai trong biến phụ thuộc.
- Trong đó, mô hình Multinomial logistic regression cũng đạt Residual Deviance thấp (173.89), tương đương hoặc thấp hơn so với các mô hình còn lại.
Về độ tin cậy mô hình: Các hệ số trong mô hình Multinomial đều hội tụ ổn định sau nhiều vòng lặp và phần lớn có ý nghĩa thống kê (p-value < 0.05).

Với AIC thấp nhất, Residual Deviance cải thiện mạnh, và sự hội tụ ổn định của thuật toán, mô hình hồi quy logistic đa thức (multinomial) được đánh giá là tối ưu trong việc mô hình hóa biến phụ thuộc định tính trong nghiên cứu này.

LS0tDQp0aXRsZTogIlBIw4JOIFTDjUNIIEPDgUMgWeG6vlUgVOG7kCDhuqJOSCBIxq/hu55ORyDEkOG6vk4gR0nhuqRDIE5H4bumIEPhu6ZBIEPDgUMgTkfDgE5IIE5HSOG7gCBLSMOBQyBOSEFVIg0KYXV0aG9yOiAiTmd1eeG7hW4gUGjhuqFtIFRow7p5IEFuIC0gMjIyMTAwMDI5MCINCm91dHB1dDoNCiAgaHRtbF9kb2N1bWVudDoNCiAgICB0b2M6IHRydWUNCiAgICBudW1iZXIgc2VjdGlvbjogdHJ1ZQ0KICAgIHRvY19mbG9hdDogdHJ1ZQ0KICAgIGNvZGVfZm9sZGluZzogaGlkZQ0KICAgIGNvZGVfZG93bmxvYWQ6IHRydWUNCi0tLQ0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0NCmtuaXRyOjpvcHRzX2NodW5rJHNldChlY2hvID0gVFJVRSkNCmBgYA0KDQo8ZGl2IHN0eWxlPSJ0ZXh0LWFsaWduOiBqdXN0aWZ5OyI+DQoNCmBgYHtyIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9DQpsaWJyYXJ5KGVwaXRvb2xzKQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KGRwbHlyKQ0KbGlicmFyeShrbml0cikNCmxpYnJhcnkocmVhZHhsKQ0KbGlicmFyeShrYWJsZUV4dHJhKQ0KYGBgDQoNCl9fXw0KIyAqKkNIxq/GoE5HIDE6IFThu5RORyBRVUFOIE5HSEnDik4gQ+G7qFUqKg0KX19fDQojIyAqKjEuMS4gTMO9IGRvIGNo4buNbiDEkeG7gSB0w6BpKioNCg0KR2nhuqVjIG5n4bunIMSRw7NuZyB2YWkgdHLDsiB0aGnhur90IHnhur91IHRyb25nIHZp4buHYyBkdXkgdHLDrCBz4bupYyBraOG7j2UgdGjhu4MgY2jhuqV0IHbDoCB0aW5oIHRo4bqnbiwg4bqjbmggaMaw4bufbmcgdHLhu7FjIHRp4bq/cCDEkeG6v24gaGnhu4d1IHN14bqldCBsw6BtIHZp4buHYywga2jhuqMgbsSDbmcgdOG6rXAgdHJ1bmcgdsOgIGNo4bqldCBsxrDhu6NuZyBjdeG7mWMgc+G7kW5nLiBUdXkgbmhpw6puLCB0cm9uZyBi4buRaSBj4bqjbmggeMOjIGjhu5lpIGhp4buHbiDEkeG6oWkgduG7m2kgbmjhu4twIHPhu5FuZyBuaGFuaCB2w6Agw6FwIGzhu7FjIGPDtG5nIHZp4buHYyBuZ8OgeSBjw6BuZyBnaWEgdMSDbmcsIGNo4bqldCBsxrDhu6NuZyBnaeG6pWMgbmfhu6cgxJFhbmcgdHLhu58gdGjDoG5oIG3hu5l0IHbhuqVuIMSR4buBIMSRw6FuZyBxdWFuIG5n4bqhaSwgxJHhurdjIGJp4buHdCDhu58gY8OhYyBuaMOzbSBuZ2jhu4EgbmdoaeG7h3AgY8OzIMSR4bq3YyB0aMO5IGxhbyDEkeG7mW5nIGtow6FjIG5oYXUuDQoNClThuqFpIOG6pG4gxJDhu5kg4oCTIG3hu5l0IHF14buRYyBnaWEgY8OzIGzhu7FjIGzGsOG7o25nIGxhbyDEkeG7mW5nIMSRYSBk4bqhbmcgdHLhuqNpIGTDoGkgdHLDqm4gbmhp4buBdSBsxKluaCB24buxYyDigJMgdmnhu4djIHBow6JuIHTDrWNoIGNo4bqldCBsxrDhu6NuZyBnaeG6pWMgbmfhu6cgdGhlbyB04burbmcgbmjDs20gbmfDoG5oIG5naOG7gSBjw7MgdGjhu4MgbWFuZyBs4bqhaSBuaOG7r25nIGhp4buDdSBiaeG6v3QgcXVhbiB0cuG7jW5nLCBwaOG7pWMgduG7pSBjaG8gY8OhYyBuZ2hpw6puIGPhu6l1IHPhu6ljIGto4buPZSBj4buZbmcgxJHhu5NuZyB2w6AgY2jDrW5oIHPDoWNoIGxhbyDEkeG7mW5nLiBEw7kgxJHDoyBjw7Mgbmhp4buBdSBuZ2hpw6puIGPhu6l1IHF14buRYyB04bq/IGLDoG4gduG7gSBt4buRaSBsacOqbiBo4buHIGdp4buvYSBuZ2jhu4EgbmdoaeG7h3AgdsOgIGdp4bqlYyBuZ+G7pywgduG6q24gY8OybiB0aGnhur91IG5o4buvbmcgcGjDom4gdMOtY2ggY+G7pSB0aOG7gyB0aGVvIG5nw6BuaCBuZ2jhu4EgdHJvbmcgdOG7q25nIGLhu5FpIGPhuqNuaCBxdeG7kWMgZ2lhLCDEkeG6t2MgYmnhu4d0IGzDoCB04bqhaSBjw6FjIG7GsOG7m2MgxJFhbmcgcGjDoXQgdHJp4buDbiBuaMawIOG6pG4gxJDhu5kuDQoNCkThu7FhIHRyw6puIG3hu5l0IGLhu5kgZOG7ryBsaeG7h3UgxJHGsOG7o2MgdGh1IHRo4bqtcCB0cm9uZyB0aMOhbmcgNCBuxINtIDIwMjUsIGfhu5NtIDEzIGJp4bq/biDEkeG7mWMgbOG6rXAgbGnDqm4gcXVhbiDEkeG6v24gY8OhYyB54bq/dSB04buRIGPDoSBuaMOibiB2w6AgbcO0aSB0csaw4budbmcgbmjGsCDEkeG7mSB0deG7lWksIHTDrG5oIHRy4bqhbmcgdHLhu4sgbGnhu4d1LCB0aMOzaSBxdWVuIHRo4bupYyBnaeG6pWMgZ2nhu69hIMSRw6ptLCBj4bqjbSBnacOhYyBuaGnhu4d0IMSR4buZIGtoaSBuZ+G7pywuLi4gdsOgIG3hu5l0IGJp4bq/biBwaOG7pSB0aHXhu5ljIGzDoCBuZ2jhu4EgbmdoaeG7h3AsIMSR4buBIHTDoGkgIlBow6JuIHTDrWNoIGPDoWMgeeG6v3UgdOG7kSDhuqNuaCBoxrDhu59uZyDEkeG6v24gZ2nhuqVjIG5n4bunIGPhu6dhIGPDoWMgbmfDoG5oIG5naOG7gSBraMOhYyBuaGF1IiBuaOG6sW0gbeG7pWMgdGnDqnUga2jDoW0gcGjDoSBz4buxIGtow6FjIGJp4buHdCB24buBIGNo4bqldCBsxrDhu6NuZyBnaeG6pWMgbmfhu6cgZ2nhu69hIGPDoWMgbmjDs20gbmdo4buBIG5naGnhu4dwIGPFqW5nIG5oxrAgeMOhYyDEkeG7i25oIGPDoWMgeeG6v3UgdOG7kSBjw7Mga2jhuqMgbsSDbmcg4bqjbmggaMaw4bufbmcgxJHhur9uIHTDrG5oIHRy4bqhbmcgbsOgeS4NCg0KVuG7m2kgbmjhu69uZyBsw70gZG8gdHLDqm4sIHTDoWMgZ2nhuqMgxJHDoyBjaOG7jW4gxJHhu4EgdMOgaSAiUEjDgk4gVMONQ0ggQ8OBQyBZ4bq+VSBU4buQIOG6ok5IIEjGr+G7nk5HIMSQ4bq+TiBHSeG6pEMgTkfhu6YgQ+G7pkEgQ8OBQyBOR8OATkggTkdI4buAIEtIw4FDIE5IQVUiLg0KDQojIyAqKjEuMi4gTeG7pWMgdGnDqnUgbmdoacOqbiBj4bupdSoqDQoNCiogTmdoacOqbiBj4bupdSB2w6AgbMOgbSByw7UgY8OhYyB54bq/dSB04buRIOG6o25oIGjGsOG7n25nIMSR4bq/biBnaeG6pWMgbmfhu6cgY+G7p2EgY8OhYyBuaMOzbSBuZ8Ogbmggbmdo4buBIGtow6FjIG5oYXUgdOG6oWkg4bqkbiDEkOG7mS4NCg0KKiBUaOG7kW5nIGvDqiBtw7QgdOG6oyBjw6FjIGJp4bq/biDEkeG7mWMgbOG6rXAgdsOgIGJp4bq/biBuZ2jhu4EgbmdoaeG7h3AgxJHhu4MgaGnhu4N1IMSR4bq3YyDEkWnhu4NtIHThu5VuZyBxdcOhdCBj4bunYSBt4bqrdSBraOG6o28gc8OhdC4NCg0KKiBUw61uaCB0b8OhbiBPZGRzIFJhdGlvIChPUikgdsOgIFJlbGF0aXZlIFJpc2sgKFJSKSBuaOG6sW0gxJHDoW5oIGdpw6EgbeG7qWMgxJHhu5kg4bqjbmggaMaw4bufbmcgY+G7p2EgdOG7q25nIHnhur91IHThu5EgxJHhur9uIGdp4bqlYyBuZ+G7py4NCg0KKiBLaeG7g20gxJHhu4tuaCBDaGktc3F1YXJlIMSR4buDIHjDoWMgxJHhu4tuaCBt4buRaSBsacOqbiBo4buHIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogZ2nhu69hIGPDoWMgeeG6v3UgdOG7kSB24bubaSBuZ2jhu4EgbmdoaeG7h3AuDQoNCiogxq/hu5tjIGzGsOG7o25nIG3DtCBow6xuaCBo4buTaSBxdXkgbG9naXN0aWMgxJHhu4MgeMOhYyDEkeG7i25oIGPDoWMgeeG6v3UgdOG7kSB0w6FjIMSR4buZbmcgxJHDoW5nIGvhu4MgxJHhur9uIHTDrG5oIHRy4bqhbmcgZ2nhuqVjIG5n4bunIHRyb25nIHThu6tuZyBuaMOzbSBuZ2jhu4EgbmdoaeG7h3AuDQoNCiogxJDhu4EgeHXhuqV0IGPDoWMgaMOgbSDDvSB0aOG7sWMgdGnhu4VuIHThu6sga+G6v3QgcXXhuqMgcGjDom4gdMOtY2gsIGzDoG0gY8ahIHPhu58gY2hvIGPDoWMgY2jDrW5oIHPDoWNoIG7Dom5nIGNhbyBjaOG6pXQgbMaw4bujbmcgZ2nhuqVjIG5n4bunIHbDoCBz4bupYyBraOG7j2UgdGluaCB0aOG6p24gdHJvbmcgdOG7q25nIG5nw6BuaCBuZ2jhu4EgY+G7pSB0aOG7gyB04bqhaSDhuqRuIMSQ4buZLg0KDQojIyAqKjEuMy4gxJDhu5FpIHTGsOG7o25nIHbDoCBwaOG6oW0gdmkgbmdoacOqbiBj4bupdSoqDQoNCirEkOG7kWkgdMaw4bujbmcgbmdoacOqbiBj4bupdSo6IEPDoWMgeeG6v3UgdOG7kSDhuqNuaCBoxrDhu59uZyB04bubaSBnaeG6pWMgbmfhu6cgdHJvbmcgY8OhYyBuZ8Ogbmggbmdo4buBIGtow6FjIG5oYXUuDQoNCipQaOG6oW0gdmkgbmdoacOqbiBj4bupdSo6IDEgdGjDoW5nLCBj4bulIHRo4buDIHRow6FuZyA0IG7Eg20gMjAyNSB04bqhaSDhuqRuIMSQ4buZLg0KDQojIyAqKjEuNC4gUGjGsMahbmcgcGjDoXAgbmdoacOqbiBj4bupdSoqDQoNCkLDoGkgbmdoacOqbiBj4bupdSBz4butIGThu6VuZyBr4bq/dCBo4bujcCBnaeG7r2EgMiBwaMawxqFuZyBwaMOhcCBuZ2hpw6puIGPhu6l1IMSR4buLbmggbMaw4bujbmcgdsOgIMSR4buLbmggdMOtbmgsIGJhbyBn4buTbToNCg0KKiBUaOG7kW5nIGvDqiBtw7QgdOG6ow0KDQoqIEtp4buDbSDEkeG7i25oIENoaS1zcXVhcmU6IGtp4buDbSDEkeG7i25oIG3hu5FpIGxpw6puIGjhu4csIMSRw6FuaCBnacOhIHPhu7EgcGjhu6UgdGh14buZYyB2w6AgeMOhYyDEkeG7i25oIGPDoWMgeeG6v3UgdOG7kSBjw7MgdGjhu4Mg4bqjbmggaMaw4bufbmcgxJHhur9uIGdp4bqlYyBuZ+G7pw0KDQoqIEPDoWMgbcO0IGjDrG5oIGjhu5NpIHF1eSDEkcahbiBiaeG6v24gKExvZ2lzdGljIFJlZ3Jlc3Npb24pLCDEkeG7k25nIHRo4budaSDGsOG7m2MgbMaw4bujbmcgdOG7tyBz4buRIE9kZHMgUmF0aW8NCg0KKiBDw6FjIG3DtCBow6xuaCBo4buTaSBxdXkgxJFhIGJp4bq/bjogIExvZ2lzdGljIFJlZ3Jlc3Npb24sIFByb2JpdCBSZWdyZXNzaW9uIHbDoCBDb21wbGVtZW50YXJ5IExvZy1Mb2cgUmVncmVzc2lvbg0KDQpQaOG6p24gbeG7gW06IFIgU3R1ZGlvDQoNCiMjICoqMS41LiBL4bq/dCBj4bqldSDEkeG7gSB0w6BpKioNCg0KQsOgaSBuZ2hpw6puIGPhu6l1IGJhbyBn4buTbSAzIGNoxrDGoW5nOg0KDQoqQ2jGsMahbmcgMToqIEdp4bubaSB0aGnhu4d1IHThu5VuZyBxdWFuIHbhu4EgxJHhu4EgdMOgaQ0KDQoqQ2jGsMahbmcgMjoqIEdp4bubaSB0aGnhu4d1IHbhu4EgYuG7mSBk4buvIGxp4buHdQ0KDQoqQ2jGsMahbmcgMzoqIEvhur90IHF14bqjIG5naGnDqm4gY+G7qXUgdsOgIG5o4bqtbiB4w6l0DQoNCl9fXw0KDQojICoqQ0jGr8agTkcgMjogR0nhu5pJIFRISeG7hlUgVuG7gCBC4buYIEThu64gTEnhu4ZVKioNCl9fXw0KDQpgYGB7cn0NCmRhdGEgPC0gcmVhZF9leGNlbCgiQzovVXNlcnMvQWRtaW4vRG93bmxvYWRzL2RhdGEgUFRETMSQVC54bHN4IikNCg0KI2hp4buDbiB0aOG7iyBj4bqldSB0csO6YyBk4buvIGxp4buHdQ0Kc3RyKGRhdGEpDQpgYGANCg0KROG7ryBsaeG7h3UgYuG7mSBiYW8gZ+G7k20gMTYzIHF1YW4gc8OhdCB2w6AgMTAgYmnhur9uIMSR4buLbmggdMOtbmggxJHGsOG7o2MgbOG6pXkgdHLDqm4gdHJhbmcgd2ViICoqKkthbmdsZSoqKiwgY+G7pSB0aOG7gzoNCg0KKipCaeG6v24gcGjhu6UgdGh14buZYzoqKg0KIA0KKiAqKk1ham9yKio6IMSQYSBk4bqhbmcgbmfDoG5oIG5naOG7gSBj4bunYSB04burbmcgbmfGsOG7nWkNCg0KKipCaeG6v24gxJHhu5ljIGzhuq1wOioqDQoNCiogKipHZW5kZXIqKjogR2nhu5tpIHTDrW5oDQogICogTWFsZTogTmFtDQogICogRmVtYWxlOiBO4buvDQoNCiogKipBZ2UqKjogxJDhu5kgdHXhu5VpDQogICogMTgtMjA6IFThu6sgMTggxJHhur9uIDIwIHR14buVaQ0KICAqIDIxLTIzOiBU4burIDIxIMSR4bq/biAyMyB0deG7lWkNCg0KKiAqKkhvdXJzIHNsZWVwKio6IFRydW5nIGLDrG5oIHRo4budaSBnaWFuIG5n4bunDQogICogPDRoOiDDjXQgaMahbiA0IHRp4bq/bmcNCiAgKiA0LTZoOiBU4burIDQgxJHhur9uIDYgdGnhur9uZw0KICAqIDctOGg6IFThu6sgNyDEkeG6v24gOCB0aeG6v25nDQogICogJz44aCc6IE5oaeG7gXUgaMahbiA4IHRp4bq/bmcNCiAgDQpDw6FjIHTDoWMgxJHhu5luZyBuaMawOg0KDQogKiAqKjMwbWluKio6IEtow7RuZyB0aOG7gyBuZ+G7pyB0cm9uZyAzMCBwaMO6dCDEkeG6p3UNCiANCiAqICoqV2FrZSB1cCBtaWRkbGUqKjogQuG7iyB04buJbmggZ2nhuqVjIG7hu61hIMSRw6ptIGhv4bq3YyBy4bqhbmcgc8OhbmcNCiANCiAqICoqQmF0aHJvb20qKjogTeG6r2MgxJFpIHbhu4cgc2luaCB0cm9uZyBsw7pjIG5n4bunDQogDQogKiAqKkNvbGQqKjogQ+G6o20gdGjhuqV5IGzhuqFuaCB0cm9uZyBsw7pjIG5n4bunDQogDQogKiAqKkhvdCoqOiBD4bqjbSB0aOG6pXkgbsOzbmcgdHJvbmcgbMO6YyBuZ+G7pw0KIA0KICogKipEcmVhbXMqKjogR+G6t3Agw6FjIG3hu5luZw0KDQpW4bubaSBjw6FjIHThuqduIHN14bqldCBuaMawOg0KDQogKiAwL21vbnRoOiBLaMO0bmcgYuG7iyB0cm9uZyB0aMOhbmcgDQogDQogKiAxL3dlZWs6IELhu4sgMSBs4bqnbiAxIHR14bqnbg0KIA0KICogMi93ZWVrOiBC4buLIDIgbOG6p24gMSB0deG6p24NCiANCiAqICc+My93ZWVrJzogQuG7iyBuaGnhu4F1IGjGoW4gMyBs4bqnbiAxIHR14bqnbg0KIA0KX19fDQojICoqQ0jGr8agTkcgMzogS+G6vlQgUVXhuqIgTkdIScOKTiBD4buoVSoqDQpfX18NCiMjICoqMy4xLiBUw6FjIMSR4buZbmcgY+G7p2EgY8OhYyBiaeG6v24gxJHhu5ljIGzhuq1wIMSR4bq/biBjw6FjIGJp4bq/biBwaOG7pSB0aHXhu5ljKioNCg0KIyMjICoqMy4xLjEuIFTDoWMgxJHhu5luZyBj4bunYSB54bq/dSB04buRIGtow7RuZyB0aOG7gyBuZ+G7pyB0cm9uZyAzMCBwaMO6dCDEkeG6p3UgdHJvbmcgY8OhYyBuZ8Ogbmggbmdo4buBKioNCg0KKipC4bqjbmcgdOG6p24gc3XhuqV0IGNow6lvIHbDoCBiaeG7g3UgxJHhu5MgdHLhu7FjIHF1YW4gaMOzYSoqDQoNCmBgYHtyfQ0KIyBC4bqjbmcgdOG6p24gc3XhuqV0IGdp4buvYSAzMCBtaW4gdsOgIE1ham9yDQp0YW5zdWF0MSA8LSB0YWJsZShkYXRhJGAzMCBtaW5gLCBkYXRhJE1ham9yKSANCnBoYW50cmFtMSA8LSBwcm9wLnRhYmxlKHRhbnN1YXQxLCBtYXJnaW4gPSAxKSAqIDEwMA0KDQpiYW5nMSA8LSBhcy5kYXRhLmZyYW1lKHRhbnN1YXQxKQ0KYmFuZzEkUGVyY2VudGFnZSA8LSByb3VuZChhcy52ZWN0b3IocGhhbnRyYW0xKSwgMikNCg0KY29sbmFtZXMoYmFuZzEpIDwtIGMoIjMwIG1pbiIsICJNYWpvciIsICJGcmVxdWVuY3kiLCAiUGVyY2VudGFnZSIpDQogICAgDQprYWJsZShiYW5nMSwgY2FwdGlvbiA9ICJC4bqjbmcgMTogTcO0IHThuqMgdGjhu5FuZyBrw6ogY+G7p2EgMiBiaeG6v24gMzAgbWluIHbDoCBNYWpvciIpDQoNCiNUcuG7sWMgcXVhbiBow7NhIGNobyBt4buRaSBxdWFuIGjhu4cgZ2nhu69hIDIgYmnhur9uDQpiYW5nMSA8LSBkYXRhICU+JQ0KICBncm91cF9ieShNYWpvciwgYDMwIG1pbmApICU+JQ0KICBzdW1tYXJpc2UoQ291bnQgPSBuKCksIC5ncm91cHMgPSAiZHJvcCIpDQoNCmdncGxvdChiYW5nMSwgYWVzKHggPSBNYWpvciwgeSA9IENvdW50LCBmaWxsID0gYDMwIG1pbmApKSArDQogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsNCiAgbGFicyh0aXRsZSA9ICJIw6xuaCAxOiBN4buRaSBxdWFuIGjhu4cgZ2nhu69hIHnhur91IHThu5Ega2jDtG5nIHRo4buDIG5n4bunIHRyb25nIDMwIHBow7p0IMSR4bqndSB2w6AgbmfDoG5oIG5naOG7gSIsDQogICAgICAgeCA9ICJOZ8Ogbmggbmdo4buBIiwgeSA9ICJLaMO0bmcgdGjhu4Mgbmfhu6cgdHJvbmcgMzAgcGjDunQgxJHhuqd1IiwgZmlsbCA9ICJU4bqnbiBzdeG6pXQiKSArDQogIHRoZW1lX21pbmltYWwoKSArDQogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIGhqdXN0ID0gMSkpDQoNCmBgYA0KDQpOaMOsbiBjaHVuZywgbmjDs20gKipDw7RuZyBuZ2jhu4cgdGjDtG5nIHRpbiAoQ29tcHV0ZXIgU2NpZW5jZSkqKiBjw7MgdOG7tyBs4buHIGfhurdwIGtow7Mga2jEg24gduG7gSBnaeG6pWMgbmfhu6cgY2FvIG5o4bqldCB0cm9uZyB04bqldCBj4bqjIGPDoWMgbeG7qWMgdOG6p24gc3XhuqV0LiBD4bulIHRo4buDLCBjw7MgKioxNiBuZ8aw4budaSoqIHRodeG7mWMgbmjDs20gbsOgeSBn4bq3cCB0w6xuaCB0cuG6oW5nIGtow7Mgbmfhu6cgKioyIGzhuqduL3R14bqnbioqLCBjaGnhur9tICoqNzYuMTklKiogc+G7kSBuZ8aw4budaSB0cm9uZyBuZ8OgbmguIE5nb8OgaSByYSwgY8OzIHThu5tpICoqMTEgbmfGsOG7nWkqKiBjaGnhur9tICoqNTUlKiogY2hvIGJp4bq/dCBn4bq3cCB24bqlbiDEkeG7gSBt4bqldCBuZ+G7pyAqKmjGoW4gMyBs4bqnbi90deG6p24qKiDigJMgxJHDonkgbMOgIHThu7cgbOG7hyBjYW8gbmjhuqV0IHNvIHbhu5tpIGPDoWMgbmjDs20gbmfDoG5oIGtow6FjLCBjaG8gdGjhuqV5IGNo4bqldCBsxrDhu6NuZyBnaeG6pWMgbmfhu6cgdHJvbmcgbmjDs20gbsOgeSDhu58gbeG7qWMgYsOhbyDEkeG7mW5nLiBOZ3V5w6puIG5ow6JuIGPDsyB0aOG7gyDEkeG6v24gdOG7qyDEkeG6t2MgdGjDuSBjw7RuZyB2aeG7h2MgY8SDbmcgdGjhurNuZywgw6FwIGzhu7FjIGNhbywgdGjhu51pIGdpYW4gbMOgbSB2aeG7h2Mga8OpbyBkw6BpIHbDoCB0aeG6v3AgeMO6YyB0aMaw4budbmcgeHV5w6puIHbhu5tpIHRoaeG6v3QgYuG7iyDEkWnhu4duIHThu60uDQoNCk5nw6BuaCAqKkPDtG5nIG5naOG7hyBzaW5oIGjhu41jIChCaW90ZWNoKSoqIGPFqW5nIGNobyB0aOG6pXkgcGjDom4gYuG7kSBt4bqldCBuZ+G7pyDhu58gbeG7qWMgY2FvIG5oxrBuZyBy4bqjaSByw6FjIOG7nyBuaGnhu4F1IG3hu6ljIMSR4buZLiBD4bulIHRo4buDLCBjw7MgKioxMiBuZ8aw4budaSoqIGtow7RuZyBn4bq3cCB24bqlbiDEkeG7gSB24buBIGdp4bqlYyBuZ+G7pyAoKiowIGzhuqduL3Row6FuZyoqKSwgY2hp4bq/bSAqKjQxLjM4JSoqIHPhu5EgbmfGsOG7nWkgdHJvbmcgbmfDoG5oOyAqKjUgbmfGsOG7nWkqKiBi4buLIG3huqV0IG5n4bunICoqMSBs4bqnbi90deG6p24qKiAoKioyOS40MSUqKikgdsOgICoqNCBuZ8aw4budaSoqIHRoxrDhu51uZyB4dXnDqm4gbeG6pXQgbmfhu6cgKipoxqFuIDMgbOG6p24vdHXhuqduKiosIGNoaeG6v20gKioyMCUqKi4gxJBp4buBdSBuw6B5IGNobyB0aOG6pXkgdMOsbmggdHLhuqFuZyBt4bqldCBuZ+G7pyB0dXkgY8OzIHh14bqldCBoaeG7h24gbmjGsG5nIGtow7RuZyB04bqtcCB0cnVuZyB2w6BvIG3hu5l0IG3hu6ljIHThuqduIHN14bqldCBj4bulIHRo4buDLCBwaOG6o24gw6FuaCB0w61uaCBjaOG6pXQgY8O0bmcgdmnhu4djIGhv4bq3YyBs4buRaSBz4buRbmcgxJFhIGThuqFuZyB0cm9uZyBuaMOzbSBuw6B5Lg0KDQpN4buZdCBz4buRIG5nw6BuaCBuaMawICoqQ8ahIGtow60gKE1lY2hhbmljYWwgRW5naW5lZXJpbmcpKiogdsOgICoqWMOieSBk4buxbmcgKENpdmlsIEVuZ2luZWVyaW5nKSoqIGPDsyBz4buRIGzGsOG7o25nIG5nxrDhu51pIHRoYW0gZ2lhIHbhu6thIHBo4bqjaSB2w6AgcGjDom4gYuG7kSB0xrDGoW5nIMSR4buRaSDEkeG7gXUgZ2nhu69hIGPDoWMgbeG7qWMgxJHhu5kgdOG6p24gc3XhuqV0LiBWw60gZOG7pSwgbmfDoG5oIEPGoSBraMOtIGPDsyAqKjQgbmfGsOG7nWkqKiBraMO0bmcgZ+G6t3AgduG6pW4gxJHhu4EgduG7gSBnaeG6pWMgbmfhu6cgdsOgICoqMiBuZ8aw4budaSoqIG3huqV0IG5n4bunIHThu6sgMeKAkzIgbOG6p24vdHXhuqduLCBraMO0bmcgY8OzIHRyxrDhu51uZyBo4bujcCBuw6BvIG3huqV0IG5n4bunIG5naGnDqm0gdHLhu41uZywgY2hvIHRo4bqleSBt4bupYyDEkeG7mSDhuqNuaCBoxrDhu59uZyB0aOG6pXAgaMahbi4NCg0KTmfGsOG7o2MgbOG6oWksIGPDoWMgbmfDoG5oIG5oxrAgKirhu6huZyBk4bulbmcgbcOheSB0w61uaCAoQ29tcHV0ZXIgQXBwbGljYXRpb24pKiosICoqVGjhu51pIHRyYW5nIChGYXNoaW9uIFRlY2hub2xvZ3kpKiosICoqTmdo4buHIHRodeG6rXQgKEZpbmUgQXJ0cykqKiB2w6AgKipUaOG7kW5nIGvDqiAoU3RhdGlzdGljcykqKiBjw7Mgc+G7kSBuZ8aw4budaSB0aGFtIGdpYSBy4bqldCDDrXQgKGTGsOG7m2kgNCBuZ8aw4budaS9uaMOzbSksIGtoaeG6v24gdmnhu4djIHLDunQgcmEgbmjhuq1uIMSR4buLbmggdOG7lW5nIHF1w6F0IGfhurdwIG5oaeG7gXUgaOG6oW4gY2jhur8uIFR1eSBuaGnDqm4sIHBo4bqnbiBs4bubbiBjw6FjIGPDoSBuaMOibiB0aHXhu5ljIG5ow7NtIG7DoHkgxJHhu4F1IGtow7RuZyBiw6FvIGPDoW8gdMOsbmggdHLhuqFuZyBt4bqldCBuZ+G7pyB0aMaw4budbmcgeHV5w6puLg0KDQpDdeG7kWkgY8O5bmcsIGPDoWMgbmfDoG5oIG5oxrAgKipUaGnhur90IGvhur8gKERlc2lnbmVyKSoqIHbDoCAqKkjDs2EgaOG7jWMgKENoZW1pY2FsKSoqIGdoaSBuaOG6rW4gcuG6pXQgw610IHRyxrDhu51uZyBo4bujcCBt4bqldCBuZ+G7pyBuZ2hpw6ptIHRy4buNbmcuIENo4buJIGPDsyAqKjEgbmfGsOG7nWkqKiB0cm9uZyBt4buXaSBuZ8OgbmggbsOgeSBjaG8gYmnhur90IGLhu4sgbeG6pXQgbmfhu6cgKipoxqFuIDMgbOG6p24vdHXhuqduKiogaG/hurdjICoqMSBs4bqnbi90deG6p24qKiwgY8OybiBs4bqhaSDEkeG7gXUga2jDtG5nIGfhurdwIHbhuqVuIMSR4buBIHbhu4EgZ2nhuqVjIG5n4bunLg0KDQpCaeG7g3UgxJHhu5MgY+G7mXQgxJFpIGvDqG0gxJHDoyB0cuG7sWMgcXVhbiBow7NhIHLDtSBz4buxIGtow6FjIGJp4buHdCBnaeG7r2EgY8OhYyBuZ8OgbmgsIMSR4bq3YyBiaeG7h3QgbuG7lWkgYuG6rXQg4bufIG5ow7NtIEPDtG5nIG5naOG7hyB0aMO0bmcgdGluIHbDoCBDw7RuZyBuZ2jhu4cgc2luaCBo4buNYyDigJMgduG7m2kgY2hp4buBdSBjYW8gY+G7mXQgbOG7m24g4bufIG3hu6ljICoqMiBs4bqnbi90deG6p24qKiB2w6AgKipoxqFuIDMgbOG6p24vdHXhuqduKiosIHBo4bqjbiDDoW5oIHPhu7EgY8SDbmcgdGjhurNuZyDEkeG6t2MgdGjDuSB0cm9uZyBjw7RuZyB2aeG7h2MgY+G7p2EgaGFpIG5ow7NtIG5nw6BuaCBuw6B5Lg0KDQoqKktp4buDbSDEkeG7i25oIENoaSBiw6xuaCBwaMawxqFuZyoqDQoNCkdp4bqjIHRodXnhur90IGtp4buDbSDEkeG7i25oOg0KDQokSF8wJDogSGFpIGJp4bq/biDEkeG7mWMgbOG6rXAga2jDtG5nIGPDsyBz4buxIHTDoWMgxJHhu5luZw0KDQokSF8xJDogSGFpIGJp4bq/biBjw7MgbeG7kWkgcXVhbiBo4buHDQoNCmBgYHtyfQ0KY2hpc3EudGVzdCh0YW5zdWF0MSkNCmBgYA0KDQpE4buxYSB0csOqbiBr4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaCBDaGktYsOsbmggcGjGsMahbmcgZ2nhu69hIGhhaSBiaeG6v24gKiozMCBtaW4qKiB2w6AgKipNYWpvcioqLCBnacOhIHRy4buLIHRo4buRbmcga8OqIENoaS1zcXVhcmVkIMSR4bqhdCAyNS42MjYgduG7m2kgMzAgYuG6rWMgdOG7sSBkbyB2w6AgZ2nDoSB0cuG7iyBwIHTGsMahbmcg4bupbmcgbMOgIDAuNjk0MS4gVsOsIHAtdmFsdWUgPiAwLjA1LCBuw6puIGtow7RuZyDEkeG7pyBi4bqxbmcgY2jhu6luZyDEkeG7gyBiw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3QgJEhfMCQsIHThu6ljIGzDoCAqKktIw5RORyBjw7MgbeG7kWkgbGnDqm4gaOG7hyBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqKiogZ2nhu69hIGLhu4sga2jDsyBuZ+G7pyB0cm9uZyAzMCBwaMO6dCDEkeG6p3UgdsOgIG5nw6BuaCBuZ2jhu4EuIE7Ds2kgY8OhY2gga2jDoWMsIGhhaSBiaeG6v24gbsOgeSBjw7MgeHUgaMaw4bubbmcgxJHhu5ljIGzhuq1wIHbhu5tpIG5oYXUgdHJvbmcgYuG7mSBk4buvIGxp4buHdS4gVHV5IG5oacOqbiwgY+G6p24gbMawdSDDvSBy4bqxbmcga2nhu4NtIMSR4buLbmggY8OzIGPhuqNuaCBiw6FvIHbhu4EgxJHhu5kgY2jDrW5oIHjDoWMgZG8gbeG7mXQgc+G7kSDDtCB0cm9uZyBi4bqjbmcgdOG6p24gc3XhuqV0IGPDsyBnacOhIHRy4buLIGvhu7MgduG7jW5nIG5o4buPLCDEkWnhu4F1IG7DoHkgY8OzIHRo4buDIOG6o25oIGjGsOG7n25nIMSR4bq/biDEkeG7mSB0aW4gY+G6rXkgY+G7p2Ega+G6v3QgcXXhuqMuDQoNCiMjIyAqKjMuMS4yLiBUw6FjIMSR4buZbmcgY+G7p2EgeeG6v3UgdOG7kSB04buJbmggZOG6rXkgbuG7rWEgxJHDqm0gaG/hurdjIHLhuqFuZyBzw6FuZyB0cm9uZyBjw6FjIG5nw6BuaCBuZ2jhu4EqKg0KDQoqKkLhuqNuZyB04bqnbiBzdeG6pXQgY2jDqW8gdsOgIGJp4buDdSDEkeG7kyB0cuG7sWMgcXVhbiBow7NhKioNCg0KYGBge3J9DQojIELhuqNuZyB04bqnbiBzdeG6pXQgZ2nhu69hIFdha2UgdXAgbWlkZGxlIHbDoCBNYWpvcg0KdGFuc3VhdDIgPC0gdGFibGUoZGF0YSRgV2FrZSB1cCBtaWRkbGVgLCBkYXRhJE1ham9yKSANCnBoYW50cmFtMiA8LSBwcm9wLnRhYmxlKHRhbnN1YXQyLCBtYXJnaW4gPSAxKSAqIDEwMA0KDQpiYW5nMiA8LSBhcy5kYXRhLmZyYW1lKHRhbnN1YXQyKQ0KYmFuZzIkUGVyY2VudGFnZSA8LSByb3VuZChhcy52ZWN0b3IocGhhbnRyYW0yKSwgMikNCg0KY29sbmFtZXMoYmFuZzIpIDwtIGMoIldha2UgdXAgMSIsICJNYWpvciIsICJGcmVxdWVuY3kiLCAiUGVyY2VudGFnZSIpDQogICAgDQprYWJsZShiYW5nMiwgY2FwdGlvbiA9ICJC4bqjbmcgMjogTcO0IHThuqMgdGjhu5FuZyBrw6ogY+G7p2EgMiBiaeG6v24gV2FrZSB1cCBtaWRkbGUgdsOgIE1ham9yIikNCg0KI1RydWMgcXVhbiBob2ENCmJhbmcyIDwtIGRhdGEgJT4lDQogIGdyb3VwX2J5KE1ham9yLCBgV2FrZSB1cCBtaWRkbGVgKSAlPiUNCiAgc3VtbWFyaXNlKENvdW50ID0gbigpLCAuZ3JvdXBzID0gImRyb3AiKQ0KDQpnZ3Bsb3QoYmFuZzIsIGFlcyh4ID0gTWFqb3IsIHkgPSBDb3VudCwgZmlsbCA9IGBXYWtlIHVwIG1pZGRsZWApKSArDQogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsNCiAgbGFicyh0aXRsZSA9ICJIw6xuaCAyOiBN4buRaSBxdWFuIGjhu4cgZ2nhu69hIHnhur91IHThu5EgdOG7iW5oIGdp4bqlYyBu4butYSDEkcOqbSBob+G6t2MgcuG6oW5nIHPDoW5nIHbDoCBuZ8Ogbmggbmdo4buBIiwNCiAgICAgICB4ID0gIk5nw6BuaCBuZ2jhu4EiLCB5ID0gIkLhu4sgdOG7iW5oIGdp4bqlYyBu4butYSDEkcOqbSBob+G6t2MgcuG6oW5nIHPDoW5nIiwgZmlsbCA9ICJU4bqnbiBzdeG6pXQiKSArDQogIHRoZW1lX21pbmltYWwoKSArDQogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIGhqdXN0ID0gMSkpDQpgYGANCg0KTmjDs20gKipDw7RuZyBuZ2jhu4cgdGjDtG5nIHRpbiAoQ29tcHV0ZXIgU2NpZW5jZSkqKiB0aeG6v3AgdOG7pWMgY2hvIHRo4bqleSB04buJIGzhu4cgcuG7kWkgbG/huqFuIGdp4bqlYyBuZ+G7pyBjYW8sIMSR4bq3YyBiaeG7h3QgbMOgIHTDrG5oIHRy4bqhbmcgKip04buJbmggZ2nhuqVjIG7hu61hIMSRw6ptIGhv4bq3YyBy4bqhbmcgc8OhbmcqKi4gQ8OzIHThu5tpICoqMzUgbmfGsOG7nWkqKiB0cm9uZyBuaMOzbSBuw6B5IGNobyBiaeG6v3QgKipraMO0bmcgZ+G6t3AgdMOsbmggdHLhuqFuZyBuw6B5KiogKDAgbOG6p24vdGjDoW5nKSwgY2hp4bq/bSAqKjQ2LjY3JSoqLiBUdXkgbmhpw6puLCBiw6puIGPhuqFuaCDEkcOzLCBjw7MgKioxMiBuZ8aw4budaSoqICgqKjI5LjI3JSoqKSBjaG8gYmnhur90IGjhu40gdOG7iW5oIGdp4bqlYyAqKjEgbOG6p24vdHXhuqduKiosIHbDoCAqKjggbmfGsOG7nWkqKiAoKio0MCUqKikgdOG7iW5oIGdp4bqlYyB04bubaSAqKjIgbOG6p24vdHXhuqduKiosIGNobyB0aOG6pXkgduG6pW4gxJHhu4EgZ2nhuqVjIG5n4bunIHbhuqtuIGtow6EgcGjhu5UgYmnhur9uLiBOZ2/DoGkgcmEsIGPDsyAqKjYgbmfGsOG7nWkqKiAoKioyNi4wOSUqKikgdOG7iW5oIGdp4bqlYyAqKmjGoW4gMyBs4bqnbi90deG6p24qKiwgY2hvIHRo4bqleSBt4bupYyDEkeG7mSBy4buRaSBsb+G6oW4gZ2nhuqVjIG5n4bunIHRydW5nIGLDrG5oIMSR4bq/biBu4bq3bmcg4bufIG5nw6BuaCBuw6B5Lg0KDQpOZ8OgbmggKipDw7RuZyBuZ2jhu4cgc2luaCBo4buNYyAoQmlvdGVjaCkqKiBjxaluZyBnaGkgbmjhuq1uIHPhu7EgcGjDom4gYuG7kSDEkcOhbmcga+G7gyDhu58gbmhp4buBdSBt4bupYyB04bqnbiBzdeG6pXQ6ICoqMjEgbmfGsOG7nWkqKiAoY2hp4bq/bSAqKjI4LjAwJSoqKSBraMO0bmcgZ+G6t3AgdMOsbmggdHLhuqFuZyB04buJbmggZ2nhuqVjIHbDoG8gYmFuIMSRw6ptLCAqKjggbmfGsOG7nWkqKiAoKioxOS41MSUqKikgYuG7iyB04buJbmggZ2nhuqVjICoqMSBs4bqnbi90deG6p24qKiwgdsOgICoqNCBuZ8aw4budaSoqICgqKjE3LjM5JSoqKSBi4buLIOG6o25oIGjGsOG7n25nICoqaMahbiAzIGzhuqduL3R14bqnbioqLCBjw7luZyB24bubaSAqKjMgbmfGsOG7nWkqKiAoKioxNS4wMCUqKikgdOG7iW5oIGdp4bqlYyAqKjIgbOG6p24vdHXhuqduKiouIE5o4buvbmcgY29uIHPhu5EgbsOgeSBjaG8gdGjhuqV5IG3hu6ljIMSR4buZIOG6o25oIGjGsOG7n25nIGtow6EgxJFhIGThuqFuZyB0cm9uZyBuaMOzbSBuw6B5LCB0xrDGoW5nIHThu7EgeHUgaMaw4bubbmcgcXVhbiBzw6F0IMSRxrDhu6NjIHRyb25nIGJp4bq/biBraMOzIG5n4bunLg0KDQpN4buZdCBz4buRIG5nw6BuaCBuaMawICoqVGhp4bq/dCBr4bq/IChEZXNpZ25lcikqKiB2w6AgKipUaOG7nWkgdHJhbmcgKEZhc2hpb24gVGVjaG5vbG9neSkqKiBjw7Mgc+G7sSBwaMOibiBi4buRIHTGsMahbmcgxJHhu5FpIMSR4buBdSBuaMawbmcgduG7m2kgcXV5IG3DtCBt4bqrdSBuaOG7jy4gVHJvbmcgbmfDoG5oIFRoaeG6v3Qga+G6vywgKio4IG5nxrDhu51pKiogKGNoaeG6v20gKioxMC42NyUqKikga2jDtG5nIGLhu4sgdOG7iW5oIGdp4bqlYywgKiozIG5nxrDhu51pKiogKCoqNy4zMiUqKikgYuG7iyB04buJbmggZ2nhuqVjICoqMSBs4bqnbi90deG6p24qKiwgdsOgIGNo4buJIGPDsyAqKjIgbmfGsOG7nWkqKiBn4bq3cCB0w6xuaCB0cuG6oW5nIG7DoHkg4bufIG3hu6ljIG5naGnDqm0gdHLhu41uZyBoxqFuLiBUxrDGoW5nIHThu7EsIG5nw6BuaCBUaOG7nWkgdHJhbmcgY8OzICoqMyBuZ8aw4budaSoqIGtow7RuZyBi4buLIOG6o25oIGjGsOG7n25nLCAqKjMgbmfGsOG7nWkqKiB04buJbmggZ2nhuqVjICoqMSBs4bqnbi90deG6p24qKiwgdsOgICoqMiBuZ8aw4budaSoqIHThu6sgbeG7qWMgKioyIGzhuqduL3R14bqnbiB0cuG7nyBsw6puKiouDQoNCk5nw6BuaCAqKk3hu7kgdGh14bqtdCAoRmluZSBBcnRzKSoqIGPDsyBxdXkgbcO0IG3huqt1IG5o4buPIG5oxrBuZyDEkeG6t2MgYmnhu4d0IMSRw6FuZyBjaMO6IMO9IHbhu5tpICoqMSB0csaw4budbmcgaOG7o3AgZHV5IG5o4bqldCoqIGLhu4sgdOG7iW5oIGdp4bqlYyAqKmjGoW4gMyBs4bqnbi90aMOhbmcqKiAoKioxMDAlKiogdHJvbmcgbmjDs20gbsOgeSksIHbDoCAqKjQgbmfGsOG7nWkqKiBraMOhYyB04buJbmggZ2nhuqVjICoqMSBs4bqnbi90deG6p24qKiwgY2hvIHRo4bqleSDEkcOieSBsw6AgbmjDs20gY8OzIHh1IGjGsOG7m25nIGdp4bqlYyBuZ+G7pyBrw6ltIG5oxrBuZyBj4bqnbiB0aMOqbSBk4buvIGxp4buHdSDEkeG7gyB4w6FjIG5o4bqtbiBjaOG6r2MgY2jhuq9uLg0KDQrEkOG7kWkgduG7m2kgY8OhYyBuZ8OgbmggbmjGsCAqKkPGoSBraMOtIChNZWNoYW5pY2FsKSoqLCAqKkvhu7kgdGh14bqtdCBjw7RuZyBuZ2hp4buHcCAoSW5kdXN0cmlhbCBFbmdpbmVlcmluZykqKiB2w6AgKipUaOG7kW5nIGvDqiAoU3RhdGlzdGljcykqKiwgdOG7tyBs4buHIHThu4luaCBnaeG6pWMgZGFvIMSR4buZbmcgY2jhu6cgeeG6v3UgcXVhbmggbeG7qWMgKiox4oCTMiBs4bqnbi90deG6p24qKiwgbmjGsG5nIGtow7RuZyBjw7MgdHLGsOG7nW5nIGjhu6NwIG7DoG8gdsaw4bujdCBxdcOhICoqMyBs4bqnbi90aMOhbmcqKiwgY2hvIHRo4bqleSBt4bupYyDEkeG7mSDhuqNuaCBoxrDhu59uZyB0aOG6pXAgaMahbiBzbyB24bubaSBjw6FjIG5nw6BuaCBjw7RuZyBuZ2jhu4cuDQoNCkJp4buDdSDEkeG7kyBtaW5oIGjhu41hIHRy4buxYyBxdWFuIMSRw6MgdGjhu4MgaGnhu4duIHLDtSBz4buxIGNow6puaCBs4buHY2gsIMSR4bq3YyBiaeG7h3QgbMOgIOG7nyBuaMOzbSAqKkNvbXB1dGVyIFNjaWVuY2UqKiwgduG7m2kgY2hp4buBdSBjYW8gY+G7mXQgdsaw4bujdCB0cuG7mWkg4bufIG3hu6ljICoqMCBs4bqnbi90aMOhbmcqKiB2w6AgKiox4oCTMiBs4bqnbi90deG6p24qKiwgcGjhuqNuIMOhbmggc+G7sSBwaOG7lSBiaeG6v24gY+G7p2EgdMOsbmggdHLhuqFuZyB04buJbmggZ2nhuqVjIGdp4buvYSDEkcOqbSDhu58gbmjDs20gbmfDoG5oIG7DoHkgc28gduG7m2kgY8OhYyBuZ8OgbmggY8OybiBs4bqhaS4NCg0KYGBge3J9DQpjaGlzcS50ZXN0KHRhbnN1YXQyKQ0KYGBgDQpE4buxYSB0csOqbiBr4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaCBDaGktYsOsbmggcGjGsMahbmcgZ2nhu69hIGhhaSBiaeG6v24gKipXYWtlIHVwIG1pZGRsZSoqIHbDoCAqKk1ham9yKiosIGdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQgxJHhuqF0IDgwLjU5IHbhu5tpIDUwIGLhuq1jIHThu7EgZG8gdsOgIGdpw6EgdHLhu4sgcCB0xrDGoW5nIOG7qW5nIGzDoCAwLjAwMzk0Ni4gVsOsIHAtdmFsdWUgPCAwLjA1LCBuw6puIGtow7RuZyDEkeG7pyBi4bqxbmcgY2jhu6luZyDEkeG7gyBiw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3QgJEhfMSQsIHThu6ljIGzDoCAqKkPDkyBt4buRaSBsacOqbiBo4buHIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6oqKiBnaeG7r2EgeeG6v3UgdOG7kSBi4buLIHThu4luaCBk4bqteSBnaeG7r2EgbuG7rWEgxJHDqm0gaG/hurdjIHLhuqFuZyBzw6FuZyB2w6AgbmfDoG5oIG5naOG7gS4gTsOzaSBjw6FjaCBraMOhYywgaGFpIGJp4bq/biBuw6B5IGPDsyB4dSBoxrDhu5tuZyBwaOG7pSB0aHXhu5ljIHbhu5tpIG5oYXUgdHJvbmcgYuG7mSBk4buvIGxp4buHdS4gVHV5IG5oacOqbiwgY+G6p24gbMawdSDDvSBy4bqxbmcga2nhu4NtIMSR4buLbmggY8OzIGPhuqNuaCBiw6FvIHbhu4EgxJHhu5kgY2jDrW5oIHjDoWMgZG8gbeG7mXQgc+G7kSDDtCB0cm9uZyBi4bqjbmcgdOG6p24gc3XhuqV0IGPDsyBnacOhIHRy4buLIGvhu7MgduG7jW5nIG5o4buPLCDEkWnhu4F1IG7DoHkgY8OzIHRo4buDIOG6o25oIGjGsOG7n25nIMSR4bq/biDEkeG7mSB0aW4gY+G6rXkgY+G7p2Ega+G6v3QgcXXhuqMuDQoNCioqVOG7tyBs4buHIFJSIHbDoCBPUioqDQoNCmBgYHtyfQ0Kcmlza3JhdGlvKHRhbnN1YXQyKQ0KYGBgDQoNCkvhur90IHF14bqjIHBow6JuIHTDrWNoIGNobyB0aOG6pXkgY2jhu4kgbmjDs20gbmfDoG5oIG5naOG7gSBi4buLIHThu4luaCBnaeG6pWMgbuG7rWEgxJHDqm0gaG/hurdjIHLDoW5nIHPDoW5nIHRyw6puIDMgbOG6p24gbeG7mXQgdHXhuqduIGPDsyBnacOhIHRy4buLIFJSID0gMSwgbmjGsG5nIGRvIGNo4buJIGPDsyAqKjEgbmfGsOG7nWkqKiB0cm9uZyBuaMOzbSBuw6B5IG7Dqm4ga2jDtG5nIMSR4bunIGPGoSBz4bufIMSR4buDIHN1eSBsdeG6rW4uIEPDoWMgbmjDs20gY8OybiBs4bqhaSDEkeG7gXUgY2hvIGvhur90IHF14bqjICoqTmFOKiogZG8gc+G7kSBsxrDhu6NuZyBxdcOhIMOtdCBob+G6t2MgcGjDom4gYuG7kSBt4bqldCBjw6JuIMSR4buRaSwga2hp4bq/biBraMO0bmcgdGjhu4MgdMOtbmggdG/DoW4gUlIgdsOgIGtob+G6o25nIHRpbiBj4bqteS4gRMO5IG5ow7NtICoqPjMvd2VlayoqIGPDsyBnacOhIHRy4buLIHAgdHJvbmcga2nhu4NtIMSR4buLbmggY2hpIGLDrG5oIHBoxrDGoW5nIGzDoCAqKjAuMDA3NioqLCBuaMawbmcgY8OhYyBraeG7g20gxJHhu4tuaCBraMOhYyBs4bqhaSBraMO0bmcgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqi4gTmjGsCB24bqteSwgKipjaMawYSB0aOG7gyBr4bq/dCBsdeG6rW4gcsO1IHLDoG5nIHbhu4EgbeG7kWkgbGnDqm4gaOG7hyBnaeG7r2EgbmfDoG5oIG5naOG7gSB2w6AgdOG6p24gc3XhuqV0IGLhu4sgdOG7iW5oIGdp4bqlYyBu4butYSDEkcOqbSBob+G6t2MgcsOhbmcgc8OhbmcqKiB0cm9uZyBuZ2hpw6puIGPhu6l1IG7DoHkuDQoNCiMjIyAqKjMuMS4zLiBUw6FjIMSR4buZbmcgY+G7p2EgeeG6v3UgdOG7kSBt4bqvYyB24buHIHNpbmggdHJvbmcgbMO6YyBuZ+G7pyB0cm9uZyBjw6FjIG5nw6BuaCBuZ2jhu4EqKg0KDQpgYGB7cn0NCiMgQuG6o25nIHThuqduIHN14bqldCBnaeG7r2EgQmF0aHJvb20gdsOgIE1ham9yDQp0YW5zdWF0MyA8LSB0YWJsZShkYXRhJGBCYXRocm9vbWAsIGRhdGEkTWFqb3IpIA0KcGhhbnRyYW0zIDwtIHByb3AudGFibGUodGFuc3VhdDMsIG1hcmdpbiA9IDEpICogMTAwDQoNCmJhbmczIDwtIGFzLmRhdGEuZnJhbWUodGFuc3VhdDMpDQpiYW5nMyRQZXJjZW50YWdlIDwtIHJvdW5kKGFzLnZlY3RvcihwaGFudHJhbTMpLCAyKQ0KDQpjb2xuYW1lcyhiYW5nMykgPC0gYygiQmF0aHJvb20iLCAiTWFqb3IiLCAiRnJlcXVlbmN5IiwgIlBlcmNlbnRhZ2UiKQ0KICAgIA0Ka2FibGUoYmFuZzMsIGNhcHRpb24gPSAiQuG6o25nIDM6IE3DtCB04bqjIHRo4buRbmcga8OqIGPhu6dhIDIgYmnhur9uIEJhdGhyb29tIHbDoCBNYWpvciIpDQoNCiNUcnVjIHF1YW4gaG9hDQpiYW5nMyA8LSBkYXRhICU+JQ0KICBncm91cF9ieShNYWpvciwgYEJhdGhyb29tYCkgJT4lDQogIHN1bW1hcmlzZShDb3VudCA9IG4oKSwgLmdyb3VwcyA9ICJkcm9wIikNCg0KZ2dwbG90KGJhbmczLCBhZXMoeCA9IE1ham9yLCB5ID0gQ291bnQsIGZpbGwgPSBgQmF0aHJvb21gKSkgKw0KICBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5IiwgcG9zaXRpb24gPSAiZG9kZ2UiKSArDQogIGxhYnModGl0bGUgPSAiSMOsbmggMzogTeG7kWkgcXVhbiBo4buHIGdp4buvYSB54bq/dSB04buRIG3huq9jIHbhu4cgc2luaCB0cm9uZyBsw7pjIG5n4bunIHbDoCBuZ8Ogbmggbmdo4buBIiwNCiAgICAgICB4ID0gIk5nYW5oIG5naGUiLCB5ID0gIm3huq9jIHbhu4cgc2luaCB0cm9uZyBsw7pjIG5n4bunIiwgZmlsbCA9ICJU4bqnbiBzdeG6pXQiKSArDQogIHRoZW1lX21pbmltYWwoKSArDQogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIGhqdXN0ID0gMSkpDQpgYGANCg0KTmjDrG4gY2h1bmcsIG5nw6BuaCAqKkPDtG5nIG5naOG7hyB0aMO0bmcgdGluIChDb21wdXRlciBTY2llbmNlKSoqIGdoaSBuaOG6rW4gdOG7tyBs4buHIG5nxrDhu51pIGfhurdwIHbhuqVuIMSR4buBIHBo4bqjaSDEkWkgduG7hyBzaW5oIHRyb25nIGzDumMgbmfhu6cgY2FvIGjGoW4gcsO1IHLhu4d0IHNvIHbhu5tpIGPDoWMgbmfDoG5oIGtow6FjLiBD4bulIHRo4buDLCBjw7MgdOG7m2kgKioxOCBuZ8aw4budaSoqIChjaGnhur9tICoqMzUuMjklKiopIGLDoW8gY8OhbyBi4buLIHTDrG5oIHRy4bqhbmcgbsOgeSDhu58gbeG7qWMgKioxIGzhuqduL3R14bqnbioqLCB2w6AgKio3IG5nxrDhu51pKiogKDM1JSkgZ+G6t3AgdMOsbmggdHLhuqFuZyAqKmjGoW4gMyBs4bqnbi90deG6p24qKi4gTmdvw6BpIHJhLCAqKjQgbmfGsOG7nWkga2jDoWMqKiAoMTYlKSBjxaluZyBn4bq3cCB0w6xuaCB0cuG6oW5nIG7DoHkgKioyIGzhuqduL3R14bqnbioqLiDEkMOieSBsw6AgbmjDs20gZHV5IG5o4bqldCBjw7Mgc+G7kSBsxrDhu6NuZyDEkcOhbmcga+G7gyDhu58gdOG6pXQgY+G6oyBjw6FjIG3hu6ljIHThuqduIHN14bqldCwgcGjhuqNuIMOhbmggcsO1IG7DqXQgduG7gSBjaOG6pXQgbMaw4bujbmcgZ2nhuqVjIG5n4bunIGvDqW0gdsOgIGto4bqjIG7Eg25nIGLhu4sgZ2nDoW4gxJFv4bqhbiBzaW5oIGzDvSBjYW8uDQoNCk5nw6BuaCAqKkPDtG5nIG5naOG7hyBzaW5oIGjhu41jIChCaW90ZWNoKSoqIGPFqW5nIGPDsyB04bu3IGzhu4cgZ+G6t3AgdMOsbmggdHLhuqFuZyBuw6B5IGtow6EgY2FvLCDEkeG6t2MgYmnhu4d0IGzDoCDhu58gbeG7qWMgKioyIGzhuqduL3R14bqnbioqIHbhu5tpICoqOSBuZ8aw4budaSoqICgzNiUpIHbDoCAqKjEgbOG6p24vdHXhuqduKiogduG7m2kgKioxMCBuZ8aw4budaSoqICgxOS42MSUpLiBUdXkgbmhpw6puLCBuaMOzbSBuw6B5IHbhuqtuIGPDsyB04bubaSAqKjE0IG5nxrDhu51pKioga2jDtG5nIGfhurdwIHTDrG5oIHRy4bqhbmcgbsOgeSwgY2hp4bq/bSAqKjIxLjg4JSoqIOKAkyBjaG8gdGjhuqV5IHPhu7EgcGjDom4gYuG7kSB0xrDGoW5nIMSR4buRaSBwaMOibiB0w6FuLg0KDQpN4buZdCBz4buRIG5nw6BuaCBuaMawICoqVGjhu51pIHRyYW5nIChGYXNoaW9uIFRlY2hub2xvZ3kpKiosICoqTmdo4buHIHRodeG6rXQgKEZpbmUgQXJ0cykqKiB2w6AgKipUaGnhur90IGvhur8gKERlc2lnbmVyKSoqIGdoaSBuaOG6rW4gY8OhYyBt4bupYyB04bqnbiBzdeG6pXQgdHJ1bmcgYsOsbmggKDHigJMyIGzhuqduL3R14bqnbiksIG5oxrBuZyBraMO0bmcgY8OzIHRyxrDhu51uZyBo4bujcCBuw6BvIG3huq9jIOG7nyBt4bupYyBuZ2hpw6ptIHRy4buNbmcgKD4zIGzhuqduL3Row6FuZykuIE5nxrDhu6NjIGzhuqFpLCBuZ8OgbmggKipYw6J5IGThu7FuZyAoQ2l2aWwpKiogY8OzIDEgbmfGsOG7nWkgZHV5IG5o4bqldCBt4bqvYyB0csOqbiAzIGzhuqduL3Row6FuZyAoY2hp4bq/bSAxMDAlIHRyb25nIG5ow7NtIG7DoHkpLCB0dXkgbmhpw6puIHPhu5EgbMaw4bujbmcgbeG6q3UgcXXDoSBuaOG7jyDEkeG7gyDEkcawYSByYSBr4bq/dCBsdeG6rW4gY2jhuq9jIGNo4bqvbi4NCg0KTeG7mXQgc+G7kSBuZ8OgbmggbmjGsCAqKlRo4buRbmcga8OqIChTdGF0aXN0aWNzKSoqIHbDoCAqKkvhu7kgdGh14bqtdCBjxqEga2jDrSAoTWVjaGFuaWNhbCkqKiBjw7Mgc+G7sSB4deG6pXQgaGnhu4duIHLhuqNpIHLDoWMg4bufIGPDoWMgbeG7qWMgdOG6p24gc3XhuqV0LCBuaMawbmcgcGjhuqduIGzhu5tuIGPDoWMgZ2nDoSB0cuG7iyDEkeG7gXUg4bufIG3hu6ljIHRo4bqlcCBob+G6t2MgdHJ1bmcgYsOsbmgsIGtow7RuZyBjaG8gdGjhuqV5IGThuqV1IGhp4buHdSBu4buVaSBi4bqtdC4gxJDDoW5nIGNow7ogw70sIGtow7RuZyBuZ8OgbmggbsOgbyBjw7MgdHLGsOG7nW5nIGjhu6NwIG7DoG8gbeG6r2Mgbmhp4buBdSBoxqFuIDMgbOG6p24gbeG7l2kgdGjDoW5nIG5nb+G6oWkgdHLhu6sgKipYw6J5IGThu7FuZyoqLCDEkWnhu4F1IG7DoHkgY2hvIHRo4bqleSB04bqnbiBzdeG6pXQgbeG6r2MgY2FvIGzDoCBraMOhIGhp4bq/bSB0cm9uZyB0b8OgbiBi4buZIG3huqt1IGto4bqjbyBzw6F0Lg0KDQpCaeG7g3UgxJHhu5MgY+G7mXQgbWluaCBo4buNYSByw7Ugc+G7sSBraMOhYyBiaeG7h3QgbuG7lWkgYuG6rXQg4bufIG5nw6BuaCBDw7RuZyBuZ2jhu4cgdGjDtG5nIHRpbiwga2hpIGNoaeG7gXUgY2FvIGPhu5l0IHThuqFpIGPDoWMgbeG7qWMgMeKAkzIgbOG6p24vdHXhuqduIHbDoCA+MyBs4bqnbi90deG6p24gdsaw4bujdCB0cuG7mWkgaMahbiBo4bqzbiBjw6FjIG5nw6BuaCBraMOhYywgY2hvIHRo4bqleSBuZ8Ogbmggbmdo4buBIG7DoHkgY8OzIHRo4buDIGxpw6puIHF1YW4gxJHhur9uIHLhu5FpIGxv4bqhbiBzaW5oIGzDvSB0cm9uZyBraGkgbmfhu6cg4oCTIGPhuqduIMSRxrDhu6NjIHhlbSB4w6l0IHPDonUgaMahbiB0cm9uZyBjw6FjIHBow6JuIHTDrWNoIHRp4bq/cCB0aGVvLg0KDQpgYGB7cn0NCmNoaXNxLnRlc3QodGFuc3VhdDMpDQpgYGANCg0KROG7sWEgdHLDqm4ga+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmggQ2hpLWLDrG5oIHBoxrDGoW5nIGdp4buvYSBoYWkgYmnhur9uICoqQmF0aHJvb20qKiB2w6AgKipNYWpvcioqLCBnacOhIHRy4buLIHRo4buRbmcga8OqIENoaS1zcXVhcmVkIMSR4bqhdCAxMTIuNjYgduG7m2kgNjAgYuG6rWMgdOG7sSBkbyB2w6AgZ2nDoSB0cuG7iyBwIHTGsMahbmcg4bupbmcgbMOgIDQuNTc3ZS0wNS4gVsOsIHAtdmFsdWUgPCAwLjA1LCBuw6puIGtow7RuZyDEkeG7pyBi4bqxbmcgY2jhu6luZyDEkeG7gyBiw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3QgJEhfMSQsIHThu6ljIGzDoCAqKkPDkyBt4buRaSBsacOqbiBo4buHIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6oqKiBnaeG7r2EgeeG6v3UgdOG7kSBt4bqvYyDEkWkgduG7hyBzaW5oIGzDumMgbmfhu6cgdsOgIG5nw6BuaCBuZ2jhu4EuIE7Ds2kgY8OhY2gga2jDoWMsIGhhaSBiaeG6v24gbsOgeSBjw7MgeHUgaMaw4bubbmcgcGjhu6UgdGh14buZYyB24bubaSBuaGF1IHRyb25nIGLhu5kgZOG7ryBsaeG7h3UuIFR1eSBuaGnDqm4sIGPhuqduIGzGsHUgw70gcuG6sW5nIGtp4buDbSDEkeG7i25oIGPDsyBj4bqjbmggYsOhbyB24buBIMSR4buZIGNow61uaCB4w6FjIGRvIG3hu5l0IHPhu5Egw7QgdHJvbmcgYuG6o25nIHThuqduIHN14bqldCBjw7MgZ2nDoSB0cuG7iyBr4buzIHbhu41uZyBuaOG7jywgxJFp4buBdSBuw6B5IGPDsyB0aOG7gyDhuqNuaCBoxrDhu59uZyDEkeG6v24gxJHhu5kgdGluIGPhuq15IGPhu6dhIGvhur90IHF14bqjLg0KDQpgYGB7cn0NCnJpc2tyYXRpbyh0YW5zdWF0MykNCmBgYA0KDQpL4bq/dCBxdeG6oyBjaG8gdGjhuqV5IGto4bqjIG7Eg25nIGLhu4sgdOG7iW5oIGdp4bqlYyBnaeG7r2EgxJHDqm0g4bufIG5ow7NtIOKAnDHigJMyIGzhuqduL3R14bqnbuKAnSBjYW8gaMahbiBraG/huqNuZyAzLjUgbOG6p24gc28gduG7m2kgbmjDs20g4oCca2jDtG5nIGJhbyBnaeG7nSBi4buLIHThu4luaOKAnSwgduG7m2kga2hv4bqjbmcgdGluIGPhuq15IDk1JSBsw6AgWzIuMTEsIDUuODNdLiDEkGnhu4F1IG7DoHkgY2hvIHRo4bqleSBt4buRaSBsacOqbiBo4buHIMSRw6FuZyBr4buDIGdp4buvYSB04bqnbiBzdeG6pXQgYuG7iyB04buJbmggZ2nhuqVjIHbDoCBt4buZdCBz4buRIG5nw6BuaCBuZ2jhu4EgbmjhuqV0IMSR4buLbmguIOG7niBjw6FjIG5ow7NtIGPDsyB04bqnbiBzdeG6pXQgY2FvIGjGoW4sIG5oxrAg4oCcbeG7l2kgxJHDqm3igJ0sIGto4bqjIG7Eg25nIGLhu4sgdOG7iW5oIGdp4bqlYyBjYW8gaMahbiBn4bqnbiA1IGzhuqduIHNvIHbhu5tpIG5ow7NtIOKAnGtow7RuZyBiYW8gZ2nhu53igJ0sIHR1eSBuaGnDqm4gbeG7mXQgc+G7kSBSUiBraMO0bmcgeMOhYyDEkeG7i25oIMSRxrDhu6NjIGRvIGThu68gbGnhu4d1IGPDsyBuaMOzbSBi4bqxbmcgMC4gTmjDrG4gY2h1bmcsIG5ndXkgY8ahIGLhu4sgdOG7iW5oIGdp4bqlYyB0xINuZyB0aGVvIHThuqduIHN14bqldCB2w6AgY8OzIGxpw6puIHF1YW4gxJHhur9uIG5nw6BuaCBuZ2jhu4EsIGPhuqduIMSRxrDhu6NjIHhlbSB4w6l0IHRow6ptIHRyb25nIG3DtCBow6xuaCDEkWEgYmnhur9uLg0KDQojIyMgKiozLjEuNC4gVMOhYyDEkeG7mW5nIGPhu6dhIHnhur91IHThu5EgY+G6o20gdGjhuqV5IGzhuqFuaCBraGkgbmfhu6cgdHJvbmcgY8OhYyBuZ8Ogbmggbmdo4buBKioNCg0KYGBge3J9DQojIELhuqNuZyB04bqnbiBzdeG6pXQgZ2nhu69hIENvbGQgdsOgIE1ham9yDQp0YW5zdWF0NCA8LSB0YWJsZShkYXRhJGBDb2xkYCwgZGF0YSRNYWpvcikgDQpwaGFudHJhbTQgPC0gcHJvcC50YWJsZSh0YW5zdWF0NCwgbWFyZ2luID0gMSkgKiAxMDANCg0KYmFuZzQgPC0gYXMuZGF0YS5mcmFtZSh0YW5zdWF0NCkNCmJhbmc0JFBlcmNlbnRhZ2UgPC0gcm91bmQoYXMudmVjdG9yKHBoYW50cmFtNCksIDIpDQoNCmNvbG5hbWVzKGJhbmc0KSA8LSBjKCJDb2xkIiwgIk1ham9yIiwgIkZyZXF1ZW5jeSIsICJQZXJjZW50YWdlIikNCiAgICANCmthYmxlKGJhbmc0LCBjYXB0aW9uID0gIkLhuqNuZyA0OiBNw7QgdOG6oyB0aOG7kW5nIGvDqiBj4bunYSAyIGJp4bq/biBDb2xkIHbDoCBNYWpvciIpDQoNCiNUcnVjIHF1YW4gaG9hDQpiYW5nNCA8LSBkYXRhICU+JQ0KICBncm91cF9ieShNYWpvciwgYENvbGRgKSAlPiUNCiAgc3VtbWFyaXNlKENvdW50ID0gbigpLCAuZ3JvdXBzID0gImRyb3AiKQ0KDQpnZ3Bsb3QoYmFuZzQsIGFlcyh4ID0gTWFqb3IsIHkgPSBDb3VudCwgZmlsbCA9IGBDb2xkYCkpICsNCiAgZ2VvbV9iYXIoc3RhdCA9ICJpZGVudGl0eSIsIHBvc2l0aW9uID0gImRvZGdlIikgKw0KICBsYWJzKHRpdGxlID0gIkjDrG5oIDQ6IE3hu5FpIHF1YW4gaOG7hyBnaeG7r2EgeeG6v3UgdOG7kSBj4bqjbSB0aOG6pXkgbOG6oW5oIGtoaSBuZ+G7pyB2w6AgbmfDoG5oIG5naOG7gSIsDQogICAgICAgeCA9ICJOZ2FuaCBuZ2hlIiwgeSA9ICJj4bqjbSB0aOG6pXkgbOG6oW5oIGtoaSBuZ+G7pyIsIGZpbGwgPSAiVOG6p24gc3XhuqV0IikgKw0KICB0aGVtZV9taW5pbWFsKCkgKw0KICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCBoanVzdCA9IDEpKQ0KYGBgDQoNCkThu68gbGnhu4d1IGNobyB0aOG6pXkgdOG6p24gc3XhuqV0IGPhuqNtIHRo4bqleSBs4bqhbmgga2hpIG5n4bunIGPDsyBz4buxIGtow6FjIGJp4buHdCDEkcOhbmcga+G7gyBnaeG7r2EgY8OhYyBuZ8Ogbmggbmdo4buBLiBD4bulIHRo4buDLCBuZ8aw4budaSBsw6BtIHRyb25nIGzEqW5oIHbhu7FjIENvbXB1dGVyIFNjaWVuY2UgY8OzIHThu7cgbOG7hyBj4bqjbSB0aOG6pXkgbOG6oW5oIGNhbyBoxqFuIHNvIHbhu5tpIGPDoWMgbmfDoG5oIGtow6FjOiBn4bqnbiAzMyUgY+G6o20gdGjhuqV5IGzhuqFuaCBraG/huqNuZyBoYWkgbOG6p24gbeG7l2kgdHXhuqduLCB2w6AgaMahbiAyNyUgYsOhbyBjw6FvIHTDrG5oIHRy4bqhbmcgbsOgeSB44bqjeSByYSB0aMaw4budbmcgeHV5w6puIGjGoW4gYmEgbOG6p24gbeG7l2kgdHXhuqduLiBUcm9uZyBraGkgxJHDsywgY8OhYyBuZ8OgbmggbmjGsCBDaXZpbCwgQ2hlbWljYWwgaGF5IE1lY2hhbmljYWwgY8OzIHh1IGjGsOG7m25nIGPhuqNtIHRo4bqleSBs4bqhbmgga2hpIG5n4bunIMOtdCBoxqFuIMSRw6FuZyBr4buDIOKAkyB0aOG6rW0gY2jDrSBt4buZdCBz4buRIHThuqduIHN14bqldCBuaMawIOKAnDIgbOG6p24vdGjDoW5n4oCdIGhv4bq3YyDigJwxIGzhuqduL3R14bqnbuKAnSDhu58gbeG7mXQgc+G7kSBuZ8OgbmggaOG6p3UgbmjGsCBraMO0bmcgeHXhuqV0IGhp4buHbi4NCg0KxJBp4buBdSBuw6B5IGNobyB0aOG6pXkgbmfDoG5oIG5naOG7gSBjw7MgdGjhu4MgbGnDqm4gcXVhbiDEkeG6v24gxJHhurdjIMSRaeG7g20gbcO0aSB0csaw4budbmcgbMOgbSB2aeG7h2MgaG/hurdjIGzhu5FpIHPhu5FuZywgdOG7qyDEkcOzIOG6o25oIGjGsOG7n25nIMSR4bq/biBj4bqjbSBnacOhYyBs4bqhbmgga2hpIG5n4bunLiBDw6FjIGvhur90IHF14bqjIG7DoHkgxJHhurd0IHJhIGdp4bqjIHRoaeG6v3QgcuG6sW5nIG5o4buvbmcgbmfGsOG7nWkgbMOgbSB2aeG7h2MgdHJvbmcgbmfDoG5oIGPDtG5nIG5naOG7hyBjw7MgdGjhu4MgdGnhur9wIHjDumMgbmhp4buBdSB24bubaSBtw7RpIHRyxrDhu51uZyDEkWnhu4F1IGjDsmEsIGzDoG0gdmnhu4djIGtodXlhIGhv4bq3YyB0aGnhur91IHbhuq1uIMSR4buZbmcsIGThuqtuIMSR4bq/biBk4buFIGLhu4sgbOG6oW5oIGjGoW4ga2hpIG5n4bunIHbDoG8gYmFuIMSRw6ptLg0KDQpgYGB7cn0NCmNoaXNxLnRlc3QodGFuc3VhdDQpDQpgYGANCg0KROG7sWEgdHLDqm4ga+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmggQ2hpLWLDrG5oIHBoxrDGoW5nIGdp4buvYSBoYWkgYmnhur9uICoqQ29sZCoqIHbDoCAqKk1ham9yKiosIGdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQgxJHhuqF0IDU0LjUwMiB24bubaSA0MCBi4bqtYyB04buxIGRvIHbDoCBnacOhIHRy4buLIHAgdMawxqFuZyDhu6luZyBsw6AgMC4wLjA2Mjg1LiBWw6wgcC12YWx1ZSA+IDAuMDUsIG7Dqm4ga2jDtG5nIMSR4bunIGLhurFuZyBjaOG7qW5nIMSR4buDIGLDoWMgYuG7jyBnaeG6oyB0aHV54bq/dCAkSF8wJCwgdOG7qWMgbMOgICoqS0jDlE5HIGPDsyBt4buRaSBsacOqbiBo4buHIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6oqKiBnaeG7r2EgYuG7iyBs4bqhbmggdHJvbmcga2hpIG5n4bunIHbDoCBuZ8Ogbmggbmdo4buBLiBOw7NpIGPDoWNoIGtow6FjLCBoYWkgYmnhur9uIG7DoHkgY8OzIHh1IGjGsOG7m25nIMSR4buZYyBs4bqtcCB24bubaSBuaGF1IHRyb25nIGLhu5kgZOG7ryBsaeG7h3UuIFR1eSBuaGnDqm4sIGPhuqduIGzGsHUgw70gcuG6sW5nIGtp4buDbSDEkeG7i25oIGPDsyBj4bqjbmggYsOhbyB24buBIMSR4buZIGNow61uaCB4w6FjIGRvIG3hu5l0IHPhu5Egw7QgdHJvbmcgYuG6o25nIHThuqduIHN14bqldCBjw7MgZ2nDoSB0cuG7iyBr4buzIHbhu41uZyBuaOG7jywgxJFp4buBdSBuw6B5IGPDsyB0aOG7gyDhuqNuaCBoxrDhu59uZyDEkeG6v24gxJHhu5kgdGluIGPhuq15IGPhu6dhIGvhur90IHF14bqjLg0KDQojIyMgKiozLjEuNS4gVMOhYyDEkeG7mW5nIGPhu6dhIHnhur91IHThu5EgY+G6o20gdGjhuqV5IG7Ds25nIGtoaSBuZ+G7pyB0cm9uZyBjw6FjIG5nw6BuaCBuZ2jhu4EqKg0KDQpgYGB7cn0NCiMgQuG6o25nIHThuqduIHN14bqldCBnaeG7r2EgSG90dsOgIE1ham9yDQp0YW5zdWF0NSA8LSB0YWJsZShkYXRhJGBIb3RgLCBkYXRhJE1ham9yKSANCnBoYW50cmFtNSA8LSBwcm9wLnRhYmxlKHRhbnN1YXQ1LCBtYXJnaW4gPSAxKSAqIDEwMA0KDQpiYW5nNSA8LSBhcy5kYXRhLmZyYW1lKHRhbnN1YXQ1KQ0KYmFuZzUkUGVyY2VudGFnZSA8LSByb3VuZChhcy52ZWN0b3IocGhhbnRyYW01KSwgMikNCg0KY29sbmFtZXMoYmFuZzUpIDwtIGMoIkhvdCIsICJNYWpvciIsICJGcmVxdWVuY3kiLCAiUGVyY2VudGFnZSIpDQogICAgDQprYWJsZShiYW5nNSwgY2FwdGlvbiA9ICJC4bqjbmcgNDogTcO0IHThuqMgdGjhu5FuZyBrw6ogY+G7p2EgMiBiaeG6v24gSG90IHbDoCBNYWpvciIpDQoNCiNUcnVjIHF1YW4gaG9hDQpiYW5nNSA8LSBkYXRhICU+JQ0KICBncm91cF9ieShNYWpvciwgYEhvdGApICU+JQ0KICBzdW1tYXJpc2UoQ291bnQgPSBuKCksIC5ncm91cHMgPSAiZHJvcCIpDQoNCmdncGxvdChiYW5nNSwgYWVzKHggPSBNYWpvciwgeSA9IENvdW50LCBmaWxsID0gYEhvdGApKSArDQogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsNCiAgbGFicyh0aXRsZSA9ICJIw6xuaCA1OiBN4buRaSBxdWFuIGjhu4cgZ2nhu69hIHnhur91IHThu5EgY+G6o20gdGjhuqV5IG7Ds25nIGtoaSBuZ+G7pyB2w6AgbmfDoG5oIG5naOG7gSIsDQogICAgICAgeCA9ICJOZ2FuaCBuZ2hlIiwgeSA9ICJj4bqjbSB0aOG6pXkgbOG6oW5oIGtoaSBuZ+G7pyIsIGZpbGwgPSAiVOG6p24gc3XhuqV0IikgKw0KICB0aGVtZV9taW5pbWFsKCkgKw0KICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCBoanVzdCA9IDEpKQ0KYGBgDQoNCkThu7FhIHRyw6puIGLhuqNuZyB04bqnbiBzdeG6pXQgY2jDqW8sIGPDsyB0aOG7gyB0aOG6pXkgbeG7qWMgxJHhu5kgZ+G6t3AgcGjhuqNpIHTDrG5oIHRy4bqhbmcgdGnhu4N1IMSRw6ptIGhv4bq3YyBzw6FuZyBz4bubbSBjw7Mgc+G7sSBraMOhYyBiaeG7h3QgxJHDoW5nIGvhu4MgZ2nhu69hIGPDoWMgbmjDs20gbmfDoG5oIG5naOG7gS4gTmjDs20gbMOgbSB2aeG7h2MgdHJvbmcgbMSpbmggduG7sWMgZ2nDoW8gZOG7pWMgY8OzIHThu7cgbOG7hyAqKmtow7RuZyBn4bq3cCoqIHbhuqVuIMSR4buBIG7DoHkgY2FvIG5o4bqldCAoNjgsOCUpLCB0cm9uZyBraGkgbmjDs20gbMOgbSB0cm9uZyBsxKluaCB24buxYyBjw7RuZyBuZ2hp4buHcCDigJMgeMOieSBk4buxbmcgbOG6oWkgY8OzIHThu7cgbOG7hyAqKmfhurdwIHBo4bqjaSoqIGNhbyBuaOG6pXQgKDU0LDUlKS4gxJBp4buBdSBuw6B5IGNobyB0aOG6pXkgbmfGsOG7nWkgbMOgbSB0cm9uZyBsxKluaCB24buxYyBjw7RuZyBuZ2hp4buHcCDigJMgeMOieSBk4buxbmcgY8OzIGto4bqjIG7Eg25nIGfhurdwIHbhuqVuIMSR4buBIHbhu4EgZ2nhuqVjIG5n4bunIChj4bulIHRo4buDIGzDoCB0aeG7g3UgxJHDqm0gaG/hurdjIHRp4buDdSByw6FuZyBzw6FuZykgKipjYW8gaMahbioqIMSRw6FuZyBr4buDIHNvIHbhu5tpIG5o4buvbmcgbmfGsOG7nWkgbMOgbSB0cm9uZyBjw6FjIG5nw6BuaCBraMOhYy4gxJDDonkgY8OzIHRo4buDIGzDoCBkbyB0w61uaCBjaOG6pXQgY8O0bmcgdmnhu4djIG7hurduZyBuaOG7jWMsIHRo4budaSBnaWFuIGzDoG0gdmnhu4djIGtow7RuZyBj4buRIMSR4buLbmggaG/hurdjIG3hu6ljIMSR4buZIGPEg25nIHRo4bqzbmcgY2FvIHRyb25nIG5nw6BuaCBjw7RuZyBuZ2hp4buHcCDigJMgeMOieSBk4buxbmcuDQoNCg0KYGBge3J9DQpjaGlzcS50ZXN0KHRhbnN1YXQ1KQ0KYGBgDQoNCkThu7FhIHRyw6puIGvhur90IHF14bqjIGtp4buDbSDEkeG7i25oIENoaS1iw6xuaCBwaMawxqFuZyBnaeG7r2EgaGFpIGJp4bq/biAqKkhvdCoqIHbDoCAqKk1ham9yKiosIGdpw6EgdHLhu4sgdGjhu5FuZyBrw6ogQ2hpLXNxdWFyZWQgxJHhuqF0IDEwMi45NSB24bubaSA2MCBi4bqtYyB04buxIGRvIHbDoCBnacOhIHRy4buLIHAgdMawxqFuZyDhu6luZyBsw6AgMC4wMDA0NzIgVsOsIHAtdmFsdWUgPCAwLjA1LCBuw6puIGtow7RuZyDEkeG7pyBi4bqxbmcgY2jhu6luZyDEkeG7gyBiw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3QgJEhfMSQsIHThu6ljIGzDoCAqKkPDkyBt4buRaSBsacOqbiBo4buHIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6oqKiBnaeG7r2EgeeG6v3UgdOG7kSBj4bqjbSB0aOG6pXkgbsOzbmcgbMO6YyBuZ+G7pyB2w6AgbmfDoG5oIG5naOG7gS4gTsOzaSBjw6FjaCBraMOhYywgaGFpIGJp4bq/biBuw6B5IGPDsyB4dSBoxrDhu5tuZyBwaOG7pSB0aHXhu5ljIHbhu5tpIG5oYXUgdHJvbmcgYuG7mSBk4buvIGxp4buHdS4gVHV5IG5oacOqbiwgY+G6p24gbMawdSDDvSBy4bqxbmcga2nhu4NtIMSR4buLbmggY8OzIGPhuqNuaCBiw6FvIHbhu4EgxJHhu5kgY2jDrW5oIHjDoWMgZG8gbeG7mXQgc+G7kSDDtCB0cm9uZyBi4bqjbmcgdOG6p24gc3XhuqV0IGPDsyBnacOhIHRy4buLIGvhu7MgduG7jW5nIG5o4buPLCDEkWnhu4F1IG7DoHkgY8OzIHRo4buDIOG6o25oIGjGsOG7n25nIMSR4bq/biDEkeG7mSB0aW4gY+G6rXkgY+G7p2Ega+G6v3QgcXXhuqMuDQoNCmBgYHtyfQ0Kcmlza3JhdGlvKHRhbnN1YXQ1KQ0KYGBgDQoNCkvhur90IHF14bqjIHBow6JuIHTDrWNoIGNobyB0aOG6pXkgbeG7kWkgbGnDqm4gaOG7hyBnaeG7r2EgbmfDoG5oIG5naOG7gSB2w6AgdOG6p24gc3XhuqV0IGPhuqNtIHRo4bqleSBuw7NuZyBraGkgbmfhu6cgbMOgIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6og4bufIG3hu5l0IHPhu5EgbeG7qWMgY+G7pSB0aOG7gy4gQ+G7pSB0aOG7gywgbmjDs20ga2jDtG5nIGPhuqNtIHRo4bqleSBuw7NuZyBs4bqnbiBuw6BvIHRyb25nIHRow6FuZyAoMC9tb250aCkgY8OzIG5ndXkgY8ahIHRo4bqlcCBoxqFuIHLDtSBy4buHdCBzbyB24bubaSBuaMOzbSBj4bqjbSB0aOG6pXkgbsOzbmcgaMahbiAzIGzhuqduIG3hu5dpIHRow6FuZyAocmlzayByYXRpbyA9IDAuMDc0LCBwLXZhbHVlIEZpc2hlciA9IDAuMDAwNCkuIFTGsMahbmcgdOG7sSwgbmjDs20gY+G6o20gdGjhuqV5IG7Ds25nIGhhaSBs4bqnbiBt4buXaSB0deG6p24gKDIvd2VlaykgY8WpbmcgY8OzIHPhu7Ega2jDoWMgYmnhu4d0IMSRw6FuZyBr4buDIChwID0gMC4wMDE0KS4gTmjhu69uZyBr4bq/dCBxdeG6oyBuw6B5IGfhu6NpIMO9IHLhurFuZyB04bqnbiBzdeG6pXQgY+G6o20gdGjhuqV5IG7Ds25nIGtoaSBuZ+G7pyBjw7MgdGjhu4Mga2jDoWMgbmhhdSB0w7l5IHRoZW8gbmfDoG5oIG5naOG7gSwgdHV5IG5oacOqbiBt4buZdCBz4buRIG5ow7NtIGtow6FjIGzhuqFpIGtow7RuZyBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqIHLDtSByw6BuZyAocCA+IDAuMDUpLiDEkOG7k25nIHRo4budaSwgZG8gY+G6o25oIGLDoW8g4oCcQ2hpLXNxdWFyZWQgYXBwcm94aW1hdGlvbiBtYXkgYmUgaW5jb3JyZWN04oCdIHh14bqldCBoaeG7h24gbmhp4buBdSBs4bqnbiwgdGEgbsOqbiDGsHUgdGnDqm4gxJHhu41jIGvhur90IHF14bqjIHRoZW8gRmlzaGVyIGV4YWN0IHRlc3QgxJHhu4MgxJHhuqNtIGLhuqNvIMSR4buZIHRpbiBj4bqteS4NCg0KIyMjICoqMy4xLjYuIFTDoWMgxJHhu5luZyBj4bunYSB54bq/dSB04buRIGfhurdwIMOhYyBt4buZbmcga2hpIG5n4bunIHRyb25nIGPDoWMgbmfDoG5oIG5naOG7gSoqDQoNCmBgYHtyfQ0KIyBC4bqjbmcgdOG6p24gc3XhuqV0IGdp4buvYSBEcmVhbXMgdsOgIE1ham9yDQp0YW5zdWF0NiA8LSB0YWJsZShkYXRhJGBEcmVhbXNgLCBkYXRhJE1ham9yKSANCnBoYW50cmFtNiA8LSBwcm9wLnRhYmxlKHRhbnN1YXQ2LCBtYXJnaW4gPSAxKSAqIDEwMA0KDQpiYW5nNiA8LSBhcy5kYXRhLmZyYW1lKHRhbnN1YXQ2KQ0KYmFuZzYkUGVyY2VudGFnZSA8LSByb3VuZChhcy52ZWN0b3IocGhhbnRyYW02KSwgMikNCg0KY29sbmFtZXMoYmFuZzUpIDwtIGMoIkRyZWFtcyIsICJNYWpvciIsICJGcmVxdWVuY3kiLCAiUGVyY2VudGFnZSIpDQogICAgDQprYWJsZShiYW5nNiwgY2FwdGlvbiA9ICJC4bqjbmcgNjogTcO0IHThuqMgdGjhu5FuZyBrw6ogY+G7p2EgMiBiaeG6v24gRHJlYW1zIHbDoCBNYWpvciIpDQoNCiNUcnVjIHF1YW4gaG9hDQpiYW5nNiA8LSBkYXRhICU+JQ0KICBncm91cF9ieShNYWpvciwgYEhvdGApICU+JQ0KICBzdW1tYXJpc2UoQ291bnQgPSBuKCksIC5ncm91cHMgPSAiZHJvcCIpDQoNCmdncGxvdChiYW5nNiwgYWVzKHggPSBNYWpvciwgeSA9IENvdW50LCBmaWxsID0gYEhvdGApKSArDQogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsNCiAgbGFicyh0aXRsZSA9ICJIw6xuaCA2OiBN4buRaSBxdWFuIGjhu4cgZ2nhu69hIHnhur91IHThu5EgZ+G6t3Agw6FjIG3hu5luZyBraGkgbmfhu6cgdsOgIG5nw6BuaCBuZ2jhu4EiLA0KICAgICAgIHggPSAiTmdhbmggbmdoZSIsIHkgPSAiR+G6t3Agw6FjIG3hu5luZyBraGkgbmfhu6ciLCBmaWxsID0gIlThuqduIHN14bqldCIpICsNCiAgdGhlbWVfbWluaW1hbCgpICsNCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgaGp1c3QgPSAxKSkNCmBgYA0KDQpL4bq/dCBxdeG6oyBjaG8gdGjhuqV5IG3hu5FpIGxpw6puIGjhu4cgcsO1IHLhu4d0IGdp4buvYSBuZ8Ogbmggbmdo4buBIHbDoCB04bqnbiBzdeG6pXQgZ+G6t3Agw6FjIG3hu5luZyBraGkgbmfhu6cuIFRyb25nIMSRw7MsbmfDoG5oIEJpb3RlY2ggdsOgIENvbXB1dGVyIFNjaWVuY2UgY8OzIHThu7cgbOG7hyBn4bq3cCDDoWMgbeG7mW5nIGNhbyBoxqFuIMSRw6FuZyBr4buDIHNvIHbhu5tpIGPDoWMgbmfDoG5oIGtow6FjLiBD4bulIHRo4buDLCBuaMOzbSBuZ8OgbmggQmlvdGVjaCBjw7MgdOG7tyBs4buHIGfhurdwIMOhYyBt4buZbmcgdOG7qyAx4oCTMiBs4bqnbi90deG6p24gdsOgIDIgbOG6p24vdHXhuqduIGNoaeG6v20gZ+G6p24gbeG7mXQgbuG7rWEgdOG7lW5nIHPhu5EsIHRyb25nIGtoaSBuZ8OgbmggQ29tcHV0ZXIgU2NpZW5jZSBjw7MgdOG7tyBs4buHIMOhYyBt4buZbmcgbeG7l2kgdHXhuqduIHLhuqV0IGNhbyAoMS93ZWVrOiA0Ni41NSUpLiBOZ8aw4bujYyBs4bqhaSwgY8OhYyBuZ8OgbmggbmjGsCBNZWNoYW5pY2FsLCBDaGVtaWNhbCB2w6AgQ2l2aWwgbOG6oWkgY8OzIHThu7cgbOG7hyBn4bq3cCDDoWMgbeG7mW5nIHRo4bqlcCBoxqFuLCBwaOG6p24gbOG7m24gcsahaSB2w6BvIG5ow7NtIGtow7RuZyBn4bq3cCDDoWMgbeG7mW5nIHRyb25nIHRow6FuZy4gxJBp4buBdSBuw6B5IGfhu6NpIMO9IHLhurFuZyBuZ8Ogbmggbmdo4buBIGPDsyB0aOG7gyDhuqNuaCBoxrDhu59uZyDEkeG6v24gY2jhuqV0IGzGsOG7o25nIGdp4bqlYyBuZ+G7pywgdHJvbmcgxJHDsyBjw6FjIG5nw6BuaCBjw7Mgw6FwIGzhu7FjIGPDtG5nIHZp4buHYyBob+G6t2MgdMOtbmggY2jhuqV0IGPDtG5nIHZp4buHYyBjYW8gY8OzIHRo4buDIGzDoG0gdMSDbmcgbmd1eSBjxqEgZ+G6t3Agw6FjIG3hu5luZy4NCg0KDQpgYGB7cn0NCmNoaXNxLnRlc3QodGFuc3VhdDYpDQpgYGANCg0KROG7sWEgdHLDqm4ga+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmggQ2hpLWLDrG5oIHBoxrDGoW5nIGdp4buvYSBoYWkgYmnhur9uICoqRHJlYW1zKiogdsOgICoqTWFqb3IqKiwgZ2nDoSB0cuG7iyB0aOG7kW5nIGvDqiBDaGktc3F1YXJlZCDEkeG6oXQgMTE3LjQxIHbhu5tpIDYwIGLhuq1jIHThu7EgZG8gdsOgIGdpw6EgdHLhu4sgcCB0xrDGoW5nIOG7qW5nIGzDoCAxLjM1N2UtMDUgVsOsIHAtdmFsdWUgPCAwLjA1LCBuw6puIGtow7RuZyDEkeG7pyBi4bqxbmcgY2jhu6luZyDEkeG7gyBiw6FjIGLhu48gZ2nhuqMgdGh1eeG6v3QgJEhfMSQsIHThu6ljIGzDoCAqKkPDkyBt4buRaSBsacOqbiBo4buHIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6oqKiBnaeG7r2EgeeG6v3UgdOG7kSBn4bq3cCDDoWMgbeG7mW5nIGzDumMgbmfhu6cgdsOgIG5nw6BuaCBuZ2jhu4EuIE7Ds2kgY8OhY2gga2jDoWMsIGhhaSBiaeG6v24gbsOgeSBjw7MgeHUgaMaw4bubbmcgcGjhu6UgdGh14buZYyB24bubaSBuaGF1IHRyb25nIGLhu5kgZOG7ryBsaeG7h3UuIFR1eSBuaGnDqm4sIGPhuqduIGzGsHUgw70gcuG6sW5nIGtp4buDbSDEkeG7i25oIGPDsyBj4bqjbmggYsOhbyB24buBIMSR4buZIGNow61uaCB4w6FjIGRvIG3hu5l0IHPhu5Egw7QgdHJvbmcgYuG6o25nIHThuqduIHN14bqldCBjw7MgZ2nDoSB0cuG7iyBr4buzIHbhu41uZyBuaOG7jywgxJFp4buBdSBuw6B5IGPDsyB0aOG7gyDhuqNuaCBoxrDhu59uZyDEkeG6v24gxJHhu5kgdGluIGPhuq15IGPhu6dhIGvhur90IHF14bqjLg0KDQoNCmBgYHtyfQ0Kcmlza3JhdGlvKHRhbnN1YXQ2KQ0KYGBgDQoNCkvhur90IHF14bqjIHBow6JuIHTDrWNoIGNobyB0aOG6pXkgY8OzIG3hu5FpIGxpw6puIGjhu4cgxJHDoW5nIGvhu4MgZ2nhu69hIG5nw6BuaCBuZ2jhu4EgdsOgIHThuqduIHN14bqldCBn4bq3cCDDoWMgbeG7mW5nIGtoaSBuZ+G7py4gQ+G7pSB0aOG7gywgbmjDs20gbmfGsOG7nWkga2hhaSBiw6FvIGtow7RuZyBn4bq3cCDDoWMgbeG7mW5nICgwIGzhuqduL3Row6FuZykgdsOgIG5ow7NtIGfhurdwIMOhYyBt4buZbmcga2hv4bqjbmcgMSBs4bqnbi90deG6p24gY8OzIHPhu7Ega2jDoWMgYmnhu4d0IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogc28gduG7m2kgbmjDs20gdGjGsOG7nW5nIHh1ecOqbiBn4bq3cCDDoWMgbeG7mW5nICg+MyBs4bqnbi90aMOhbmcpLCB24bubaSBwLXZhbHVlIGzhuqduIGzGsOG7o3QgbMOgIDAuMDAwMDM2IHbDoCAwLjAwMDQ2LiBUdXkgbmhpw6puLCBkbyBt4buZdCBz4buRIMO0IHRyb25nIGLhuqNuZyB04bqnbiBzdeG6pXQgY8OzIGdpw6EgdHLhu4sgcXXDoSBuaOG7jyBob+G6t2MgYuG6sW5nIDAsIHThu7cgbOG7hyBuZ3V5IGPGoSAoUlIpIGPhu6dhIGPDoWMgbmjDs20gY8OybiBs4bqhaSBraMO0bmcgxJHGsOG7o2MgxrDhu5tjIGzGsOG7o25nIHbDoCBoaeG7g24gdGjhu4sgxJHhuqd5IMSR4bunIChi4buLIHRy4bqjIHbhu4EgTmFOKSwgbsOqbiBraMO0bmcgdGjhu4Mga+G6v3QgbHXhuq1uIGNow61uaCB4w6FjIG3hu6ljIMSR4buZIGNow6puaCBs4buHY2ggbmd1eSBjxqEgZ2nhu69hIGPDoWMgbmjDs20uIETDuSB24bqteSwga+G6v3QgcXXhuqMgduG6q24gY2hvIHRo4bqleSBuZ8Ogbmggbmdo4buBIGPDsyB0aOG7gyBsw6AgbeG7mXQgeeG6v3UgdOG7kSBsacOqbiBxdWFuIMSR4bq/biB04bqnbiBzdeG6pXQgZ+G6t3Agw6FjIG3hu5luZyBj4bunYSBuZ8aw4budaSB0aGFtIGdpYSBraOG6o28gc8OhdC4NCl9fXw0KDQojIyAqKjMuMi4gUGjDom4gdMOtY2ggaOG7k2kgcXV5IMSRYSBiaeG6v24gdsOgIGzhu7FhIGNo4buNbiBtw7QgaMOsbmggdOG7kWkgxrB1KioNCg0KQ8OhYyBiaeG6v24gxJHhu5ljIGzhuq1wIGPDsyB0w6FjIMSR4buZbmcgxJHhur9uIGJp4bq/biBuZ2jhu4EgbmdoaeG7h3AgKE1ham9yKSBiYW8gZ+G7k206IFdha2UgdXAgbWlkZGxlLCBCYXRocm9vbSwgSG90IHbDoCBEcmVhbXMuDQoNCioqTcO0IGjDrG5oIGjhu5NpIHF1eSBNdWx0aW5vbWlhbCBsb2dpc3RpYyByZWdyZXNzaW9uKioNCg0KYGBge3J9DQpsaWJyYXJ5KG5uZXQpDQpkYXRhJE1ham9yIDwtIGZhY3RvcihkYXRhJE1ham9yKSAgIyBiaeG6v24gxJHhu4tuaCB0w61uaCBuaGnhu4F1IG5ow7NtDQptb2RlbF9tdWx0aW5vbSA8LSBtdWx0aW5vbShNYWpvciB+IGBXYWtlIHVwIG1pZGRsZWAgKyBCYXRocm9vbSArIEhvdCArIERyZWFtcywgZGF0YSA9IGRhdGEpDQpgYGANCg0KS+G6v3QgcXXhuqMgdHLDqm4gdGjhu4MgaGnhu4duIHF1w6EgdHLDrG5oIGjhu5lpIHThu6UgY+G7p2EgbcO0IGjDrG5oIGjhu5NpIHF1eSBQcm9iaXQgduG7m2kgMjc1IHF1YW4gc8OhdCB2w6AgMjQwIGJp4bq/biDEkeG7mWMgbOG6rXAgc2F1IGtoaSBiaeG6v24gcGjDom4gbG/huqFpIMSRxrDhu6NjIG3DoyBow7NhLiBHacOhIHRy4buLIGjDoG0gbG9nLWxpa2VsaWhvb2QgYmFuIMSR4bqndSBsw6AgMzkwLjg1NjkyOSB2w6AgZ2nhuqNtIGThuqduIHF1YSB04burbmcgdsOybmcgbOG6t3AsIGNobyB0aOG6pXkgdGh14bqtdCB0b8OhbiB04buRaSDGsHUgaMOzYSDEkWFuZyBk4bqnbiB0aeG7h20gY+G6rW4gbmdoaeG7h20gdOG7kXQuIFNhdSA5MCB2w7JuZyBs4bq3cCwgbcO0IGjDrG5oIMSR4bqhdCBnacOhIHRy4buLIGxvZy1saWtlbGlob29kIGN14buRaSBjw7luZyBsw6AgMTczLjg4NzEwNSB2w6AgxJHDoyBo4buZaSB04bulIHRow6BuaCBjw7RuZy4gxJBp4buBdSBuw6B5IGNo4bupbmcgdOG7jyBtw7QgaMOsbmggxJHDoyDGsOG7m2MgbMaw4bujbmcg4buVbiDEkeG7i25oIGPDoWMgaOG7hyBz4buRLCBz4bq1biBzw6BuZyDEkeG7gyB0aeG6v24gaMOgbmggcGjDom4gdMOtY2ggw70gbmdoxKlhIHRo4buRbmcga8OqIGPhu6dhIGPDoWMgYmnhur9uIHbDoCDEkcOhbmggZ2nDoSDEkeG7mSBwaMO5IGjhu6NwIGPhu6dhIG3DtCBow6xuaC4NCg0KKipNw7QgaMOsbmggaOG7k2kgcXV5IENvbXBsZW1lbnRhcnkgTG9nLUxvZyAoQ2xvZ2xvZykqKg0KYGBge3J9DQptb2RlbF9sb2cgPC0gZ2xtKE1ham9yIH4gYFdha2UgdXAgbWlkZGxlYCArIA0KICAgICAgICAgICAgICAgICAgICAgICAgIGBCYXRocm9vbWAgKw0KICAgICAgICAgICAgICAgICAgICAgICAgIGBIb3RgICsgDQogICAgICAgICAgICAgICAgICAgICAgICAgYERyZWFtc2AsDQogICAgICAgICAgICAgICAgICAgIGRhdGEgPSBkYXRhLCBmYW1pbHkgPSBiaW5vbWlhbChsaW5rID0gImNsb2dsb2ciKSkNCnN1bW1hcnkobW9kZWxfbG9nKQ0KYGBgDQoNCk3DtCBow6xuaCBo4buTaSBxdXkgbmjhu4sgcGjDom4gduG7m2kgbGnDqm4ga+G6v3QgY29tcGxlbWVudGFyeSBsb2ctbG9nIChjbG9nbG9nKSDEkcaw4bujYyDGsOG7m2MgbMaw4bujbmcgbmjhurFtIHBow6JuIHTDrWNoIGPDoWMgeeG6v3UgdOG7kSDhuqNuaCBoxrDhu59uZyDEkeG6v24ga2jhuqMgbsSDbmcgeOG6o3kgcmEgYmnhur9uIHBo4bulIHRodeG7mWMgxJHhu4tuaCB0w61uaCAqTWFqb3IqLiBUdXkgbmhpw6puLCBr4bq/dCBxdeG6oyBjaG8gdGjhuqV5IG3DtCBow6xuaCBn4bq3cCB24bqlbiDEkeG7gSBuZ2hpw6ptIHRy4buNbmcgduG7gSDEkeG7mSDhu5VuIMSR4buLbmggdsOgIMO9IG5naMSpYSB0aOG7kW5nIGvDqi4gQ+G7pSB0aOG7gywgY+G6o25oIGLDoW8gYCJmaXR0ZWQgcHJvYmFiaWxpdGllcyBudW1lcmljYWxseSAwIG9yIDEgb2NjdXJyZWQiYCBjaG8gdGjhuqV5IGhp4buHbiB0xrDhu6NuZyBwaMOibiB0w6FjaCBob8OgbiB0b8OgbiAoY29tcGxldGUgc2VwYXJhdGlvbiksIGtoaSBt4buZdCBz4buRIHThu5UgaOG7o3AgYmnhur9uIGdp4bqjaSB0aMOtY2ggZOG7sSDEkW/DoW4gaG/DoG4gdG/DoG4gxJHDum5nIGJp4bq/biBwaOG7pSB0aHXhu5ljLCBk4bqrbiDEkeG6v24gc2FpIHPhu5EgY2h14bqpbiBy4bqldCBs4bubbi4gSOG6p3UgaOG6v3QgY8OhYyBo4buHIHPhu5EgxJHhu4F1IGPDsyBnacOhIHRy4buLIHAgY2FvIChQcig+fHp8KSA+IDAuOSksIGNobyB0aOG6pXkgY2jDum5nIGtow7RuZyBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqIHRyb25nIG3DtCBow6xuaC4gTeG6t2MgZMO5IMSR4buZIGzhu4djaCBkxrAgKHJlc2lkdWFsIGRldmlhbmNlKSBnaeG6o20gc28gduG7m2kgxJHhu5kgbOG7h2NoIGJhbiDEkeG6p3UgKG51bGwgZGV2aWFuY2UpLCDEkWnhu4F1IG7DoHkga2jDtG5nIMSR4bunIMSR4buDIGto4bqzbmcgxJHhu4tuaCBtw7QgaMOsbmggcGjDuSBo4bujcCB04buRdCwgxJHhurdjIGJp4buHdCB0cm9uZyBi4buRaSBj4bqjbmggc2FpIHPhu5EgY2h14bqpbiBwaMOzbmcgxJHhuqFpLiBEbyDEkcOzLCBr4bq/dCBxdeG6oyBtw7QgaMOsbmggbsOgeSBraMO0bmcgxJHDoW5nIHRpbiBj4bqteSDEkeG7gyBkaeG7hW4gZ2nhuqNpIHbDoCBj4bqnbiDEkcaw4bujYyB4ZW0geMOpdCBs4bqhaSBi4bqxbmcgY8OhY2ggeOG7rSBsw70gduG6pW4gxJHhu4EgcGjDom4gdMOhY2ggZOG7ryBsaeG7h3UsIGxv4bqhaSBi4buPIGhv4bq3YyBn4buZcCBiaeG6v24sIGhv4bq3YyBz4butIGThu6VuZyBwaMawxqFuZyBwaMOhcCDGsOG7m2MgbMaw4bujbmcgdGhheSB0aOG6vyBuaMawIEZpcnRoIGxvZ2lzdGljIHJlZ3Jlc3Npb24uDQoNCioqTcO0IGjDrG5oIGjhu5NpIHF1eSBMb2dpc3RpYyAoTG9naXQpKioNCg0KYGBge3J9DQptb2RlbF9sb2dpdCA8LSBnbG0oTWFqb3IgfiBgV2FrZSB1cCBtaWRkbGVgICsgDQogICAgICAgICAgICAgICAgICAgICAgICAgICBgQmF0aHJvb21gICsgDQogICAgICAgICAgICAgICAgICAgICAgICAgICBgSG90YCArDQogICAgICAgICAgICAgICAgICAgICAgICAgICBgRHJlYW1zYCwNCiAgICAgICAgICAgICAgZGF0YSA9IGRhdGEsIGZhbWlseSA9IGJpbm9taWFsKGxpbmsgPSAibG9naXQiKSkNCnN1bW1hcnkobW9kZWxfbG9naXQpDQpgYGANCg0KTcO0IGjDrG5oIGjhu5NpIHF1eSBsb2dpc3RpYyDEkcaw4bujYyB4w6J5IGThu7FuZyBuaOG6sW0gxJHDoW5oIGdpw6EgbeG7kWkgbGnDqm4gaOG7hyBnaeG7r2EgY8OhYyB54bq/dSB04buRIGdp4bqlYyBuZ+G7pyAoKldha2UgdXAgbWlkZGxlKiwgKkJhdGhyb29tKiwgKkhvdCosICpEcmVhbXMqKSB2w6AgYmnhur9uIHBo4bulIHRodeG7mWMgxJHhu4tuaCB0w61uaCAqTWFqb3IqIChuZ8Ogbmggbmdo4buBKS4gVHV5IG5oacOqbiwga+G6v3QgcXXhuqMgbcO0IGjDrG5oIGNobyB0aOG6pXkgaGnhu4duIHTGsOG7o25nICoqcGjDom4gdMOhY2ggaG/DoG4gdG/DoG4gKGNvbXBsZXRlIHNlcGFyYXRpb24pKiogaG/hurdjIGfhuqduIG5oxrAgaG/DoG4gdG/DoG4sIMSRxrDhu6NjIHRo4buDIGhp4buHbiBxdWEgY+G6o25oIGLDoW8gYCJmaXR0ZWQgcHJvYmFiaWxpdGllcyBudW1lcmljYWxseSAwIG9yIDEgb2NjdXJyZWQiYCB2w6AgZ2nDoSB0cuG7iyBo4buHIHPhu5EgxrDhu5tjIGzGsOG7o25nIGPDuW5nIHNhaSBz4buRIGNodeG6qW4gcuG6pXQgbOG7m24gKG5oaeG7gXUgYmnhur9uIGPDsyBTdGQuIEVycm9yID4gNywwMDApLiBOZ2/DoGkgcmEsIGjhuqd1IGjhur90IGdpw6EgdHLhu4sgcCDEkeG7gXUgeOG6pXAgeOG7iSAxLjAwMCwgY2hvIHRo4bqleSAqKmtow7RuZyBjw7MgYmnhur9uIG7DoG8gxJHhuqF0IMO9IG5naMSpYSB0aOG7kW5nIGvDqioqIHRyb25nIG3DtCBow6xuaC4gTeG6t2MgZMO5IMSR4buZIGzhu4djaCBjw7JuIGzhuqFpIChyZXNpZHVhbCBkZXZpYW5jZSA9IDE0NC42MikgZ2nhuqNtIHNvIHbhu5tpIMSR4buZIGzhu4djaCBiYW4gxJHhuqd1IChudWxsIGRldmlhbmNlID0gMTcyLjEzKSwgbeG7qWMgY+G6o2kgdGhp4buHbiBuw6B5IGzDoCBraMO0bmcgxJHDoW5nIGvhu4Mga2hpIHjDqXQgxJHhur9uIG3hu6ljIMSR4buZIGtow7RuZyDhu5VuIMSR4buLbmggY+G7p2EgY8OhYyBo4buHIHPhu5EuIE5o4buvbmcgZOG6pXUgaGnhu4d1IG7DoHkgY2hvIHRo4bqleSAqKm3DtCBow6xuaCBraMO0bmcgxJHDoW5nIHRpbiBj4bqteSDEkeG7gyBkaeG7hW4gZ2nhuqNpKiogdsOgIGPhuqduIMSRxrDhu6NjIGhp4buHdSBjaOG7iW5oIGLhurFuZyBjw6FjaCB44butIGzDvSBs4bqhaSBk4buvIGxp4buHdSwgcsO6dCBn4buNbiBiaeG6v24sIGhv4bq3YyDDoXAgZOG7pW5nIHBoxrDGoW5nIHBow6FwIGjhu5NpIHF1eSBwZW5hbGl6ZWQgbmjGsCAqRmlydGggbG9naXN0aWMgcmVncmVzc2lvbiouDQoNCioqTOG7sWEgY2jhu41uIG3DtCBow6xuaCB04buRaSDGsHUqKg0KDQrEkOG7gyBs4buxYSBjaOG7jW4gbcO0IGjDrG5oIGjhu5NpIHF1eSB04buRaSDGsHUgdHJvbmcgYmEgbcO0IGjDrG5oOiAqKkxvZ2lzdGljKiosICoqQ29tcGxlbWVudGFyeSBMb2ctbG9nIChjbG9nbG9nKSoqIHbDoCAqKk11bHRpbm9taWFsIGxvZ2lzdGljIHJlZ3Jlc3Npb24qKiwgdGEgZOG7sWEgdsOgbyBjw6FjIHRpw6p1IGNow60gc28gc8OhbmggxJHhu4tuaCBsxrDhu6NuZyBj4bulIHRo4buDIG5oxrAgKipBSUMqKiwgKipOdWxsIERldmlhbmNlKiosICoqUmVzaWR1YWwgRGV2aWFuY2UqKiwgdsOgIG3hu6ljIMO9IG5naMSpYSB0aOG7kW5nIGvDqiBj4bunYSBjw6FjIGjhu4cgc+G7kS4NCg0KfCBNw7QgaMOsbmggICAgICAgICAgICAgICAgICAgICAgICAgfCBBSUMgICAgfCBOdWxsIERldmlhbmNlIHwgUmVzaWR1YWwgRGV2aWFuY2UgfA0KfCAtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tIHwgLS0tLS0tIHwgLS0tLS0tLS0tLS0tLSB8IC0tLS0tLS0tLS0tLS0tLS0tIHwgDQp8IExvZ2lzdGljICAgICAgICAgICAgICAgICAgICAgICAgfCAxOTUuMjQgfCAyNDAuNTMgICAgICAgIHwgMTczLjg5ICAgICAgICAgICAgfCANCnwgQ29tcGxlbWVudGFyeSBsb2ctbG9nIChjbG9nbG9nKSB8IDE5OC40NyB8IDI0MC41MyAgICAgICAgfCAxNzYuODcgICAgICAgICAgICB8IA0KfCBNdWx0aW5vbWlhbCBsb2dpc3RpYyByZWdyZXNzaW9uIHwgMTgyLjMxIHwgMzkwLjg2ICAgICAgICB8IDE3My44OSAgICAgICAgICAgIHwgDQoNCg0KKiAqKkFJQyoqOiBNw7QgaMOsbmggKipNdWx0aW5vbWlhbCBsb2dpc3RpYyByZWdyZXNzaW9uKiogY8OzIEFJQyAqKnRo4bqlcCBuaOG6pXQgKDE4Mi4zMSkqKiwgcGjhuqNuIMOhbmggbcO0IGjDrG5oIG7DoHkgduG7q2EgY8OzIMSR4buZIHBow7kgaOG7o3AgdOG7kXQsIHbhu6thIHRyw6FuaCDEkcaw4bujYyB2aeG7h2MgcXXDoSBwaOG7qWMgdOG6oXAuDQoqICoqRGV2aWFuY2UgKMSR4buZIGzhu4djaCkqKjoNCg0KICAqIFThuqV0IGPhuqMgY8OhYyBtw7QgaMOsbmggxJHhu4F1IGdp4bqjbSDEkcOhbmcga+G7gyB04burICoqTnVsbCBEZXZpYW5jZSoqIChtw7QgaMOsbmgga2jDtG5nIGPDsyBiaeG6v24gxJHhu5ljIGzhuq1wKSB4deG7kW5nICoqUmVzaWR1YWwgRGV2aWFuY2UqKiAobcO0IGjDrG5oIGPDsyBiaeG6v24pLCB0aOG7gyBoaeG7h24gY8OhYyBiaeG6v24gxJHhu5ljIGzhuq1wIMSRw6MgZ2nhuqNpIHRow61jaCB04buRdCBwaMawxqFuZyBzYWkgdHJvbmcgYmnhur9uIHBo4bulIHRodeG7mWMuDQogICogVHJvbmcgxJHDsywgbcO0IGjDrG5oICoqTXVsdGlub21pYWwgbG9naXN0aWMgcmVncmVzc2lvbioqIGPFqW5nIMSR4bqhdCAqKlJlc2lkdWFsIERldmlhbmNlKiogdGjhuqVwICgxNzMuODkpLCB0xrDGoW5nIMSRxrDGoW5nIGhv4bq3YyB0aOG6pXAgaMahbiBzbyB24bubaSBjw6FjIG3DtCBow6xuaCBjw7JuIGzhuqFpLg0KKiAqKlbhu4EgxJHhu5kgdGluIGPhuq15IG3DtCBow6xuaCoqOiBDw6FjIGjhu4cgc+G7kSB0cm9uZyBtw7QgaMOsbmggTXVsdGlub21pYWwgxJHhu4F1ICoqaOG7mWkgdOG7pSDhu5VuIMSR4buLbmgqKiBzYXUgbmhp4buBdSB2w7JuZyBs4bq3cCB2w6AgcGjhuqduIGzhu5tuIGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogKHAtdmFsdWUgPCAwLjA1KS4NCg0KVuG7m2kgKipBSUMgdGjhuqVwIG5o4bqldCoqLCAqKlJlc2lkdWFsIERldmlhbmNlIGPhuqNpIHRoaeG7h24gbeG6oW5oKiosIHbDoCAqKnPhu7EgaOG7mWkgdOG7pSDhu5VuIMSR4buLbmggY+G7p2EgdGh14bqtdCB0b8OhbioqLCAqKm3DtCBow6xuaCBo4buTaSBxdXkgbG9naXN0aWMgxJFhIHRo4bupYyAobXVsdGlub21pYWwpKiogxJHGsOG7o2MgxJHDoW5oIGdpw6EgbMOgICoqdOG7kWkgxrB1KiogdHJvbmcgdmnhu4djIG3DtCBow6xuaCBow7NhIGJp4bq/biBwaOG7pSB0aHXhu5ljIMSR4buLbmggdMOtbmggdHJvbmcgbmdoacOqbiBj4bupdSBuw6B5Lg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQo=

PHÂN TÍCH CÁC YẾU TỐ ẢNH HƯỞNG ĐẾN GIẤC NGỦ CỦA CÁC NGÀNH NGHỀ KHÁC NHAU

Nguyễn Phạm Thúy An - 2221000290