Tổng quan nghiên cứu

Kể từ khi mô hình nhân tố chung ra đời hơn 100 năm trước từ các công trình nền tảng năm 1947 và phát triển hoàn thiện vào những năm 1969 - 1971, phân tích nhân tố đã trở thành phương pháp thống kê nhiều chiều thiết yếu trong khoa học xã hội, kinh tế và môi trường. Vấn đề nghiên cứu trọng tâm của luận văn là giải quyết những thách thức trong việc phân tích các biến tiềm ẩn từ dữ liệu có thứ bậc, nơi các thang đo thứ tự thường xuyên chịu ảnh hưởng bởi sai số đo lường mà phương pháp hồi quy cổ điển không thể xử lý triệt để.

Mục tiêu cụ thể của công trình gồm hai phần: hệ thống hóa cơ sở toán học của phân tích nhân tố khám phá (EFA) và phân tích nhân tố khẳng định (CFA); đồng thời ứng dụng quy trình phân tích nâng cao này vào bộ dữ liệu thực tế về nhận thức và hành vi thích ứng biến đổi khí hậu. Nghiên cứu được thực hiện hoàn tất vào tháng 6 năm 2017 tại Trường Đại học Khoa học Tự nhiên thuộc Đại học Quốc gia Hà Nội, với dữ liệu khảo sát thực địa thu thập từ 243 người dân tại tỉnh Yên Bái.

Về mặt ý nghĩa thực tiễn, luận văn đã chuyển hóa thành công 34 biến quan sát ban đầu thành 5 nhân tố cấu trúc rõ ràng, giải thích trên 50% tổng phương sai trích xuất. Quy trình này nâng cao độ chính xác trong việc ước lượng tương quan giữa các biến tiềm ẩn, loại bỏ độ lệch do sai số ngẫu nhiên gây ra, đồng thời thiết lập hệ thống chỉ số kiểm định đạt chuẩn với hệ số tải chuẩn hóa vượt ngưỡng 0.50 và hệ số tin cậy trên 0.70.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu vận dụng mô hình nhân tố chung do Thurstone khởi xướng năm 1947 và được mở rộng bởi Jöreskog giai đoạn 1969 - 1971. Mô hình này phân rã tổng phương sai của mỗi biến quan sát thành hai thành phần độc lập: phương sai chung do các biến tiềm ẩn giải thích và phương sai riêng kết hợp giữa phương sai tin cậy với sai số đo lường ngẫu nhiên.

Khung phân tích được đặt trong cấu trúc mô hình phương trình cấu trúc (SEM), phân biệt rạch ròi giữa mô hình đo lường (CFA) và mô hình cấu trúc. Năm khái niệm then chốt được chuẩn hóa bao gồm:

  • Biến tiềm ẩn (Latent construct): Cấu trúc phản ánh không thể quan sát trực tiếp.
  • Biến chỉ số quan sát (Observed indicators): Các biến đo lường thực nghiệm mang tính chất phản ánh.
  • Hệ số tải nhân tố (Factor loading): Hệ số phản ánh độ liên kết giữa biến tiềm ẩn và biến quan sát.
  • Phương sai trích xuất trung bình (AVE): Thước đo giá trị tương đồng của cấu trúc.
  • Chỉ số điều chỉnh (Modification indices): Đại lượng chẩn đoán mức độ tương quan sai số chưa được giải thích trong mô hình.

Phương pháp nghiên cứu

Nguồn dữ liệu của luận văn bao gồm bộ dữ liệu điều tra xã hội học gồm 243 mẫu hợp lệ tại tỉnh Yên Bái được khảo sát trực tiếp vào tháng 6 năm 2017, sử dụng thang đo thứ bậc Likert từ 1 đến 5 điểm. Bên cạnh đó, luận văn còn đối chiếu với các tập dữ liệu thực nghiệm chuẩn gồm 300 mẫu bệnh học và 1.050 quan sát hành vi để kiểm chứng thuật toán.

Phương pháp phân tích được tiến hành qua hai bước liên hoàn:

  • Phân tích nhân tố khám phá (EFA) trên phần mềm SPSS: Áp dụng phương pháp phân tích thành phần chính (PCA) và phân tích trục chính (Principal Axis Factoring), kết hợp phép xoay vuông góc Varimax và xoay xiên Promax.
  • Phân tích nhân tố khẳng định (CFA) trên phần mềm AMOS: Sử dụng phương pháp ước lượng hợp lý cực đại (Maximum Likelihood - ML) để tối thiểu hóa hàm mục tiêu sai biệt ma trận phương sai - hiệp phương sai giữa dữ liệu mẫu và mô hình lý thuyết.

