Tổng quan luận án

Luận án tiến sĩ chuyên ngành Hệ thống thông tin (mã ngành 62480104) của nghiên cứu sinh Huỳnh Phước Hải, do PGS. Đỗ Thanh Nghị và TS. Nguyễn Văn Hòa hướng dẫn khoa học tại Trường Đại học Cần Thơ (năm 2019), mang tiêu đề "Mô hình xử lý hiệu quả dữ liệu biểu hiện gen". Luận án giải quyết bài toán phân lớp dữ liệu biểu hiện gen nhằm hỗ trợ chẩn đoán và điều trị bệnh ung thư ở người.

Tính cấp thiết và khoảng trống nghiên cứu

Bệnh ung thư là một trong những nguyên nhân hàng đầu gây tử vong trên thế giới, với thống kê từ Tổ chức Y tế Thế giới (WHO) ghi nhận 18,1 triệu ca mắc mới và 9,6 triệu ca tử vong trong năm 2018; đồng thời được dự báo sẽ gia tăng 20% đến 40% trong giai đoạn 2020–2030 (Rahib et al., 2014). Việc áp dụng các phác đồ điều trị chung thường mang lại hiệu quả khác nhau trên từng bệnh nhân do sự phức tạp của các biến dị di truyền. Do đó, hướng tiếp cận y học cá thể hóa (personalized medicine) dựa trên hồ sơ di truyền và mức độ biểu hiện gen đóng vai trò thiết yếu trong việc nâng cao độ chính xác của công tác chẩn đoán và điều trị.

Tuy nhiên, dữ liệu biểu hiện gen thu được từ các công nghệ đo đạc hiện đại (như DNA Microarray và RNA-Seq) có đặc thù là số chiều rất lớn (hàng chục nghìn đến hàng trăm nghìn gen/chiều) trong khi số lượng mẫu bệnh phẩm thu thập được lại rất nhỏ (thường chỉ vài chục đến vài trăm mẫu do chi phí xét nghiệm cao). Hiện tượng này dẫn đến thách thức "số chiều lớn, số mẫu nhỏ" (large p, small n) và "lời nguyền số chiều" (curse of dimensionality - Bellman, 1961). Dữ liệu bị phân tán thưa thớt trong không gian số chiều lớn khiến các mô hình học máy thông thường dễ rơi vào tình trạng quá khớp (overfitting), đạt kết quả phân loại rất cao trên tập huấn luyện nhưng suy giảm nghiêm trọng trên tập kiểm thử độc lập.

Mục tiêu nghiên cứu

Mục tiêu tổng quát của luận án là đề xuất các phương pháp tiếp cận mới cho bài toán phân lớp dữ liệu biểu hiện gen nhằm nâng cao độ chính xác phân loại bệnh ung thư. Các mục tiêu cụ thể được luận án xác định bao gồm:

  1. Nghiên cứu xây dựng mô hình rút trích đặc trưng hiệu quả cho dữ liệu biểu hiện gen để nâng cao độ chính xác phân lớp.
  2. Nghiên cứu xây dựng mô hình tăng cường dữ liệu cho dữ liệu biểu hiện gen để nâng cao độ chính xác phân lớp.
  3. Nghiên cứu xây dựng mô hình phân lớp hiệu quả dữ liệu biểu hiện gen.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các mô hình rút trích đặc trưng phi tuyến, mô hình sinh dữ liệu nhân tạo và mô hình phân lớp dữ liệu biểu hiện gen bệnh ung thư ở người.
  • Phạm vi nghiên cứu: Bài toán phân lớp nhị phân và đa lớp trên dữ liệu biểu hiện gen thu được từ các thí nghiệm y sinh trên người thuộc hai nền tảng công nghệ chính: DNA Microarray (các chip Affymetrix Oligonucleotide Arrays) và giải trình tự RNA-Seq thế hệ mới (Next Generation Sequencing), sử dụng nguồn dữ liệu công khai từ các ngân hàng dữ liệu Kent Ridge, ArrayExpress và The Cancer Genome Atlas (TCGA).

Tổng quan tài liệu và vị trí của luận án

Luận án đã tổng hợp và phân tích các công trình nghiên cứu liên quan đến bài toán phân lớp dữ liệu biểu hiện gen theo ba nhóm tiếp cận kỹ thuật chính:

