GVHD: ThS. Trần Mạnh Tường

Nhóm sinh viên thực hiện:

1. Mai Thị Hoài Thương - 2121001546

2. Nguyễn Thị Ngọc Hà - 2121012474

PHÂN TÍCH NHỮNG YẾU TỐ ẢNH HƯỞNG ĐẾN MỨC ĐỘ TIÊU THỤ RƯỢU CỦA SINH VIÊN

Bộ dữ liệu “Student Alcohol Consumption EDA” từ Kaggle cung cấp thông tin chi tiết về việc sinh viên tiêu thụ rượu tại các trường trung học ở Bồ Đào Nha. Bộ dữ liệu này không chỉ tập trung vào hành vi tiêu thụ rượu mà còn bao gồm các thông tin về môi trường, gia đình và xã hội của sinh viên. Mục đích của nghiên cứu là để hiểu rõ hơn về mối liên hệ giữa các yếu tố này và thói quen tiêu thụ rượu của sinh viên trong 1 ngày. Dữ liệu được thu thập từ các cuộc khảo sát trong cộng đồng học đường, mang lại cơ hội để nghiên cứu tác động của môi trường và gia đình đối với hành vi tiêu thụ rượu ở độ tuổi sinh viên.

LỜI CẢM ƠN

Lời đầu tiên, chúng em xin gửi lời cảm ơn chân thành đến Giảng viên - ThS. Trần Mạnh Tường. Trong quá trình học tập và tìm hiểu môn “Phân tích dữ liệu định tính”, chúng em đã nhận được sự quan tâm giúp đỡ, hướng dẫn tận tình và cống hiến tâm huyết của thầy. Thầy đã giúp chúng em tích lũy thêm nhiều kiến thức để có cái nhìn hoàn thiện và tổng quan về bộ môn này. Từ những kiến thức được thầy truyền tải, chúng em dần hiểu được tầm quan trọng của môn học mà bấy lâu nay còn nhiều thắc mắc. Thông qua bài tiểu luận này, chúng em xin trình bày những gì chúng em đã tìm hiểu về môn học, và cụ thể là chúng em lựa chọn đề tài “Phân tích những yếu tố ảnh hưởng đến mức độ tiêu thụ rượu của sinh viên”

Kiến thức là vô hạn nhưng sự tiếp nhận kiến thức của bản thân mỗi chúng em luôn hạn chế một cách nhất định. Do đó, trong quá trình hoàn thành bài tiểu luận, chắc chắn chúng em không tránh khỏi những thiếu sót, vậy nên chúng em rất mong nhận được những đóng góp đến từ thầy để bài tiểu luận của nhóm được hoàn thiện hơn.

Chúng em kính chúc thầy sức khỏe, hạnh phúc và thành công trên con đường sự nghiệp của mình.

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

1.1 Lí do chọn đề tài

Hiện nay, hành vi tiêu thụ rượu là một hiện tượng phổ biến trong đời sống sinh viên. Nhiều sinh viên tham gia các hoạt động xã hội như tiệc tùng, dã ngoại và các sự kiến của trường, nơi rượu thường được phục vụ và tiêu thụ. Nhiều sinh viên không chỉ uống rượu mà còn lạm dụng nó, uống đến mức say xỉn thường xuyên. Tình trạng này đặc biệt phổ biến trong các dịp lễ, kỳ nghỉ và các buổi tiệc lớn. Việc uống rượu say xỉn có thể dẫn đến mất kiểm soát hành vi và gây ra nhiều hậu quả tiêu cực như tai nạn, bạo lực và các vấn đề về sức khỏe. Rượu không chỉ là dùng cho những cuộc vui chơi, tiệc tùng, giao lưu bạn bè mà sinh viên hiện nay còn dùng rượu để giải tỏa nỗi buồn, những áp lực trong học tập, bạn bè, gia đình, tình cảm nên việc tiêu thụ rượu của sinh viên hiện nay rất đáng báo động. Rượu rất dễ dàng được tiếp cận và những chỗ lưu thông rượu không được giám sát chặt chẽ nên ngày càng làm gia tăng mức độ tiêu thụ rượu.Thực trạng tiêu thụ rượu của sinh viên hiện nay là một vấn đề nghiêm trọng cần được quan tâm. Việc nghiên cứu và hiểu rõ hơn về hành vi tiêu thụ rượu của sinh viên trong một ngôi trường ở Bồ Đào Nha sẽ giúp đưa ra được một vài biện pháp giáo dục, phòng ngừa và hỗ trợ thích hợp, nhằm giảm thiểu những hậu quả tiêu cực và bảo vệ sức khỏe và tương lai của sinh viên.

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

Mục tiêu của chúng tôi là phân tích các yếu tố ảnh hưởng đến mức độ tiêu thụ rượu của sinh viên là thấp hay cao. Nghiên cứu nhằm tìm ra các thông tin và đưa ra một biện pháp giúp gia đình và nhà trường quản lý con cái, sinh viên của mình được tốt hơn.

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

Tập dữ liệu chứa thông tin từ cuộc khảo sát về sinh viên khóa học Toán ở 2 trường trung học ở Bồ Đào Nha. Bộ dữ liệu chưa các thông tin liên quan đến các sinh viên bao gồm 395 quan sát trên tổng 9 biến.

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

Bài nghiên cứu sử dụng phương pháp thống kê và các mô hình hồi quy để phân tích dữ liệu cho thấy tác động của các yếu tố ảnh hưởng đến mức độ tiêu thụ rượu của sinh viên.

1.5 Kết cấu bài nghiên cứu

Chương 1: Tổng quan nghiên cứu.
Chương 2: Dữ liệu nghiên cứu.
Chương 3: Phân tích dữ liệu.
Chương 4: Kết quả mô hình nghiên cứu.
Chương 5: Kết luận và khuyến nghị.

Chương 2: DỮ LIỆU NGHIÊN CỨU

2.1 Giới thiệu bộ dữ liệu

Bộ dữ liệu bao gồm 395 quan sát trên tổng 9 biến của sinh viên khóa học Toán của 2 trường trung học ở Bồ Đào Nha. Dưới đây là bảng chi tiết từng biến trên từng quan sát khác nhau.

std <- read.csv("student-math.csv", header = T)
library(data.table)

## Warning: package 'data.table' was built under R version 4.3.1

data.table(std)

2.2 Mô tả các biến trong bộ dữ liệu

school: ngôi trường sinh viên đang theo học (GP: Gabriel Pereira; MS: Mousinho da Silveria)
sex: giới tính của sinh viên (F: nữ; M: nam)
address: nơi ở của sinh viên (U: thành thị; R: nông thôn)
schoolsup: hỗ trợ giáo dục bổ sung (yes: có; no: không)
famsup: hỗ trợ giáo dục gia đình (yes: có; no: không)
activities: hoạt động ngoại khóa (yes: có; no: không)
internet: truy cập mạng tại nhà (yes: có; no: không)
romantic: mối quan hệ tình cảm (yes: có; no: không)
Dalc: mức tiêu thụ rượu trong 1 ngày (1: rất thấp; 2: thấp; 3: trung bình; 4: cao; 5: rất cao)

Chương 3: PHÂN TÍCH DỮ LIỆU

Trong bài luận này, biến phụ thuộc được chọn là biến Dalc (mức độ tiêu thụ rượu của sinh viên trong 1 ngày), do biến Dalc là dữ liệu định lượng với các mức độ 1, 2, 3, 4, 5 nên nhóm biến nó thành dữ liệu định tính là Dalc_2 với 2 biểu hiện cụ thể là thấp và cao. Ở các mức độ 1, 2 và 3 sẽ biến đổi thành biểu hiện tiêu thụ rượu thấp, còn mức độ 4 và 5 sẽ biến đổi thành biểu hiện tiêu thụ rượu cao như bảng dưới đây:

std <- read.csv("student-math.csv", header = T)
std$Dalc_2 <- ifelse(std$Dalc <= 2, "thấp", "cao")
print(std)

