XÁC SUẤT THỐNG kê (MT2013) báo cáo bài tập lớnđề tài 6

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA  XÁC SUẤT THỐNG KÊ (MT2013) BÁO CÁO BÀI TẬP LỚN ĐỀ TÀI GVHD: Nguyễn Kiều Dung Nhóm: MT24 Sinh viên thực hiện: Họ tên MSSV Lớp Nguyễn Thị Vân Anh 2012606 L13 Võ Nguyễn Phương Long Vĩnh 2015081 L13 Nguyễn Xuân Thắng 2014545 L13 Nguyễn Thanh Tùng 2014999 L02 Nguyễn Lâm Tùng 2012388 L15 Thành phố Hồ Chí Minh, tháng năm 2022 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat MỤC LỤC HOẠT ĐỘNG 1.1 Cơ sở lí thuyết 1.1.1 Phân tích phương sai – ANOVA nhân tố 1.1.2 Mơ hình hồi quy tuyến tính bội 1.2 Thực 1.2.1 Đọc làm liệu, thực thống kê mô tả .6 1.2.2 Chia liệu .11 1.2.3 Chọn mơ hình 12 1.2.4 Kiểm tra giả thiết mơ hình 16 1.2.5 Dự báo (Prediction) .18 HOẠT ĐỘNG 21 2.1 Tổng quan đề hoạt động 21 2.2 Cơ sở lí thuyết 22 2.2.1 Mơ hình hồi quy tuyến tính bội 22 2.2.2 Đánh giá phù hợp mơ hình 22 2.2.3 Phương pháp bình phương nhỏ 25 2.2.4 Các giả định mơ hình hồi quy 26 2.2 Thực 27 2.3.1 Đọc liệu (Import data): .27 2.3.2 Làm liệu (Data cleaning): 28 2.3.3 Làm rõ liệu (Data visualization): .30 2.3.4 Mơ hình hồi quy tuyến tính: 36 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat HOẠT ĐỘNG 1.1 Cơ sở lí thuyết 1.1.1 Phân tích phương sai – ANOVA nhân tố - Phân tích phương sai mơ hình dùng để xem xét biến động biến ngẫu nhiên định lượng X chịu tác động trực tiếp hay nhiều yếu tố ngun nhân - Trong mơ hình phân tích phương sai yếu tố, ta kiểm định so sánh trung bình biến ngẫu nhiên X tổng thể khác dựa vào mẫu quan sát lấy từ tổng thể Các tổng thể phân biệt mức độ khác yếu tố xem xét  Giả thiết toán ANOVA o Các tổng thể có phân phối chuẩn �(�� ; �2� ) với i = 1,2,…,k với k số tổng thể (� ≥ 3) o Phương sai tổng thể �21 = �22 = … = �2� o Các mẫu quan sát (từ tổng thể) lấy độc lập  Đặt giả thuyết kiểm định o Giả thuyết không �0 (null hypothesis): �1 = �2 = … = ��  o Giả thuyết đối �1 (alternative hypothesis): ∃�� ≠ �� với � ≠ � Tính giá trị kiểm định thống kê Source of Tổng bình phương chênh variation lệch Between Groups Within Groups ��� = ��� = � �=1 � Phương sai Bậc tự BPCL) �� ∗ (�� − �)2 k-1 (��� − �� )2 N-k �� �=1 �=1 (Trung bình ��� = ��� = ��� �−1 ��� �−� Tiêu chuẩn kiểm định F �= ��� ��� TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Total ��� = � �� �=1 �=1 (��� − �)2 N-1 SST = SSB + SSW 1.1.2 Mơ hình hồi quy tuyến tính bội - Hồi quy tuyến tính bội phần mở rộng hồi quy tuyến tính đơn Nó sử dụng ta muốn dự đoán giá trị biến phản hồi dựa giá trị hai nhiều biến giải thích khác Biến mà muốn dự đốn gọi biến phản hồi (biến phụ thuộc) Các biến mà ta sử dụng để dự đoán giá trị biến phản hồi gọi biến giải thích (biến dự báo, biến phụ thuộc) Ví dụ: sử dụng hồi quy bội số để dự đoán kết kỳ thi XSTK dựa thời gian ôn tập, niên khóa, giới tính sinh viên - Hồi quy bội cho phép xác định phù hợp tổng thể mơ hình đóng góp tương đối yếu tố dự báo tổng phương sai giải thích - Vì biến dự báo đơn lẻ khơng cung cấp mơ tả đầy đủ số lượng biến dự báo chìa khóa tác động đến biến đáp ứng theo cách đặc biệt quan trọng Các dự báo biến đáp ứng dựa vào mơ hình có biến dự báo riêng lẻ khơng xác Vì mơ hình hồi quy tuyến tính bội đưa ra: oXét trường hợp có � − biến dự báo �1 , …, ��−1 Mơ hình hồi quy:  �� = �0 + �1 ��1 + �2 ��2 + … + ��−1 ���−1 + �� gọi mơ hình bậc với � − biến dự báo  Hay: �� = �0 +  Trong đó: �−1 � � �=1 � �� + �� �0 , �1 , …, ��−1 : �à �á� �ℎ�� �ố ��1 , ��2 , …, ���−1 : �à �á� ℎằ�� �ố ��ế� �� ~� 0; �2 � = 1,2, …, � oKhi � − = mơ hình hồi quy là: �� = �0 + �1 ��1 + �� mơ hình hồi quy tuyến tính đơn oGiả sử: � �� = 0, hàm đáp ứng với mơ hình: XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6  � � = �0 + �1 ��1 + �2 ��2 + … + ��−1 ���−1 oTham số �� thay đổi trung bình đáp ứng E(Y) với đơn vị tăng biến dự báo �� tất biến dự báo lại coi số oẢnh hưởng biến dự báo trung bình đáp ứng biến dự báo khác cố định Do đó, mơ hình hồi quy bậc thiết kế cho biến dự báo mà ảnh hưởng trung bình đáp ứng cộng tính hay khơng có tương tác oVậy mơ hình hồi quy tuyến tính tổng qt với điều kiện sai số chuẩn có quan sát �� biến chuẩn độc lập, với trung bình � � phương sai không đổi �2  Các biến dự báo định tính Mơ hình hồi quy tuyến tính tổng quát: �� = �0 + �1 ��1 + �2 ��2 + … + ��−1 ���−1 + �� bao gồm biến dự báo định lượng biến dự báo định tính Nên ta sử dụng biến số nhận giá trị để định nghĩa lớp giá trị biến định tính  Ước lượng hệ số hồi quy Tiêu chuẩn bình phương cực tiểu tổng qt hóa cho mơ hình hồi quy tuyến tính tổng quát sau: �= � �=1 �2� = � �=1 �� − �0 − �1 ��1 − … − ��−1 ���−1 Các ước lượng bình phương cực tiểu giá trị �0 , �1 , …, ��−1 làm cực tiểu hóa Q Ta biểu diễn vector ước lượng hệ số hồi quy �0 , �1 , …, ��−1 b: �0 �1 �= ��−1 Các phương trình chuẩn bình phương cực tiểu cho mơ hình hồi quy tuyến tính tổng quát là: �' � � = �' � ước lượng bình phương cực tiểu là: � = (�' �)−1 �'� XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6  Bảng ANOVA cho mơ hình hồi quy tuyến tính tổng quát Nguồn SS biến đổi Hồi quy Sai số Tổng số df ��� = �' �' � − �'�� � ���� = �' � − �'�� � �−1 ��� = �' � − �' �' � �−� �−1  Kiểm định F cho quan hệ hồi quy MS ��� = ��� = ��� �−1 ��� �−� Để kiểm định liệu có hay khơng quan hệ hồi quy biến đáp ứng biến �: �1 , …, ��−1 tức lựa chọn giả thuyết: �0 : �1 = �2 = … = ��−1 = �� : �ℎô�� �ℎả� �ấ� �ả �� , � = 1, …, � − đề� �ằ�� �ℎô�� Ta dùng thống kê kiểm định: � ∗= ��� ��� Nếu � ∗≤ � −∝ ; � − 1; � − � chấp nhận �0 Nếu � ∗> � −∝ ; � − 1; � − � chấp nhận ��  Hệ số xác định bội �2 ��� ��� =1− ���� ���� �2 = Hệ số xác định bội cho biết biến dự báo mơ hình giải thích bao nhiều phần tram thay đổi biến đáp ứng Vì ta có: ≤ �2 ≤ �2 = tất giá trị �� = (� = 1, …� − 1) �2 = tất quan sát nằm mặt đáp ứng, �ứ� �� = �� , �ớ� �ọ� � Thêm nhiều biến dự báo X vào mơ hình làm tăng thêm �2 Hệ số xác định bội hiệu chỉnh, ký hiệu �2 � , điều chỉnh �2 cách chia tổng bình phương cho bậc tự nó: XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 ��� � − ��� �−� �2� = − =1− ���� � − � ���� �−1 Hệ số xác định bội hiệu chỉnh thực nhỏ biến X khác đưa vào mơ hình  Hệ số tương quan bội Hệ số tương quan bội � = �2 Khi có biến X mơ hình hồi quy (khi � − = 1), hệ số tương quan bội R trị tuyệt đối hệ số tương quan r tương quan đơn 1.2 Thực 1.2.1 Đọc làm liệu, thực thống kê mô tả Đọc dư liệu (Import data) Đọc tệp tin "auto_mpg.csv" gán với tên grade Làm sach dư liệu (Data cleaning) Từ liệu mydata , trích liệu bao gồm biến đề đặt tên "new_DF" Kiêm tra dư liệu khuyết tệp tin "new_DF" XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Nhận xét: Ta nhận thấy có liệu khuyết biến Horsepower(cụ thể dòng 33 127 331 337 355 375), số lượng liệu khuyết chiếm 0.01507538% so với tổng số quan sát liệu Do liệu bị khuyết hoàn toàn ngẫu nhiên (Missing at Random) tỉ lệ liệu khuyết thấp Nên ta chọn cách xóa quan sát có chứa liệu khuyết Xóa liệu khuyết tệp tin "new_DF" XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Nhận xét: Sau xóa, ta thấy dịng tệp tin "new_DF" chứa liệu khuyết đa biến khỏi liệu tệp tin Các thống kê mô tả: * Thống kê số lượng cho biến horsepower Nhận xét: Dựa vào bảng thống kê ta biết được: Cơng suất động số lượng loại xe có chung cơng suất * Thống kê số lượng cho biến cylinders Nhận xét: Dựa vào bảng thống kê ta biết được: Phần lớn xe sản xuất với lượng xi lanh 4, (lần lượt 199, 83, 103 loại) với loại xe có xi lanh loại có xi lanh XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 * Thống kê số lượng cho biến displacement Nhận xét: Dựa vào bảng thống kê ta biết được: Các loại dung tích xy lanh số lượng loại xe tương ứng * Thống kê số lượng cho biến weight XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Mục tiêu việc xây dựng mơ hình hồi quy để xác định giá trị biến phụ thuộc biết trước giá trị cụ thể biến độc lập Một số thống kê cho thấy mơ hình hồi quy thực mục tiêu tốt đến đâu lệch chuẩn mơ hình hồi quy (còn gọi tên Sai số chuẩn ước lượng) Giá trị ước lượng từ thông tin mẫu độ lệch chuẩn mơ hình hồi quy (sai số chuẩn ước lượng) tính tốn sau đây: ��/� = ��� �−�−1 Trong n cỡ mẫu, k biến độc lập mơ hình Sai số chuẩn ước lượng đo lường phân tán giá trị thực tế đo lường biến phụ thuộc quanh giá trị biến phụ thuộc dự đoán đường hồi quy Đánh giá ý nghĩa của biến độc lập riêng biệt : Ở kiểm định F, giả sử H1 chấp nhận ta kết luận mơ hình tồn diện có ý nghĩa Điều có ý nghĩa có biến độc lập mơ hình giải thích cách có ý nghĩa cho biến thiên phụ thuộc Tuy nhiên điều khơng có ý nghĩa tất biến độc lập đưa vào mơ hình có ý nghĩa, để xác định biến độc lập có ý nghĩa kiểm định giả thuyết sau: − H 0: β = − H 1: β j = Chúng ta dùng kiểm định t để kiểm định nghĩa hệ số hồi quy với độ tin cậy chọn trước, t xác định công thức: �= �� − ��� Trong bj hệ số dốc mơ hình hồi quy mẫu cho biến độc lập thứ j, ��� sai số chuẩn ước lượng lượng hệ số độ dốc biến độc lập thứ j Giá trị t tính tốn so sánh với giá trị t tới hạn tra từ bảng phân phối student với � (n – k - 1) bậc tự mức ý nghĩa 2.2.3 Phương pháp bình phương nhỏ 25 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Phương pháp bình phương nhỏ đưa nhà Toán học Đức Carl Friedrich Gauss - phương pháp ước lượng hồi quy tuyến tính phổ biến Với tổng thể, sai số (error) kí hiệu e, cịn mẫu nghiên cứu lúc gọi phần dư kí hiệu � Biến thiên phần dư tính tổng bình phương tất phần dư cộng lại Nguyên tắc phương pháp hồi quy OLS làm cho biến thiên phần dư phép hồi quy nhỏ Khi biểu diễn mặt phăng Oxy, đường hồi quy đường thăng qua đám đơng điểm liệu mà đó, khoảng cách từ điểm liệu (tuyệt đối � đến đường hồi quy ngắn nhất) Từ đồ thị scatter biểu diễn mối quan hệ biến độc lập biến phụ thuộc, điểm liệu nằm phân tán có xu hướng chung tạo thành đường thăng Chúng ta có nhiều đường thăng hồi quy qua đám đông điểm liệu đường nhất, vấn đề ta phải chọn đường thăng mô tả sát xu hướng liệu Bình phương nhỏ OLS tìm đường thăng dựa nguyên tác cực tiểu hóa khoảng cách từ điểm liệu đến đường thăng Trong hình đường màu đỏ đường hồi quy OLS 2.2.4 Các giả định mơ hình hồi quy a Hàm hồi quy tuyến tính theo tham số Điều có nghĩa q trình thực hành hồi quy thực tế miêu tả mối quan hệ dạng: y = β1 + β2x2 + β3x3 + β4x4 + + βkxk + � mối quan hệ thực viết lại ví dụ dạng lấy loga hai 26 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 vế b E(�� ) = 0: Kỳ vọng yếu tố ngẫu nhiên Trung bình tổng thể sai số Điều có nghĩa có số giá trị sai số mang dấu dương số mang dấu âm Do hàm xem đường trung bình nên giả định sai số ngẫu nhiên bị loại trừ nhau, mức trung bình tổng thể c Cov (�� , �� ) = 0: Khơng có tương quan �� Khơng có tương quan quan sát yếu tố sai số Nếu ta xem xét chuỗi số liệu thời gian (dữ liệu thu nhập từ nguồn nhiều khoảng thời gian khác nhau), yếu tố sai số �� khoảng thời gian khơng có tương quan với yếu tố sai số khoảng thời gian trước d.Cov (�� , xi) = : � X khơng có tương quan với Khi biến giải thích lớn hay nhỏ yếu tố sai số khơng thay đổi theo e Var (�� = σ2): Phương sai với � � σ2 Tất giá trị �� phân phối giống với σ2, cho: Var(�� ) = E(�2 ) = f �� phân phối chuẩn Điều quan trọng phát sinh khoảng tin cậy thực kiểm định giả thuyết phạm vi mẫu nhỏ Nhưng phạm vi mẫu lớn hơn, điều trở nên không quan trọng g Giữa x2, x3, , xk khơng có quan hệ tuyến tính Nếu x2, x3, , xk có quan hệ tuyến tính người ta nói có tượng đa cộng tuyến Hay không tồn λi ≡ : λ1x1i + λ2x2i + λ3x3i + + λkx ki + vi = 2.2 Thực 2.3.1 Đọc liệu (Import data): Đọc tệp Range-Queries-Aggregates.csv lưu liệu vào biến query Code: > query head(query,3) # Xem dong dau tien cua du lieu Kết thực thi: Hình 1: Kết xem dịng tệp tin "query" 2.3.2 Làm liệu (Data cleaning): a) Từ liệu query, trích liệu bao gồm biến đề đặt tên "new_DF": Code: > new_DF head(new_DF,3) Kết thực thi: Hình 2: Kết xem dòng tệp tin "new_DF" b) Kiểm tra liệu khuyết tệp tin "new_DF": Code: > apply(is.na(new_DF),2,which) #kiem tra NA va xuat vi tri dong chua NA > apply(is.na(new_DF),2,sum) #kiem tra va dem gia tri NA > apply(is.na(new_DF),2,mean) #tinh ti le NA Kết thực thi: 28 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Hinh 3: Kết kiểm tra liệu khuyết tệp tin "new_DF" Nhận xét: Ta thấy có 157 liệu khuyết biến avg, số lượng liệu khuyết chiếm 0.0785% so với tổng số quan sát liệu Do liệu bị khuyết có tỉ lệ thấp nên ta chọn cách xóa quan sát có liệu khuyết c) Xóa liệu khuyết tệp tin “new_DF”: Code: > new_DF apply(is.na(new_DF),2,sum) #kiem tra lai NA Kết thực thi: Hinh 4: Kết kiểm tra lại liệu khuyết tệp tin "new_DF" d) Kiểm tra tính xác liệu (Dữ liệu có dạng số khơng?) Code: 29 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 > is.numeric(new_DF$x) > is.numeric(new_DF$y) > is.numeric(new_DF$x_range) > is.numeric(new_DF$y_range) > is.numeric(new_DF$count) > is.numeric(new_DF$sum_) > is.numeric(new_DF$avg) Kết thực thi: 2.3.3 Làm rõ liệu (Data visualization): Thống kê mô tả: Dùng thống kê mẫu: Các biến liệu biến liên tục a) Tính giá trị thống kê mơ tả bao gồm: trung bình (mean), trung vị (median), độ lệch chuẩn (sd), giá trị lớn (max) giá trị nhỏ (min) biến liên tục Code: > mean sd Q1 median Q3 max data.frame(mean, sd, Q1, median, Q3, min, max) Kết thực thi: b) Vẽ biểu đồ Histogram thể phân phối biến AVG: Code: > hist(new_DF[,"avg"], xlab = "avg", main = "Histogram of avg", label = T, col = "lightgreen", border = "seagreen") Kết thực thi: 31 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 c) Dùng lệnh pairs() vẽ phân phối biến AVG theo biến x, y, x_range, y_range, count sum_: Code: > pairs(new_DF$avg ~ new_DF$x, main = "Pairs of AVG for each category of x", col = "blue") > pairs(new_DF$avg ~ new_DF$y, main = "Pairs of AVG for each category of y", col = "blue") > pairs(new_DF$avg ~ new_DF$x_range, main = "Pairs of AVG for each category of x_range", col = "blue") > pairs(new_DF$avg ~ new_DF$y_range, main = "Pairs of AVG for each category of y_range", col = "blue") > pairs(new_DF$avg ~ new_DF$count, main = "Pairs of AVG for each category of count", col = "blue") > pairs(new_DF$avg ~ new_DF$sum_, main = "Pairs of AVG for each category of sum_", col = "blue") Kết thực thi: 32 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 33 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 34 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Nhận xét: Nhận thấy rằng, ta sử dụng lệnh pair() rút đồ thị phân tán, từ đồ thị ta thấy biến x, y, x_range, y_range, count sum_ có quan hệ tuyến 35 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 tính với biến AVG Nghĩa biến x, y, x_range, y_range, count sum_ tăng lên biến AVG có xu hướng tăng theo 2.3.4 Mơ hình hồi quy tuyến tính: a) Ta xây dựng mơ hình hồi quy bội (gọi mơ hình 1) bao gồm: - Biến phụ thuộc: avg - Biến độc lập: x, y, x_range, y_range, count sum_ - Mơ hình biểu diễn sau: avg = �0 + �1 × x + �2 × y x_range Code: + �4 × y_range + �5 × count + �6 × sum_ + �3 × Ta thực ước lượng hệ số �� , i = 0, …, > model_1 = lm(avg ~ x + y + x_range + y_range + count + sum_, new_DF) #Xay dung mo hinh va luu voi ten model_1 > summary(model_1) #Ket qua mo hinh Kết thực thi: Nhận xét: Từ kết phân tích, ta thu được: 36 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Hệ số R2 hiệu chỉnh 0.7528 có nghĩa 75.28% biến thiên AVG giải thích mơ hình với biến độc lập x, y, x_range, y_range, count, sum_ Dựa vào kết mơ hình hồi quy, ta thấy giá trị Pr (> t ) hệ số biến độc lập nhỏ 2.2×10-16 nhỏ so với mức ý nghĩa 5% Do ta bác bỏ giả thiết H0 chấp nhận giả thiết H1 tất hệ số ứng với biến độc lập có ý nghĩa thống kê Vậy ta khơng loại bỏ biến khỏi mơ hình hồi quy b) Kiểm tra giả định mơ hình: Nhắc lại giả định mơ hình hồi quy: �� = �0 + �1 × �1 + + �� × �� + ϵ� , i = n − Tính tuyến tính liệu: mối quan hệ biến dự báo X biến phụ thuộc Y giả sử tuyến tính − Sai số có phân phối chuẩn − Phương sai sai số số − Các sai số ϵ có kỳ vọng = − Các sai số ϵ1 , , ϵ� độc lập với Ta thực phân tích thặng dư để kiểm tra giả định mơ hình Code: > par(mfrow = c(2,2)) #xep bieu hang cot > plot(model_3) #ve thi phan tich thang du Kết thực thi: 37 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 Nhận xét: − Đồ thị (Residuals vs Fitted) vẽ giá trị sai số tương ứng với giá trị dự báo, kiểm tra giả định tuyến tính liệu giả định sai số có kỳ vọng * Ta nhận thấy đường màu đỏ gần đường thăng nên giả định tuyến tính liệu thỏa man * Các sai số tập trung xung quanh đường thăng y = (ngoại trừ số ngoại lai) nên giả định sai số có kỳ vọng thỏa man − Đồ thị (Normal Q - Q) vẽ sai số đa chuẩn hóa kiểm tra giả định phân phối chuẩn Ta nhận thấy sai số đa chuẩn hóa chưa nằm đường thăng nên giả định chưa thỏa man − Đồ thị (Scale - Location) vẽ bậc hai sai số, kiểm tra giả định phương sai sai số số Ta nhận thấy bậc hai sai số tập trung quanh đường thăng màu đỏ, có ngoại lai tạm chấp nhận − Đồ thị (Residuals vs Leverage) xác định điểm có ảnh hưởng cao diện liệu Các điểm 123740, 153412, 154245 điểm có 38 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 TIEU LUAN MOI download : skknchat123@gmail.com moi nhat XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6 XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6XAC.SUAT.THONG.ke.(MT2013).bao.cao.bai.tap.londe.tai.6