Đầu trang
Bỏ qua đến nội dung chính

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.

Bắt đầu đọc

Biến kiểm soát cần lý do và kiểm tra độ nhạy

  1. Cơ sở lý thuyết

    Giải thích vì sao cần kiểm soát

  2. Mã hóa

    Khai báo loại biến và nhóm gốc

  3. Đặc tả

    Đưa vào đúng phương trình

  4. Độ 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 chọn control bằng p-value

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ò của biến và quyết định điều chỉnh
Vai trò dự kiếnVí dụ cấu trúcXử lý khi ước lượng total effectRủi ro
ConfounderZ gây X và YThường cần điều chỉnhBỏ sót có thể tạo confounding
MediatorX gây Z, Z gây YKhông điều chỉnh nếu cần total effectBiến total effect thành direct effect
ColliderX gây Z và U gây ZKhông điều chỉnhMở đường liên hệ giả
Hậu quả của XX gây ZThường không dùng làm controlPost-treatment bias
Predictor chỉ của YZ gây Y, không gây XCó 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ả.

Nhật ký lựa chọn biến kiểm soát
BiếnVai trò dự kiếnThời điểmCách đo và mã hóaQuyết định
TuổiConfounder theo lý thuyết lĩnh vựcTrước XSố năm, giữ liên tụcĐưa vào
Thâm niênProxy gần với tuổiTrước XSố nămKiểm tra đa cộng tuyến
Hài lòngMediator giả thuyếtSau XLatent 4 itemKhông control khi ước lượng total effect
Nhóm ngànhKhác biệt nềnTrước XK 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.

Mô hình AMOS có các đường từ biến kiểm soát đến biến nội sinh trong bài gốc. Hình không thể hiện đầy đủ cách mã hóa biến phân loại; cần khai báo dummy phù hợp trước khi ước lượng.
Ảnh thao tác 1. Mô hình AMOS có các đường từ biến kiểm soát đến biến nội sinh trong bài gốc. Hình không thể hiện đầy đủ cách mã hóa biến phân loại; cần khai báo dummy phù hợp trước khi ước lượng. Nguồn: amosleminh.com · Xem ảnh đầy đủ

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 và moderator không cùng câu hỏi

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.

Bảng đối chiếu tối thiểu giữa các mô hình
Thành phầnMô hình cơ sởMô hình điều chỉnhCâu hỏi
b của X đến YGiá trị và CIGiá trị và CIĐộ lớn hoặc dấu có đổi không?
N phân tíchSố trường hợpSố trường hợpMẫu có giữ nguyên không?
R bình phương của YGiá trịGiá trịKhả năng giải thích thay đổi bao nhiêu?
Fit và residualsHồ sơ fitHồ sơ fitCó sai đặc tả mới không?
EstimandTotal hoặc direct effectTotal hoặc direct effectHai hệ số còn trả lời cùng câu hỏi không?
Bảng Regression Weights trong ví dụ nguồn có Estimate, S.E., C.R. và P cho các control. Đây là hệ số chưa chuẩn hóa; không quyết định giữ hay bỏ control chỉ từ p-value.
Ảnh thao tác 2. Bảng Regression Weights trong ví dụ nguồn có Estimate, S.E., C.R. và P cho các control. Đây là hệ số chưa chuẩn hóa; không quyết định giữ hay bỏ control chỉ từ p-value. Nguồn: amosleminh.com · Xem ảnh đầy đủ

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.

Diễn giải ba đặc tả minh họa
Mô hìnhĐặc tảb X đến YCI 95%Estimand
M0Không control0,42[0,28; 0,56]Liên hệ chưa điều chỉnh
M1Điều chỉnh tuổi0,31[0,15; 0,47]Total effect có adjustment set giả định
M2Điều chỉnh tuổi và mediator M0,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.

Không dùng specification search cho control

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

Lỗi và cách sửa
Cách làm yếuVấn đềCách xử lý
Thêm mọi biến nhân khẩu họcKhô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 effectChặn đường cần ước lượngTách total và direct effect
Điều chỉnh colliderCó 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 controlKhông thấy mức nhạyBáo mô hình cơ sở và điều chỉnh
So beta chuẩn hóa duy nhấtThang 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

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.