##     school sex address schoolsup famsup activities internet romantic Dalc
## 1       GP   F       U       yes     no         no       no       no    1
## 2       GP   F       U        no    yes         no      yes       no    1
## 3       GP   F       U       yes     no         no      yes       no    2
## 4       GP   F       U        no    yes        yes      yes      yes    1
## 5       GP   F       U        no    yes         no       no       no    1
## 6       GP   M       U        no    yes        yes      yes       no    1
## 7       GP   M       U        no     no         no      yes       no    1
## 8       GP   F       U       yes    yes         no       no       no    1
## 9       GP   M       U        no    yes         no      yes       no    1
## 10      GP   M       U        no    yes        yes      yes       no    1
## 11      GP   F       U        no    yes         no      yes       no    1
## 12      GP   F       U        no    yes        yes      yes       no    1
## 13      GP   M       U        no    yes        yes      yes       no    1
## 14      GP   M       U        no    yes         no      yes       no    1
## 15      GP   M       U        no    yes         no      yes      yes    1
## 16      GP   F       U        no    yes         no      yes       no    1
## 17      GP   F       U        no    yes        yes      yes       no    1
## 18      GP   F       U       yes    yes        yes       no       no    1
## 19      GP   M       U        no    yes        yes      yes       no    2
## 20      GP   M       U        no     no        yes      yes       no    1
## 21      GP   M       U        no     no         no      yes       no    1
## 22      GP   M       U        no    yes         no      yes       no    1
## 23      GP   M       U        no     no        yes      yes       no    1
## 24      GP   M       U        no    yes        yes      yes       no    2
## 25      GP   F       R       yes    yes        yes      yes       no    1
## 26      GP   F       U        no    yes         no      yes       no    1
## 27      GP   M       U        no    yes         no      yes       no    1
## 28      GP   M       U        no     no         no      yes       no    2
## 29      GP   M       U       yes    yes        yes      yes       no    1
## 30      GP   M       U        no    yes        yes      yes      yes    5
## 31      GP   M       U        no    yes         no      yes       no    3
## 32      GP   M       U        no    yes        yes      yes       no    1
## 33      GP   M       R        no    yes        yes      yes      yes    1
## 34      GP   M       U        no     no        yes      yes       no    1
## 35      GP   M       U        no    yes         no      yes       no    1
## 36      GP   F       U        no    yes        yes       no       no    1
## 37      GP   M       U        no    yes        yes      yes       no    1
## 38      GP   M       R        no    yes        yes      yes      yes    1
## 39      GP   F       R       yes    yes        yes      yes       no    1
## 40      GP   F       R       yes    yes        yes       no       no    1
## 41      GP   F       U        no    yes        yes      yes      yes    1
## 42      GP   M       U        no    yes         no      yes      yes    2
## 43      GP   M       U        no    yes        yes      yes       no    1
## 44      GP   M       U       yes    yes         no      yes       no    1
## 45      GP   F       U       yes     no        yes      yes       no    2
## 46      GP   F       U       yes    yes        yes      yes      yes    1
## 47      GP   F       U        no    yes         no      yes       no    1
## 48      GP   M       U        no     no        yes      yes       no    1
## 49      GP   M       U        no    yes         no       no       no    2
## 50      GP   F       U       yes    yes        yes      yes       no    1
## 51      GP   F       U        no    yes         no      yes       no    2
## 52      GP   F       U        no    yes         no      yes       no    1
## 53      GP   M       U        no     no         no       no       no    3
## 54      GP   F       U       yes    yes         no      yes       no    2
## 55      GP   F       U        no     no         no      yes       no    4
## 56      GP   F       U        no     no        yes      yes      yes    1
## 57      GP   F       U        no    yes        yes      yes       no    1
## 58      GP   M       U        no    yes        yes       no       no    1
## 59      GP   M       U       yes    yes        yes      yes       no    1
## 60      GP   F       U        no    yes         no      yes       no    1
## 61      GP   F       R        no    yes        yes       no       no    2
## 62      GP   F       U       yes    yes        yes      yes      yes    5
## 63      GP   F       U       yes     no        yes      yes       no    1
## 64      GP   F       U       yes    yes        yes      yes       no    2
## 65      GP   F       U       yes     no        yes      yes      yes    2
## 66      GP   F       U        no    yes        yes      yes       no    1
## 67      GP   M       U        no    yes        yes      yes      yes    5
## 68      GP   F       U       yes    yes         no      yes       no    1
## 69      GP   F       R       yes    yes         no      yes       no    1
## 70      GP   F       R        no    yes         no      yes       no    2
## 71      GP   M       U        no    yes         no      yes       no    1
## 72      GP   M       U        no     no         no      yes       no    1
## 73      GP   F       R       yes    yes         no      yes      yes    2
## 74      GP   M       U        no     no        yes       no       no    2
## 75      GP   F       U       yes    yes        yes      yes       no    2
## 76      GP   M       U        no    yes        yes      yes       no    2
## 77      GP   M       U        no     no        yes      yes       no    1
## 78      GP   F       U        no     no         no      yes      yes    1
## 79      GP   M       U       yes    yes        yes      yes       no    1
## 80      GP   F       U        no    yes         no      yes       no    1
## 81      GP   M       U       yes    yes        yes      yes      yes    1
## 82      GP   M       U       yes     no         no      yes       no    1
## 83      GP   F       U        no    yes         no      yes       no    1
## 84      GP   M       U        no     no        yes      yes       no    1
## 85      GP   F       U        no    yes        yes      yes       no    2
## 86      GP   F       U        no     no         no      yes      yes    2
## 87      GP   F       U        no    yes         no       no       no    1
## 88      GP   F       U        no    yes        yes      yes       no    1
## 89      GP   M       U        no     no        yes      yes       no    1
## 90      GP   M       U        no    yes         no       no       no    3
## 91      GP   F       U        no    yes         no      yes      yes    1
## 92      GP   F       U        no     no        yes      yes       no    1
## 93      GP   F       U       yes    yes         no       no       no    2
## 94      GP   F       U        no    yes        yes      yes       no    1
## 95      GP   M       U        no    yes        yes      yes       no    1
## 96      GP   F       R       yes    yes        yes      yes       no    1
## 97      GP   M       R       yes    yes        yes      yes       no    1
## 98      GP   F       U        no    yes         no       no      yes    1
## 99      GP   F       U        no     no        yes      yes       no    1
## 100     GP   F       U       yes    yes         no      yes       no    1
## 101     GP   M       U       yes    yes        yes      yes       no    5
## 102     GP   M       U        no    yes        yes      yes      yes    1
## 103     GP   M       U        no    yes        yes      yes       no    1
## 104     GP   F       U       yes    yes         no      yes       no    1
## 105     GP   M       U        no    yes        yes      yes       no    1
## 106     GP   F       U       yes     no         no       no       no    1
## 107     GP   F       U       yes    yes         no      yes       no    1
## 108     GP   M       U        no    yes        yes      yes       no    1
## 109     GP   M       R        no    yes        yes      yes      yes    3
## 110     GP   F       U        no    yes        yes      yes      yes    1
## 111     GP   M       U        no     no        yes      yes       no    1
## 112     GP   F       R       yes    yes        yes      yes       no    1
## 113     GP   F       U       yes     no        yes      yes       no    1
## 114     GP   M       U        no     no         no      yes       no    1
## 115     GP   M       R        no     no        yes      yes      yes    1
## 116     GP   M       U        no    yes        yes      yes       no    1
## 117     GP   M       U        no    yes        yes       no       no    1
## 118     GP   M       U        no     no        yes      yes       no    1
## 119     GP   M       R        no    yes        yes      yes       no    1
## 120     GP   M       U        no     no         no      yes       no    1
## 121     GP   F       U        no     no         no      yes       no    1
## 122     GP   M       U        no    yes        yes      yes       no    1
## 123     GP   F       U        no    yes        yes      yes      yes    1
## 124     GP   M       U        no    yes        yes      yes       no    1
## 125     GP   F       U        no     no         no      yes      yes    1
## 126     GP   M       U       yes     no         no      yes       no    3
## 127     GP   F       U       yes     no        yes      yes      yes    1
## 128     GP   F       U        no    yes         no       no       no    1
## 129     GP   M       R        no    yes        yes      yes       no    1
## 130     GP   M       R        no     no        yes      yes       no    2
## 131     GP   F       R        no    yes         no      yes      yes    2
## 132     GP   F       U        no    yes        yes      yes      yes    1
## 133     GP   F       U        no    yes         no      yes      yes    1
## 134     GP   F       U        no     no         no      yes       no    1
## 135     GP   M       R        no    yes         no       no      yes    1
## 136     GP   F       U        no    yes        yes      yes      yes    1
## 137     GP   M       R        no     no         no       no       no    2
## 138     GP   F       U        no    yes        yes      yes      yes    1
## 139     GP   M       U        no     no         no       no      yes    1
## 140     GP   F       U        no     no        yes      yes       no    1
## 141     GP   M       U       yes    yes         no      yes       no    1
## 142     GP   M       U        no    yes        yes      yes       no    2
## 143     GP   F       U        no    yes        yes      yes       no    1
## 144     GP   F       U        no     no         no      yes       no    3
## 145     GP   M       U        no    yes         no      yes       no    1
## 146     GP   F       U        no    yes         no      yes       no    1
## 147     GP   F       U        no    yes         no      yes       no    1
## 148     GP   F       U        no    yes         no      yes       no    1
## 149     GP   M       U        no    yes         no      yes      yes    2
## 150     GP   M       U        no     no         no      yes       no    2
## 151     GP   M       U        no     no         no      yes      yes    2
## 152     GP   M       U        no     no        yes       no      yes    3
## 153     GP   F       R        no    yes        yes      yes      yes    2
## 154     GP   M       U        no    yes         no      yes      yes    1
## 155     GP   F       U       yes    yes         no       no      yes    1
## 156     GP   M       R       yes     no        yes       no       no    1
## 157     GP   M       R        no     no         no       no       no    3
## 158     GP   F       R        no    yes        yes       no       no    1
## 159     GP   M       R        no     no         no       no       no    1
## 160     GP   M       U        no    yes         no      yes      yes    4
## 161     GP   M       R        no     no        yes      yes      yes    2
## 162     GP   M       R       yes    yes         no      yes      yes    1
## 163     GP   M       U        no     no        yes       no       no    2
## 164     GP   M       U        no     no         no      yes       no    1
## 165     GP   M       R        no     no        yes       no      yes    1
## 166     GP   M       U        no    yes        yes       no       no    1
## 167     GP   M       U        no     no         no      yes       no    2
## 168     GP   F       U        no     no         no      yes      yes    1
## 169     GP   F       U        no    yes         no      yes       no    1
## 170     GP   F       U        no    yes         no      yes      yes    1
## 171     GP   M       U        no    yes        yes      yes       no    2
## 172     GP   M       U        no    yes        yes      yes      yes    1
## 173     GP   M       U        no    yes        yes      yes       no    1
## 174     GP   F       U        no     no        yes      yes      yes    1
## 175     GP   F       U        no    yes        yes      yes       no    1
## 176     GP   M       U        no     no        yes      yes       no    4
## 177     GP   F       U        no     no        yes      yes       no    1
## 178     GP   M       U        no     no        yes      yes       no    1
## 179     GP   M       R        no    yes        yes      yes      yes    3
## 180     GP   M       U        no    yes        yes      yes      yes    1
## 181     GP   M       U        no    yes        yes      yes       no    2
## 182     GP   M       U        no     no        yes      yes      yes    1
## 183     GP   F       U        no    yes        yes       no       no    2
## 184     GP   F       U        no    yes        yes      yes      yes    2
## 185     GP   F       U        no    yes         no      yes       no    1
## 186     GP   M       U        no    yes        yes      yes      yes    2
## 187     GP   M       U        no    yes        yes      yes      yes    1
## 188     GP   M       U        no     no        yes      yes      yes    1
## 189     GP   F       U        no    yes         no      yes      yes    1
## 190     GP   M       R        no     no         no       no       no    1
## 191     GP   F       U        no     no         no      yes       no    1
## 192     GP   F       U        no     no        yes      yes       no    1
## 193     GP   M       U        no     no        yes      yes       no    4
## 194     GP   M       R        no    yes        yes      yes       no    3
## 195     GP   M       U        no     no         no      yes       no    1
## 196     GP   F       U        no     no        yes      yes      yes    1
## 197     GP   M       U        no     no         no      yes       no    1
## 198     GP   M       R        no    yes        yes      yes       no    3
## 199     GP   F       U        no    yes         no      yes       no    2
## 200     GP   F       U        no    yes         no      yes       no    1
## 201     GP   F       U        no    yes        yes      yes       no    1
## 202     GP   F       U       yes    yes        yes       no       no    1
## 203     GP   F       U        no    yes         no       no       no    1
## 204     GP   F       R        no    yes         no      yes       no    1
## 205     GP   F       R        no    yes        yes      yes       no    1
## 206     GP   F       U        no    yes         no      yes      yes    3
## 207     GP   F       U        no    yes         no      yes       no    2
## 208     GP   F       U        no     no        yes      yes      yes    1
## 209     GP   F       U        no    yes         no       no       no    1
## 210     GP   F       R        no    yes        yes      yes      yes    1
## 211     GP   F       U        no    yes        yes      yes       no    1
## 212     GP   M       U        no    yes         no      yes      yes    4
## 213     GP   F       U       yes    yes         no      yes       no    1
## 214     GP   M       U        no    yes        yes      yes       no    2
## 215     GP   F       R        no    yes         no       no       no    1
## 216     GP   F       U        no     no         no      yes       no    1
## 217     GP   F       U        no     no         no      yes      yes    2
## 218     GP   M       U        no    yes         no      yes       no    2
## 219     GP   F       U        no    yes         no       no       no    1
## 220     GP   F       U        no    yes        yes      yes       no    1
## 221     GP   F       R        no    yes        yes      yes       no    1
## 222     GP   F       U        no    yes        yes       no      yes    1
## 223     GP   F       U       yes     no         no      yes       no    1
## 224     GP   M       U        no    yes         no      yes       no    5
## 225     GP   F       U        no    yes        yes      yes       no    1
## 226     GP   F       R        no     no        yes      yes      yes    1
## 227     GP   F       U        no     no        yes      yes       no    1
## 228     GP   M       U        no    yes         no      yes       no    1
## 229     GP   M       U       yes    yes        yes      yes      yes    4
## 230     GP   F       U        no     no        yes      yes      yes    1
## 231     GP   F       U        no     no        yes      yes      yes    1
## 232     GP   M       R        no    yes        yes      yes       no    1
## 233     GP   M       U       yes    yes        yes      yes      yes    1
## 234     GP   M       U        no    yes        yes      yes       no    2
## 235     GP   M       U        no    yes         no      yes       no    1
## 236     GP   M       U        no     no        yes      yes      yes    1
## 237     GP   M       U        no     no        yes      yes      yes    5
## 238     GP   F       U        no     no         no      yes      yes    1
## 239     GP   F       R        no     no        yes       no       no    1
## 240     GP   M       U        no     no         no      yes       no    3
## 241     GP   M       U        no     no        yes      yes      yes    1
## 242     GP   M       R        no    yes         no      yes       no    2
## 243     GP   M       U        no     no        yes      yes       no    1
## 244     GP   M       U        no     no        yes      yes       no    1
## 245     GP   F       U        no    yes         no      yes      yes    1
## 246     GP   M       U        no     no         no      yes       no    1
## 247     GP   M       U        no     no         no      yes       no    1
## 248     GP   M       U        no     no         no      yes      yes    5
## 249     GP   M       R        no    yes         no      yes      yes    1
## 250     GP   M       U        no     no         no      yes       no    2
## 251     GP   M       U        no     no         no      yes       no    2
## 252     GP   M       U       yes    yes         no      yes       no    1
## 253     GP   M       U        no     no         no      yes       no    2
## 254     GP   M       R        no     no        yes       no       no    1
## 255     GP   M       R        no     no         no      yes       no    2
## 256     GP   M       U        no    yes        yes      yes       no    1
## 257     GP   F       U        no    yes        yes      yes       no    1
## 258     GP   M       U        no    yes         no      yes       no    1
## 259     GP   M       U        no     no        yes      yes       no    1
## 260     GP   F       U        no     no        yes      yes      yes    1
## 261     GP   F       U        no    yes         no      yes      yes    1
## 262     GP   M       U        no    yes         no      yes       no    1
## 263     GP   M       R        no     no        yes       no       no    1
## 264     GP   F       U        no     no        yes       no       no    1
## 265     GP   F       U        no    yes        yes      yes      yes    1
## 266     GP   M       R        no    yes        yes      yes       no    3
## 267     GP   M       U        no     no        yes      yes      yes    3
## 268     GP   F       R        no     no        yes      yes       no    2
## 269     GP   M       U        no    yes        yes      yes      yes    1
## 270     GP   F       R        no    yes         no      yes      yes    1
## 271     GP   F       U        no    yes        yes      yes       no    3
## 272     GP   F       U        no    yes        yes      yes      yes    1
## 273     GP   F       U        no    yes         no       no       no    1
## 274     GP   M       R        no    yes        yes       no      yes    2
## 275     GP   F       U        no    yes         no      yes      yes    1
## 276     GP   F       U       yes    yes         no      yes      yes    2
## 277     GP   F       R        no     no         no      yes      yes    1
## 278     GP   M       U        no     no        yes      yes       no    1
## 279     GP   F       U       yes    yes        yes      yes      yes    1
## 280     GP   M       U        no     no        yes      yes       no    1
## 281     GP   M       U        no     no        yes      yes      yes    2
## 282     GP   M       U        no     no         no      yes       no    3
## 283     GP   F       R        no    yes        yes       no       no    1
## 284     GP   F       U       yes     no        yes      yes       no    1
## 285     GP   F       U        no    yes         no      yes       no    1
## 286     GP   M       U        no     no         no      yes       no    1
## 287     GP   F       U        no    yes         no      yes       no    1
## 288     GP   F       U        no    yes         no      yes       no    1
## 289     GP   M       U        no     no        yes      yes       no    1
## 290     GP   M       U        no    yes        yes      yes       no    1
## 291     GP   M       U        no    yes        yes      yes      yes    1
## 292     GP   F       U        no    yes         no      yes       no    1
## 293     GP   F       U        no     no         no      yes      yes    1
## 294     GP   F       R        no    yes         no       no       no    1
## 295     GP   M       R        no    yes        yes      yes       no    1
## 296     GP   M       U        no    yes         no      yes       no    1
## 297     GP   F       U        no    yes        yes      yes       no    2
## 298     GP   F       U        no    yes         no      yes      yes    1
## 299     GP   F       U        no    yes         no      yes       no    1
## 300     GP   M       U        no    yes         no      yes      yes    2
## 301     GP   F       U        no    yes         no      yes      yes    1
## 302     GP   M       U        no     no         no      yes       no    2
## 303     GP   F       U        no    yes         no      yes       no    1
## 304     GP   F       U        no    yes        yes      yes       no    1
## 305     GP   M       U        no    yes        yes      yes      yes    1
## 306     GP   F       U        no    yes        yes      yes       no    1
## 307     GP   M       U        no     no        yes       no       no    1
## 308     GP   M       U        no    yes         no      yes      yes    1
## 309     GP   M       R        no     no        yes       no      yes    1
## 310     GP   F       U       yes    yes        yes      yes       no    1
## 311     GP   F       U        no     no        yes       no      yes    2
## 312     GP   F       U        no    yes         no      yes      yes    1
## 313     GP   M       U        no     no         no      yes       no    2
## 314     GP   F       U        no    yes         no      yes      yes    1
## 315     GP   F       U        no     no         no      yes      yes    1
## 316     GP   F       R        no     no         no      yes      yes    1
## 317     GP   F       U        no    yes        yes      yes       no    1
## 318     GP   F       U        no    yes        yes      yes      yes    1
## 319     GP   F       R        no    yes        yes      yes       no    2
## 320     GP   F       U        no    yes         no      yes       no    3
## 321     GP   F       U        no    yes         no      yes      yes    1
## 322     GP   F       U        no    yes         no       no      yes    1
## 323     GP   F       R        no    yes        yes      yes       no    2
## 324     GP   F       U        no    yes         no      yes       no    2
## 325     GP   F       U        no     no         no      yes       no    2
## 326     GP   M       U        no     no        yes      yes       no    2
## 327     GP   M       U        no     no        yes      yes       no    3
## 328     GP   M       R        no    yes         no      yes       no    5
## 329     GP   F       U        no    yes        yes      yes       no    1
## 330     GP   F       U        no    yes         no      yes      yes    1
## 331     GP   M       U        no    yes        yes      yes       no    2
## 332     GP   F       R        no    yes         no      yes      yes    1
## 333     GP   F       U        no     no        yes      yes       no    1
## 334     GP   F       U        no     no        yes      yes      yes    1
## 335     GP   F       R        no     no        yes       no       no    1
## 336     GP   F       U        no     no         no      yes       no    1
## 337     GP   F       R        no     no         no       no       no    1
## 338     GP   F       U        no    yes         no      yes      yes    2
## 339     GP   F       U        no    yes         no      yes       no    1
## 340     GP   F       R        no    yes         no      yes       no    2
## 341     GP   F       U        no     no        yes      yes      yes    1
## 342     GP   M       U        no    yes        yes      yes       no    2
## 343     GP   M       U        no     no        yes      yes      yes    1
## 344     GP   F       U        no    yes         no      yes      yes    1
## 345     GP   F       U        no    yes         no      yes       no    1
## 346     GP   F       U        no     no         no      yes      yes    2
## 347     GP   M       R        no     no         no      yes      yes    1
## 348     GP   M       U        no    yes         no      yes      yes    2
## 349     GP   F       U        no    yes        yes      yes      yes    1
## 350     MS   M       R        no    yes         no      yes       no    5
## 351     MS   M       R        no     no         no      yes       no    3
## 352     MS   M       U        no    yes         no      yes       no    2
## 353     MS   M       U        no     no         no      yes      yes    2
## 354     MS   M       R        no    yes         no      yes       no    3
## 355     MS   M       R        no    yes        yes      yes      yes    1
## 356     MS   F       U        no    yes         no       no      yes    1
## 357     MS   F       R        no    yes        yes      yes       no    1
## 358     MS   F       U        no     no         no       no      yes    1
## 359     MS   M       U        no     no         no      yes      yes    1
## 360     MS   F       U        no     no         no      yes       no    1
## 361     MS   F       R        no     no         no       no      yes    1
## 362     MS   M       R        no     no        yes       no       no    2
## 363     MS   F       U        no    yes         no      yes      yes    1
## 364     MS   F       U        no    yes        yes      yes      yes    1
## 365     MS   F       R        no     no         no       no       no    1
## 366     MS   M       R        no    yes         no       no       no    2
## 367     MS   M       U        no     no         no      yes      yes    2
## 368     MS   F       R        no    yes         no      yes      yes    1
## 369     MS   F       U        no    yes         no      yes      yes    1
## 370     MS   F       R        no    yes         no      yes      yes    4
## 371     MS   F       U        no     no        yes       no      yes    1
## 372     MS   M       R        no    yes        yes      yes      yes    2
## 373     MS   F       U        no     no        yes       no      yes    1
## 374     MS   F       R        no     no        yes      yes       no    1
## 375     MS   F       R        no     no         no      yes       no    1
## 376     MS   F       R        no     no         no      yes       no    1
## 377     MS   F       U        no    yes         no      yes      yes    1
## 378     MS   F       R        no     no        yes      yes       no    3
## 379     MS   F       U        no     no         no      yes      yes    1
## 380     MS   F       R        no    yes        yes      yes       no    2
## 381     MS   M       U        no     no        yes      yes       no    1
## 382     MS   M       R        no     no        yes      yes      yes    1
## 383     MS   M       U        no     no        yes      yes       no    1
## 384     MS   M       R        no     no         no       no       no    1
## 385     MS   M       R        no     no         no       no       no    4
## 386     MS   F       R        no     no         no       no       no    1
## 387     MS   F       R        no    yes        yes      yes      yes    2
## 388     MS   F       R        no     no        yes      yes       no    1
## 389     MS   F       U        no    yes         no      yes       no    1
## 390     MS   F       U        no     no        yes       no       no    1
## 391     MS   M       U        no    yes         no       no       no    4
## 392     MS   M       U        no     no         no      yes       no    3
## 393     MS   M       R        no     no         no       no       no    3
## 394     MS   M       R        no     no         no      yes       no    3
## 395     MS   M       U        no     no         no      yes       no    3
##     Dalc_2
## 1     thấp
## 2     thấp
## 3     thấp
## 4     thấp
## 5     thấp
## 6     thấp
## 7     thấp
## 8     thấp
## 9     thấp
## 10    thấp
## 11    thấp
## 12    thấp
## 13    thấp
## 14    thấp
## 15    thấp
## 16    thấp
## 17    thấp
## 18    thấp
## 19    thấp
## 20    thấp
## 21    thấp
## 22    thấp
## 23    thấp
## 24    thấp
## 25    thấp
## 26    thấp
## 27    thấp
## 28    thấp
## 29    thấp
## 30     cao
## 31     cao
## 32    thấp
## 33    thấp
## 34    thấp
## 35    thấp
## 36    thấp
## 37    thấp
## 38    thấp
## 39    thấp
## 40    thấp
## 41    thấp
## 42    thấp
## 43    thấp
## 44    thấp
## 45    thấp
## 46    thấp
## 47    thấp
## 48    thấp
## 49    thấp
## 50    thấp
## 51    thấp
## 52    thấp
## 53     cao
## 54    thấp
## 55     cao
## 56    thấp
## 57    thấp
## 58    thấp
## 59    thấp
## 60    thấp
## 61    thấp
## 62     cao
## 63    thấp
## 64    thấp
## 65    thấp
## 66    thấp
## 67     cao
## 68    thấp
## 69    thấp
## 70    thấp
## 71    thấp
## 72    thấp
## 73    thấp
## 74    thấp
## 75    thấp
## 76    thấp
## 77    thấp
## 78    thấp
## 79    thấp
## 80    thấp
## 81    thấp
## 82    thấp
## 83    thấp
## 84    thấp
## 85    thấp
## 86    thấp
## 87    thấp
## 88    thấp
## 89    thấp
## 90     cao
## 91    thấp
## 92    thấp
## 93    thấp
## 94    thấp
## 95    thấp
## 96    thấp
## 97    thấp
## 98    thấp
## 99    thấp
## 100   thấp
## 101    cao
## 102   thấp
## 103   thấp
## 104   thấp
## 105   thấp
## 106   thấp
## 107   thấp
## 108   thấp
## 109    cao
## 110   thấp
## 111   thấp
## 112   thấp
## 113   thấp
## 114   thấp
## 115   thấp
## 116   thấp
## 117   thấp
## 118   thấp
## 119   thấp
## 120   thấp
## 121   thấp
## 122   thấp
## 123   thấp
## 124   thấp
## 125   thấp
## 126    cao
## 127   thấp
## 128   thấp
## 129   thấp
## 130   thấp
## 131   thấp
## 132   thấp
## 133   thấp
## 134   thấp
## 135   thấp
## 136   thấp
## 137   thấp
## 138   thấp
## 139   thấp
## 140   thấp
## 141   thấp
## 142   thấp
## 143   thấp
## 144    cao
## 145   thấp
## 146   thấp
## 147   thấp
## 148   thấp
## 149   thấp
## 150   thấp
## 151   thấp
## 152    cao
## 153   thấp
## 154   thấp
## 155   thấp
## 156   thấp
## 157    cao
## 158   thấp
## 159   thấp
## 160    cao
## 161   thấp
## 162   thấp
## 163   thấp
## 164   thấp
## 165   thấp
## 166   thấp
## 167   thấp
## 168   thấp
## 169   thấp
## 170   thấp
## 171   thấp
## 172   thấp
## 173   thấp
## 174   thấp
## 175   thấp
## 176    cao
## 177   thấp
## 178   thấp
## 179    cao
## 180   thấp
## 181   thấp
## 182   thấp
## 183   thấp
## 184   thấp
## 185   thấp
## 186   thấp
## 187   thấp
## 188   thấp
## 189   thấp
## 190   thấp
## 191   thấp
## 192   thấp
## 193    cao
## 194    cao
## 195   thấp
## 196   thấp
## 197   thấp
## 198    cao
## 199   thấp
## 200   thấp
## 201   thấp
## 202   thấp
## 203   thấp
## 204   thấp
## 205   thấp
## 206    cao
## 207   thấp
## 208   thấp
## 209   thấp
## 210   thấp
## 211   thấp
## 212    cao
## 213   thấp
## 214   thấp
## 215   thấp
## 216   thấp
## 217   thấp
## 218   thấp
## 219   thấp
## 220   thấp
## 221   thấp
## 222   thấp
## 223   thấp
## 224    cao
## 225   thấp
## 226   thấp
## 227   thấp
## 228   thấp
## 229    cao
## 230   thấp
## 231   thấp
## 232   thấp
## 233   thấp
## 234   thấp
## 235   thấp
## 236   thấp
## 237    cao
## 238   thấp
## 239   thấp
## 240    cao
## 241   thấp
## 242   thấp
## 243   thấp
## 244   thấp
## 245   thấp
## 246   thấp
## 247   thấp
## 248    cao
## 249   thấp
## 250   thấp
## 251   thấp
## 252   thấp
## 253   thấp
## 254   thấp
## 255   thấp
## 256   thấp
## 257   thấp
## 258   thấp
## 259   thấp
## 260   thấp
## 261   thấp
## 262   thấp
## 263   thấp
## 264   thấp
## 265   thấp
## 266    cao
## 267    cao
## 268   thấp
## 269   thấp
## 270   thấp
## 271    cao
## 272   thấp
## 273   thấp
## 274   thấp
## 275   thấp
## 276   thấp
## 277   thấp
## 278   thấp
## 279   thấp
## 280   thấp
## 281   thấp
## 282    cao
## 283   thấp
## 284   thấp
## 285   thấp
## 286   thấp
## 287   thấp
## 288   thấp
## 289   thấp
## 290   thấp
## 291   thấp
## 292   thấp
## 293   thấp
## 294   thấp
## 295   thấp
## 296   thấp
## 297   thấp
## 298   thấp
## 299   thấp
## 300   thấp
## 301   thấp
## 302   thấp
## 303   thấp
## 304   thấp
## 305   thấp
## 306   thấp
## 307   thấp
## 308   thấp
## 309   thấp
## 310   thấp
## 311   thấp
## 312   thấp
## 313   thấp
## 314   thấp
## 315   thấp
## 316   thấp
## 317   thấp
## 318   thấp
## 319   thấp
## 320    cao
## 321   thấp
## 322   thấp
## 323   thấp
## 324   thấp
## 325   thấp
## 326   thấp
## 327    cao
## 328    cao
## 329   thấp
## 330   thấp
## 331   thấp
## 332   thấp
## 333   thấp
## 334   thấp
## 335   thấp
## 336   thấp
## 337   thấp
## 338   thấp
## 339   thấp
## 340   thấp
## 341   thấp
## 342   thấp
## 343   thấp
## 344   thấp
## 345   thấp
## 346   thấp
## 347   thấp
## 348   thấp
## 349   thấp
## 350    cao
## 351    cao
## 352   thấp
## 353   thấp
## 354    cao
## 355   thấp
## 356   thấp
## 357   thấp
## 358   thấp
## 359   thấp
## 360   thấp
## 361   thấp
## 362   thấp
## 363   thấp
## 364   thấp
## 365   thấp
## 366   thấp
## 367   thấp
## 368   thấp
## 369   thấp
## 370    cao
## 371   thấp
## 372   thấp
## 373   thấp
## 374   thấp
## 375   thấp
## 376   thấp
## 377   thấp
## 378    cao
## 379   thấp
## 380   thấp
## 381   thấp
## 382   thấp
## 383   thấp
## 384   thấp
## 385    cao
## 386   thấp
## 387   thấp
## 388   thấp
## 389   thấp
## 390   thấp
## 391    cao
## 392    cao
## 393    cao
## 394    cao
## 395    cao

