Giao diện
Trong bài học trước, ta đã làm quen với mô hình tuyến tính đơn: Một phương trình đường thẳng mô tả mối liên hệ giữa một biến giải thích duy nhất và một biến phản hồi. Tuy nhiên, các hiện tượng kinh tế, xã hội và khoa học trong thực tế hiếm khi vận hành đơn giản như vậy. Một hiện tượng thường chịu sự tác động đồng thời của nhiều yếu tố đa dạng: Giá nhà không chỉ phụ thuộc vào diện tích mà còn vào vị trí, số phòng ngủ và năm xây dựng. Tương tự, mức lương của một kỹ sư không chỉ dựa vào số năm kinh nghiệm mà còn phụ thuộc vào trình độ học vấn, chứng chỉ nghề nghiệp và quy mô doanh nghiệp.
Khi dữ liệu mở rộng sang nhiều chiều, hệ thống thị giác của con người gặp phải rào cản sinh học: Ta không thể vẽ đồ thị trực quan vượt quá 3 chiều không gian. Để phân tích và định lượng các mối quan hệ đa chiều phức tạp đó, thống kê học cung cấp một công cụ mạnh mẽ: Mô hình Hồi quy tuyến tính bội (Multiple Linear Regression).
1. Khung toán học của Hồi quy tuyến tính bội
Khái niệm: Hồi quy tuyến tính bội (Multiple Linear Regression)
Hồi quy tuyến tính bội là phương pháp mô hình hóa một biến phản hồi định lượng liên tục
Khai triển tường minh bằng ký hiệu tổng dồn:
trong đó:
: Giá trị dự báo (Fitted value) của biến phản hồi. : Hệ số chặn (Intercept), là giá trị kỳ vọng của khi toàn bộ các biến giải thích đều đồng thời nhận giá trị bằng 0. ( ): Các hệ số hồi quy riêng phần (Partial regression coefficients). : Các biến giải thích, có thể là biến định lượng hoặc biến định tính đã được mã hóa.
Trong môi trường R, phương pháp bình phương tối thiểu cho hồi quy bội được thực hiện bằng hàm lm(), trong đó các biến giải thích được kết nối với nhau bằng dấu cộng +:
r
lm(y ~ x1 + x2 + ... + xp, data = my_data)Để thấu hiểu sâu sắc bản chất hình học và ý nghĩa thống kê của mô hình hồi quy bội, ta sẽ lần lượt phân tích hai kịch bản phân tích kinh điển:
- Kịch bản 1: Mô hình kết hợp giữa các biến định lượng và biến định tính.
- Kịch bản 2: Mô hình gồm nhiều biến định lượng liên tục (Mặt phẳng hồi quy).
2. Kịch bản 1: Kết hợp Biến định lượng & Biến định tính (Biến giả Dummy)
Để có một bối cảnh thực nghiệm sinh động, ta cùng nghiên cứu bảng dữ liệu đánh giá ẩm thực nổi tiếng từ cẩm nang Zagat (Zagat Survey). Zagat đánh giá các nhà hàng dựa trên điểm số chuyên nghiệp từ thang điểm 30 cho ba tiêu chí: Chất lượng đồ ăn (food), phong cách bài trí (decor), chất lượng phục vụ (service), cùng với mức giá trung bình của một bữa ăn (price, USD) và khu vực địa lý (geo).

Hình 8.1: Đánh giá ẩm thực của Zagat đối với nhà hàng Marea tại New York: Điểm số đồ ăn đạt 27/30, khu vực Midtown, đi kèm các mô tả chi tiết về chi phí và trải nghiệm.
Xét tập dữ liệu gồm
text
# Bảng dữ liệu zagat: 168 nhà hàng x 6 biến số
restaurant price food decor service geo
<chr> <dbl> <dbl> <dbl> <dbl> <chr>
1 Daniella Ristorante 43 22 18 20 west
2 Tello's Ristorante 32 20 19 19 west
3 Biricchino 34 21 13 18 west
4 Bottino 41 20 20 17 west
5 Da Umberto 54 24 19 21 west
6 Le Madri 52 22 22 21 west
7 Le Zie 34 22 16 21 west
8 Pasticcio 34 20 18 21 east
9 Belluno 39 22 19 22 east
10 Cinque Terre 44 21 17 19 east
# ... 158 nhà hàng tiếp theoTrong bảng dữ liệu này, biến vị trí địa lý geo là một biến định tính danh nghĩa (Nominal Categorical) nhận hai mức giá trị: Bờ Đông (east) hoặc Bờ Tây (west) Manhattan.

