CHƯƠNG VI PHÂN TÍCH HỒI QUI LOGISTIC (LOGISTIC REGRESSION)

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC Thí dụ 6.1 : Một điều tra về hội chứng hô hấp sinh sản của heo nái PRRS liên quan đến việc tiêm ngừa vaccin phòng bệnh này c

Title and Content Layout with List

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC

Thí dụ 6.1 : Một điều tra về hội chứng hô hấp sinh sản của heo

nái (PRRS) liên quan đến việc tiêm ngừa vaccin phòng bệnh này cho 2000 heo nái được chọn ngẫu nhiên tại một trang trại chăn nuôi của thành phố H với kết qủa qua bảng sau:

Làm sao chúng ta biết có mối liên quan giữa việc có chích ngừa và bệnh PRRS?

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC

(1) Chỉ số nguy cơ tương đối (Relative risk – RR)

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC

(1) Chỉ số nguy cơ tương đối (Relative risk – RR)

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC

(2) Chỉ số chênh (ODD)

- Chỉ số chênh là tỉ số của 2 xác suất

- Nếu P là xác suất bị bệnh thì 1- P là xác suất không

bị bệnh Chỉ số chênh được tính như sau:

ODD = P/(1-P)

- Nếu ODD > 1: khả năng bị bệnh cao hơn khả năng không bị

bệnh

- Nếu ODD < 1: khả năng bị bệnh thấp khả năng không bị bệnh

- Nếu ODD = 1: khả năng bị bệnh bằng khả năng không bị

bệnh

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC

(2) Chỉ số chênh (ODD)

Trở lại với thí dụ trên, ta có:

› Chỉ số chênh trong nhóm có chích ngừa là:

ODD1 = P1/(1-P1) = 0,15/(1-0,15)

= 0,15/0,85 =0,1767 = 17,64%

Nghĩa là heo nái có chích ngừa vẫn có nguy cơ bịPRRS với khoảng 0,1764 lần (17,64%) so với heo nái cóchích ngừa mà không bị PRRS.

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC

6.1-MỘT SỐ KHÁI NIỆM LIÊN QUAN ĐẾN PHÂN TÍCH HỒI QUI LOGISTIC

(3) Tỉ số chênh (Odds Ratio – OR )

- Tỉ số chênh (OR) là tỉ số của 2 ODD

4,41%

0,0441 (0,85)0,80

0,15(0,20) P1)P2

(1

P2)

P1(1 P2)

P2/(1

P1)

P1/(1 ODD

ODD OR

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂNXét 2 biến X và Y:

X: biến độc lập có dạng định tính nhị phân như sốnghoặc chết, bệnh hoặc không, sạch hoặc dơ, tốt hoặcxấu…hay dạng định lượng như các chỉ tiêu về sinh lý sinhhóa máu, sinh trưởng, sản xuất của thú…

Y: biến phụ thuộc chỉ có duy nhất dạng nhị phân

 phương trình hồi qui giữa Y và một hay nhiều biến Xđược gọi là phương trình hồi qui logistic nhị phân

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Nếu X là biến độc lập, Y là biến phụ thuộc Y và P là xácsuất xãy ra sự kiện của biến Y (ví dụ chết) và 1-P là xácxuất để sự kiện đó không xãy ra (ví dụ sống), thì phươngtrình hồi qui logistic nhị phân có dạng như sau:

(1)

k

=

X b

+

… +

X b

+ X

b + b

(

P 1

X b

+

… +

X b

+ X

b + b

) (

P

k

=

X b

+

… +

X b

+ X

b + b

) (

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

ODDP

AODD

1

ODDP









Muốn biết xác suất P của biến Y: cần có các kết quả của biến Y

và các X1, X2, Xn  xác định và thiết lập được phương trình Logit (P)  tìm được các hệ số b0, b1, b2…bn của phương trình,

từ đó tính ra được A và thế vào pt trên sẽ tính được xác suất P

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

› Để xác định và thiết lập được phương trình Logit (P),người ta phải dùng phương pháp tương đồng tối đa(The maximum likelihook) là một thuật toán rất phức tạpkhi tính toán thủ công nên cần sự hổ trợ giải quyết bằngcác phần mềm thống kê chuyên dụng như Minitab 16.2,PSSR, SAS 9.1, R

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

› * Thí dụ 6.2 :Một điều tra về bệnh Viêm phổi địa

phương (biến phụ thuộc Y) trên heo liên quan đến một

số yếu tố như chích ngừa, tự sản xuất giống và điềukiện vệ sinh chuồng trại (theo thứ tự biến độc lập X1,X2 và X3) Tìm phương trình hồi qui logistic giữa cácbiến này và xác suất dự đoán với kết qủa ghi nhận ởbảng sau?

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Sử dụng phần mềm Minitab 16.2 để tính toán cho ta kết

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Sử dụng phần mềm Minitab 16.2 để tính toán cho ta kết qủa các hệ số b0, b1, b2…

Phương trình hồi qui logistic giữa các biến còn lại :

Log e (ODD) = -0,8311 + 0,6954X 2 + 0,8654X 3.

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Phương trình hồi qui logistic giữa các biến còn lại :

Log e (ODD) = -0,8311 + 0,6954X 2 + 0,8654X 3.

phổi địa phương là bao nhiêu?

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Phương trình hồi qui logistic giữa các biến còn lại :

Log e (ODD) = -0,8311 + 0,6954X 2 + 0,8654X 3.

Viêm phổi địa phương là bao nhiêu?

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Phương trình hồi qui logistic giữa các biến còn lại :

Log e (ODD) = -0,8311 + 0,6954X 2 + 0,8654X 3.

bệnh Viêm phổi địa phương là bao nhiêu?

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Phương trình hồi qui logistic giữa các biến còn lại :

Log e (ODD) = -0,8311 + 0,6954X 2 + 0,8654X 3.

bệnh Viêm phổi địa phương là bao nhiêu?

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

6.2- ĐỊNH NGHĨA HỒI QUI LOGISTIC NHỊ PHÂN

Phương trình hồi qui logistic giữa các biến còn lại :

Log e (ODD) = -0,8311 + 0,6954X 2 + 0,8654X 3.

nhiêu?