3.1 Thống kê mô tả biến phụ thuộc

3.1.1 Lập bảng tần số và bảng tần suất

table(std$Dalc_2)

## 
##  cao thấp 
##   44  351

table(std$Dalc_2)/length(std$Dalc_2)

## 
##       cao      thấp 
## 0.1113924 0.8886076

3.1.2 Lập biểu đồ

library(dplyr)

## Warning: package 'dplyr' was built under R version 4.3.1

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:data.table':
## 
##     between, first, last

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(tidyverse)

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

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

## Warning: package 'tibble' was built under R version 4.3.1

## Warning: package 'tidyr' was built under R version 4.3.1

## Warning: package 'readr' was built under R version 4.3.1

## Warning: package 'purrr' was built under R version 4.3.1

## Warning: package 'forcats' was built under R version 4.3.1

## Warning: package 'lubridate' was built under R version 4.3.1

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats   1.0.0     ✔ readr     2.1.4
## ✔ ggplot2   3.5.1     ✔ stringr   1.5.0
## ✔ lubridate 1.9.2     ✔ tibble    3.2.1
## ✔ purrr     1.0.1     ✔ tidyr     1.3.0

## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::between()     masks data.table::between()
## ✖ dplyr::filter()      masks stats::filter()
## ✖ dplyr::first()       masks data.table::first()
## ✖ lubridate::hour()    masks data.table::hour()
## ✖ lubridate::isoweek() masks data.table::isoweek()
## ✖ dplyr::lag()         masks stats::lag()
## ✖ dplyr::last()        masks data.table::last()
## ✖ lubridate::mday()    masks data.table::mday()
## ✖ lubridate::minute()  masks data.table::minute()
## ✖ lubridate::month()   masks data.table::month()
## ✖ lubridate::quarter() masks data.table::quarter()
## ✖ lubridate::second()  masks data.table::second()
## ✖ purrr::transpose()   masks data.table::transpose()
## ✖ lubridate::wday()    masks data.table::wday()
## ✖ lubridate::week()    masks data.table::week()
## ✖ lubridate::yday()    masks data.table::yday()
## ✖ lubridate::year()    masks data.table::year()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

std %>% group_by(Dalc_2) %>%  summarise(n = n()) %>% 
    ggplot(aes(x=Dalc_2,y=n)) + 
    geom_col(fill = "#89CFF0")+
    geom_text(aes(label =n),vjust = 1, color = "white")+
    theme_bw()

Nhận xét: Biểu đồ cột và bảng ở trên thể hiện mức độ tiêu thụ rượu của sinh viên trong 1 ngày. Theo khảo sát thì mức tiêu thụ rượu của sinh viên trong 1 ngày với mức độ tiêu thụ rượu cao chỉ có 44 sinh viên chiếm tỷ lệ 11.14%, còn lại 351 sinh viên với mức độ tiêu thụ rượu thấp chiếm tỷ lệ 88.86%. Như kết quả cho thấy, tại 2 ngôi trường này thì sinh viên tiêu thụ rượu trong 1 ngày rất là ít.

Biểu đồ tròn dưới đây cho thấy rõ hơn về mức độ tiêu thụ rượu của sinh viên trong 1 ngày.

ap <- prop.table(table(std$Dalc_2)) * 100
ggplot(mapping = aes(x = "", y = ap, fill = factor(names(ap)))) +
  geom_bar(width = 1, stat = "identity") +
  coord_polar("y") +
  labs(title = "Phân bố mức độ tiêu thụ rượu trong 1 ngày",
       fill = "Mức tiêu thụ") +
  theme_void() +
  geom_text(aes(label = paste0(round(ap, 1), "%")),
            position = position_stack(vjust = 0.5))

## Don't know how to automatically pick scale for object of type <table>.
## Defaulting to continuous.

3.2 Thống kê suy diễn

3.2.1 Kiểm định tính độc lập

Phương pháp Chi bình phương

\[\chi^2 = \sum_{i,j} \frac{(n_{ij} -\hat{\mu}_{ij})^2}{\hat{\mu}_{ij}}\]

Với \(n_{i,j}\) là giá trị của ô \(i,j,\hat{\mu}_{ij}=\frac{n_{i+} \cdot n_{+j}}{n}\)

3.2.1.1 Biến Dalc_2 và school

Giả thuyết:

\(H_0\): Dalc_2 và school độc lập

\(H_1\): Dalc_2 và school không độc lập

mtht <- table(std$Dalc_2, std$school)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 7.1839, df = 1, p-value = 0.007356

Nhận xét: Qua kết quả kiểm định cho ta p−value = 0.007356 < 0.05, bác bỏ \(H_0\), nghĩa là giữa việc sinh viên học ở trường nào và mức độ tiêu thụ rượu là có liên quan với nhau.

3.2.1.2 Biến Dalc_2 và sex

Giả thuyết:

\(H_0\): Dalc_2 và sex độc lập

\(H_1\): Dalc_2 và sex không độc lập

mtht <- table(std$Dalc_2, std$sex)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 22.078, df = 1, p-value = 2.618e-06

Nhận xét: Qua kết quả kiểm định cho ta p−value = 2.618e-06 < 0.05, bác bỏ \(H_0\), nghĩa là giữa giới tính của sinh viên và mức độ tiêu thụ rượu là có liên quan với nhau.

3.2.1.3 Biến Dalc_2 và address

Giả thuyết:

\(H_0\): Dalc_2 và address độc lập

\(H_1\): Dalc_2 và address không độc lập

mtht <- table(std$Dalc_2, std$address)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 3.2594, df = 1, p-value = 0.07101

Nhận xét: Qua kết quả kiểm định cho ta p−value = 0.07101 > 0.05, chưa có cơ sở bác bỏ \(H_0\), nghĩa là giữa nơi ở của sinh viên và mức độ tiêu thụ rượu không có liên quan với nhau.

3.2.1.4 Biến Dalc_2 và schoolsup

Giả thuyết:

\(H_0\): Dalc_2 và schoolsup độc lập

\(H_1\): Dalc_2 và schoolsup không độc lập

mtht <- table(std$Dalc_2, std$schoolsup)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 0.31726, df = 1, p-value = 0.5733

Nhận xét: Qua kết quả kiểm định cho ta p−value = 0.5733 > 0.05, chưa có cơ sở bác bỏ \(H_0\), nghĩa là giữa việc hỗ trợ giáo dục bổ sung cho sinh viên và mức độ tiêu thụ rượu không có liên quan với nhau.

3.2.1.5 Biến Dalc_2 và famsup

Giả thuyết:

\(H_0\): Dalc_2 và famsup độc lập

\(H_1\): Dalc_2 và famsup không độc lập

mtht <- table(std$Dalc_2, std$famsup)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 1.288, df = 1, p-value = 0.2564

Nhận xét: Qua kết quả kiểm định cho ta p−value = 0.2564 > 0.05, chưa có cơ sở bác bỏ \(H_0\), nghĩa là giữa việc hỗ trợ giáo dục gia đình và mức độ tiêu thụ rượu không có liên quan với nhau.

3.2.1.6 Biến Dalc_2 và activities

Giả thuyết:

\(H_0\): Dalc_2 và activities độc lập

\(H_1\): Dalc_2 và activities không độc lập

mtht <- table(std$Dalc_2, std$activities)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 1.5485, df = 1, p-value = 0.2134

