Phân tích biến kiểm soát trong AMOS: chọn đúng, đặc tả đúng và báo cáo minh bạch
Biến kiểm soát không phải danh sách nhân khẩu học thêm vào để mô hình trông chặt chẽ hơn. Mỗi biến được điều chỉnh làm thay đổi đại lượng đang ước lượng và có thể giảm hoặc tạo thêm sai lệch.
Quy trình an toàn bắt đầu bằng xác định estimand, vẽ quan hệ nhân quả dự kiến, chọn adjustment set trước khi xem p-value, mã hóa đúng trong AMOS và so hệ số chưa chuẩn hóa cùng khoảng tin cậy giữa mô hình cơ sở với mô hình điều chỉnh.
Biến kiểm soát cần lý do và kiểm tra độ nhạy
- Cơ sở lý thuyết
Giải thích vì sao cần kiểm soát
- Mã hóa
Khai báo loại biến và nhóm gốc
- Đặc tả
Đưa vào đúng phương trình
- Độ nhạy
So sánh kết quả có và không control
Sơ đồ khái niệm. Thêm nhiều biến kiểm soát không tự làm suy luận đáng tin hơn.
Nguồn: Lê Minh biên soạn · Minh họa khái niệm, không phải đầu ra phần mềm.1. Mục tiêu học tập và estimand phải xác định trước
Sau chuyên đề này, bạn cần phân biệt confounder với mediator và collider, quyết định biến nào nên điều chỉnh, đưa biến liên tục hoặc phân loại vào AMOS, rồi đánh giá kết luận chính có nhạy với tập control hay không.
Trước khi chọn control, phải nêu rõ đại lượng cần ước lượng. Total effect của X lên Y bao gồm các đường trung gian hợp lệ. Direct effect giữ mediator cố định trả lời câu hỏi khác. Nếu mục tiêu thay đổi giữa hai mô hình, sự thay đổi của hệ số X đến Y không thể được gọi đơn giản là mô hình đã kiểm soát tốt hơn (Cinelli, Forney và Pearl, 2024).
- Đã nêu rõ X, Y và hiệu ứng cần ước lượng.
- Đã xác định thời điểm đo của từng biến.
- Đã phân biệt total effect với direct effect.
- Tập control được chọn trước khi xem kết quả chính.
- Mỗi control có căn cứ lý thuyết hoặc thiết kế.
- Dữ liệu chứa phép đo đủ tin cậy cho control quan trọng.
2. Thêm control là một giả định nhân quả
Trong dữ liệu quan sát, một control thường được dùng để chặn đường back-door giữa X và Y. Trường hợp dễ hiểu nhất là Z gây ra cả X và Y. Điều chỉnh Z có thể giúp tách phần liên hệ không đi qua đường nhân quả X đến Y, nhưng kết luận này phụ thuộc cấu trúc giả định chứ không phụ thuộc riêng tương quan hoặc p-value (Rohrer, 2018).
Spector và Brannick (2011) cảnh báo việc thêm control một cách tự động với kỳ vọng làm sạch quan hệ. Một biến có tương quan với Y chưa chắc là confounder. Một biến không có ý nghĩa thống kê vẫn có thể cần điều chỉnh nếu thiết kế và lý thuyết cho thấy nó chặn đường gây nhiễu.
AMOS ước lượng mô hình covariance đã đặc tả. Phần mềm không thể xác định một mũi tên control là nguyên nhân, mediator hay collider. Trách nhiệm nhận dạng đại lượng vẫn thuộc về thiết kế nghiên cứu và lập luận của người phân tích.
Không loại một confounder chỉ vì p lớn hơn 0,05 và không thêm một biến chỉ vì nó dự báo Y. Quyết định phải dựa trên vai trò của biến trong cấu trúc nhân quả dự kiến.
3. Confounder, mediator, collider và hậu quả của X
Confounder là nguyên nhân chung của X và Y hoặc là biến giúp chặn một đường gây nhiễu thích hợp. Mediator nằm trên đường X đến Y. Collider là hậu quả chung của hai biến. Điều chỉnh mediator khi mục tiêu là total effect sẽ chặn một phần hiệu ứng cần ước lượng; điều chỉnh collider có thể mở một đường liên hệ giả (Cinelli và cộng sự, 2024).
Biến được đo trước X thường an toàn hơn biến hậu xử lý, nhưng pre-treatment không tự động có nghĩa là control tốt. Vai trò của một biến còn phụ thuộc các biến khác trong adjustment set. Cùng một Z có thể hữu ích trong một cấu trúc và gây bias trong cấu trúc khác.
Nếu biến giới tính, tuổi hoặc thâm niên chỉ là proxy cho cơ chế thật sự quan trọng, cần nói rõ giới hạn. Spector và Brannick (2011) khuyến nghị tránh dùng nhân khẩu học thay cho khái niệm lý thuyết mà nghiên cứu thực sự quan tâm.
| Vai trò dự kiến | Ví dụ cấu trúc | Xử lý khi ước lượng total effect | Rủi ro |
|---|---|---|---|
| Confounder | Z gây X và Y | Thường cần điều chỉnh | Bỏ sót có thể tạo confounding |
| Mediator | X gây Z, Z gây Y | Không điều chỉnh nếu cần total effect | Biến total effect thành direct effect |
| Collider | X gây Z và U gây Z | Không điều chỉnh | Mở đường liên hệ giả |
| Hậu quả của X | X gây Z | Thường không dùng làm control | Post-treatment bias |
| Predictor chỉ của Y | Z gây Y, không gây X | Có thể tăng precision nếu hợp lệ | Không tự động giải quyết confounding |
4. Chọn và đo lường control trước khi chạy AMOS
Becker (2005) và Bernerth và Aguinis (2016) yêu cầu giải thích biến nào được xem xét, vì sao được đưa vào hoặc loại ra, cách đo và kết quả có phụ thuộc lựa chọn đó hay không. Câu vì nghiên cứu trước đã kiểm soát là chưa đủ nếu nghiên cứu trước không có cùng câu hỏi nhân quả.
Control đo kém có thể không chặn đầy đủ confounding. Nếu control là một khái niệm tiềm ẩn quan trọng, dùng một chỉ báo đơn không đáng tin có thể để lại residual confounding. Khi có thang đo nhiều item và cỡ mẫu phù hợp, có thể mô hình hóa control như latent variable; khi dùng điểm tổng, phải công khai cách tạo điểm và độ tin cậy.
Nên lập một bảng quyết định trước phân tích gồm tên biến, thời điểm đo, vai trò dự kiến, đường cần điều chỉnh, cách mã hóa và phân tích độ nhạy. Bảng này giảm nguy cơ thêm hoặc bớt control sau khi nhìn thấy kết quả.
| Biến | Vai trò dự kiến | Thời điểm | Cách đo và mã hóa | Quyết định |
|---|---|---|---|---|
| Tuổi | Confounder theo lý thuyết lĩnh vực | Trước X | Số năm, giữ liên tục | Đưa vào |
| Thâm niên | Proxy gần với tuổi | Trước X | Số năm | Kiểm tra đa cộng tuyến |
| Hài lòng | Mediator giả thuyết | Sau X | Latent 4 item | Không control khi ước lượng total effect |
| Nhóm ngành | Khác biệt nền | Trước X | K nhóm, tạo K - 1 dummy | Đưa vào nếu phù hợp DAG |
5. Đặc tả biến kiểm soát trong AMOS
Trong IBM SPSS Amos 29, control quan sát được đưa vào như biến ngoại sinh nếu nó không có phương trình giải thích trong mô hình. Vẽ đường một chiều từ control đến đúng biến nội sinh cần điều chỉnh. Cho control hiệp phương sai với các biến ngoại sinh liên quan khi cấu trúc lý thuyết cho phép. Không nối control vào mọi biến nội sinh theo thói quen (IBM Corp., 2022).
Nếu control là latent variable, cần đặc tả mô hình đo lường và kiểm tra nghiệm như các cấu trúc khác. Nếu control là điểm tổng quan sát, sai số đo không được AMOS tự sửa. Hai cách này không tương đương.
Với dữ liệu thiếu, estimator và cơ chế missing phải được giữ nhất quán giữa mô hình cơ sở với mô hình điều chỉnh. So sánh hai mô hình trên hai tập trường hợp khác nhau có thể làm thay đổi hệ số vì mẫu đổi, không chỉ vì control.
Bước 1
Chốt estimand, adjustment set và mô hình cơ sở trước khi mở AMOS.
Bước 2
Kiểm tra phạm vi, missing, coding và thời điểm đo của từng control.
Bước 3
Đưa control quan sát hoặc latent vào đúng dạng và đặt nhóm tham chiếu nếu dùng dummy.
Bước 4
Vẽ đường từ control đến các biến nội sinh theo causal graph, không theo p-value.
Bước 5
Cho covariance với predictor ngoại sinh khi quan hệ nền được giả định và cần ước lượng.
Bước 6
Ước lượng mô hình cơ sở và mô hình điều chỉnh trên cùng tập dữ liệu.
Bước 7
So hệ số chưa chuẩn hóa, SE, CI, dấu và kết luận của đường trọng tâm.
Bước 8
Lưu file mô hình, output và nhật ký lý do cho từng control.