Hình 8.2: Biểu đồ phân tán khảo sát mối liên hệ giữa chất lượng đồ ăn (food) và giá bữa ăn (price), được tô màu phân biệt theo vị trí địa lý (geo: Bờ Đông màu đỏ, bờ Tây màu xanh lam).
Quan sát Hình 8.2, ta nhận thấy:
- Có một mối tương quan dương khá rõ giữa chất lượng đồ ăn và giá tiền: Đồ ăn càng ngon thì giá bữa ăn càng đắt.
- Khi quan sát kỹ theo màu sắc, các nhà hàng ở bờ Tây dường như có mức giá nằm thấp hơn một chút so với các nhà hàng ở bờ Đông có cùng mức điểm đồ ăn.

Hình 8.3: Ước lượng đồ họa bằng hai đường thẳng song song: Một đường cho bờ Đông và một đường cho bờ Tây.
Cơ chế hoạt động của Biến giả (Dummy / Indicator Variable)
Theo hệ thống phân loại dữ liệu, các phép toán đại số như phép nhân và phép cộng chỉ được định nghĩa trên các biến định lượng. Làm thế nào để đưa một biến định tính như geo vào phương trình hồi quy?
Giải pháp toán học mẫu mực là sử dụng Biến giả (Dummy Variable hay Indicator Variable):
Khái niệm: Biến giả (Dummy Variable)
Biến giả là một biến số nhị phân nhận giá trị
Hệ thống thống kê tự động chọn một mức làm Mức cơ sở (Reference level) và tạo biến giả cho các mức còn lại. Vì chữ cái east đứng trước west trong bảng chữ cái tiếng Anh, R mặc định chọn east làm mức cơ sở và tạo ra biến giả geowest:
Khớp mô hình hồi quy trong R:
r
lm(price ~ food + geo, data = zagat)Kết quả ước lượng các hệ số:
text
Coefficients:
(Intercept) food geowest
-15.970 2.875 -1.459Phương trình hồi quy cụ thể:
Bản chất hình học: Hai đường thẳng song song
Để giải mã hình học của phương trình trên, ta thay giá trị của biến giả
- Với các nhà hàng ở bờ Tây (
):
\widehat{\text{price}} &= -15.97 + 2.87 \times \text{food} - 1.46 \times 1 \ &= (-15.97 - 1.46) + 2.87 \times \text{food} \ &= -17.43 + 2.87 \times \text{food} \end
\begin{aligned} \widehat{\text{price}} &= -15.97 + 2.87 \times \text{food} - 1.46 \times 0 \ &= -15.97 + 2.87 \times \text{food} \end
Độ dốc của chất lượng đồ ăn hoàn toàn không đổi (
3. Kịch bản 2: Ba biến định lượng & Mặt phẳng hồi quy
Khi mô hình có từ hai biến giải thích định lượng trở lên, quỹ đạo hình học không còn là các đường thẳng song song nữa.
Xét mô hình giải thích giá bữa ăn (price) đồng thời bằng cả điểm chất lượng đồ ăn (food) và điểm phong cách bài trí (decor):
r
lm(price ~ food + decor, data = zagat)Kết quả ước lượng tham số:
text
Coefficients:
(Intercept) food decor
-24.500 1.646 1.882Phương trình mô hình hồi quy:
Hình học của Mặt phẳng hồi quy trong không gian 3 chiều
Vì cả food và decor đều là các biến định lượng liên tục, mô hình phương trình trên mô tả một Mặt phẳng hồi quy (Regression Plane) trong không gian tọa độ 3 chiều
- Điểm cắt trục tung trong không gian 3 chiều: Khi cả điểm đồ ăn và bài trí đều bằng 0, mặt phẳng cắt trục giá tại
. - Độ dốc theo trục
decor() lớn hơn độ dốc theo trục food(), biểu thị rằng mặt phẳng nghiêng dốc hơn về phía biến phong cách bài trí.
4. Nghệ thuật giải thích hệ số hồi quy: Điều kiện "Giữ các yếu tố khác không đổi"
Điểm mấu chốt sâu sắc nhất khi bước từ hồi quy đơn sang hồi quy bội nằm ở cách diễn giải các hệ số
Trong hồi quy bội, mỗi hệ số hồi quy
Diễn giải các hệ số trong mô hình hai biến định lượng
Xét lại phương trình mô hình:
Hệ số
: - Cách diễn giải chuẩn mực: Khi kiểm soát phong cách bài trí (
decorgiữ không đổi), cứ mỗi điểm tăng thêm trong chất lượng đồ ăn, giá bữa ăn dự kiến sẽ tăng thêm $1.65. - Nói cách khác: Nếu so sánh hai nhà hàng có không gian bài trí đẹp như nhau, nhà hàng có điểm đồ ăn cao hơn 1 điểm sẽ có giá đắt hơn trung bình $1.65.
- Cách diễn giải chuẩn mực: Khi kiểm soát phong cách bài trí (
Hệ số
: - Cách diễn giải chuẩn mực: Khi kiểm soát chất lượng đồ ăn (
foodgiữ không đổi), cứ mỗi điểm tăng thêm trong phong cách bài trí, giá bữa ăn dự kiến sẽ tăng thêm $1.88. - Bài học kinh tế thú vị: Tại Manhattan, việc đầu tư thêm một điểm bài trí không gian thậm chí còn đem lại mức tăng giá cao hơn (+$1.88) so với việc cải thiện chất lượng món ăn (+$1.65).
- Cách diễn giải chuẩn mực: Khi kiểm soát chất lượng đồ ăn (
So sánh với hồi quy đơn biến & Hiện tượng Biến gây nhiễu
Hãy chú ý sự biến thiên của hệ số food:
- Trong mô hình hồi quy đơn chỉ gồm biến
food:(hệ số là ). - Trong mô hình hồi quy bội có thêm biến
decor: Hệ số củafoodgiảm xuống còn.
Tại sao lại có sự sụt giảm này? Lý do là trong thực tế, các nhà hàng có đồ ăn ngon thường cũng đầu tư không gian bài trí đẹp mắt. Trong mô hình đơn biến, biến food đã vô tình "gánh" luôn cả hiệu ứng tăng giá do không gian bài trí đem lại. Khi ta đưa decor vào mô hình bội, hiệu ứng của bài trí được bóc tách riêng ra, giúp ta đo lường được tác động thuần túy của chất lượng đồ ăn khi không gian bài trí đã được cố định.
5. Tóm tắt chương
- Hồi quy tuyến tính bội: Mở rộng mô hình tuyến tính lên
biến giải thích dưới dạng , cho phép nghiên cứu mối quan hệ đa chiều vượt qua giới hạn thị giác của con người. - Biến giả (Dummy variable): Mã hóa biến định tính gồm
mức thành biến nhị phân . Mô hình kết hợp giữa một biến định lượng và một biến giả hai mức sẽ tạo ra hai đường thẳng hồi quy song song có chung độ dốc nhưng khác hệ số chặn. - Mặt phẳng hồi quy: Khi mô hình gồm hai biến giải thích định lượng, quỹ đạo hình học là một mặt phẳng trong không gian 3 chiều, với độ nghiêng theo từng trục được quyết định bởi hệ số tương ứng.
- Nguyên lý kiểm soát biến số (Ceteris Paribus): Hệ số hồi quy riêng phần
phản ánh mức thay đổi kỳ vọng của trên mỗi đơn vị tăng thêm của khi tất cả các biến giải thích còn lại trong mô hình được giữ nguyên cố định.
Tài liệu tham khảo
- James, G., Witten, D., Hastie, T., & Tibshirani, R. (2021). An Introduction to Statistical Learning: with Applications in R (2nd ed.). Springer.
- Wickham, H., & Grolemund, G. (2017). R for Data Science: Import, Tidy, Transform, Visualize, and Model Data. O'Reilly Media.
- UC Berkeley STAT 20: Introduction to Probability and Statistics, Department of Statistics, UC Berkeley.
Bài tập 1: Giải mã Biến giả Đa mức & Hệ phương trình Hồi quy song song
Một nhà nghiên cứu sinh thái học khảo sát thời gian ngủ trung bình mỗi ngày (carni, ăn cỏ herbi, ăn tạp omni).
Chọn nhóm ăn thịt (carni) làm mức cơ sở. Nhà nghiên cứu định nghĩa hai biến giả:
nếu loài ăn cỏ, và cho các trường hợp khác. nếu loài ăn tạp, và cho các trường hợp khác.
Kết quả khớp mô hình hồi quy tuyến tính bội thu được:
- Hãy viết phương trình hồi quy dự báo thời gian ngủ theo
cho từng nhóm động vật riêng biệt ( carni,herbi,omni). - Nhận xét về mối quan hệ hình học giữa ba đường hồi quy trên mặt phẳng tọa độ
. - Khi hai loài động vật có cùng khối lượng cơ thể (
bằng nhau), thời gian ngủ của loài ăn cỏ dự kiến chênh lệch bao nhiêu giờ so với loài ăn thịt? So với loài ăn tạp?
Gợi ý: Gợi ý Bài tập 1
- Với nhóm ăn thịt (
carni), cả hai biến giảvà đều nhận giá trị . - Với nhóm ăn cỏ (
herbi),và . - Các đường thẳng có cùng hệ số trước
sẽ song song với nhau.
Lời giải chi tiết Bài tập 1
- Viết phương trình hồi quy cho từng nhóm động vật:
- Nhóm ăn thịt (
carni): Đây là mức cơ sở,và :
- Nhóm ăn thịt (
\hat{y}_{\text{carni}} &= 15.2 - 0.45 x_1 - 1.80(0) + 0.60(0) \ &= 15.2 - 0.45 x_1 \end
\begin{aligned} \hat{y}_{\text{herbi}} &= 15.2 - 0.45 x_1 - 1.80(1) + 0.60(0) \ &= (15.2 - 1.80) - 0.45 x_1 \ &= 13.4 - 0.45 x_1 \end
\begin{aligned} \hat{y}_{\text{omni}} &= 15.2 - 0.45 x_1 - 1.80(0) + 0.60(1) \ &= (15.2 + 0.60) - 0.45 x_1 \ &= 15.8 - 0.45 x_1 \end
Khi có cùng khối lượng cơ thể, loài ăn cỏ dự kiến ngủ **ít hơn $1.80$ giờ mỗi ngày** so với loài ăn thịt.
- Giữa loài ăn cỏ và loài ăn tạp:Khi có cùng khối lượng cơ thể, loài ăn cỏ dự kiến ngủ ít hơn
giờ mỗi ngày so với loài ăn tạp.
Bài tập 2: Phân tích Cơ chế Biến gây nhiễu & Bóc tách Tác động
Một chuyên gia kinh tế lao động nghiên cứu mức lương khởi điểm (
- Mô hình A (Hồi quy đơn):trong đó
là điểm trung bình tích lũy GPA đại học (thang điểm 4). - Mô hình B (Hồi quy bội):trong đó
là số lượng dự án thực tế sinh viên đã hoàn thành trong thời gian học.
- Hãy diễn giải ý nghĩa của hệ số
trong Mô hình A và hệ số trong Mô hình B. - Tại sao khi bổ sung biến số dự án thực tế
vào mô hình, hệ số của biến GPA lại giảm mạnh từ xuống còn ? - Nếu một sinh viên có GPA bằng 3.5 và đã hoàn thành 4 dự án thực tế, mức lương dự báo theo Mô hình B là bao nhiêu? Nếu mức lương thực tế sinh viên này nhận được là
triệu đồng/tháng thì phần dư bằng bao nhiêu? Sinh viên này đạt mức lương trên hay dưới kỳ vọng của mô hình?
Gợi ý: Gợi ý Bài tập 2
- Chú ý cụm từ điều kiện bắt buộc "Khi kiểm soát số lượng dự án thực tế" khi diễn giải Mô hình B.
- Sinh viên có GPA cao thường có xu hướng hoàn thành nhiều dự án thực tế hơn (mối tương quan giữa
và ). - Tính giá trị dự báo:
, sau đó tính phần dư .
Lời giải chi tiết Bài tập 2
Diễn giải các hệ số:
- Trong Mô hình A (
): Mỗi điểm GPA tăng thêm đi kèm với mức tăng lương khởi điểm kỳ vọng là triệu đồng/tháng. - Trong Mô hình B (
): Khi kiểm soát số lượng dự án thực tế ( giữ không đổi), mỗi điểm GPA tăng thêm đi kèm với mức tăng lương khởi điểm kỳ vọng là triệu đồng/tháng.
- Trong Mô hình A (
Giải thích sự sụt giảm của hệ số GPA:
- Trong thực tế, các sinh viên chăm chỉ có điểm GPA cao thường cũng tích cực tham gia và hoàn thành nhiều dự án thực nghiệm hơn (
và có tương quan dương mạnh). - Trong Mô hình A, biến số dự án
bị bỏ sót (Omitted variable). Hệ số của GPA đã vô tình hấp thu cả phần tác động tăng lương đến từ kỹ năng thực chiến của các dự án. - Khi đưa
vào Mô hình B, tác động của kinh nghiệm dự án được bóc tách độc lập (+ triệu đồng/dự án). Hệ số phản ánh chính xác tác động thuần túy của kiến thức lý thuyết học thuật khi kinh nghiệm dự án đã được cố định ngang nhau.
- Trong thực tế, các sinh viên chăm chỉ có điểm GPA cao thường cũng tích cực tham gia và hoàn thành nhiều dự án thực nghiệm hơn (
Tính giá trị dự báo và phần dư:
- Thay
và vào Mô hình B:
- Thay
\hat{y} &= 6.0 + 1.1(3.5) + 3.2(4) \ &= 6.0 + 3.85 + 12.8 \ &= 22.65 \text{ million VND/month} \end
\begin{aligned} \hat{e} &= y - \hat{y} \ &= 28.0 - 22.65 \ &= +5.35 \text{ million VND/month} \end
Bài tập 3: Mô hình Hồi quy Doanh thu & Ứng dụng Quản trị Kinh doanh
Một chuỗi cửa hàng bán lẻ xây dựng mô hình hồi quy bội để dự báo doanh thu hàng tháng (
: Diện tích mặt sàn kinh doanh ( , đơn vị: ). : Ngân sách tiếp thị quảng cáo trực tuyến (triệu đồng/tháng). : Vị trí cửa hàng ( location, gồm hai mức: Trung tâm thương mạimallvà Mặt phố độc lậpstreet). Chọn mặt phố độc lậpstreetlàm mức cơ sở (nếu nằm trong TTTM, nếu mặt phố).
Phương trình hồi quy ước lượng được:
- Hãy giải thích ý nghĩa thực tế của hệ số
, , và . - Cửa hàng A có diện tích
( ), ngân sách tiếp thị triệu đồng/tháng ( ), đặt tại trung tâm thương mại. Doanh thu dự báo của cửa hàng A là bao nhiêu? - Nếu chủ chuỗi cửa hàng muốn tăng doanh thu dự báo của một cửa hàng mặt phố thêm
trăm triệu đồng (170 triệu đồng) mà không thay đổi ngân sách tiếp thị, cửa hàng đó cần mở rộng diện tích thêm bao nhiêu mét vuông?
Gợi ý: Gợi ý Bài tập 3
- Hệ số
trăm triệu đồng tương đương với bao nhiêu triệu đồng doanh thu trên mỗi triệu đồng tiếp thị? - Với câu 3: Sử dụng mối liên hệ
.
Lời giải chi tiết Bài tập 3
Giải thích các hệ số mô hình:
: Khi kiểm soát ngân sách tiếp thị và vị trí cửa hàng, cứ mỗi diện tích mặt sàn mở rộng thêm, doanh thu hàng tháng của cửa hàng dự kiến tăng thêm trăm triệu đồng (85 triệu đồng). : Khi kiểm soát diện tích và vị trí, cứ mỗi 1 triệu đồng chi thêm cho tiếp thị quảng cáo, doanh thu dự kiến tăng thêm trăm triệu đồng (4 triệu đồng). Tỷ suất doanh thu trên chi phí tiếp thị đạt . : Khi hai cửa hàng có cùng diện tích mặt sàn và cùng ngân sách tiếp thị, cửa hàng đặt tại trung tâm thương mại dự kiến đạt doanh thu hàng tháng cao hơn trăm triệu đồng (210 triệu đồng) so với cửa hàng đặt tại mặt phố độc lập.
Dự báo doanh thu cho cửa hàng A:
- Cửa hàng A có
, , :
- Cửa hàng A có
\hat{y}_A &= 1.2 + 0.85(2) + 0.04(50) + 2.10(1) \ &= 1.2 + 1.70 + 2.00 + 2.10 \ &= 7.00 \text{ hundred million VND} = 700 \text{ million VND/month} \end
- Để
trăm triệu đồng: - Vì đơn vị của
là , mức tăng diện tích cần thiết là: - Kết luận: Cửa hàng cần mở rộng thêm đúng
diện tích sàn kinh doanh.