Nhận xét: Qua kết quả kiểm định cho ta p−value = 0.2134 > 0.05, chưa đủ cơ sở bác bỏ \(H_0\), nghĩa là giữa việc hoạt động ngoại khóa của sinh viên và mức độ tiêu thụ rượu không có liên quan với nhau.

3.2.1.7 Biến Dalc_2 và internet

Giả thuyết:

\(H_0\): Dalc_2 và internet độc lập

\(H_1\): Dalc_2 và internet không độc lập

mtht <- table(std$Dalc_2, std$internet)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 1.7168e-30, df = 1, p-value = 1

Nhận xét: Qua kết quả kiểm định cho ta p−value = 1 > 0.05, chưa đủ cơ sở bác bỏ \(H_0\), nghĩa là giữa việc sử dụng mạng tại nhà của sinh viên và mức độ tiêu thụ rượu không có liên quan với nhau.

3.2.1.6 Biến Dalc_2 và romantics

Giả thuyết:

\(H_0\): Dalc_2 và romantic độc lập

\(H_1\): Dalc_2 và sex romantic độc lập

mtht <- table(std$Dalc_2, std$romantic)
chisq.test(mtht)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  mtht
## X-squared = 0.0047742, df = 1, p-value = 0.9449

Nhận xét: Qua kết quả kiểm định cho ta p−value = 0.9449 > 0.05, chưa đủ cơ sở bác bỏ \(H_0\), nghĩa là giữa mối quan hệ tình cảm của sinh viên và mức độ tiêu thụ rượu không có liên quan với nhau.

Kết luận: Chỉ có biến school (trường học) và sex (giới tính) là có ảnh hưởng tới biến Dalc_2. Tiếp theo, nhóm tiến hành mô tả 2 biến, từ đó đưa ra cái nhìn khái quát hơn về mức độ ảnh hưởng lên biến Dalc_2.

3.3 Thống kê mô tả 2 biến

3.3.1 Biến Dalc_2 và sex

3.3.1.1 Lập bảng tần số và bảng tần suất

htnh <- table(std$Dalc_2,std$school)
addmargins(htnh)

##       
##         GP  MS Sum
##   cao   33  11  44
##   thấp 316  35 351
##   Sum  349  46 395

htnh <- prop.table(htnh)
addmargins(htnh)

##       
##                GP         MS        Sum
##   cao  0.08354430 0.02784810 0.11139241
##   thấp 0.80000000 0.08860759 0.88860759
##   Sum  0.88354430 0.11645570 1.00000000

3.3.1.2 Lập biểu đồ

std %>% group_by(Dalc_2,school) %>%  summarise(n = n()) %>% 
     ggplot(aes(x=Dalc_2,y=n, fill = school)) + 
     geom_col(position = 'dodge')+
     theme_bw()

## `summarise()` has grouped output by 'Dalc_2'. You can override using the
## `.groups` argument.

Nhận xét: Bảng và biểu đồ trên cho thấy ở trường GP có 33 sinh viên có mức tiêu thụ rượu cao chiếm tỷ lệ 8.35%, còn lại 316 sinh viên có mức tiêu thụ rượu thấp chiếm tỷ lệ 80%; ở trường MS có 11 sinh viên có mức tiêu thụ rượu cao chiếm tỷ lệ 2.78%, còn lại 35 sinh viên có mức tiêu thụ rượu thấp chiếm 8.86%. Tóm lại thì cả hai trường đều có số lượng ít sinh viên có mức tiêu thụ rượu cao, đa số đều tiêu thụ rượu thấp hoặc hầu như là không tiêu thụ.

3.3.2 Biến Dalc_2 và sex

3.3.2.1 Lập bảng tần số và bảng tần suất

htnh <- table(std$Dalc_2,std$sex)
addmargins(htnh)

##       
##          F   M Sum
##   cao    8  36  44
##   thấp 200 151 351
##   Sum  208 187 395

htnh <- prop.table(htnh)
addmargins(htnh)

##       
##                 F          M        Sum
##   cao  0.02025316 0.09113924 0.11139241
##   thấp 0.50632911 0.38227848 0.88860759
##   Sum  0.52658228 0.47341772 1.00000000

3.3.2.2 Lập biểu đồ

std %>% group_by(Dalc_2,sex) %>%  summarise(n = n()) %>% 
     ggplot(aes(x=Dalc_2,y=n, fill = sex)) + 
     geom_col(position = 'dodge')+
     theme_bw()

## `summarise()` has grouped output by 'Dalc_2'. You can override using the
## `.groups` argument.

Nhận xét: Bảng và biểu đồ trên cho thấy ở giới tính nữ có 8 sinh viên có mức tiêu thụ rượu cao chiếm tỷ lệ 2.03%, còn lại 200 sinh viên nữ có mức tiêu thụ rượu thấp chiếm tỷ lệ 50.63%%; ở giới tính nam có 36 sinh viên có mức tiêu thụ rượu cao chiếm tỷ lệ 9.11%, còn lại 151 sinh viên nam có mức tiêu thụ rượu thấp chiếm 38.23%. Tóm lại thì cả hai giới tính đều có số lượng ít sinh viên có mức tiêu thụ rượu cao, đa số đều tiêu thụ rượu thấp hoặc hầu như là không tiêu thụ.

3.4 Relative Risk và khoảng ước lượng cho Relative Risk

3.4.1 Biến Dalc_2 và school

library(DescTools)

## Warning: package 'DescTools' was built under R version 4.3.3

## 
## Attaching package: 'DescTools'

## The following object is masked from 'package:data.table':
## 
##     %like%

table(std$Dalc_2,std$school)

##       
##         GP  MS
##   cao   33  11
##   thấp 316  35

mtht <- table(std$Dalc_2,std$school)
RelRisk(mtht, conf.level = .95)

## rel. risk    lwr.ci    upr.ci 
## 0.8330696 0.6711519 0.9549721

Nhận xét:

Relative Risk của Dalc_2 và school là 0.8330696 nhỏ hơn 1. Điều này có nghĩa là sinh viên trường GP có mức độ tiêu thụ rượu cao nhiều hơn so với trường MS (mức độ tiêu thụ rượu cao ở trường GP bằng 0.833 lần so với trường MS).
Khoảng tin cậy này cho thấy rằng, với mức độ tin cậy (thường là 95%), mức độ tiêu thụ rượu thực sự nằm trong khoảng từ 0.6711519 đến 0.9549721. Vì khoảng tin cậy này không bao gồm giá trị 1, kết quả này cho thấy có ý nghĩa thống kê giữa mức độ tiêu thụ rượu cao ở hai trường.

3.4.2 Biến Dalc_2 và sex

library(DescTools)
table(std$Dalc_2,std$sex)

##       
##          F   M
##   cao    8  36
##   thấp 200 151

mtht <- table(std$Dalc_2,std$sex)
RelRisk(mtht, conf.level = .95)

## rel. risk    lwr.ci    upr.ci 
## 0.3190909 0.1659926 0.5656633

Nhận xét:

Relative Risk của Dalc_2 và sex là 0.3190909 nhỏ hơn 1. Điều này có nghĩa là sinh viên nữ có mức độ tiêu thụ rượu cao ít hơn so với sinh viên nam (mức độ tiêu thụ rượu cao ở sinh viên nữ bằng 0.319 lần so với sinh viên nam).
Khoảng tin cậy này cho thấy rằng, với mức độ tin cậy (thường là 95%), mức độ tiêu thụ rượu thực sự nằm trong khoảng từ 0.1659926 đến 0.5656633. Vì khoảng tin cậy này không bao gồm giá trị 1, kết quả này cho thấy có ý nghĩa thống kê giữa mức độ tiêu thụ rượu cao ở giới tính nữ và nam.

3.5 Odds ratio và khoảng ước lượng cho Odds ratio

3.5.1 Biến Dalc_2 và school

table(std$Dalc_2,std$school)

##       
##         GP  MS
##   cao   33  11
##   thấp 316  35

ntnh <- table(std$Dalc_2,std$school)
OddsRatio(ntnh,conf.level = .95)

## odds ratio     lwr.ci     upr.ci 
##  0.3322785  0.1543865  0.7151466

Nhận xét:

Odds ratio của Dalc_2 và school là 0.3322785 nhỏ hơn 1. Điều này có nghĩa là tỷ lệ chênh của sinh viên tiêu thụ rượu cao tại trường GP so với trường MS là 0.332, tức là sinh viên tại trường GP có tỷ lệ tiêu thụ rượu cao thấp hơn so với sinh viên tại trường MS.
Khoảng tin cậy này cho thấy rằng, với mức độ tin cậy (thường là 95%), odds ratio thực sự nằm trong khoảng từ 0.1543865 đến 0.7151466. Vì khoảng tin cậy này không bao gồm giá trị 1, kết quả này có ý nghĩa thống kê.

3.5.2 Biến Dalc_2 và sex

table(std$Dalc_2,std$sex)

##       
##          F   M
##   cao    8  36
##   thấp 200 151

ntnh <- table(std$Dalc_2,std$sex)
OddsRatio(ntnh,conf.level = .95)

## odds ratio     lwr.ci     upr.ci 
## 0.16777778 0.07578856 0.37141995

Nhận xét:

Odds ratio của Dalc_2 và sex là 0.16777778 nhỏ hơn 1. Điều này có nghĩa là tỷ lệ chênh của sinh viên tiêu thụ rượu cao là nữ so với nam là 0.167, tức là sinh viên nữ có tỷ lệ tiêu thụ rượu cao thấp hơn so với sinh viên nam.
Khoảng tin cậy này cho thấy rằng, với mức độ tin cậy (thường là 95%), odds ratio thực sự nằm trong khoảng từ 0.07578856 đến 0.37141995. Vì khoảng tin cậy này không bao gồm giá trị 1, kết quả này có ý nghĩa thống kê.

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

library(tidyverse)
tn <- std %>% select(Dalc_2, school, sex)
tn <- tn %>% mutate(Dalc_2 = if_else(Dalc_2 == "thấp", 0, 1))
tn <- tn %>% mutate(school = case_when(
  school == "GP" ~ 0,
  school == "MS" ~ 1
))
tn <- tn %>% mutate(sex = case_when(
  sex == "F" ~ 0,
  sex == "M" ~ 1
))
head(tn)

3.6.1 Hồi quy giữa Dalc_2 và school

lpm_school <- lm(Dalc_2 ~ school, data = tn)
summary(lpm_school)

## 
## Call:
## lm(formula = Dalc_2 ~ school, data = tn)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.23913 -0.09456 -0.09456 -0.09456  0.90544 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.09456    0.01670   5.662 2.89e-08 ***
## school       0.14457    0.04894   2.954  0.00332 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.312 on 393 degrees of freedom
## Multiple R-squared:  0.02173,    Adjusted R-squared:  0.01924 
## F-statistic: 8.729 on 1 and 393 DF,  p-value: 0.003321

3.6.2 Hồi quy giữa Dalc_2 và sex

lpm_sex <- lm(Dalc_2 ~ sex, data = tn)
summary(lpm_sex)

## 
## Call:
## lm(formula = Dalc_2 ~ sex, data = tn)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.19251 -0.19251 -0.03846 -0.03846  0.96154 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.03846    0.02121   1.814   0.0705 .  
## sex          0.15405    0.03082   4.998 8.73e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3058 on 393 degrees of freedom
## Multiple R-squared:  0.05977,    Adjusted R-squared:  0.05738 
## F-statistic: 24.98 on 1 and 393 DF,  p-value: 8.728e-07

3.7 Mô hình hồi quy Logistic

3.7.1 Mô hình Logistic cho biến school

lpm_school <- glm(Dalc_2 ~ school, data = tn, family = binomial)
summary(lpm_school)

## 
## Call:
## glm(formula = Dalc_2 ~ school, family = binomial, data = tn)
## 
## Coefficients:
##             Estimate Std. Error z value Pr(>|z|)    
## (Intercept)  -2.2592     0.1829 -12.349  < 2e-16 ***
## school        1.1018     0.3911   2.817  0.00484 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 276.04  on 394  degrees of freedom
## Residual deviance: 269.05  on 393  degrees of freedom
## AIC: 273.05
## 
## Number of Fisher Scoring iterations: 5

3.7.2 Mô hình Logistic cho biến sex

lpm_sex <- glm(Dalc_2 ~ sex, data = tn, family = binomial)
summary(lpm_sex)

## 
## Call:
## glm(formula = Dalc_2 ~ sex, family = binomial, data = tn)
## 
## Coefficients:
##             Estimate Std. Error z value Pr(>|z|)    
## (Intercept)  -3.2189     0.3606  -8.928  < 2e-16 ***
## sex           1.7851     0.4055   4.403 1.07e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 276.04  on 394  degrees of freedom
## Residual deviance: 251.02  on 393  degrees of freedom
## AIC: 255.02
## 
## Number of Fisher Scoring iterations: 6

3.8 Mô hình Logit

3.8.1 Mô hình Logit cho biến school

reglogit <- glm (Dalc_2 ~ school, data = tn, family = binomial(link = 'logit'))
print(reglogit)

## 
## Call:  glm(formula = Dalc_2 ~ school, family = binomial(link = "logit"), 
##     data = tn)
## 
## Coefficients:
## (Intercept)       school  
##      -2.259        1.102  
## 
## Degrees of Freedom: 394 Total (i.e. Null);  393 Residual
## Null Deviance:       276 
## Residual Deviance: 269   AIC: 273

3.8.2 Mô hình Logit cho biến sex

3.9 Mô hình Probit

3.9.1 Mô hình Probit cho biến school

3.9.2 Mô hình Probit cho biến sex

regprobit <- glm (Dalc_2 ~ sex, data = tn, family = binomial(link = 'probit'))
print(regprobit)

## 
## Call:  glm(formula = Dalc_2 ~ sex, family = binomial(link = "probit"), 
##     data = tn)
## 
## Coefficients:
## (Intercept)          sex  
##     -1.7688       0.9002  
## 
## Degrees of Freedom: 394 Total (i.e. Null);  393 Residual
## Null Deviance:       276 
## Residual Deviance: 251   AIC: 255

(mục lớn nha) Đánh giá mô hình Chỉ số AIC