1. Hướng tiếp cận giảm chiều và rút trích đặc trưng

  • Chọn tập con đặc trưng (Feature Selection): Các công trình của Golub et al. (1999) với kỹ thuật Signal-to-Noise Ratio (S2N), Li et al. (2001) với giải thuật di truyền, Dudoit et al. (2002) với Linear Discriminant Analysis (LDA), hay Guyon et al. (2002) với SVM-RFE. Hướng tiếp cận này chọn lọc các gen có trọng số cao dựa trên các kiểm định thống kê hoặc mô hình bao đóng. Tuy nhiên, việc loại bỏ bớt các gen có thể gây thất thoát thông tin quan trọng và độ chính xác phân lớp bị phụ thuộc chặt chẽ vào từng tập dữ liệu cụ thể.
  • Rút trích đặc trưng tuyến tính (Feature Extraction): Các phương pháp như Phân tích thành phần chính (PCA), Phân tích thành phần độc lập (ICA), Bình phương tối thiểu từng phần (PLS), và Nhúng tuyến tính cục bộ (LLE) (Bicciato et al., 2003; Dai et al., 2006). Nhược điểm của nhóm phương pháp này là chỉ nắm bắt được mối quan hệ tuyến tính, chưa khai thác được các quan hệ phi tuyến tiềm ẩn giữa các gen trong không gian biểu hiện phức tạp.

2. Hướng tiếp cận sử dụng các thuật toán phân lớp cơ sở và mô hình tập hợp

  • Máy học véc-tơ hỗ trợ (SVM): Được áp dụng rộng rãi bởi Vapnik (1995), Guyon et al. (2002), Yeoh et al. (2002), Furey et al. (2000), Statnikov et al. (2005). SVM thể hiện khả năng chịu nhiễu tốt và kiểm soát lề phân hoạch trong không gian nhiều chiều.
  • Thuật toán k láng giềng (k-NN) và Cây quyết định (C4.5/CART): Cover & Hart (1967), Quinlan (1993), Breiman et al. (1984), Tan et al. (2005). Các mô hình này gặp hạn chế lớn khi chạy trực tiếp trên không gian số chiều gốc do sự phân tán dữ liệu và việc phân tách đơn thuộc tính song song với trục tọa độ làm mất tính liên kết đa chiều giữa các gen.
  • Mô hình tập hợp (Ensemble Learning): Rừng ngẫu nhiên (Random Forests - Breiman, 2001; Diaz-Uriarte & Alvarez de Andres, 2006), Bagging (Breiman, 1996; Dudoit & Fridlyand, 2003), Boosting / AdaBoost (Freund & Schapire, 1997; Dettling & Buhlmann, 2003). Các nghiên cứu cải tiến như rừng xiên phân đa thuộc tính (Do et al., 2010) cho thấy tiềm năng nâng cao độ chính xác phân loại.

3. Hướng tiếp cận tăng cường dữ liệu và học sâu

  • Kỹ thuật SMOTE (Chawla et al., 2002) ban đầu được thiết kế để giải quyết mất cân bằng lớp trên dữ liệu thông thường.
  • Mạng nơ-ron tích chập sâu (DCNN) và Mạng đối kháng sinh mẫu (GAN - Goodfellow et al., 2014) được ứng dụng mạnh trong thị giác máy tính và xử lý ngôn ngữ tự nhiên, nhưng việc thiết kế cấu trúc DCNN và GAN chuyên biệt cho ma trận biểu hiện gen dạng bảng còn là khoảng trống nghiên cứu chưa được khai thác hoàn chỉnh.

Khoảng trống luận án lựa chọn giải quyết

Luận án tập trung giải quyết đồng thời hai thách thức căn bản của dữ liệu biểu hiện gen: số chiều cực lớn ($p$) và số mẫu rất nhỏ ($n$). Luận án xây dựng mô hình DCNN để rút trích đặc trưng phi tuyến bậc cao; phát triển mô hình GAN chuyên biệt để tăng cường mẫu dữ liệu tổng hợp; kết hợp DCNN với SMOTE; và đề xuất mô hình tập hợp cây xiên phân ngẫu nhiên đơn giản (Bag-RODS và Boost-RODS) dựa trên siêu phẳng tối ưu của SVM tuyến tính để phân lớp trực tiếp hiệu quả trên không gian số chiều gốc cũng như trên không gian đặc trưng mới.

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