6. Mã hóa biến phân loại và tránh đổi câu hỏi nghiên cứu
Biến danh nghĩa K nhóm cần K - 1 biến dummy khi được đưa vào như predictor quan sát. Nhóm bị bỏ là nhóm tham chiếu, nên mọi hệ số phải được diễn giải so với nhóm này. Đổi nhóm tham chiếu thay đổi cách viết hệ số, không thay đổi các giá trị dự báo của mô hình.
Không mã hóa các nhóm không có thứ tự thành 1, 2, 3 rồi dùng như biến liên tục. Cách đó áp một khoảng cách và xu hướng tuyến tính không được thiết kế hỗ trợ. Với ordinal control, quyết định giữ liên tục hay tạo dummy phải dựa trên số mức, hình dạng quan hệ và mục tiêu.
Nếu mục tiêu là kiểm tra hệ số X đến Y có khác giữa các nhóm, biến nhóm đang đóng vai trò moderator. Khi đó, multigroup SEM và measurement invariance thường phù hợp hơn việc gọi nhóm là control.
Control điều chỉnh mức trung bình hoặc quan hệ nền theo đặc tả. Moderator hỏi hiệu ứng X đến Y có thay đổi theo Z hay không. Một biến có thể được dùng trong hai vai trò, nhưng mô hình và cách diễn giải phải tách biệt.
7. So sánh mô hình cơ sở và mô hình điều chỉnh
Không chỉ nhìn p-value của control. Bằng chứng chính là hệ số trọng tâm thay đổi ra sao khi adjustment set hợp lý được áp dụng. Nên ưu tiên hệ số chưa chuẩn hóa và khoảng tin cậy vì hệ số chuẩn hóa có thể đổi do phương sai được tái phân bổ giữa các mô hình.
Fit tốt hơn không chứng minh control đúng. Thêm đường thường giảm chi-square trong mẫu. Kline (2023) nhấn mạnh phải đọc fit cùng đặc tả, tham số và mục tiêu suy luận. Một control không có căn cứ vẫn có thể cải thiện fit.
Khi các control tương quan mạnh, sai số chuẩn có thể tăng và dấu hệ số dao động. Cần xem correlation, variance inflation trong mô hình hồi quy tương ứng hoặc condition diagnostics, đồng thời tránh biến nhiều proxy gần nhau thành adjustment set mặc định.
| Thành phần | Mô hình cơ sở | Mô hình điều chỉnh | Câu hỏi |
|---|---|---|---|
| b của X đến Y | Giá trị và CI | Giá trị và CI | Độ lớn hoặc dấu có đổi không? |
| N phân tích | Số trường hợp | Số trường hợp | Mẫu có giữ nguyên không? |
| R bình phương của Y | Giá trị | Giá trị | Khả năng giải thích thay đổi bao nhiêu? |
| Fit và residuals | Hồ sơ fit | Hồ sơ fit | Có sai đặc tả mới không? |
| Estimand | Total hoặc direct effect | Total hoặc direct effect | Hai hệ số còn trả lời cùng câu hỏi không? |