LS0tDQp0aXRsZTogIlRJ4buCVSBMVeG6rE4gTcOUTiBQSMOCTiBUw41DSCBE4buuIExJ4buGVSDEkOG7ik5IIFTDjU5IIg0Kb3V0cHV0Og0KICBodG1sX2RvY3VtZW50Og0KICAgIHRvYzogVFJVRQ0KICAgIHRvY19mbG9hdDogVFJVRQ0KICAgIGRmX3ByaW50OiBwYWdlZA0KICAgIGNvZGVfZG93bmxvYWQ6IHRydWUNCiAgICBjb2RlX2ZvbGRpbmc6IGhpZGUNCiAgcGRmX2RvY3VtZW50Og0KICAgIGV4dHJhX2RlcGVuZGVuY2llczoNCiAgICAgIHZpZXRuYW06IHV0ZjgNCiAgICB0b2M6IHllcw0KICAgIG51bWJlcl9zZWN0aW9uczogeWVzDQogIHdvcmRfZG9jdW1lbnQ6DQogICAgdG9jOiB5ZXMNCiAgICBudW1iZXJfc2VjdGlvbnM6IHllcw0KZ2VvbWV0cnk6DQogICAgICAtIGlubmVyPTNjbQ0KICAgICAgLSBvdXRlcj00Y20NCiAgICAgIC0gdG9wPTNjbQ0KICAgICAgLSBib3R0b209NGNtDQogICAgICAtIGhlYWRzZXA9MjJwdA0KICAgICAgLSBoZWFkaGVpZ2h0PTExcHQNCiAgICAgIC0gZm9vdHNraXA9MzNwdA0KICAgICAgLSBpZ25vcmVoZWFkDQogICAgICAtIGlnbm9yZWZvb3QNCiAgICAgIC0gaGVpZ2h0cm91bmRlZA0KLS0tDQoqKkdWSEQ6IFRoUy4gVHLhuqduIE3huqFuaCBUxrDhu51uZyoqDQoNCioqTmjDs20gc2luaCB2acOqbiB0aOG7sWMgaGnhu4duOioqIA0KDQpfMS4gTWFpIFRo4buLIEhvw6BpIFRoxrDGoW5nIC0gMjEyMTAwMTU0Nl8NCg0KXzIuIE5ndXnhu4VuIFRo4buLIE5n4buNYyBIw6AgLSAyMTIxMDEyNDc0Xw0KDQojIFBIw4JOIFTDjUNIIE5I4buuTkcgWeG6vlUgVOG7kCDhuqJOSCBIxq/hu55ORyDEkOG6vk4gTeG7qEMgxJDhu5ggVEnDilUgVEjhu6QgUsav4buiVSBD4bumQSBTSU5IIFZJw4pODQoNCkLhu5kgZOG7ryBsaeG7h3UgIlN0dWRlbnQgQWxjb2hvbCBDb25zdW1wdGlvbiBFREEiIHThu6sgS2FnZ2xlIGN1bmcgY+G6pXAgdGjDtG5nIHRpbiBjaGkgdGnhur90IHbhu4Egdmnhu4djIHNpbmggdmnDqm4gdGnDqnUgdGjhu6Ugcsaw4bujdSB04bqhaSBjw6FjIHRyxrDhu51uZyB0cnVuZyBo4buNYyDhu58gQuG7kyDEkMOgbyBOaGEuIELhu5kgZOG7ryBsaeG7h3UgbsOgeSBraMO0bmcgY2jhu4kgdOG6rXAgdHJ1bmcgdsOgbyBow6BuaCB2aSB0acOqdSB0aOG7pSByxrDhu6N1IG3DoCBjw7JuIGJhbyBn4buTbSBjw6FjIHRow7RuZyB0aW4gduG7gSBtw7RpIHRyxrDhu51uZywgZ2lhIMSRw6xuaCB2w6AgeMOjIGjhu5lpIGPhu6dhIHNpbmggdmnDqm4uIE3hu6VjIMSRw61jaCBj4bunYSBuZ2hpw6puIGPhu6l1IGzDoCDEkeG7gyBoaeG7g3UgcsO1IGjGoW4gduG7gSBt4buRaSBsacOqbiBo4buHIGdp4buvYSBjw6FjIHnhur91IHThu5EgbsOgeSB2w6AgdGjDs2kgcXVlbiB0acOqdSB0aOG7pSByxrDhu6N1IGPhu6dhIHNpbmggdmnDqm4gdHJvbmcgMSBuZ8OgeS4gROG7ryBsaeG7h3UgxJHGsOG7o2MgdGh1IHRo4bqtcCB04burIGPDoWMgY3Xhu5ljIGto4bqjbyBzw6F0IHRyb25nIGPhu5luZyDEkeG7k25nIGjhu41jIMSRxrDhu51uZywgbWFuZyBs4bqhaSBjxqEgaOG7mWkgxJHhu4MgbmdoacOqbiBj4bupdSB0w6FjIMSR4buZbmcgY+G7p2EgbcO0aSB0csaw4budbmcgdsOgIGdpYSDEkcOsbmggxJHhu5FpIHbhu5tpIGjDoG5oIHZpIHRpw6p1IHRo4bulIHLGsOG7o3Ug4bufIMSR4buZIHR14buVaSBzaW5oIHZpw6puLg0KDQojIEzhu5xJIEPhuqJNIMagTg0KDQpM4budaSDEkeG6p3UgdGnDqm4sIGNow7puZyBlbSB4aW4gZ+G7rWkgbOG7nWkgY+G6o20gxqFuIGNow6JuIHRow6BuaCDEkeG6v24gR2nhuqNuZyB2acOqbiAtIFRoUy4gVHLhuqduIE3huqFuaCBUxrDhu51uZy4gVHJvbmcgcXXDoSB0csOsbmggaOG7jWMgdOG6rXAgdsOgIHTDrG0gaGnhu4N1IG3DtG4gIlBow6JuIHTDrWNoIGThu68gbGnhu4d1IMSR4buLbmggdMOtbmgiLCBjaMO6bmcgZW0gxJHDoyBuaOG6rW4gxJHGsOG7o2Mgc+G7sSBxdWFuIHTDom0gZ2nDunAgxJHhu6EsIGjGsOG7m25nIGThuqtuIHThuq1uIHTDrG5oIHbDoCBj4buRbmcgaGnhur9uIHTDom0gaHV54bq/dCBj4bunYSB0aOG6p3kuIFRo4bqneSDEkcOjIGdpw7pwIGNow7puZyBlbSB0w61jaCBsxal5IHRow6ptIG5oaeG7gXUga2nhur9uIHRo4bupYyDEkeG7gyBjw7MgY8OhaSBuaMOsbiBob8OgbiB0aGnhu4duIHbDoCB04buVbmcgcXVhbiB24buBIGLhu5kgbcO0biBuw6B5LiBU4burIG5o4buvbmcga2nhur9uIHRo4bupYyDEkcaw4bujYyB0aOG6p3kgdHJ1eeG7gW4gdOG6o2ksIGNow7puZyBlbSBk4bqnbiBoaeG7g3UgxJHGsOG7o2MgdOG6p20gcXVhbiB0cuG7jW5nIGPhu6dhIG3DtG4gaOG7jWMgbcOgIGLhuqV5IGzDonUgbmF5IGPDsm4gbmhp4buBdSB0aOG6r2MgbeG6r2MuIFRow7RuZyBxdWEgYsOgaSB0aeG7g3UgbHXhuq1uIG7DoHksIGNow7puZyBlbSB4aW4gdHLDrG5oIGLDoHkgbmjhu69uZyBnw6wgY2jDum5nIGVtIMSRw6MgdMOsbSBoaeG7g3UgduG7gSBtw7RuIGjhu41jLCB2w6AgY+G7pSB0aOG7gyBsw6AgY2jDum5nIGVtIGzhu7FhIGNo4buNbiDEkeG7gSB0w6BpICJQaMOibiB0w61jaCBuaOG7r25nIHnhur91IHThu5Eg4bqjbmggaMaw4bufbmcgxJHhur9uIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgY+G7p2Egc2luaCB2acOqbiINCg0KS2nhur9uIHRo4bupYyBsw6AgdsO0IGjhuqFuIG5oxrBuZyBz4buxIHRp4bq/cCBuaOG6rW4ga2nhur9uIHRo4bupYyBj4bunYSBi4bqjbiB0aMOibiBt4buXaSBjaMO6bmcgZW0gbHXDtG4gaOG6oW4gY2jhur8gbeG7mXQgY8OhY2ggbmjhuqV0IMSR4buLbmguIERvIMSRw7MsIHRyb25nIHF1w6EgdHLDrG5oIGhvw6BuIHRow6BuaCBiw6BpIHRp4buDdSBsdeG6rW4sIGNo4bqvYyBjaOG6r24gY2jDum5nIGVtIGtow7RuZyB0csOhbmgga2jhu49pIG5o4buvbmcgdGhp4bq/dSBzw7N0LCB24bqteSBuw6puIGNow7puZyBlbSBy4bqldCBtb25nIG5o4bqtbiDEkcaw4bujYyBuaOG7r25nIMSRw7NuZyBnw7NwIMSR4bq/biB04burIHRo4bqneSDEkeG7gyBiw6BpIHRp4buDdSBsdeG6rW4gY+G7p2EgbmjDs20gxJHGsOG7o2MgaG/DoG4gdGhp4buHbiBoxqFuLg0KDQpDaMO6bmcgZW0ga8OtbmggY2jDumMgdGjhuqd5IHPhu6ljIGto4buPZSwgaOG6oW5oIHBow7pjIHbDoCB0aMOgbmggY8O0bmcgdHLDqm4gY29uIMSRxrDhu51uZyBz4buxIG5naGnhu4dwIGPhu6dhIG3DrG5oLg0KDQojIENIxq/GoE5HIDE6IFThu5RORyBRVUFOIE5HSEnDik4gQ+G7qFUNCg0KIyMgMS4xIEzDrSBkbyBjaOG7jW4gxJHhu4EgdMOgaQ0KDQpIaeG7h24gbmF5LCBow6BuaCB2aSB0acOqdSB0aOG7pSByxrDhu6N1IGzDoCBt4buZdCBoaeG7h24gdMaw4bujbmcgcGjhu5UgYmnhur9uIHRyb25nIMSR4budaSBz4buRbmcgc2luaCB2acOqbi4gTmhp4buBdSBzaW5oIHZpw6puIHRoYW0gZ2lhIGPDoWMgaG/huqF0IMSR4buZbmcgeMOjIGjhu5lpIG5oxrAgdGnhu4djIHTDuW5nLCBkw6Mgbmdv4bqhaSB2w6AgY8OhYyBz4buxIGtp4bq/biBj4bunYSB0csaw4budbmcsIG7GoWkgcsaw4bujdSB0aMaw4budbmcgxJHGsOG7o2MgcGjhu6VjIHbhu6UgdsOgIHRpw6p1IHRo4bulLiBOaGnhu4F1IHNpbmggdmnDqm4ga2jDtG5nIGNo4buJIHXhu5FuZyByxrDhu6N1IG3DoCBjw7JuIGzhuqFtIGThu6VuZyBuw7MsIHXhu5FuZyDEkeG6v24gbeG7qWMgc2F5IHjhu4luIHRoxrDhu51uZyB4dXnDqm4uIFTDrG5oIHRy4bqhbmcgbsOgeSDEkeG6t2MgYmnhu4d0IHBo4buVIGJp4bq/biB0cm9uZyBjw6FjIGThu4twIGzhu4UsIGvhu7Mgbmdo4buJIHbDoCBjw6FjIGJ14buVaSB0aeG7h2MgbOG7m24uIFZp4buHYyB14buRbmcgcsaw4bujdSBzYXkgeOG7iW4gY8OzIHRo4buDIGThuqtuIMSR4bq/biBt4bqldCBraeG7g20gc2/DoXQgaMOgbmggdmkgdsOgIGfDonkgcmEgbmhp4buBdSBo4bqtdSBxdeG6oyB0acOqdSBj4buxYyBuaMawIHRhaSBu4bqhbiwgYuG6oW8gbOG7sWMgdsOgIGPDoWMgduG6pW4gxJHhu4EgduG7gSBz4bupYyBraOG7j2UuIFLGsOG7o3Uga2jDtG5nIGNo4buJIGzDoCBkw7luZyBjaG8gbmjhu69uZyBjdeG7mWMgdnVpIGNoxqFpLCB0aeG7h2MgdMO5bmcsIGdpYW8gbMawdSBi4bqhbiBiw6ggbcOgIHNpbmggdmnDqm4gaGnhu4duIG5heSBjw7JuIGTDuW5nIHLGsOG7o3UgxJHhu4MgZ2nhuqNpIHThu49hIG7hu5dpIGJ14buTbiwgbmjhu69uZyDDoXAgbOG7sWMgdHJvbmcgaOG7jWMgdOG6rXAsIGLhuqFuIGLDqCwgZ2lhIMSRw6xuaCwgdMOsbmggY+G6o20gbsOqbiB2aeG7h2MgdGnDqnUgdGjhu6Ugcsaw4bujdSBj4bunYSBzaW5oIHZpw6puIGhp4buHbiBuYXkgcuG6pXQgxJHDoW5nIGLDoW8gxJHhu5luZy4gUsaw4bujdSBy4bqldCBk4buFIGTDoG5nIMSRxrDhu6NjIHRp4bq/cCBj4bqtbiB2w6Agbmjhu69uZyBjaOG7lyBsxrB1IHRow7RuZyByxrDhu6N1IGtow7RuZyDEkcaw4bujYyBnacOhbSBzw6F0IGNo4bq3dCBjaOG6vSBuw6puIG5nw6B5IGPDoG5nIGzDoG0gZ2lhIHTEg25nIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UuVGjhu7FjIHRy4bqhbmcgdGnDqnUgdGjhu6Ugcsaw4bujdSBj4bunYSBzaW5oIHZpw6puIGhp4buHbiBuYXkgbMOgIG3hu5l0IHbhuqVuIMSR4buBIG5naGnDqm0gdHLhu41uZyBj4bqnbiDEkcaw4bujYyBxdWFuIHTDom0uIFZp4buHYyBuZ2hpw6puIGPhu6l1IHbDoCBoaeG7g3UgcsO1IGjGoW4gduG7gSBow6BuaCB2aSB0acOqdSB0aOG7pSByxrDhu6N1IGPhu6dhIHNpbmggdmnDqm4gdHJvbmcgbeG7mXQgbmfDtGkgdHLGsOG7nW5nIOG7nyBC4buTIMSQw6BvIE5oYSBz4bq9IGdpw7pwIMSRxrBhIHJhIMSRxrDhu6NjIG3hu5l0IHbDoGkgYmnhu4duIHBow6FwIGdpw6FvIGThu6VjLCBwaMOybmcgbmfhu6thIHbDoCBo4buXIHRy4bujIHRow61jaCBo4bujcCwgbmjhurFtIGdp4bqjbSB0aGnhu4N1IG5o4buvbmcgaOG6rXUgcXXhuqMgdGnDqnUgY+G7sWMgdsOgIGLhuqNvIHbhu4cgc+G7qWMga2jhu49lIHbDoCB0xrDGoW5nIGxhaSBj4bunYSBzaW5oIHZpw6puLg0KDQojIyAxLjIgTeG7pWMgdGnDqnUgbmdoacOqbiBj4bupdQ0KDQpN4bulYyB0acOqdSBj4bunYSBjaMO6bmcgdMO0aSBsw6AgcGjDom4gdMOtY2ggY8OhYyB54bq/dSB04buRIOG6o25oIGjGsOG7n25nIMSR4bq/biBt4bupYyDEkeG7mSB0acOqdSB0aOG7pSByxrDhu6N1IGPhu6dhIHNpbmggdmnDqm4gbMOgIHRo4bqlcCBoYXkgY2FvLiBOZ2hpw6puIGPhu6l1IG5o4bqxbSB0w6xtIHJhIGPDoWMgdGjDtG5nIHRpbiB2w6AgxJHGsGEgcmEgbeG7mXQgYmnhu4duIHBow6FwIGdpw7pwIGdpYSDEkcOsbmggdsOgIG5ow6AgdHLGsOG7nW5nIHF14bqjbiBsw70gY29uIGPDoWksIHNpbmggdmnDqm4gY+G7p2EgbcOsbmggxJHGsOG7o2MgdOG7kXQgaMahbi4NCg0KIyMgMS4zIMSQ4buRaSB0xrDhu6NuZyB2w6AgcGjhuqFtIHZpIG5naGnDqm4gY+G7qXUNCg0KVOG6rXAgZOG7ryBsaeG7h3UgY2jhu6lhIHRow7RuZyB0aW4gdOG7qyBjdeG7mWMga2jhuqNvIHPDoXQgduG7gSBzaW5oIHZpw6puIGtow7NhIGjhu41jIFRvw6FuIOG7nyAyIHRyxrDhu51uZyB0cnVuZyBo4buNYyDhu58gQuG7kyDEkMOgbyBOaGEuIELhu5kgZOG7ryBsaeG7h3UgY2jGsGEgY8OhYyB0aMO0bmcgdGluIGxpw6puIHF1YW4gxJHhur9uIGPDoWMgc2luaCB2acOqbiBiYW8gZ+G7k20gMzk1IHF1YW4gc8OhdCB0csOqbiB04buVbmcgOSBiaeG6v24uDQoNCiMjIDEuNCBQaMawxqFuZyBwaMOhcCBuZ2hpw6puIGPhu6l1DQoNCkLDoGkgbmdoacOqbiBj4bupdSBz4butIGThu6VuZyBwaMawxqFuZyBwaMOhcCB0aOG7kW5nIGvDqiB2w6AgY8OhYyBtw7QgaMOsbmggaOG7k2kgcXV5IMSR4buDIHBow6JuIHTDrWNoIGThu68gbGnhu4d1IGNobyB0aOG6pXkgdMOhYyDEkeG7mW5nIGPhu6dhIGPDoWMgeeG6v3UgdOG7kSDhuqNuaCBoxrDhu59uZyDEkeG6v24gbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBj4bunYSBzaW5oIHZpw6puLg0KDQojIyAxLjUgS+G6v3QgY+G6pXUgYsOgaSBuZ2hpw6puIGPhu6l1DQoNCi0gQ2jGsMahbmcgMTogVOG7lW5nIHF1YW4gbmdoacOqbiBj4bupdS4NCg0KLSBDaMawxqFuZyAyOiBE4buvIGxp4buHdSBuZ2hpw6puIGPhu6l1Lg0KDQotIENoxrDGoW5nIDM6IFBow6JuIHTDrWNoIGThu68gbGnhu4d1Lg0KDQotIENoxrDGoW5nIDQ6IEvhur90IHF14bqjIG3DtCBow6xuaCBuZ2hpw6puIGPhu6l1Lg0KDQotIENoxrDGoW5nIDU6IEvhur90IGx14bqtbiB2w6Aga2h1eeG6v24gbmdo4buLLg0KDQojIENoxrDGoW5nIDI6IEThu64gTEnhu4ZVIE5HSEnDik4gQ+G7qFUNCg0KIyMgMi4xIEdp4bubaSB0aGnhu4d1IGLhu5kgZOG7ryBsaeG7h3UNCg0KQuG7mSBk4buvIGxp4buHdSBiYW8gZ+G7k20gMzk1IHF1YW4gc8OhdCB0csOqbiB04buVbmcgOSBiaeG6v24gY+G7p2Egc2luaCB2acOqbiBraMOzYSBo4buNYyBUb8OhbiBj4bunYSAyIHRyxrDhu51uZyB0cnVuZyBo4buNYyDhu58gQuG7kyDEkMOgbyBOaGEuIETGsOG7m2kgxJHDonkgbMOgIGLhuqNuZyBjaGkgdGnhur90IHThu6tuZyBiaeG6v24gdHLDqm4gdOG7q25nIHF1YW4gc8OhdCBraMOhYyBuaGF1Lg0KYGBge3J9DQpzdGQgPC0gcmVhZC5jc3YoInN0dWRlbnQtbWF0aC5jc3YiLCBoZWFkZXIgPSBUKQ0KbGlicmFyeShkYXRhLnRhYmxlKQ0KZGF0YS50YWJsZShzdGQpDQpgYGANCg0KIyMgMi4yIE3DtCB04bqjIGPDoWMgYmnhur9uIHRyb25nIGLhu5kgZOG7ryBsaeG7h3UNCg0KLSBzY2hvb2w6IG5nw7RpIHRyxrDhu51uZyBzaW5oIHZpw6puIMSRYW5nIHRoZW8gaOG7jWMgKEdQOiBHYWJyaWVsIFBlcmVpcmE7IE1TOiBNb3VzaW5obyBkYSBTaWx2ZXJpYSkNCg0KLSBzZXg6IGdp4bubaSB0w61uaCBj4bunYSBzaW5oIHZpw6puIChGOiBu4buvOyBNOiBuYW0pDQoNCi0gYWRkcmVzczogbsahaSDhu58gY+G7p2Egc2luaCB2acOqbiAoVTogdGjDoG5oIHRo4buLOyBSOiBuw7RuZyB0aMO0bikNCg0KLSBzY2hvb2xzdXA6IGjhu5cgdHLhu6MgZ2nDoW8gZOG7pWMgYuG7lSBzdW5nICh5ZXM6IGPDszsgbm86IGtow7RuZykNCg0KLSBmYW1zdXA6IGjhu5cgdHLhu6MgZ2nDoW8gZOG7pWMgZ2lhIMSRw6xuaCAoeWVzOiBjw7M7IG5vOiBraMO0bmcpDQoNCi0gYWN0aXZpdGllczogaG/huqF0IMSR4buZbmcgbmdv4bqhaSBraMOzYSAoeWVzOiBjw7M7IG5vOiBraMO0bmcpDQoNCi0gaW50ZXJuZXQ6IHRydXkgY+G6rXAgbeG6oW5nIHThuqFpIG5ow6AgKHllczogY8OzOyBubzoga2jDtG5nKQ0KDQotIHJvbWFudGljOiBt4buRaSBxdWFuIGjhu4cgdMOsbmggY+G6o20gKHllczogY8OzOyBubzoga2jDtG5nKQ0KDQotIERhbGM6IG3hu6ljIHRpw6p1IHRo4bulIHLGsOG7o3UgdHJvbmcgMSBuZ8OgeSAoMTogcuG6pXQgdGjhuqVwOyAyOiB0aOG6pXA7IDM6IHRydW5nIGLDrG5oOyA0OiBjYW87IDU6IHLhuqV0IGNhbykNCg0KIyBDaMawxqFuZyAzOiBQSMOCTiBUw41DSCBE4buuIExJ4buGVQ0KDQpUcm9uZyBiw6BpIGx14bqtbiBuw6B5LCBiaeG6v24gcGjhu6UgdGh14buZYyDEkcaw4bujYyBjaOG7jW4gbMOgIGJp4bq/biBEYWxjICht4bupYyDEkeG7mSB0acOqdSB0aOG7pSByxrDhu6N1IGPhu6dhIHNpbmggdmnDqm4gdHJvbmcgMSBuZ8OgeSksIGRvIGJp4bq/biBEYWxjIGzDoCBk4buvIGxp4buHdSDEkeG7i25oIGzGsOG7o25nIHbhu5tpIGPDoWMgbeG7qWMgxJHhu5kgMSwgMiwgMywgNCwgNSBuw6puIG5ow7NtIGJp4bq/biBuw7MgdGjDoG5oIGThu68gbGnhu4d1IMSR4buLbmggdMOtbmggbMOgIERhbGNfMiB24bubaSAyIGJp4buDdSBoaeG7h24gY+G7pSB0aOG7gyBsw6AgdGjhuqVwIHbDoCBjYW8uIOG7niBjw6FjIG3hu6ljIMSR4buZIDEsIDIgdsOgIDMgc+G6vSBiaeG6v24gxJHhu5VpIHRow6BuaCBiaeG7g3UgaGnhu4duIHRpw6p1IHRo4bulIHLGsOG7o3UgdGjhuqVwLCBjw7JuIG3hu6ljIMSR4buZIDQgdsOgIDUgc+G6vSBiaeG6v24gxJHhu5VpIHRow6BuaCBiaeG7g3UgaGnhu4duIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIG5oxrAgYuG6o25nIGTGsOG7m2kgxJHDonk6DQpgYGB7cn0NCnN0ZCA8LSByZWFkLmNzdigic3R1ZGVudC1tYXRoLmNzdiIsIGhlYWRlciA9IFQpDQpzdGQkRGFsY18yIDwtIGlmZWxzZShzdGQkRGFsYyA8PSAyLCAidGjhuqVwIiwgImNhbyIpDQpwcmludChzdGQpDQpgYGANCg0KIyMgMy4xIFRo4buRbmcga8OqIG3DtCB04bqjIGJp4bq/biBwaOG7pSB0aHXhu5ljDQoNCiMjIyAzLjEuMSBM4bqtcCBi4bqjbmcgdOG6p24gc+G7kSB2w6AgYuG6o25nIHThuqduIHN14bqldA0KYGBge3J9DQp0YWJsZShzdGQkRGFsY18yKQ0KdGFibGUoc3RkJERhbGNfMikvbGVuZ3RoKHN0ZCREYWxjXzIpDQpgYGANCg0KIyMjIDMuMS4yIEzhuq1wIGJp4buDdSDEkeG7kw0KYGBge3J9DQpsaWJyYXJ5KGRwbHlyKQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpzdGQgJT4lIGdyb3VwX2J5KERhbGNfMikgJT4lICBzdW1tYXJpc2UobiA9IG4oKSkgJT4lIA0KICAgIGdncGxvdChhZXMoeD1EYWxjXzIseT1uKSkgKyANCiAgICBnZW9tX2NvbChmaWxsID0gIiM4OUNGRjAiKSsNCiAgICBnZW9tX3RleHQoYWVzKGxhYmVsID1uKSx2anVzdCA9IDEsIGNvbG9yID0gIndoaXRlIikrDQogICAgdGhlbWVfYncoKQ0KYGBgDQoNCl9OaOG6rW4geMOpdDpfIEJp4buDdSDEkeG7kyBj4buZdCB2w6AgYuG6o25nIOG7nyB0csOqbiB0aOG7gyBoaeG7h24gbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBj4bunYSBzaW5oIHZpw6puIHRyb25nIDEgbmfDoHkuIFRoZW8ga2jhuqNvIHPDoXQgdGjDrCBt4bupYyB0acOqdSB0aOG7pSByxrDhu6N1IGPhu6dhIHNpbmggdmnDqm4gdHJvbmcgMSBuZ8OgeSB24bubaSBt4bupYyDEkeG7mSB0acOqdSB0aOG7pSByxrDhu6N1IGNhbyBjaOG7iSBjw7MgNDQgc2luaCB2acOqbiBjaGnhur9tIHThu7cgbOG7hyAxMS4xNCUsIGPDsm4gbOG6oWkgMzUxIHNpbmggdmnDqm4gduG7m2kgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSB0aOG6pXAgY2hp4bq/bSB04bu3IGzhu4cgODguODYlLiBOaMawIGvhur90IHF14bqjIGNobyB0aOG6pXksIHThuqFpIDIgbmfDtGkgdHLGsOG7nW5nIG7DoHkgdGjDrCBzaW5oIHZpw6puIHRpw6p1IHRo4bulIHLGsOG7o3UgdHJvbmcgMSBuZ8OgeSBy4bqldCBsw6Agw610LiANCg0KQmnhu4N1IMSR4buTIHRyw7JuIGTGsOG7m2kgxJHDonkgY2hvIHRo4bqleSByw7UgaMahbiB24buBIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgY+G7p2Egc2luaCB2acOqbiB0cm9uZyAxIG5nw6B5Lg0KYGBge3J9DQphcCA8LSBwcm9wLnRhYmxlKHRhYmxlKHN0ZCREYWxjXzIpKSAqIDEwMA0KZ2dwbG90KG1hcHBpbmcgPSBhZXMoeCA9ICIiLCB5ID0gYXAsIGZpbGwgPSBmYWN0b3IobmFtZXMoYXApKSkpICsNCiAgZ2VvbV9iYXIod2lkdGggPSAxLCBzdGF0ID0gImlkZW50aXR5IikgKw0KICBjb29yZF9wb2xhcigieSIpICsNCiAgbGFicyh0aXRsZSA9ICJQaMOibiBi4buRIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgdHJvbmcgMSBuZ8OgeSIsDQogICAgICAgZmlsbCA9ICJN4bupYyB0acOqdSB0aOG7pSIpICsNCiAgdGhlbWVfdm9pZCgpICsNCiAgZ2VvbV90ZXh0KGFlcyhsYWJlbCA9IHBhc3RlMChyb3VuZChhcCwgMSksICIlIikpLA0KICAgICAgICAgICAgcG9zaXRpb24gPSBwb3NpdGlvbl9zdGFjayh2anVzdCA9IDAuNSkpDQoNCmBgYA0KDQojIyAzLjIgVGjhu5FuZyBrw6ogc3V5IGRp4buFbg0KDQojIyMgMy4yLjEgS2nhu4NtIMSR4buLbmggdMOtbmggxJHhu5ljIGzhuq1wIA0KDQpQaMawxqFuZyBwaMOhcCBDaGkgYsOsbmggcGjGsMahbmcNCg0KJCRcY2hpXjIgPSBcc3VtX3tpLGp9IFxmcmFjeyhuX3tpan0gLVxoYXR7XG11fV97aWp9KV4yfXtcaGF0e1xtdX1fe2lqfX0kJA0KDQpW4bubaSAkbl97aSxqfSQgbMOgIGdpw6EgdHLhu4sgY+G7p2Egw7QgJGksaixcaGF0e1xtdX1fe2lqfT1cZnJhY3tuX3tpK30gXGNkb3Qgbl97K2p9fXtufSQNCg0KIyMjIyAzLjIuMS4xIEJp4bq/biBEYWxjXzIgdsOgIHNjaG9vbA0KDQpHaeG6oyB0aHV54bq/dDoNCg0KJEhfMCQ6IERhbGNfMiB2w6Agc2Nob29sIMSR4buZYyBs4bqtcCANCg0KJEhfMSQ6IERhbGNfMiB2w6Agc2Nob29sIGtow7RuZyDEkeG7mWMgbOG6rXANCmBgYHtyfQ0KbXRodCA8LSB0YWJsZShzdGQkRGFsY18yLCBzdGQkc2Nob29sKQ0KY2hpc3EudGVzdChtdGh0KQ0KYGBgDQpfTmjhuq1uIHjDqXQ6XyBRdWEga+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmggY2hvIHRhIHDiiJJ2YWx1ZSA9IDAuMDA3MzU2IDwgMC4wNSwgYsOhYyBi4buPICRIXzAkLCBuZ2jEqWEgbMOgIGdp4buvYSB2aeG7h2Mgc2luaCB2acOqbiBo4buNYyDhu58gdHLGsOG7nW5nIG7DoG8gdsOgIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgbMOgIGPDsyBsacOqbiBxdWFuIHbhu5tpIG5oYXUuDQoNCiMjIyMgMy4yLjEuMiBCaeG6v24gRGFsY18yIHbDoCBzZXgNCg0KR2nhuqMgdGh1eeG6v3Q6DQoNCiRIXzAkOiBEYWxjXzIgdsOgIHNleCDEkeG7mWMgbOG6rXAgDQoNCiRIXzEkOiBEYWxjXzIgdsOgIHNleCBraMO0bmcgxJHhu5ljIGzhuq1wDQpgYGB7cn0NCm10aHQgPC0gdGFibGUoc3RkJERhbGNfMiwgc3RkJHNleCkNCmNoaXNxLnRlc3QobXRodCkNCmBgYA0KDQpfTmjhuq1uIHjDqXQ6XyBRdWEga+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmggY2hvIHRhIHDiiJJ2YWx1ZSA9IDIuNjE4ZS0wNiA8IDAuMDUsIGLDoWMgYuG7jyAkSF8wJCwgbmdoxKlhIGzDoCBnaeG7r2EgZ2nhu5tpIHTDrW5oIGPhu6dhIHNpbmggdmnDqm4gdsOgIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgbMOgIGPDsyBsacOqbiBxdWFuIHbhu5tpIG5oYXUuDQoNCiMjIyMgMy4yLjEuMyBCaeG6v24gRGFsY18yIHbDoCBhZGRyZXNzDQoNCkdp4bqjIHRodXnhur90Og0KDQokSF8wJDogRGFsY18yIHbDoCBhZGRyZXNzIMSR4buZYyBs4bqtcCANCg0KJEhfMSQ6IERhbGNfMiB2w6AgYWRkcmVzcyBraMO0bmcgxJHhu5ljIGzhuq1wDQpgYGB7cn0NCm10aHQgPC0gdGFibGUoc3RkJERhbGNfMiwgc3RkJGFkZHJlc3MpDQpjaGlzcS50ZXN0KG10aHQpDQpgYGANCg0KX05o4bqtbiB4w6l0Ol8gUXVhIGvhur90IHF14bqjIGtp4buDbSDEkeG7i25oIGNobyB0YSBw4oiSdmFsdWUgPSAwLjA3MTAxID4gMC4wNSwgY2jGsGEgY8OzIGPGoSBz4bufIGLDoWMgYuG7jyAkSF8wJCwgbmdoxKlhIGzDoCBnaeG7r2EgbsahaSDhu58gY+G7p2Egc2luaCB2acOqbiB2w6AgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBraMO0bmcgY8OzIGxpw6puIHF1YW4gduG7m2kgbmhhdS4NCg0KIyMjIyAzLjIuMS40IEJp4bq/biBEYWxjXzIgdsOgIHNjaG9vbHN1cA0KDQpHaeG6oyB0aHV54bq/dDoNCg0KJEhfMCQ6IERhbGNfMiB2w6Agc2Nob29sc3VwIMSR4buZYyBs4bqtcCANCg0KJEhfMSQ6IERhbGNfMiB2w6Agc2Nob29sc3VwIGtow7RuZyDEkeG7mWMgbOG6rXANCmBgYHtyfQ0KbXRodCA8LSB0YWJsZShzdGQkRGFsY18yLCBzdGQkc2Nob29sc3VwKQ0KY2hpc3EudGVzdChtdGh0KQ0KYGBgDQoNCl9OaOG6rW4geMOpdDpfIFF1YSBr4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaCBjaG8gdGEgcOKIknZhbHVlID0gMC41NzMzID4gMC4wNSwgY2jGsGEgY8OzIGPGoSBz4bufIGLDoWMgYuG7jyAkSF8wJCwgbmdoxKlhIGzDoCBnaeG7r2Egdmnhu4djIGjhu5cgdHLhu6MgZ2nDoW8gZOG7pWMgYuG7lSBzdW5nIGNobyBzaW5oIHZpw6puIHbDoCBt4bupYyDEkeG7mSB0acOqdSB0aOG7pSByxrDhu6N1IGtow7RuZyBjw7MgbGnDqm4gcXVhbiB24bubaSBuaGF1Lg0KDQojIyMjIDMuMi4xLjUgQmnhur9uIERhbGNfMiB2w6AgZmFtc3VwDQoNCkdp4bqjIHRodXnhur90Og0KDQokSF8wJDogRGFsY18yIHbDoCBmYW1zdXAgxJHhu5ljIGzhuq1wIA0KDQokSF8xJDogRGFsY18yIHbDoCBmYW1zdXAga2jDtG5nIMSR4buZYyBs4bqtcA0KYGBge3J9DQptdGh0IDwtIHRhYmxlKHN0ZCREYWxjXzIsIHN0ZCRmYW1zdXApDQpjaGlzcS50ZXN0KG10aHQpDQpgYGANCg0KX05o4bqtbiB4w6l0Ol8gUXVhIGvhur90IHF14bqjIGtp4buDbSDEkeG7i25oIGNobyB0YSBw4oiSdmFsdWUgPSAwLjI1NjQgPiAwLjA1LCBjaMawYSBjw7MgY8ahIHPhu58gYsOhYyBi4buPICRIXzAkLCBuZ2jEqWEgbMOgIGdp4buvYSB2aeG7h2MgaOG7lyB0cuG7oyBnacOhbyBk4bulYyBnaWEgxJHDrG5oICB2w6AgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBraMO0bmcgY8OzIGxpw6puIHF1YW4gduG7m2kgbmhhdS4NCg0KIyMjIyAzLjIuMS42IEJp4bq/biBEYWxjXzIgdsOgIGFjdGl2aXRpZXMNCg0KR2nhuqMgdGh1eeG6v3Q6DQoNCiRIXzAkOiBEYWxjXzIgdsOgIGFjdGl2aXRpZXMgxJHhu5ljIGzhuq1wIA0KDQokSF8xJDogRGFsY18yIHbDoCBhY3Rpdml0aWVzIGtow7RuZyDEkeG7mWMgbOG6rXANCmBgYHtyfQ0KbXRodCA8LSB0YWJsZShzdGQkRGFsY18yLCBzdGQkYWN0aXZpdGllcykNCmNoaXNxLnRlc3QobXRodCkNCmBgYA0KDQpfTmjhuq1uIHjDqXQ6XyBRdWEga+G6v3QgcXXhuqMga2nhu4NtIMSR4buLbmggY2hvIHRhIHDiiJJ2YWx1ZSA9IDAuMjEzNCA+IDAuMDUsIGNoxrBhIMSR4bunIGPGoSBz4bufIGLDoWMgYuG7jyAkSF8wJCwgbmdoxKlhIGzDoCBnaeG7r2Egdmnhu4djIGhv4bqhdCDEkeG7mW5nIG5nb+G6oWkga2jDs2EgY+G7p2Egc2luaCB2acOqbiB2w6AgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBraMO0bmcgY8OzIGxpw6puIHF1YW4gduG7m2kgbmhhdS4NCg0KIyMjIyAzLjIuMS43IEJp4bq/biBEYWxjXzIgdsOgIGludGVybmV0DQoNCkdp4bqjIHRodXnhur90Og0KDQokSF8wJDogRGFsY18yIHbDoCBpbnRlcm5ldCDEkeG7mWMgbOG6rXAgDQoNCiRIXzEkOiBEYWxjXzIgdsOgIGludGVybmV0IGtow7RuZyDEkeG7mWMgbOG6rXANCmBgYHtyfQ0KbXRodCA8LSB0YWJsZShzdGQkRGFsY18yLCBzdGQkaW50ZXJuZXQpDQpjaGlzcS50ZXN0KG10aHQpDQpgYGANCg0KX05o4bqtbiB4w6l0Ol8gUXVhIGvhur90IHF14bqjIGtp4buDbSDEkeG7i25oIGNobyB0YSBw4oiSdmFsdWUgPSAxID4gMC4wNSwgY2jGsGEgxJHhu6cgY8ahIHPhu58gYsOhYyBi4buPICRIXzAkLCBuZ2jEqWEgbMOgIGdp4buvYSB2aeG7h2Mgc+G7rSBk4bulbmcgbeG6oW5nIHThuqFpIG5ow6AgY+G7p2Egc2luaCB2acOqbiB2w6AgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBraMO0bmcgY8OzIGxpw6puIHF1YW4gduG7m2kgbmhhdS4NCg0KIyMjIyAzLjIuMS42IEJp4bq/biBEYWxjXzIgdsOgIHJvbWFudGljcw0KDQpHaeG6oyB0aHV54bq/dDoNCg0KJEhfMCQ6IERhbGNfMiB2w6Agcm9tYW50aWMgxJHhu5ljIGzhuq1wIA0KDQokSF8xJDogRGFsY18yIHbDoCBzZXggcm9tYW50aWMgxJHhu5ljIGzhuq1wDQpgYGB7cn0NCm10aHQgPC0gdGFibGUoc3RkJERhbGNfMiwgc3RkJHJvbWFudGljKQ0KY2hpc3EudGVzdChtdGh0KQ0KYGBgDQoNCl9OaOG6rW4geMOpdDpfIFF1YSBr4bq/dCBxdeG6oyBraeG7g20gxJHhu4tuaCBjaG8gdGEgcOKIknZhbHVlID0gMC45NDQ5ID4gMC4wNSwgY2jGsGEgxJHhu6cgY8ahIHPhu58gYsOhYyBi4buPICRIXzAkLCBuZ2jEqWEgbMOgIGdp4buvYSBt4buRaSBxdWFuIGjhu4cgdMOsbmggY+G6o20gY+G7p2Egc2luaCB2acOqbiB2w6AgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBraMO0bmcgY8OzIGxpw6puIHF1YW4gduG7m2kgbmhhdS4NCg0KS+G6v3QgbHXhuq1uOiBDaOG7iSBjw7MgYmnhur9uIHNjaG9vbCAodHLGsOG7nW5nIGjhu41jKSB2w6Agc2V4IChnaeG7m2kgdMOtbmgpIGzDoCBjw7Mg4bqjbmggaMaw4bufbmcgdOG7m2kgYmnhur9uIERhbGNfMi4gVGnhur9wIHRoZW8sIG5ow7NtIHRp4bq/biBow6BuaCBtw7QgdOG6oyAyIGJp4bq/biwgdOG7qyDEkcOzIMSRxrBhIHJhIGPDoWkgbmjDrG4ga2jDoWkgcXXDoXQgaMahbiB24buBIG3hu6ljIMSR4buZIOG6o25oIGjGsOG7n25nIGzDqm4gYmnhur9uIERhbGNfMi4NCg0KIyMgMy4zIFRo4buRbmcga8OqIG3DtCB04bqjIDIgYmnhur9uDQoNCiMjIyAzLjMuMSBCaeG6v24gRGFsY18yIHbDoCBzZXgNCg0KIyMjIyAzLjMuMS4xIEzhuq1wIGLhuqNuZyB04bqnbiBz4buRIHbDoCBi4bqjbmcgdOG6p24gc3XhuqV0DQpgYGB7cn0NCmh0bmggPC0gdGFibGUoc3RkJERhbGNfMixzdGQkc2Nob29sKQ0KYWRkbWFyZ2lucyhodG5oKQ0KaHRuaCA8LSBwcm9wLnRhYmxlKGh0bmgpDQphZGRtYXJnaW5zKGh0bmgpDQpgYGANCg0KIyMjIyAzLjMuMS4yIEzhuq1wIGJp4buDdSDEkeG7kw0KYGBge3J9DQpzdGQgJT4lIGdyb3VwX2J5KERhbGNfMixzY2hvb2wpICU+JSAgc3VtbWFyaXNlKG4gPSBuKCkpICU+JSANCiAgICAgZ2dwbG90KGFlcyh4PURhbGNfMix5PW4sIGZpbGwgPSBzY2hvb2wpKSArIA0KICAgICBnZW9tX2NvbChwb3NpdGlvbiA9ICdkb2RnZScpKw0KICAgICB0aGVtZV9idygpDQpgYGANCg0KX05o4bqtbiB4w6l0Ol8gQuG6o25nIHbDoCBiaeG7g3UgxJHhu5MgdHLDqm4gY2hvIHRo4bqleSDhu58gdHLGsOG7nW5nIEdQIGPDsyAzMyBzaW5oIHZpw6puIGPDsyBt4bupYyB0acOqdSB0aOG7pSByxrDhu6N1IGNhbyBjaGnhur9tIHThu7cgbOG7hyA4LjM1JSwgY8OybiBs4bqhaSAzMTYgc2luaCB2acOqbiBjw7MgbeG7qWMgdGnDqnUgdGjhu6Ugcsaw4bujdSB0aOG6pXAgY2hp4bq/bSB04bu3IGzhu4cgODAlOyDhu58gdHLGsOG7nW5nIE1TIGPDsyAxMSBzaW5oIHZpw6puIGPDsyBt4bupYyB0acOqdSB0aOG7pSByxrDhu6N1IGNhbyBjaGnhur9tIHThu7cgbOG7hyAyLjc4JSwgY8OybiBs4bqhaSAzNSBzaW5oIHZpw6puIGPDsyBt4bupYyB0acOqdSB0aOG7pSByxrDhu6N1IHRo4bqlcCBjaGnhur9tIDguODYlLiBUw7NtIGzhuqFpIHRow6wgY+G6oyBoYWkgdHLGsOG7nW5nIMSR4buBdSBjw7Mgc+G7kSBsxrDhu6NuZyDDrXQgc2luaCB2acOqbiBjw7MgbeG7qWMgdGnDqnUgdGjhu6Ugcsaw4bujdSBjYW8sIMSRYSBz4buRIMSR4buBdSB0acOqdSB0aOG7pSByxrDhu6N1IHRo4bqlcCBob+G6t2MgaOG6p3UgbmjGsCBsw6Aga2jDtG5nIHRpw6p1IHRo4bulLg0KDQojIyMgMy4zLjIgQmnhur9uIERhbGNfMiB2w6Agc2V4DQoNCiMjIyMgMy4zLjIuMSBM4bqtcCBi4bqjbmcgdOG6p24gc+G7kSB2w6AgYuG6o25nIHThuqduIHN14bqldA0KYGBge3J9DQpodG5oIDwtIHRhYmxlKHN0ZCREYWxjXzIsc3RkJHNleCkNCmFkZG1hcmdpbnMoaHRuaCkNCmh0bmggPC0gcHJvcC50YWJsZShodG5oKQ0KYWRkbWFyZ2lucyhodG5oKQ0KYGBgDQoNCiMjIyMgMy4zLjIuMiBM4bqtcCBiaeG7g3UgxJHhu5MNCmBgYHtyfQ0Kc3RkICU+JSBncm91cF9ieShEYWxjXzIsc2V4KSAlPiUgIHN1bW1hcmlzZShuID0gbigpKSAlPiUgDQogICAgIGdncGxvdChhZXMoeD1EYWxjXzIseT1uLCBmaWxsID0gc2V4KSkgKyANCiAgICAgZ2VvbV9jb2wocG9zaXRpb24gPSAnZG9kZ2UnKSsNCiAgICAgdGhlbWVfYncoKQ0KYGBgDQoNCl9OaOG6rW4geMOpdDpfIELhuqNuZyB2w6AgYmnhu4N1IMSR4buTIHRyw6puIGNobyB0aOG6pXkg4bufIGdp4bubaSB0w61uaCBu4buvIGPDsyA4IHNpbmggdmnDqm4gY8OzIG3hu6ljIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIGNoaeG6v20gdOG7tyBs4buHIDIuMDMlLCBjw7JuIGzhuqFpIDIwMCBzaW5oIHZpw6puIG7hu68gY8OzIG3hu6ljIHRpw6p1IHRo4bulIHLGsOG7o3UgdGjhuqVwIGNoaeG6v20gdOG7tyBs4buHIDUwLjYzJSU7IOG7nyBnaeG7m2kgdMOtbmggbmFtIGPDsyAzNiBzaW5oIHZpw6puIGPDsyBt4bupYyB0acOqdSB0aOG7pSByxrDhu6N1IGNhbyBjaGnhur9tIHThu7cgbOG7hyA5LjExJSwgY8OybiBs4bqhaSAxNTEgc2luaCB2acOqbiBuYW0gY8OzIG3hu6ljIHRpw6p1IHRo4bulIHLGsOG7o3UgdGjhuqVwIGNoaeG6v20gMzguMjMlLiBUw7NtIGzhuqFpIHRow6wgY+G6oyBoYWkgZ2nhu5tpIHTDrW5oIMSR4buBdSBjw7Mgc+G7kSBsxrDhu6NuZyDDrXQgc2luaCB2acOqbiBjw7MgbeG7qWMgdGnDqnUgdGjhu6Ugcsaw4bujdSBjYW8sIMSRYSBz4buRIMSR4buBdSB0acOqdSB0aOG7pSByxrDhu6N1IHRo4bqlcCBob+G6t2MgaOG6p3UgbmjGsCBsw6Aga2jDtG5nIHRpw6p1IHRo4bulLg0KDQojIyAzLjQgUmVsYXRpdmUgUmlzayB2w6Aga2hv4bqjbmcgxrDhu5tjIGzGsOG7o25nIGNobyBSZWxhdGl2ZSBSaXNrDQoNCiMjIyAzLjQuMSBCaeG6v24gRGFsY18yIHbDoCBzY2hvb2wNCmBgYHtyfQ0KbGlicmFyeShEZXNjVG9vbHMpDQp0YWJsZShzdGQkRGFsY18yLHN0ZCRzY2hvb2wpDQptdGh0IDwtIHRhYmxlKHN0ZCREYWxjXzIsc3RkJHNjaG9vbCkNClJlbFJpc2sobXRodCwgY29uZi5sZXZlbCA9IC45NSkNCmBgYA0KDQpfTmjhuq1uIHjDqXQ6XyANCg0KLSBSZWxhdGl2ZSBSaXNrIGPhu6dhIERhbGNfMiB2w6Agc2Nob29sIGzDoCAwLjgzMzA2OTYgbmjhu48gaMahbiAxLiDEkGnhu4F1IG7DoHkgY8OzIG5naMSpYSBsw6Agc2luaCB2acOqbiB0csaw4budbmcgR1AgY8OzIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIG5oaeG7gXUgaMahbiBzbyB24bubaSB0csaw4budbmcgTVMgKG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIOG7nyB0csaw4budbmcgR1AgYuG6sW5nIDAuODMzIGzhuqduIHNvIHbhu5tpIHRyxrDhu51uZyBNUykuDQoNCi0gS2hv4bqjbmcgdGluIGPhuq15IG7DoHkgY2hvIHRo4bqleSBy4bqxbmcsIHbhu5tpIG3hu6ljIMSR4buZIHRpbiBj4bqteSAodGjGsOG7nW5nIGzDoCA5NSUpLCBt4bupYyDEkeG7mSB0acOqdSB0aOG7pSByxrDhu6N1IHRo4buxYyBz4buxIG7hurFtIHRyb25nIGtob+G6o25nIHThu6sgMC42NzExNTE5IMSR4bq/biAwLjk1NDk3MjEuIFbDrCBraG/huqNuZyB0aW4gY+G6rXkgbsOgeSBraMO0bmcgYmFvIGfhu5NtIGdpw6EgdHLhu4sgMSwga+G6v3QgcXXhuqMgbsOgeSBjaG8gdGjhuqV5IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ogZ2nhu69hIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIOG7nyBoYWkgdHLGsOG7nW5nLiANCg0KIyMjIDMuNC4yIEJp4bq/biBEYWxjXzIgdsOgIHNleA0KYGBge3J9DQpsaWJyYXJ5KERlc2NUb29scykNCnRhYmxlKHN0ZCREYWxjXzIsc3RkJHNleCkNCm10aHQgPC0gdGFibGUoc3RkJERhbGNfMixzdGQkc2V4KQ0KUmVsUmlzayhtdGh0LCBjb25mLmxldmVsID0gLjk1KQ0KYGBgDQoNCl9OaOG6rW4geMOpdDpfDQoNCi0gUmVsYXRpdmUgUmlzayBj4bunYSBEYWxjXzIgdsOgIHNleCBsw6AgMC4zMTkwOTA5IG5o4buPIGjGoW4gMS4gxJBp4buBdSBuw6B5IGPDsyBuZ2jEqWEgbMOgIHNpbmggdmnDqm4gbuG7ryBjw7MgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBjYW8gw610IGjGoW4gc28gduG7m2kgc2luaCB2acOqbiBuYW0gKG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIOG7nyBzaW5oIHZpw6puIG7hu68gYuG6sW5nIDAuMzE5IGzhuqduIHNvIHbhu5tpIHNpbmggdmnDqm4gbmFtKS4NCg0KLSBLaG/huqNuZyB0aW4gY+G6rXkgbsOgeSBjaG8gdGjhuqV5IHLhurFuZywgduG7m2kgbeG7qWMgxJHhu5kgdGluIGPhuq15ICh0aMaw4budbmcgbMOgIDk1JSksIG3hu6ljIMSR4buZIHRpw6p1IHRo4bulIHLGsOG7o3UgdGjhu7FjIHPhu7EgbuG6sW0gdHJvbmcga2hv4bqjbmcgdOG7qyAwLjE2NTk5MjYgxJHhur9uIDAuNTY1NjYzMy4gVsOsIGtob+G6o25nIHRpbiBj4bqteSBuw6B5IGtow7RuZyBiYW8gZ+G7k20gZ2nDoSB0cuG7iyAxLCBr4bq/dCBxdeG6oyBuw6B5IGNobyB0aOG6pXkgY8OzIMO9IG5naMSpYSB0aOG7kW5nIGvDqiBnaeG7r2EgbeG7qWMgxJHhu5kgdGnDqnUgdGjhu6Ugcsaw4bujdSBjYW8g4bufIGdp4bubaSB0w61uaCBu4buvIHbDoCBuYW0uDQoNCiMjIDMuNSBPZGRzIHJhdGlvIHbDoCBraG/huqNuZyDGsOG7m2MgbMaw4bujbmcgY2hvIE9kZHMgcmF0aW8NCg0KIyMjIDMuNS4xIEJp4bq/biBEYWxjXzIgdsOgIHNjaG9vbA0KYGBge3J9DQp0YWJsZShzdGQkRGFsY18yLHN0ZCRzY2hvb2wpDQpudG5oIDwtIHRhYmxlKHN0ZCREYWxjXzIsc3RkJHNjaG9vbCkNCk9kZHNSYXRpbyhudG5oLGNvbmYubGV2ZWwgPSAuOTUpDQpgYGANCg0KX05o4bqtbiB4w6l0Ol8NCg0KLSBPZGRzIHJhdGlvIGPhu6dhIERhbGNfMiB2w6Agc2Nob29sIGzDoCAwLjMzMjI3ODUgbmjhu48gaMahbiAxLiDEkGnhu4F1IG7DoHkgY8OzIG5naMSpYSBsw6AgdOG7tyBs4buHIGNow6puaCBj4bunYSBzaW5oIHZpw6puIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIHThuqFpIHRyxrDhu51uZyBHUCBzbyB24bubaSB0csaw4budbmcgTVMgbMOgIDAuMzMyLCB04bupYyBsw6Agc2luaCB2acOqbiB04bqhaSB0csaw4budbmcgR1AgY8OzIHThu7cgbOG7hyB0acOqdSB0aOG7pSByxrDhu6N1IGNhbyB0aOG6pXAgaMahbiBzbyB24bubaSBzaW5oIHZpw6puIHThuqFpIHRyxrDhu51uZyBNUy4NCg0KLSBLaG/huqNuZyB0aW4gY+G6rXkgbsOgeSBjaG8gdGjhuqV5IHLhurFuZywgduG7m2kgbeG7qWMgxJHhu5kgdGluIGPhuq15ICh0aMaw4budbmcgbMOgIDk1JSksIG9kZHMgcmF0aW8gdGjhu7FjIHPhu7EgbuG6sW0gdHJvbmcga2hv4bqjbmcgdOG7qyAwLjE1NDM4NjUgxJHhur9uIDAuNzE1MTQ2Ni4gVsOsIGtob+G6o25nIHRpbiBj4bqteSBuw6B5IGtow7RuZyBiYW8gZ+G7k20gZ2nDoSB0cuG7iyAxLCBr4bq/dCBxdeG6oyBuw6B5IGPDsyDDvSBuZ2jEqWEgdGjhu5FuZyBrw6ouDQoNCiMjIyAzLjUuMiBCaeG6v24gRGFsY18yIHbDoCBzZXgNCmBgYHtyfQ0KdGFibGUoc3RkJERhbGNfMixzdGQkc2V4KQ0KbnRuaCA8LSB0YWJsZShzdGQkRGFsY18yLHN0ZCRzZXgpDQpPZGRzUmF0aW8obnRuaCxjb25mLmxldmVsID0gLjk1KQ0KYGBgDQoNCl9OaOG6rW4geMOpdDpfDQoNCi0gT2RkcyByYXRpbyBj4bunYSBEYWxjXzIgdsOgIHNleCBsw6AgMC4xNjc3Nzc3OCBuaOG7jyBoxqFuIDEuIMSQaeG7gXUgbsOgeSBjw7MgbmdoxKlhIGzDoCB04bu3IGzhu4cgY2jDqm5oIGPhu6dhIHNpbmggdmnDqm4gdGnDqnUgdGjhu6Ugcsaw4bujdSBjYW8gbMOgIG7hu68gc28gduG7m2kgbmFtIGzDoCAwLjE2NywgdOG7qWMgbMOgIHNpbmggdmnDqm4gbuG7ryBjw7MgdOG7tyBs4buHIHRpw6p1IHRo4bulIHLGsOG7o3UgY2FvIHRo4bqlcCBoxqFuIHNvIHbhu5tpIHNpbmggdmnDqm4gbmFtLg0KDQotIEtob+G6o25nIHRpbiBj4bqteSBuw6B5IGNobyB0aOG6pXkgcuG6sW5nLCB24bubaSBt4bupYyDEkeG7mSB0aW4gY+G6rXkgKHRoxrDhu51uZyBsw6AgOTUlKSwgb2RkcyByYXRpbyB0aOG7sWMgc+G7sSBu4bqxbSB0cm9uZyBraG/huqNuZyB04burIDAuMDc1Nzg4NTYgxJHhur9uIDAuMzcxNDE5OTUuIFbDrCBraG/huqNuZyB0aW4gY+G6rXkgbsOgeSBraMO0bmcgYmFvIGfhu5NtIGdpw6EgdHLhu4sgMSwga+G6v3QgcXXhuqMgbsOgeSBjw7Mgw70gbmdoxKlhIHRo4buRbmcga8OqLg0KDQojIyAzLjYgTcO0IGjDrG5oIHjDoWMgc3XhuqV0IHR1eeG6v24gdMOtbmgNCmBgYHtyfQ0KbGlicmFyeSh0aWR5dmVyc2UpDQp0biA8LSBzdGQgJT4lIHNlbGVjdChEYWxjXzIsIHNjaG9vbCwgc2V4KQ0KdG4gPC0gdG4gJT4lIG11dGF0ZShEYWxjXzIgPSBpZl9lbHNlKERhbGNfMiA9PSAidGjhuqVwIiwgMCwgMSkpDQp0biA8LSB0biAlPiUgbXV0YXRlKHNjaG9vbCA9IGNhc2Vfd2hlbigNCiAgc2Nob29sID09ICJHUCIgfiAwLA0KICBzY2hvb2wgPT0gIk1TIiB+IDENCikpDQp0biA8LSB0biAlPiUgbXV0YXRlKHNleCA9IGNhc2Vfd2hlbigNCiAgc2V4ID09ICJGIiB+IDAsDQogIHNleCA9PSAiTSIgfiAxDQopKQ0KaGVhZCh0bikNCmBgYA0KDQojIyMgMy42LjEgSOG7k2kgcXV5IGdp4buvYSBEYWxjXzIgdsOgIHNjaG9vbA0KYGBge3J9DQpscG1fc2Nob29sIDwtIGxtKERhbGNfMiB+IHNjaG9vbCwgZGF0YSA9IHRuKQ0Kc3VtbWFyeShscG1fc2Nob29sKQ0KYGBgDQoNCiMjIyAzLjYuMiBI4buTaSBxdXkgZ2nhu69hIERhbGNfMiB2w6Agc2V4DQpgYGB7cn0NCmxwbV9zZXggPC0gbG0oRGFsY18yIH4gc2V4LCBkYXRhID0gdG4pDQpzdW1tYXJ5KGxwbV9zZXgpDQpgYGANCg0KIyMgMy43IE3DtCBow6xuaCBo4buTaSBxdXkgTG9naXN0aWMNCg0KIyMjIDMuNy4xIE3DtCBow6xuaCBMb2dpc3RpYyBjaG8gYmnhur9uIHNjaG9vbA0KYGBge3J9DQpscG1fc2Nob29sIDwtIGdsbShEYWxjXzIgfiBzY2hvb2wsIGRhdGEgPSB0biwgZmFtaWx5ID0gYmlub21pYWwpDQpzdW1tYXJ5KGxwbV9zY2hvb2wpDQpgYGANCg0KIyMjIDMuNy4yIE3DtCBow6xuaCBMb2dpc3RpYyBjaG8gYmnhur9uIHNleA0KYGBge3J9DQpscG1fc2V4IDwtIGdsbShEYWxjXzIgfiBzZXgsIGRhdGEgPSB0biwgZmFtaWx5ID0gYmlub21pYWwpDQpzdW1tYXJ5KGxwbV9zZXgpDQpgYGANCg0KIyMgMy44IE3DtCBow6xuaCBMb2dpdA0KDQojIyMgMy44LjEgTcO0IGjDrG5oIExvZ2l0IGNobyBiaeG6v24gc2Nob29sDQoNCmBgYHtyfQ0KcmVnbG9naXQgPC0gZ2xtIChEYWxjXzIgfiBzY2hvb2wsIGRhdGEgPSB0biwgZmFtaWx5ID0gYmlub21pYWwobGluayA9ICdsb2dpdCcpKQ0KcHJpbnQocmVnbG9naXQpDQpgYGANCg0KIyMjIDMuOC4yIE3DtCBow6xuaCBMb2dpdCBjaG8gYmnhur9uIHNleA0KDQoNCg0KIyMgMy45IE3DtCBow6xuaCBQcm9iaXQNCg0KIyMjIDMuOS4xIE3DtCBow6xuaCBQcm9iaXQgY2hvIGJp4bq/biBzY2hvb2wNCg0KDQojIyMgMy45LjIgTcO0IGjDrG5oIFByb2JpdCBjaG8gYmnhur9uIHNleA0KYGBge3J9DQpyZWdwcm9iaXQgPC0gZ2xtIChEYWxjXzIgfiBzZXgsIGRhdGEgPSB0biwgZmFtaWx5ID0gYmlub21pYWwobGluayA9ICdwcm9iaXQnKSkNCnByaW50KHJlZ3Byb2JpdCkNCmBgYA0KDQoobeG7pWMgbOG7m24gbmhhKSAgxJDDoW5oIGdpw6EgbcO0IGjDrG5oDQpDaOG7iSBz4buRIEFJQw0KDQoNCg0KDQoNCg0K