Cơ sở lý thuyết và khung phân tích

  • Biểu hiện gen và hệ phiên mã: Dữ liệu biểu hiện gen phản ánh số lượng bản phiên mã mRNA được tạo ra trong tế bào qua hai giai đoạn phiên mã (transcription) và dịch mã (translation). Ma trận dữ liệu sau chuẩn hóa có dạng bảng $m \times n$, trong đó các dòng đại diện cho mẫu bệnh nhân ($m$) và các cột đại diện cho mức biểu hiện của các gen ($n$).
  • Lý thuyết học thống kê và Máy học véc-tơ hỗ trợ (SVM): Tìm siêu phẳng phân cách tối ưu $w \cdot x - b = 0$ bằng cách cực đại hóa khoảng cách lề $2/|w|$ và cực tiểu hóa sai số huấn luyện qua bài toán quy hoạch toàn phương: $$\min_{\alpha} \frac{1}{2} \sum_{i=1}^m \sum_{j=1}^m y_i y_j \alpha_i \alpha_j K(x_i, x_j) - \sum_{i=1}^m \alpha_i$$ ràng buộc $\sum_{i=1}^m y_i \alpha_i = 0$ và $0 \le \alpha_i \le C, \forall i=1,\dots,m$. Hàm nhân tuyến tính $K(x_i, x_j) = x_i \cdot x_j$ và hàm nhân phi tuyến RBF $K(x_i, x_j) = \exp(-\gamma |x_i - x_j|^2)$. Đối với bài toán đa lớp, luận án áp dụng phương pháp 1-vs-1 và 1-vs-all.
  • Phân rã sai số Bias - Variance: Làm cơ sở lý thuyết cho việc kết hợp các bộ phân loại yếu thành bộ phân loại mạnh bằng kỹ thuật lấy mẫu Bootstrap (Bagging) và điều chỉnh phân phối trọng số mẫu (Boosting).
  • Mạng nơ-ron tích chập sâu (DCNN) & Mạng đối kháng sinh mẫu (GAN): Chuyển đổi dữ liệu 1 chiều của gen thành cấu trúc ma trận vuông phù hợp với các tầng tích chập (Convolutional Layers) và gộp (Pooling Layers); triển khai mạng Generator ($G$) và Discriminator ($D$) trong GAN với hàm mục tiêu minimax để học phân phối dữ liệu biểu hiện gen gốc và sinh mẫu nhân tạo.

Phương pháp nghiên cứu và quy trình thực nghiệm

Luận án sử dụng phương pháp nghiên cứu thực nghiệm định lượng trên các tập dữ liệu thực nghiệm tiêu chuẩn. Quy trình gồm các bước: thu thập dữ liệu chuẩn hóa, thiết kế kiến trúc mô hình đề xuất, tinh chỉnh siêu tham số, phân lớp dữ liệu và kiểm định ý nghĩa thống kê.

+-----------------------------------------------------------------------------------+
|                            DỮ LIỆU BIỂU HIỆN GEN GỐC                              |
|         (DNA Microarray: 50 tập; RNA-Seq: 25 tập 2 lớp & 1 tập lớn 36 lớp)        |
+-----------------------------------------+-----------------------------------------+
                                          |
          +-------------------------------+-------------------------------+
          |                               |                               |
          v                               v                               v
+-------------------+           +-------------------+           +-------------------+
|  HƯỚNG TIẾP CẬN 1 |           |  HƯỚNG TIẾP CẬN 2 |           |  HƯỚNG TIẾP CẬN 3 |
| Rút trích đặc     |           | Tăng cường dữ     |           | Phân lớp trực tiếp|
| trưng bằng DCNN & |           | liệu bằng mạng    |           | bằng Bag-RODS và  |
| kết hợp SMOTE     |           | GAN (gán nhãn SVM)|           | Boost-RODS        |
+---------+---------+           +---------+---------+           +---------+---------+
          |                               |                               |
          +-------------------------------+-------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
|                               ĐÁNH GIÁ VÀ SO SÁNH                                 |
|       - Giao thức: Train-Test độc lập, Leave-One-Out (LOO), 10-Fold CV            |
|       - Tiêu chí: Độ chính xác tổng thể (Accuracy %), Thời gian thực thi          |
|       - Kiểm định thống kê: Mức ý nghĩa p-value < 0,05                            |
+-----------------------------------------------------------------------------------+

Giao thức đánh giá và kiểm định

Theo khuyến nghị của dự án Kiểm chuẩn Chất lượng Microarray (MAQC-II):

  1. Chia tách có sẵn (trn-tst): Sử dụng tập huấn luyện và kiểm tra độc lập theo công bố gốc.
  2. Leave-One-Out (LOO): Áp dụng cho các tập dữ liệu có số mẫu $m < 300$.
  3. 10-Fold Cross-Validation: Áp dụng cho các tập dữ liệu có số mẫu $m \ge 300$.
  4. Tiêu chí: Độ chính xác phân loại tổng thể $\text{Accuracy} = \frac{TP + TN}{TP + TN + FP + FN} \times 100%$ và kiểm định thống kê $p$-value ($p < 0,05$).