Phương pháp chọn mẫu thuận tiện có phân tầng theo khu vực địa lý được áp dụng. Cỡ mẫu 243 hoàn toàn đáp ứng tiêu chuẩn khắt khe của Hair, vượt mức tối thiểu 200 mẫu cần thiết cho mô hình có 5 cấu trúc với hệ số tương quan chung trung bình dao động từ 0.55 đến 0.60. Toàn bộ tiến trình xử lý, kiểm định ma trận thông tin với công thức xác định bậc tự do toán học được hoàn thành trong giai đoạn 2016 - 2017.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích thực nghiệm trên dữ liệu biến đổi khí hậu tại Yên Bái đã mang lại 4 phát hiện quan trọng:

Thứ nhất, kiểm định tính thích hợp của dữ liệu cho thấy hệ số KMO đạt 0.844, nằm trong khoảng rất tốt từ 0.80 đến 0.90. Kiểm định Bartlett về độ cầu đạt mức ý nghĩa thống kê cao với giá trị sai số sig nhỏ hơn 0.001 (chi-square đạt 1.463 với 210 bậc tự do), bác bỏ giả thuyết các biến không tương quan trong tổng thể.

Thứ hai, kết quả EFA đã rút gọn thành công tập dữ liệu phức tạp thành 5 nhân tố tiềm ẩn có ý nghĩa thực tiễn cao: Ý định thích ứng (F1), Niềm tin vào biến đổi khí hậu (F2), Suy nghĩ chủ quan (F3), Đánh giá khả năng thích ứng (F4), và Nhận thức rủi ro (F5). Tổng phương sai trích đạt trên 55%, trong đó nhân tố F1 giải thích tỷ lệ phương sai lớn nhất với giá trị riêng khởi tạo ban đầu vượt mức 6.0.

Thứ ba, việc so sánh giữa phép xoay Promax và Varimax chứng minh rằng các biến tiềm ẩn trong hành vi con người có tương quan chặt chẽ với nhau. Hệ số tương quan giữa các nhân tố sau phép xoay xiên đạt mức 0.62, phản ánh chân thực hơn bản chất tự nhiên của dữ liệu so với giả định trực giao 90 độ.

Thứ tư, khi chuyển sang kiểm định CFA trên AMOS, việc áp dụng các điều kiện nhận dạng thống kê và cố định tham số tải đơn vị (ULI) bằng 1.0 đã giúp mô hình đạt độ tương thích cao. Chỉ số điều chỉnh (MI) đã phát hiện các cặp sai số tương quan mạnh, tiêu biểu như giữa sai số e1 và e5 với giá trị MI bằng 22.667, giúp định hướng chuẩn xác việc hiệu chỉnh mô hình.

Thảo luận kết quả

Nguyên nhân dẫn đến sự thành công của cấu trúc 5 nhân tố xuất phát từ việc các biến quan sát trong thang đo đã được thiết kế đồng nhất theo lý thuyết phản ánh. Điểm số đánh giá của người dân Yên Bái về các thiệt hại sản xuất và thu nhập có sự đồng biến cao với niềm tin thực tế về biến đổi khí hậu.

Khi đối chiếu với các nghiên cứu kinh điển, kết quả phân tích khẳng định ưu thế vượt trội của phương pháp ước lượng hợp lý cực đại (ML) trong CFA so với phương pháp bình phương nhỏ nhất (OLS). OLS thường giả định dữ liệu không có sai số đo lường, dẫn đến việc ước lượng thấp hơn khoảng 15% đến 25% độ lớn tác động thực tế giữa các cấu trúc. Ngược lại, CFA hiệu chỉnh toàn bộ phương sai riêng, mang lại hệ số đường dẫn chuẩn xác.

Dữ liệu nghiên cứu được trực quan hóa sinh động thông qua biểu đồ lắng đọng (Scree plot) kết hợp phân tích song song 50 lần lặp, cho thấy độ dốc thay đổi rõ rệt tại điểm cắt nhân tố thứ năm trước khi đường giá trị riêng giảm xuống dưới 1.0. Sơ đồ đường dẫn (Path diagram) trên AMOS hiển thị rõ các mũi tên định hướng và hệ số tương quan hiệp phương sai hai chiều, cho phép quan sát trực diện cấu trúc đo lường tổng thể.