TIỂU LUẬN MÔN PHÂN TÍCH DỮ LIỆU ĐỊNH TÍNH

PHÂN TÍCH NHỮNG YẾU TỐ ẢNH HƯỞNG ĐẾN MỨC ĐỘ TIÊU THỤ RƯỢU CỦA SINH VIÊN

LỜI CẢM ƠN

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

1.1 Lí do chọn đề tài

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

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

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

1.5 Kết cấu bài nghiên cứu

Chương 2: DỮ LIỆU NGHIÊN CỨU

2.1 Giới thiệu bộ dữ liệu

2.2 Mô tả các biến trong bộ dữ liệu

Chương 3: PHÂN TÍCH DỮ LIỆU

3.1 Thống kê mô tả biến phụ thuộc

3.1.1 Lập bảng tần số và bảng tần suất

3.1.2 Lập biểu đồ

3.2 Thống kê suy diễn

3.2.1 Kiểm định tính độc lập

3.2.1.1 Biến Dalc_2 và school

3.2.1.2 Biến Dalc_2 và sex

3.2.1.3 Biến Dalc_2 và address

3.2.1.4 Biến Dalc_2 và schoolsup

3.2.1.5 Biến Dalc_2 và famsup

3.2.1.6 Biến Dalc_2 và activities