Nguồn tư liệu và tập dữ liệu thực nghiệm

Luận án thu thập và sử dụng tổng cộng 4 nhóm dữ liệu thực nghiệm từ các ngân hàng dữ liệu quốc tế:

  • 50 tập dữ liệu DNA Microarray: Từ Kent Ridge và ArrayExpress, thuộc 3 dòng chip Affymetrix phổ biến: Human Gene 1.0 ST (A-AFFY-141, 33.297 chiều), Human Genome U133 Plus 2.0 (A-AFFY-44, 54.675 chiều), Human Genome U133A (A-AFFY-33, 22.283 chiều). Quy mô từ 31 đến 677 mẫu, từ 2.000 đến 54.755 chiều, số lớp từ 2 đến 16 lớp.
  • 20 tập dữ liệu DNA Microarray mẫu nhỏ: Số lượng mẫu dưới 130 bệnh nhân để đánh giá chuyên biệt mô hình GAN.
  • 25 tập dữ liệu RNA-Seq nhị phân: Lấy từ kho dữ liệu TCGA, mỗi tập có 20.531 gen/chiều, số mẫu từ 66 đến 1.100 mẫu.
  • 1 tập dữ liệu RNA-Seq lớn: Tổng hợp từ TCGA/GDAC gồm 36 loại bệnh ung thư và mẫu bình thường, tổng cộng 12.181 mẫu bệnh nhân.

Bảng tổng hợp đặc tính các nhóm dữ liệu thực nghiệm trong luận án:

Nhóm dữ liệu Nguồn Số tập Số mẫu (min - max) Số chiều/gen Số lớp Giao thức đánh giá
DNA Microarray chuẩn Kent Ridge, ArrayExpress 50 31 – 677 2.000 – 54.755 2 – 16 LOO, 10-fold, trn-tst
DNA Microarray mẫu nhỏ ArrayExpress 20 47 – 130 2.000 – 54.675 2 – 14 LOO
RNA-Seq 2 lớp TCGA 25 66 – 1.100 20.531 2 LOO ($<300$), 10-fold ($\ge 300$)
RNA-Seq đa lớp quy mô lớn TCGA / GDAC 1 12.181 20.531 37 (36 u + 1 chuẩn) 10-fold CV

Nội dung chính theo từng chương

Chương 1: Giới thiệu

Chương 1 thiết lập cơ sở khoa học và bối cảnh nghiên cứu. Tác giả chỉ ra tính cấp thiết của việc xử lý dữ liệu biểu hiện gen phục vụ y học cá thể hóa trong chẩn đoán ung thư; phân tích bản chất của thách thức "large p, small n" và nguy cơ quá khớp của các mô hình học máy. Tác giả xác định 3 nhiệm vụ nghiên cứu tương ứng với 3 hướng tiếp cận: rút trích đặc trưng, tăng cường dữ liệu và xây dựng mô hình phân lớp mới trên không gian chiều gốc. Chương này cũng trình bày 3 đóng góp khoa học chính và cấu trúc 7 chương của luận án.

Chương 2: Cơ sở lý thuyết và các công trình liên quan

Chương 2 trình bày chi tiết các khái niệm sinh học phân tử (DNA, mRNA, protein, phiên mã, dịch mã, hệ phiên mã) và nguyên lý hoạt động của hai nền tảng công nghệ DNA Microarray và RNA-Seq. Chương này chuẩn hóa công thức toán học cho bài toán phân lớp dữ liệu biểu hiện gen và thiết lập giao thức đánh giá mô hình (LOO, 10-fold, độ chính xác tổng thể, kiểm định $p < 0,05$). Phần trọng tâm của chương là khảo sát có hệ thống các thuật toán học máy: k-NN, cây quyết định (C4.5, CART), SVM (tuyến tính LSVM và RBF), các phương pháp Ensemble (Random Forests, Bagging, AdaBoost) và các mô hình học sâu. Tác giả tổng hợp chi tiết thông số của 50 tập dữ liệu DNA Microarray, 20 tập DNA Microarray mẫu nhỏ, 25 tập RNA-Seq 2 lớp và 1 tập RNA-Seq đa lớp lớn.

Chương 3: Mô hình rút trích đặc trưng cho dữ liệu biểu hiện gen

