Tổng quan nghiên cứu
Quản lý và duy trì tỷ lệ sinh viên theo học là bài toán mang tính sống còn đối với các cơ sở giáo dục đại học. Dữ liệu thực nghiệm từ 1.598 sinh viên khóa 48 (niên khóa 2006-2011) tại Trường Đại học Thủy Lợi Hà Nội cho thấy sự biến động học tập của người học chịu ảnh hưởng sâu sắc bởi các yếu tố nhân khẩu học và hoàn cảnh xã hội. Trong phương thức quản trị truyền thống, các trường đại học thường chỉ lưu trữ thông tin hành chính tĩnh như họ tên, quê quán, mã số sinh viên mà chưa đi sâu vào phân tích định lượng chuỗi dữ liệu theo trục thời gian.
Vấn đề cốt lõi mà nghiên cứu giải quyết là sự thiếu hụt các công cụ toán học dự báo biến động học tập. Việc chỉ đánh giá kết cục tốt nghiệp hoặc thôi học dưới dạng biến nhị phân đơn thuần sẽ làm mất đi thông tin quý giá về thời điểm sinh viên gặp khó khăn. Do đó, mục tiêu cụ thể của luận văn là ứng dụng lý thuyết xác suất và thống kê toán học, trọng tâm là phân tích sống sót và mô hình nguy cơ tỷ lệ Cox, nhằm xác định hàm sống sót, hàm nguy cơ và các yếu tố ảnh hưởng trực tiếp đến thời gian theo học của sinh viên.
Phạm vi nghiên cứu được thực hiện trọn vẹn trong chu kỳ đào tạo 5 năm (tương đương 60 tháng) đối với toàn bộ sinh viên khóa 48 phân bố tại 24 lớp chuyên ngành của Trường Đại học Thủy Lợi. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp cơ sở khoa học định lượng giúp nhà trường nâng cao tỷ lệ sinh viên tốt nghiệp đúng hạn từ mức 78% lên trên 85%, đồng thời giảm thiểu tỷ lệ sinh viên bỏ học sớm trong 24 tháng đầu khóa học.
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 được xây dựng trên nền tảng lý thuyết phân tích sống sót (Survival Analysis) – một nhánh quan trọng của thống kê toán học chuyên nghiên cứu khoảng thời gian dẫn đến một sự kiện quan tâm. Hai hàm toán học trung tâm của khung lý thuyết này bao gồm: hàm sống sót S(t) biểu thị xác suất một cá nhân duy trì trạng thái học tập vượt qua thời điểm t, và hàm nguy cơ h(t) đo lường cường độ tức thời xảy ra sự kiện thôi học tại thời điểm t với điều kiện đối tượng vẫn còn theo học ngay trước mốc thời gian đó.
Khung lý thuyết tích hợp ba mô hình phân tích chủ đạo: mô hình phi tham số Kaplan-Meier (1958) dùng để ước tính xác suất sống sót tích lũy từ dữ liệu quan sát; mô hình tham số Weibull với tham số hình dạng m giúp xác định xu hướng nguy cơ tăng, giảm hay không đổi theo thời gian; và mô hình bán tham số Cox (Cox Proportional Hazards Model - 1972). Mô hình Cox được định nghĩa qua phương trình h(t, X) = h0(t) * exp(beta^T * X), trong đó h0(t) là hàm nguy cơ cơ bản tại thời điểm t và vector hệ số beta phản ánh mức độ tác động của các biến giải thích X.
Bốn khái niệm cốt lõi được chuẩn hóa trong nghiên cứu bao gồm:
- Sự kiện (Event): Thời điểm sinh viên thôi học hoặc dừng chương trình đào tạo vì bất kỳ lý do nào.
- Thời gian sống sót (Survival Time): Số tháng tính từ khi sinh viên nhập học (tháng thứ 0) đến khi sự kiện xảy ra hoặc khi kết thúc nghiên cứu (tháng thứ 60).
- Kiểm duyệt phải (Right Censoring): Áp dụng cho những sinh viên tiếp tục theo học bình thường hoặc đã tốt nghiệp tại thời điểm kết thúc 60 tháng theo dõi.
- Tỷ số nguy cơ (Hazard Ratio - HR): Đo lường mức độ chênh lệch về nguy cơ thôi học giữa các nhóm đối tượng nghiên cứu khi thay đổi một đơn vị của biến giải thích.
Phương pháp nghiên cứu
Nghiên cứu sử dụng nguồn dữ liệu toàn bộ gồm 1.598 sinh viên hệ chính quy khóa 48 được cung cấp bởi Phòng Quản lý và Công tác sinh viên Trường Đại học Thủy Lợi Hà Nội. Phương pháp chọn mẫu là chọn mẫu toàn thể (toàn bộ sinh viên nhập học khóa 48), bảo đảm tính đại diện tuyệt đối và triệt tiêu sai số chọn mẫu. Mỗi cá nhân được mã hóa qua vector 6 biến giải thích: lớp chuyên ngành, độ tuổi khi nhập học, giới tính, tỉnh thành cư trú, khu vực ưu tiên và nhóm đối tượng tuyển sinh.
Lý do lựa chọn phân tích sống sót và mô hình Cox thay vì phương pháp hồi quy logistic thông thường là vì phân tích sống sót xử lý hoàn hảo dữ liệu kiểm duyệt phải, không làm biến dạng thông tin của những sinh viên chưa xảy ra sự kiện và tích hợp chặt chẽ yếu tố thời gian liên tục từ tháng thứ 0 đến tháng thứ 60.
Phương pháp phân tích bao gồm: tính toán ước lượng Kaplan-Meier cho từng biến số phân loại; sử dụng kiểm định Log-rank để kiểm tra sự khác biệt có ý nghĩa thống kê giữa các đường cong sống sót đơn biến; và xây dựng mô hình hồi quy đa biến Cox. Các tham số beta trong mô hình Cox được ước lượng thông qua hàm hợp lý từng phần Efron và tối ưu hóa bằng thuật toán lặp Newton-Raphson. Toàn bộ quy trình tính toán, kiểm định giả thuyết Wald và vẽ đồ thị hàm sống sót được thực thi trên môi trường ngôn ngữ thống kê R thông qua gói phần mềm chuyên dụng survival.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Thứ nhất, về cấu trúc nhân khẩu học, sinh viên nam chiếm tỷ lệ áp đảo với 80,91% (1.293 sinh viên), trong khi sinh viên nữ chỉ chiếm 19,09% (305 sinh viên). Độ tuổi nhập học trải dài từ 17 đến 29 tuổi nhưng tập trung chủ yếu ở độ tuổi 18 (686 sinh viên) và 19 tuổi (555 sinh viên), chiếm 77,66% tổng quy mô khóa học. Phân bố địa lý trải rộng trên 53 tỉnh thành, tập trung đông nhất tại Nam Định với 212 sinh viên, Thanh Hóa với 187 sinh viên và Thái Bình với 153 sinh viên.
Thứ hai, về phân bố khu vực tuyển sinh, Khu vực 4 (khu vực 2 - nông thôn) chiếm tỷ trọng vượt trội với 57,51% (919 sinh viên), tiếp đến là Khu vực 2 với 18,34% (293 sinh viên) và Khu vực 1 (vùng cao, miền núi) với 17,46% (279 sinh viên). Khu vực 3 (nội thành đô thị lớn) chiếm tỷ lệ thấp nhất với 6,51% (104 sinh viên). Về diện đối tượng, nhóm đối tượng 10 (lao động phổ thông, không ưu tiên đặc biệt) chiếm đa số với 81,35% (1.300 sinh viên), nhóm đối tượng 1 (dân tộc thiểu số) chiếm 8,20% (131 sinh viên).
Thứ ba, kết quả ước lượng Kaplan-Meier chỉ ra sự phân hóa rõ rệt về xác suất duy trì học tập đến tháng thứ 60 giữa các nhóm đối tượng. Sinh viên thuộc đối tượng 1 có xác suất hoàn thành khóa học đạt 85% (khoảng tin cậy 95%: 78,8% - 91,1%). Nhóm đối tượng 4 và đối tượng 10 đạt xác suất khoảng 79%, đối tượng 6 đạt 78%, đối tượng 3 đạt 67%, và thấp nhất là đối tượng 5 chỉ đạt mức 50% (khoảng tin cậy 95%: 12,5% - 100%).
Thứ tư, nghiên cứu phát hiện hai điểm rơi biến động học tập quan trọng nhất: tháng thứ 11 và tháng thứ 23. Đây là các mốc thời gian ghi nhận số lượng sinh viên thôi học cao nhất trong toàn bộ chu kỳ 60 tháng, trùng khớp với thời điểm kết thúc năm thứ nhất và năm thứ hai đại học.
Thảo luận kết quả
Khi so sánh giữa các vùng miền, sinh viên thuộc Khu vực 3 (nội thành đô thị) có xác suất duy trì học tập thấp nhất, chỉ đạt khoảng 70% đến 75%, thấp hơn đáng kể so với sinh viên Khu vực 1 và Khu vực 4 (cùng đạt xấp xỉ 80%). Nguyên nhân chính của hiện tượng này là sinh viên xuất thân từ khu vực nông thôn và miền núi thường có tính kiên trì cao hơn, định hướng nghề nghiệp gắn bó mật thiết với ngành tài nguyên nước và kỹ thuật công trình. Ngược lại, sinh viên đô thị có nhiều lựa chọn nghề nghiệp thay thế và dễ chuyển hướng sang các ngành kinh tế, công nghệ khác khi gặp khó khăn ban đầu.
Hiện tượng đột biến thôi học tại tháng thứ 11 và tháng thứ 23 được giải thích bởi rào cản học thuật từ các môn đại cương nặng tính toán như Giải tích, Vật lý đại cương, Cơ học lý thuyết và Sức bền vật liệu. Sinh viên không tích lũy đủ số tín chỉ tối thiểu hoặc có điểm trung bình tích lũy dưới chuẩn sẽ bị cảnh báo học vụ và tự ý bỏ cuộc.
Về mặt trình bày dữ liệu, các đường cong sống sót Kaplan-Meier phân tách rõ ràng theo từng bậc thang giảm dần sau mỗi học kỳ. Đồ thị hàm sống sót của Khu vực 3 dốc nhanh hơn hẳn so với Khu vực 1 và 4, phản ánh nguy cơ thôi học tích lũy cao hơn. Bảng hồi quy Cox đa biến lượng hóa chi tiết tỷ số nguy cơ HR của từng nhóm biến, chứng minh vai trò kiểm soát của biến đối tượng và khu vực đối với sự ổn định học tập của sinh viên.
Đề xuất và khuyến nghị
-
Thiết lập hệ thống cảnh báo học vụ sớm ứng dụng mô hình phân tích sống sót:
- Hành động: Xây dựng phần mềm theo dõi biến động học tập tự động tích hợp thuật toán phân tích sống sót, tự động quét dữ liệu chuyên cần và điểm kiểm tra định kỳ của sinh viên.
- Mục tiêu: Giảm 25% tỷ lệ sinh viên thôi học tại các mốc tháng thứ 11 và tháng thứ 23, nâng tỷ lệ hoàn thành khóa học chung lên trên 85%.
- Thời gian thực hiện: Triển khai trong 12 tháng của năm học mới.
- Chủ thể thực hiện: Phòng Quản lý và Công tác sinh viên chủ trì, phối hợp cùng Phòng Đào tạo và Trung tâm Công nghệ thông tin.
-
Triển khai chương trình phụ đạo và cố vấn học tập chuyên sâu cho năm thứ nhất và năm thứ hai:
- Hành động: Tổ chức các lớp bổ trợ miễn phí cho các môn khoa học cơ bản trước các kỳ thi học kỳ, đồng thời phân công giảng viên cố vấn hỗ trợ sinh viên có nguy cơ rủi ro cao.
- Mục tiêu: Nâng tỷ lệ đạt chuẩn học phần đại cương của sinh viên diện ưu tiên lên mức 90%.
- Thời gian thực hiện: Thực hiện liên tục từ tháng thứ 1 đến tháng thứ 24 của mỗi khóa tuyển sinh.
- Chủ thể thực hiện: Ban Chủ nhiệm các Khoa chuyên môn và Đoàn Thanh niên nhà trường.
-
Mở rộng chính sách học bổng và quỹ hỗ trợ tài chính theo đặc thù khu vực:
- Hành động: Huy động nguồn tài trợ từ các doanh nghiệp ngành thủy lợi và trích quỹ phúc lợi để cấp học bổng sinh hoạt phí cho sinh viên Khu vực 1 và Khu vực 4 có hoàn cảnh khó khăn.
- Mục tiêu: Đảm bảo 100% sinh viên nghèo vượt khó đạt học lực khá trở lên không phải bỏ học vì lý do tài chính, duy trì tỷ lệ sống sót của nhóm này trên 88%.
- Thời gian thực hiện: Xét duyệt và giải ngân định kỳ 2 lần mỗi năm học.
- Chủ thể thực hiện: Phòng Quản lý và Công tác sinh viên phối hợp cùng Hội Cựu sinh viên.
-
Đổi mới công tác hướng nghiệp và tư vấn tâm lý học đường:
- Hành động: Tổ chức tuần sinh hoạt công dân đầu khóa gắn với trải nghiệm thực tế công trình thủy lợi, giải tỏa áp lực tâm lý cho sinh viên khu vực đô thị.
- Mục tiêu: Giảm tỷ lệ tự ý chuyển ngành hoặc thôi học của sinh viên Khu vực 3 từ mức 25% xuống dưới 10%.
- Thời gian thực hiện: Thực hiện ngay trong 3 tháng đầu nhập học của tân sinh viên.
- Chủ thể thực hiện: Trung tâm Hướng nghiệp và Truyền thông kết hợp với các chuyên gia tâm lý giáo dục.
Đối tượng nên tham khảo luận văn
- Ban Giám hiệu và các nhà quản lý giáo dục đại học: Cung cấp góc nhìn chiến lược dựa trên dữ liệu định lượng để tối ưu hóa quy chế đào tạo tín chỉ, hoạch định chỉ tiêu tuyển sinh vùng miền và nâng cao hiệu quả quản trị trên quy mô từ 1.500 đến 2.000 sinh viên mỗi khóa.
- Cán bộ phòng công tác sinh viên và cố vấn học tập: Nắm bắt chính xác các nhóm đối tượng có tỷ lệ rủi ro thôi học cao (như đối tượng 3, đối tượng 5 và sinh viên đô thị) cùng các mốc thời gian nhạy cảm (tháng 11 và tháng 23) để thiết kế chương trình can thiệp kịp thời.
- Giảng viên và nhà nghiên cứu chuyên ngành Thống kê toán học, Khoa học dữ liệu: Tiếp cận tài liệu tham khảo giá trị về phương pháp chuyển hóa mô hình phân tích sống sót từ lĩnh vực y sinh sang phân tích các biến cố xã hội và quản lý giáo dục.
- Học viên cao học và sinh viên ngành Toán ứng dụng, Tin học: Học hỏi quy trình xử lý dữ liệu kiểm duyệt phải, cách thức ước lượng tham số bằng thuật toán Newton-Raphson và kỹ năng lập trình chuyên sâu với gói survival trên phần mềm R.
Câu hỏi thường gặp
-
Phân tích sống sót trong giáo dục khác biệt như thế nào so với phân tích hồi quy logistic truyền thống?
Hồi quy logistic chỉ phân loại nhị phân trạng thái tốt nghiệp hay bỏ học mà bỏ qua hoàn toàn trục thời gian. Phân tích sống sót xử lý trọn vẹn dữ liệu kiểm duyệt phải của 1.598 sinh viên trong 60 tháng, giúp xác định chính xác cường độ nguy cơ thôi học thay đổi tại từng thời điểm cụ thể.
-
Tại sao sinh viên Khu vực 1 miền núi lại có tỷ lệ duy trì học tập cao hơn sinh viên Khu vực 3 nội thành?
Số liệu thực nghiệm cho thấy sinh viên Khu vực 1 đạt tỷ lệ duy trì học tập 80%, cao hơn mức 70% đến 75% của Khu vực 3. Sinh viên nông thôn và miền núi có tính kiên định cao, ít bị phân tán bởi các cơ hội việc làm ngoài ngành so với sinh viên đô thị.
-
Thuật toán Newton-Raphson đóng vai trò gì trong việc ước lượng mô hình Cox?
Thuật toán lặp Newton-Raphson được sử dụng để tối đa hóa hàm log hợp lý từng phần Efron trong môi trường R. Thuật toán cập nhật vector hệ số beta qua từng bước lặp ma trận thông tin, đảm bảo tìm ra nghiệm hội tụ chính xác nhất để đánh giá tác động của 6 biến giải thích.
-
Mốc thời gian nào sinh viên có nguy cơ thôi học cao nhất trong toàn bộ 5 năm đào tạo?
Dữ liệu chỉ ra hai mốc nguy cơ cao nhất là tháng thứ 11 (kết thúc năm thứ nhất) và tháng thứ 23 (kết thúc năm thứ hai). Đây là giai đoạn chuyển giao kiến thức từ đại cương sang cơ sở ngành, tạo ra áp lực đào thải học thuật lớn nhất đối với sinh viên.
-
Gói thư viện survival trong ngôn ngữ R hỗ trợ những công cụ phân tích then chốt nào?
Gói survival cung cấp hàm Surv để thiết lập đối tượng thời gian và trạng thái kiểm duyệt, hàm survfit để tính toán đường cong sống sót Kaplan-Meier kèm khoảng tin cậy 95%, và hàm coxph để ước lượng mô hình nguy cơ tỷ lệ Cox đa biến nhanh chóng và chuẩn xác.
Kết luận
- Luận văn đã ứng dụng thành công phương pháp phân tích sống sót và mô hình nguy cơ tỷ lệ Cox trên bộ dữ liệu thực tế gồm 1.598 sinh viên khóa 48 Trường Đại học Thủy Lợi Hà Nội.
- Xác lập bức tranh toàn diện về hàm sống sót và chỉ ra hai mốc thời gian có nguy cơ thôi học cao nhất là tháng thứ 11 và tháng thứ 23 trong vòng đời đào tạo 60 tháng.
- Lượng hóa rõ nét sự phân hóa xác suất duy trì học tập theo khu vực và đối tượng ưu tiên, trong đó sinh viên nông thôn và miền núi đạt tỷ lệ duy trì khoảng 80%, vượt trội so với sinh viên thành thị.
- Đóng góp quy trình phân tích dữ liệu giáo dục chuẩn hóa bằng ngôn ngữ R, khắc phục triệt để hạn chế của các phương pháp thống kê mô tả truyền thống.
- Cung cấp hệ thống 4 giải pháp khả thi nhằm nâng cao tỷ lệ sinh viên tốt nghiệp đúng hạn lên trên 85% trong chu kỳ 5 năm tiếp theo.
Nghiên cứu mở ra hướng tiếp cận hiện đại trong chuyển đổi số giáo dục đại học, kết hợp chặt chẽ giữa toán học ứng dụng và quản trị đào tạo. Các cơ sở giáo dục, nhà nghiên cứu và bạn đọc quan tâm có thể khai thác và mở rộng mô hình này để nâng cao chất lượng quản lý sinh viên tại đơn vị của mình.