3.2.1.7 Biến Dalc_2 và internet

3.2.1.6 Biến Dalc_2 và romantics

3.3 Thống kê mô tả 2 biến

3.3.1 Biến Dalc_2 và sex

3.3.1.1 Lập bảng tần số và bảng tần suất

3.3.1.2 Lập biểu đồ

3.3.2 Biến Dalc_2 và sex

3.3.2.1 Lập bảng tần số và bảng tần suất

3.3.2.2 Lập biểu đồ

3.4 Relative Risk và khoảng ước lượng cho Relative Risk

3.4.1 Biến Dalc_2 và school

3.4.2 Biến Dalc_2 và sex

3.5 Odds ratio và khoảng ước lượng cho Odds ratio

3.5.1 Biến Dalc_2 và school

3.5.2 Biến Dalc_2 và sex

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

3.6.1 Hồi quy giữa Dalc_2 và school

3.6.2 Hồi quy giữa Dalc_2 và sex

3.7 Mô hình hồi quy Logistic

3.7.1 Mô hình Logistic cho biến school

3.7.2 Mô hình Logistic cho biến sex

3.8 Mô hình Logit

3.8.1 Mô hình Logit cho biến school

3.8.2 Mô hình Logit cho biến sex

3.9 Mô hình Probit

3.9.1 Mô hình Probit cho biến school

3.9.2 Mô hình Probit cho biến sex