Chương 3 đề xuất kiến trúc mạng nơ-ron tích chập sâu (DCNN) chuyên biệt để rút trích các đặc trưng tiềm ẩn phi tuyến từ dữ liệu biểu hiện gen.

  • Kiến trúc DCNN đề xuất: Véc-tơ biểu hiện gen $1 \times N$ được chuyển đổi thành ma trận vuông kích thước $K \times K$ (với $K = \lceil \sqrt{N} \rceil$, các ô thiếu được bổ sung giá trị 0). Mạng bao gồm 2 tầng tích chập (CONV1, CONV2) kết hợp hàm kích hoạt ReLU, xen kẽ với 2 tầng lấy mẫu cực đại (POOLING1, POOLING2), và kết thúc bằng tầng kết nối đầy đủ (Fully Connected) để trích xuất véc-tơ đặc trưng mới có số chiều rút giảm đáng kể nhưng mang tính đại diện phân lớp cao.
  • Kết quả thực nghiệm:
    • Trên 50 tập DNA Microarray: DCNN kết hợp SVM (DCNN-SVM) và DCNN kết hợp LSVM (DCNN-LSVM) đạt độ chính xác trung bình vượt trội so với khi phân lớp trực tiếp trên không gian số chiều gốc. DCNN-SVM cải thiện độ chính xác trên 42/50 tập dữ liệu với độ tin cậy thống kê $p < 0,05$.
    • Trên 25 tập RNA-Seq: DCNN-SVM, DCNN-LSVM và DCNN-kNN đều cải thiện độ chính xác phân lớp so với các mô hình SVM, LSVM và k-NN ban đầu.
    • Trên tập RNA-Seq lớn (12.181 mẫu, 36 loại ung thư): Mô hình DCNN kết hợp bộ phân loại cho kết quả phân loại đa lớp vượt trội so với các thuật toán học máy truyền thống chạy trên toàn bộ 20.531 đặc trưng gốc.

Chương 4: Mô hình tăng cường mẫu đặc trưng rút trích bằng SMOTE

Chương 4 kết hợp mô hình DCNN rút trích đặc trưng (Chương 3) với kỹ thuật sinh mẫu tổng hợp SMOTE nhằm giải quyết đồng thời cả hai thách thức: giảm số chiều và tăng số lượng mẫu huấn luyện.

  • Quy trình: Dữ liệu biểu hiện gen gốc đầu tiên đi qua mạng DCNN để rút trích véc-tơ đặc trưng chiều thấp. Sau đó, thuật toán SMOTE được áp dụng trên không gian đặc trưng mới này để sinh thêm các mẫu tổng hợp cho tập huấn luyện dựa trên khoảng cách giữa các $k$ láng giềng gần nhất. Dữ liệu sau khi tăng cường được dùng để huấn luyện các bộ phân loại SVM, LSVM, k-NN và Random Forest (RF).
  • Kết quả thực nghiệm trên 50 tập DNA Microarray:
    • Mô hình kết hợp DCNN-SMOTE-SVM đạt độ chính xác trung bình cao nhất trong số 15 mô hình so sánh.
    • DCNN-SMOTE-LSVM và DCNN-SMOTE-RF đều cho độ chính xác cao hơn rõ rệt so với các mô hình chỉ áp dụng DCNN đơn thuần hoặc chỉ phân lớp trên số chiều gốc ($p < 0,05$).

Chương 5: Mô hình tăng cường dữ liệu cho dữ liệu biểu hiện gen

Chương 5 đề xuất mô hình tăng cường dữ liệu trực tiếp trên không gian biểu hiện gen bằng Mạng đối kháng sinh mẫu (GAN).

  • Kiến trúc GAN cho dữ liệu biểu hiện gen: Mạng sinh $G$ nhận véc-tơ nhiễu ngẫu nhiên và sinh ra mẫu biểu hiện gen nhân tạo có phân phối tương đồng với dữ liệu gốc; mạng phân biệt $D$ đánh giá tính chân thực của mẫu. Dữ liệu nhân tạo sinh ra sau khi mạng GAN hội tụ sẽ được gán nhãn tự động bằng mô hình máy học véc-tơ hỗ trợ (SVM) đã huấn luyện trên tập gốc, sau đó bổ sung vào tập huấn luyện ban đầu.
  • Kết quả thực nghiệm trên 20 tập DNA Microarray mẫu nhỏ ($m < 130$):
    • GAN-SVM (kết hợp GAN với SVM nhân RBF) đạt độ chính xác trung bình cao nhất, cải thiện rõ rệt so với SVM gốc.
    • GAN-LSVM, GAN-RF, GAN-kNN và GAN-C4.5 đều ghi nhận sự gia tăng độ chính xác phân loại so với các phiên bản không có tăng cường dữ liệu.

Chương 6: Mô hình tập hợp cây xiên phân ngẫu nhiên đơn giản