8. Ví dụ số giả định: control tốt và mediator bị dùng sai
Ví dụ dưới đây là minh họa, không phải dữ liệu nghiên cứu thật. Mục tiêu là total effect của X lên Y. Mô hình cơ sở trên N = 420 cho b = 0,42, SE = 0,07, CI 95% [0,28; 0,56]. Lý thuyết xác định tuổi là nguyên nhân chung của X và Y. Sau khi điều chỉnh tuổi, b = 0,31, SE = 0,08, CI 95% [0,15; 0,47]. Kết luận về dấu và độ bất định vẫn ổn định, nhưng độ lớn giảm.
Một mô hình khác thêm M, trong khi M được giả thuyết là mediator do X tác động. Hệ số trực tiếp X đến Y còn b = 0,08, SE = 0,07, CI 95% [-0,05; 0,21]. Kết quả này không chứng minh X hết tác động sau kiểm soát. Nó chuyển câu hỏi từ total effect sang direct effect giữ M cố định.
Nếu bài báo cần cả hai câu hỏi, hãy báo mô hình total effect và mô hình mediation riêng. Không xếp M chung với tuổi trong một danh sách control vì hai biến có vai trò nhân quả khác nhau.
| Mô hình | Đặc tả | b X đến Y | CI 95% | Estimand |
|---|---|---|---|---|
| M0 | Không control | 0,42 | [0,28; 0,56] | Liên hệ chưa điều chỉnh |
| M1 | Điều chỉnh tuổi | 0,31 | [0,15; 0,47] | Total effect có adjustment set giả định |
| M2 | Điều chỉnh tuổi và mediator M | 0,08 | [-0,05; 0,21] | Direct effect giữ M cố định |
9. Phân tích độ nhạy cho lựa chọn control
Một bộ kết quả có thể ổn định với adjustment set chính nhưng nhạy với các control còn tranh luận. Hãy xác định trước mô hình chính, sau đó chạy một số đặc tả có lý do: bỏ proxy trùng lặp, thay cách mã hóa phi tuyến hoặc thêm một control nền còn bất định.
Bernerth và Aguinis (2016) khuyến nghị minh bạch ảnh hưởng của control lên kết luận trọng tâm. Phân tích độ nhạy không phải thử mọi tổ hợp rồi chọn mô hình thuận lợi. Mỗi đặc tả phải đại diện cho một giả định có thể giải thích.
Nếu một confounder quan trọng không được đo, thêm nhiều nhân khẩu học không bảo đảm hết omitted variable bias. Cần nêu giới hạn và, khi phù hợp, dùng công cụ sensitivity analysis ngoài AMOS để định lượng mức confounding chưa quan sát cần thiết để đảo kết luận.
Thử nhiều tập control rồi chỉ báo mô hình cho hệ số đẹp là lựa chọn phụ thuộc kết quả. Hãy công khai mô hình chính, các đặc tả độ nhạy và lý do của từng thay đổi.
10. Các lỗi thường gặp khi dùng biến kiểm soát
| Cách làm yếu | Vấn đề | Cách xử lý |
|---|---|---|
| Thêm mọi biến nhân khẩu học | Không có vai trò nhân quả rõ | Chọn control theo estimand và DAG |
| Giữ control khi p nhỏ | Selection phụ thuộc kết quả | Chốt trước khi chạy |
| Điều chỉnh mediator cho total effect | Chặn đường cần ước lượng | Tách total và direct effect |
| Điều chỉnh collider | Có thể mở đường liên hệ giả | Rà cấu trúc nguyên nhân chung |
| Mã hóa ba nhóm thành 1, 2, 3 | Áp xu hướng tuyến tính giả | Dùng K - 1 dummy hoặc mô hình phù hợp |
| Chỉ báo hệ số sau control | Không thấy mức nhạy | Báo mô hình cơ sở và điều chỉnh |
| So beta chuẩn hóa duy nhất | Thang chuẩn hóa đổi giữa mô hình | Ưu tiên b chưa chuẩn hóa và CI |
11. Mẫu báo cáo APA 7 và bảng kiểm hoàn tất
Mẫu sau dùng số giả định ở mục 8. Khi áp dụng, thay số, tên biến và vai trò nhân quả bằng thông tin của nghiên cứu thật.
Mẫu báo cáo tiếng Việt
Biến tuổi được chọn làm biến kiểm soát trước phân tích vì được giả thuyết là nguyên nhân chung của X và Y. Tuổi được giữ ở đơn vị năm và được cho phép hiệp phương sai với X; một đường trực tiếp từ tuổi đến Y được ước lượng trong IBM SPSS Amos [phiên bản]. Trong mô hình cơ sở, hiệu ứng chưa chuẩn hóa của X lên Y là b = .42, SE = .07, CI 95% [.28, .56]. Sau điều chỉnh, hệ số là b = .31, SE = .08, CI 95% [.15, .47]. Dấu và kết luận suy luận không đổi, nhưng độ lớn giảm. Mô hình có và không control được ước lượng trên cùng N = 420. Mediator M không được đưa vào adjustment set vì mục tiêu chính là total effect.
- Đã nêu estimand là total hay direct effect.
- Mỗi control có vai trò và thời điểm đo rõ.
- Không chọn control bằng p-value.
- Không điều chỉnh mediator hoặc collider ngoài mục tiêu.
- Biến phân loại dùng coding và nhóm tham chiếu rõ ràng.
- Mô hình cơ sở và điều chỉnh dùng cùng dữ liệu.
- Đã báo b chưa chuẩn hóa, SE và CI.
- Đã xem đa cộng tuyến giữa các control.
- Đã báo ảnh hưởng của control lên kết luận chính.
- Đã lưu nhật ký các control được cân nhắc nhưng không dùng.
- Đã thực hiện độ nhạy cho control còn tranh luận.
- Đã nêu giới hạn của confounding chưa quan sát.
Kết luận
Biến kiểm soát là một phần của câu hỏi nhân quả, không phải phụ kiện cho mô hình. Thêm sai biến có thể thay đổi estimand hoặc tạo bias, dù output AMOS vẫn hội tụ và fit đẹp.
Báo cáo đáng tin cậy cho biết vì sao từng control được chọn, được mã hóa và nối như thế nào, hệ số trọng tâm thay đổi bao nhiêu, và kết luận có ổn định dưới các đặc tả hợp lý khác hay không.
Tài liệu tham khảo
- Becker, T. E. (2005). Potential problems in the statistical control of variables in organizational research: A qualitative analysis with recommendations. Organizational Research Methods, 8(3), 274-289. https://doi.org/10.1177/1094428105278021
- Bernerth, J. B., & Aguinis, H. (2016). A critical review and best-practice recommendations for control variable usage. Personnel Psychology, 69(1), 229-283. https://doi.org/10.1111/peps.12103
- Cinelli, C., Forney, A., & Pearl, J. (2024). A crash course in good and bad controls. Sociological Methods & Research, 53(3), 1071-1104. https://doi.org/10.1177/00491241221099552
- Rohrer, J. M. (2018). Thinking clearly about correlations and causation: Graphical causal models for observational data. Advances in Methods and Practices in Psychological Science, 1(1), 27-42. https://doi.org/10.1177/2515245917745629
- Spector, P. E., & Brannick, M. T. (2011). Methodological urban legends: The misuse of statistical control variables. Organizational Research Methods, 14(2), 287-305. https://doi.org/10.1177/1094428110369842
- Kline, R. B. (2023). Principles and practice of structural equation modeling (5th ed.). Guilford Press. https://www.guilford.com/books/Principles-and-Practice-of-Structural-Equation-Modeling/Rex-Kline/9781462551910
- IBM Corp. (2022). IBM SPSS Amos 29 user's guide. https://www.ibm.com/docs/en/SSLVMB_29.0.0/pdf/IBM_SPSS_Amos_User_Guide.pdf
Ghi chú học thuật: Causal graph làm rõ giả định nhưng không chứng minh giả định đúng. Với dữ liệu quan sát cắt ngang, cần giới hạn ngôn ngữ nhân quả và công khai nguy cơ confounding chưa đo.