Đề xuất và khuyến nghị

Dựa trên kết luận toán học và kết quả thực nghiệm, luận văn đưa ra 4 nhóm khuyến nghị hành động cụ thể:

  • Chuẩn hóa quy trình xây dựng thang đo trong nghiên cứu định lượng: Thiết kế tối thiểu từ 3 đến 4 biến quan sát cho mỗi cấu trúc tiềm ẩn nhằm đảm bảo mô hình luôn ở trạng thái xác định trên. Đảm bảo hệ số tải chuẩn hóa của các biến đạt từ 0.50 trở lên và độ tin cậy tổng hợp đạt trên 0.70. Mục tiêu hoàn thiện bộ công cụ khảo sát chuẩn trong vòng 3 đến 6 tháng cho các đề tài nghiên cứu xã hội. Chủ thể thực hiện là các chuyên gia nghiên cứu và học viên cao học.

  • Mở rộng quy mô mẫu khảo sát theo độ phức tạp của mô hình: Thiết lập quy chuẩn quy mô mẫu tối thiểu từ 200 đến 300 quan sát đối với các mô hình SEM có từ 5 đến 6 biến tiềm ẩn nhằm triệt tiêu hiện tượng Heywood case (phương sai sai số âm). Lộ trình áp dụng liên tục trong các kỳ thu thập dữ liệu hàng năm, do các viện nghiên cứu và trường đại học chủ trì.

  • Tối ưu hóa mô hình đo lường thông qua chẩn đoán sai số chuẩn hóa: Triển khai kiểm tra định kỳ các chỉ số điều chỉnh có MI lớn hơn 4.0 và loại bỏ các biến có số dư chuẩn hóa vượt ngưỡng trị tuyệt đối 2.5 hoặc 4.0. Mục tiêu duy trì chỉ số CMIN/DF dưới ngưỡng 2.0 và chỉ số GFI vượt 0.90 trong mọi báo cáo phân tích, thực hiện bởi các chuyên viên xử lý dữ liệu thống kê.

  • Tích hợp kết quả phân tích nhận thức vào chính sách thích ứng biến đổi khí hậu: Ứng dụng mô hình 5 nhân tố để nâng cao tỷ lệ truyền thông cảnh báo thiên tai và hỗ trợ sinh kế cho 100% các xã vùng cao tại Yên Bái trong giai đoạn 2017 - 2020. Cơ quan thực hiện là Ủy ban nhân dân tỉnh Yên Bái phối hợp cùng Sở Nông nghiệp và Phát triển nông thôn.

Đối tượng nên tham khảo luận văn

Công trình luận văn thạc sĩ này mang lại giá trị học thuật và ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  • Học viên cao học và nghiên cứu sinh khối ngành Kinh tế, Toán học, Xã hội học: Nắm vững nền tảng toán học về ma trận phương sai - hiệp phương sai, công thức tính bậc tự do và quy trình chuyển tiếp chuẩn xác từ EFA sang CFA để hoàn thiện luận văn tốt nghiệp đạt chuẩn quốc tế.

  • Giảng viên và nhà nghiên cứu phương pháp định lượng: Sử dụng làm giáo trình tham khảo chuyên sâu về xử lý dữ liệu thứ bậc, phương pháp ước lượng ML, kỹ thuật phân tích song song và xử lý triệt để các lỗi mô hình như Heywood case trên SPSS và AMOS.

  • Cán bộ quản lý và hoạch định chính sách môi trường tại Yên Bái: Ứng dụng thang đo 5 thành phần đã được kiểm định độ tin cậy để đánh giá mức độ sẵn sàng thích ứng của 243 cụm dân cư, từ đó phân bổ ngân sách phòng chống thiên tai hiệu quả.

  • Chuyên viên nghiên cứu thị trường và phân tích dữ liệu hành vi: Khai thác quy trình tối ưu hóa bảng hỏi khảo sát, loại trừ các biến tải chéo không đạt yêu cầu nhằm xây dựng các chỉ số đo lường trải nghiệm khách hàng với độ chính xác trên 90%.

Câu hỏi thường gặp

Sự khác biệt căn bản giữa phân tích EFA và CFA là gì?