Chương 6 đề xuất hai thuật toán phân lớp mới: Bagging Random Oblique Decision Stumps (Bag-RODS) và Boosting Random Oblique Decision Stumps (Boost-RODS) nhằm phân lớp hiệu quả trực tiếp trên số chiều gốc cũng như trên không gian dữ liệu đã qua xử lý.

  • Cấu trúc cây xiên phân ngẫu nhiên đơn giản (RODS): Thay vì phân chia đơn thuộc tính song song với trục tọa độ như cây quyết định C4.5, RODS chọn ngẫu nhiên một tập con các đặc trưng ($n' \ll n$), sau đó áp dụng SVM tuyến tính để tìm một siêu phẳng tối ưu làm điều kiện phân tách tại gốc (tạo thành một stump xiên phân đa thuộc tính có khả năng chịu nhiễu cao).
  • Mô hình Bag-RODS & Boost-RODS: Bag-RODS kết hợp nhiều RODS độc lập huấn luyện trên các mẫu Bootstrap theo cơ chế bỏ phiếu số đông. Boost-RODS huấn luyện tuần tự các RODS bằng cách điều chỉnh phân phối trọng số các mẫu bị phân loại sai ở các vòng lặp trước.
  • Kết quả thực nghiệm:
    • Trên 50 tập DNA Microarray (phân lớp trực tiếp): Bag-RODS và Boost-RODS đạt độ chính xác trung bình cao hơn rõ rệt so với SVM, LSVM, k-NN, C4.5, RF, Bag-C4.5 và AdaBoost ($p < 0,05$).
    • Khi kết hợp với GAN và DCNN: Mô hình GAN-Bag-RODS, GAN-Boost-RODS trên 20 tập mẫu nhỏ và DCNN-Bag-RODS, DCNN-Boost-RODS trên 50 tập Microarray tiếp tục nâng cao hiệu năng phân loại tổng thể.

Chương 7: Kết luận và hướng phát triển

Tổng kết toàn diện các kết quả lý thuyết và thực nghiệm đạt được trong 6 chương; khẳng định tính khả thi của 3 hướng giải pháp; và định hướng các nghiên cứu tiếp theo.

Kết quả và những đóng góp mới

Luận án đã đóng góp 3 kết quả khoa học chính trong lĩnh vực xử lý và phân lớp dữ liệu biểu hiện gen:

  1. Đóng góp 1: Đề xuất mô hình rút trích đặc trưng mới bằng mạng nơ-ron tích chập sâu (DCNN) cho dữ liệu biểu hiện gen DNA Microarray và RNA-Seq; đồng thời phát triển giải pháp tích hợp DCNN với kỹ thuật tăng cường mẫu SMOTE (DCNN-SMOTE). Giải pháp giải quyết đồng thời hai thách thức: giảm số chiều dữ liệu phi tuyến và bù đắp số lượng mẫu nhỏ, cải thiện độ chính xác phân loại của các bộ phân loại SVM, LSVM, k-NN, RF trên 50 tập dữ liệu chuẩn. (Công bố tại các bài báo: CT1 [56], CT3 [57], CT4 [58], CT6 [59], CT7 [60]).
  2. Đóng góp 2: Đề xuất mô hình tăng cường dữ liệu biểu hiện gen bằng mạng đối kháng sinh mẫu (GAN) kết hợp gán nhãn tự động bằng SVM. Mô hình sinh ra các mẫu dữ liệu tổng hợp có phân phối tương tự dữ liệu sinh học gốc, giúp cải thiện độ chính xác phân loại của các thuật toán SVM, LSVM, k-NN, RF và C4.5 trên 20 tập dữ liệu Microarray có số mẫu nhỏ. (Công bố tại các bài báo: CT2 [62], CT7 [60]).
  3. Đóng góp 3: Đề xuất mô hình tập hợp cây xiên phân ngẫu nhiên đơn giản dựa trên SVM tuyến tính với hai thuật toán Bag-RODS và Boost-RODS. Mô hình tận dụng khả năng phân hoạch tối ưu của siêu phẳng SVM trên các tập con thuộc tính ngẫu nhiên, cho độ chính xác phân loại vượt trội so với các thuật toán SVM, LSVM, RF, k-NN, C4.5, Bagging và AdaBoost khi chạy trực tiếp trên không gian số chiều gốc, đồng thời tương thích tốt khi kết hợp với DCNN và GAN. (Công bố tại bài báo: CT5 [65]).

Kiến nghị và giải pháp đề xuất

Luận án đề xuất quy trình ứng dụng thực tế các mô hình học máy vào hệ thống hỗ trợ chẩn đoán y tế:

  • Đối với dữ liệu có số chiều cực lớn và quan hệ gen phức tạp: Ưu tiên sử dụng mạng DCNN để rút trích đặc trưng phi tuyến trước khi phân loại.
  • Đối với các thí nghiệm có số lượng bệnh nhân hạn chế ($m < 130$): Áp dụng mạng GAN hoặc kỹ thuật DCNN-SMOTE để sinh mẫu dữ liệu huấn luyện bổ sung.
  • Đối với bài toán phân lớp trực tiếp không qua tiền xử lý giảm chiều: Sử dụng giải thuật Bag-RODS hoặc Boost-RODS để đảm bảo độ chính xác cao và khả năng chống quá khớp.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu

  • Các mô hình học sâu (DCNN, GAN) đòi hỏi thời gian huấn luyện lớn hơn và cần tài nguyên phần cứng (GPU) so với các mô hình học máy truyền thống.
  • Các đặc trưng phi tuyến sau khi rút trích bằng DCNN hoặc các mẫu sinh bằng GAN mang tính trừu tượng toán học cao, khó diễn giải trực tiếp ý nghĩa sinh học phân tử cụ thể của từng gen riêng lẻ so với phương pháp chọn lọc gen (Feature Selection) truyền thống.
  • Phạm vi dữ liệu thực nghiệm chủ yếu tập trung vào các thí nghiệm microarray bulk và RNA-Seq khối u mô đồng nhất, chưa mở rộng sang dữ liệu giải trình tự đơn tế bào (single-cell RNA-Seq).

Hướng nghiên cứu tiếp theo

Theo định hướng tại Mục 7.2 của luận án:

  1. Nghiên cứu cơ chế tự động tối ưu hóa kiến trúc và siêu tham số cho mạng DCNN và GAN đối với từng loại bệnh ung thư cụ thể.
  2. Mở rộng thử nghiệm mô hình trên dữ liệu biểu hiện gen đơn tế bào (single-cell RNA-Seq) và tích hợp đa dữ liệu sinh học phân tử (Multi-omics: DNA methylation, microRNA, biểu hiện gen).
  3. Phát triển các kỹ thuật giải thích mô hình (Explainable AI) kết hợp tri thức y sinh học để diễn giải cơ chế sinh học của các siêu phẳng trong mô hình Bag-RODS và Boost-RODS.

Giá trị tham khảo

Luận án mang lại giá trị học thuật và ứng dụng thực tiễn cho nhiều nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Hệ thống thông tin, Khoa học dữ liệu và Tin sinh học: Cung cấp tài liệu tham khảo có hệ thống về cơ sở toán học của SVM, Bias-Variance, cấu trúc DCNN/GAN trên dữ liệu ma trận dạng bảng, và các giải thuật Ensemble Learning nâng cao.
  • Các nhà nghiên cứu y sinh học tính toán: Cung cấp bộ benchmark thực nghiệm toàn diện gồm thông tin chi tiết của 50 tập dữ liệu DNA Microarray, 20 tập DNA Microarray mẫu nhỏ, 25 tập RNA-Seq và 1 tập RNA-Seq quy mô lớn (12.181 mẫu) đã qua chuẩn hóa.
  • Kỹ sư phát triển hệ thống y tế thông minh: Cung cấp các khung thuật toán (Bag-RODS, Boost-RODS, DCNN-SMOTE, GAN-SVM) có thể tích hợp vào các phần mềm hỗ trợ chẩn đoán tự động giai đoạn bệnh, phân nhóm phân tử khối u và hỗ trợ ra quyết định điều trị ung thư.

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

1. Thách thức "large p, small n" trong phân lớp dữ liệu biểu hiện gen gây ra vấn đề gì cho các mô hình học máy?

Thách thức "large p, small n" biểu thị tình trạng dữ liệu có số chiều/số gen rất lớn ($p$ từ hàng chục nghìn đến hàng trăm nghìn) nhưng số lượng mẫu bệnh nhân lại rất nhỏ ($n$ từ vài chục đến vài trăm). Do số chiều quá cao, các điểm dữ liệu phân tán thưa thớt trong không gian, khiến không gian giả thuyết của mô hình học máy trở nên quá rộng. Hậu quả là mô hình rất dễ bị "quá khớp" (overfitting) — đạt độ chính xác gần như tuyệt đối trên tập huấn luyện nhưng dự đoán sai lệch lớn trên tập dữ liệu kiểm thử độc lập.

2. Mô hình DCNN chuyển đổi véc-tơ biểu hiện gen 1 chiều thành đầu vào mạng tích chập như thế nào?

Để áp dụng được phép toán tích chập 2 chiều của mạng DCNN, luận án chuyển đổi véc-tơ biểu hiện gen 1 chiều có kích thước $1 \times N$ thành một ma trận vuông 2 chiều kích thước $K \times K$, trong đó $K = \lceil \sqrt{N} \rceil$. Các vị trí còn thiếu trong ma trận vuông (nếu $K^2 > N$) được điền giá trị 0. Ma trận vuông này đóng vai trò như một bản đồ đặc trưng đầu vào để các tầng tích chập (CONV) và tầng gộp (POOLING) học các mối quan hệ phi tuyến cục bộ và toàn cục giữa các mức độ biểu hiện gen.

3. Mạng GAN tăng cường dữ liệu biểu hiện gen theo cơ chế nào và mẫu sinh ra được gán nhãn ra sao?

Mạng GAN gồm mạng sinh $G$ (học phân phối của dữ liệu gen để sinh ra các mẫu nhân tạo từ véc-tơ nhiễu) và mạng phân biệt $D$ (đánh giá phân biệt mẫu thực và mẫu sinh). Sau khi quá trình huấn luyện đối kháng minimax đạt điểm cân bằng, mạng $G$ sinh ra các mẫu biểu hiện gen tổng hợp. Các mẫu nhân tạo này sau đó được đưa qua bộ phân loại SVM (đã huấn luyện tối ưu trên tập dữ liệu gốc) để gán nhãn tự động trước khi được gộp vào tập huấn luyện mở rộng.

4. Khác biệt kỹ thuật cốt lõi giữa cây xiên phân ngẫu nhiên đơn giản (RODS) và cây quyết định truyền thống (như C4.5) là gì?

Cây quyết định truyền thống như C4.5 phân hoạch không gian dữ liệu bằng các siêu phẳng song song với trục tọa độ (phân tách đơn thuộc tính), điều này làm mất đi sự liên kết đồng thời giữa nhiều gen trong không gian chiều cao. Ngược lại, RODS chọn ngẫu nhiên một tập con các thuộc tính ($n' \ll n$) và sử dụng thuật toán SVM tuyến tính để tìm một siêu phẳng xiên tối ưu (phân tách đa thuộc tính) tại gốc. Điều này giúp RODS đạt khả năng phân tách dữ liệu hiệu quả hơn, chịu đựng nhiễu tốt hơn và giảm thiểu sai số bias lẫn variance khi kết hợp trong mô hình Bagging (Bag-RODS) và Boosting (Boost-RODS).

5. Giao thức đánh giá thực nghiệm nào được tác giả áp dụng cho các tập dữ liệu có quy mô mẫu khác nhau?

Tác giả tuân thủ khuyến nghị của dự án kiểm chuẩn MAQC-II:

  • Đối với các tập dữ liệu có sẵn phân chia huấn luyện/kiểm thử từ công bố gốc: Sử dụng trực tiếp phân chia độc lập (trn-tst).
  • Đối với các tập dữ liệu không có phân chia sẵn và có số lượng mẫu $m < 300$: Sử dụng giao thức Leave-One-Out (LOO) nhằm tận dụng tối đa số mẫu để huấn luyện.
  • Đối với các tập dữ liệu có số lượng mẫu $m \ge 300$: Sử dụng giao thức kiểm chuẩn chéo 10 lần (10-fold cross-validation). Mọi kết quả so sánh đều được kiểm định với mức ý nghĩa $p < 0,05$.

Kết luận

Luận án tiến sĩ của Huỳnh Phước Hải đã giải quyết có hệ thống bài toán phân lớp dữ liệu biểu hiện gen phục vụ chẩn đoán bệnh ung thư trong điều kiện số chiều lớn và số mẫu nhỏ. Bằng việc đề xuất thành công mô hình rút trích đặc trưng phi tuyến DCNN kết hợp SMOTE, mô hình tăng cường dữ liệu bằng mạng GAN và hai thuật toán phân lớp mới Bag-RODS, Boost-RODS, luận án đã nâng cao rõ rệt độ chính xác phân loại trên cả hai nền tảng DNA Microarray và RNA-Seq. Toàn bộ các mô hình đề xuất đều được kiểm chứng chặt chẽ thông qua hệ thống thực nghiệm quy mô lớn trên 50 tập dữ liệu chuẩn và kiểm định thống kê $p < 0,05$. Kết quả nghiên cứu là tài liệu tham khảo học thuật có giá trị cao cho các nghiên cứu thuộc chuyên ngành Hệ thống thông tin, Khoa học dữ liệu và Tin sinh học ứng dụng.