EFA là phương pháp khám phá dữ liệu không áp đặt trước số lượng nhân tố hay cấu trúc tải trọng, phù hợp cho giai đoạn đầu phát triển thang đo. Ngược lại, CFA đòi hỏi cơ sở lý thuyết vững chắc để xác định trước số lượng biến tiềm ẩn, gán các tải trọng chéo bằng 0 và kiểm định sự phù hợp của toàn bộ mô hình bằng phương pháp ML với các chỉ số như CMIN/DF và GFI.

Làm thế nào để xác định số lượng nhân tố tối ưu trong EFA?

Nhà nghiên cứu cần kết hợp đồng thời ba tiêu chuẩn: quy tắc Kaiser-Guttman giữ lại các nhân tố có giá trị riêng khởi tạo lớn hơn 1.0; kiểm tra điểm gãy trên đồ thị Scree plot; và phân tích song song Horn với 50 lần lặp dữ liệu ngẫu nhiên. Số nhân tố được chọn khi giá trị riêng thực tế lớn hơn giá trị riêng trung bình của dữ liệu ngẫu nhiên tương ứng.

Hiện tượng Heywood case trong phân tích CFA là gì và cách xử lý?

Heywood case xảy ra khi kết quả ước lượng phương sai của sai số mang giá trị âm, hoặc hệ số tương quan chuẩn hóa vượt quá trị tuyệt đối 1.0. Nguyên nhân chính là do cỡ mẫu nhỏ dưới 100 hoặc vi phạm quy tắc tối thiểu 3 biến quan sát cho một cấu trúc. Cách xử lý là bổ sung biến quan sát, tăng kích thước mẫu lên trên 200 hoặc loại bỏ các biến đo lường có hệ số tải quá thấp dưới 0.30.

Cỡ mẫu bao nhiêu là đạt chuẩn cho phân tích mô hình CFA và SEM?

Theo tiêu chuẩn thực nghiệm, mô hình có 5 cấu trúc với mỗi cấu trúc chứa hơn 3 biến và hệ số tương quan chung trên 0.60 đòi hỏi cỡ mẫu từ 100 đến 150. Khi hệ số tương quan giảm xuống khoảng 0.55 hoặc có cấu trúc chứa dưới 3 biến, cỡ mẫu bắt buộc phải đạt từ 200 đến 300. Đối với mô hình phức tạp trên 6 cấu trúc, quy mô mẫu cần mở rộng vượt 500 quan sát.

Khi nào nên ưu tiên phép xoay xiên Promax thay vì xoay trực giao Varimax?

Phép xoay trực giao Varimax cố định góc trục nhân tố 90 độ, chỉ phù hợp khi các biến tiềm ẩn hoàn toàn độc lập. Trong khoa học xã hội và hành vi, các nhân tố tâm lý luôn có mối liên hệ nội tại, do đó phép xoay xiên Promax phản ánh đúng thực tế hơn. Nếu tương quan thực tế bằng 0, Promax vẫn cho kết quả tương đương Varimax nhưng giúp chuyển tiếp sang mô hình CFA mượt mà hơn.

Kết luận

  • Luận văn hệ thống hóa toàn diện cơ sở toán thống kê từ mô hình nhân tố chung, thuật toán phân tích thành phần chính PCA đến phương pháp ước lượng hợp lý cực đại ML trong CFA.
  • Xây dựng thành công quy trình xử lý dữ liệu có thứ bậc, phân định rõ ràng vai trò bổ trợ giữa EFA trong giai đoạn dò tìm cấu trúc và CFA trong giai đoạn khẳng định lý thuyết.
  • Kiểm định thực chứng xuất sắc trên bộ dữ liệu 243 mẫu tại Yên Bái, rút gọn 34 biến quan sát thành 5 nhân tố nhận thức và ý định thích ứng biến đổi khí hậu với tổng phương sai giải thích trên 55%.
  • Đóng góp bộ chỉ dẫn phương pháp luận chuẩn mực giúp các nhà nghiên cứu nhận diện điều kiện biên, kiểm soát bậc tự do và khắc phục các sai lệch thống kê thường gặp trên SPSS và AMOS.
  • Định hướng nghiên cứu mở rộng trong giai đoạn 2018 - 2025 hướng tới việc ứng dụng phân tích đa nhóm để so sánh tính bất biến của thang đo đo lường qua nhiều vùng sinh thái khác nhau.

Các nhà nghiên cứu, giảng viên và học viên cao học hãy áp dụng ngay quy trình kiểm định thang đo liên hoàn EFA - CFA chuẩn mực này vào các đề tài nghiên cứu định lượng để đảm bảo tính chính xác và nâng cao chất lượng công bố khoa học.