Radiomics: Feature selection 2 기본 콘텐츠로 건너뛰기

Radiomics: Feature selection 2

 앞서 Radiomics에서 많이 사용되고 있는 Feature selection 방법에 대해서 이야기 하였다. 이번에는 조금 더 세분화하여 설명해보도록 하겠다.

14 feature selection methods & 12 classification methods in terms of predictive performance and stability.


Methods

❗ Radiomic Features

A total of 440 radiomic features were used and divided into 4 feature groups.

1) tumor intensity

    - intensity of histogram

2) shape

    - 3D geometric properties of the tumor

3) texture

    - GLCM: gray level co-occurrence matrices

    - GLRLM: gray level run length matrices

     ⇨ quantified the intra-tumor heterogeneity

4) wavelet features

    - transformed domain representations of the intensity and textural features.


❗ Datasets

    • survival time > 2 years ⇨ 1

    • survival time < 2 years ⇨ 0

   - 310 lung cancer patients in training cohort, and 154 patients in validation cohort.

   - All features were normalised using Z-score normalisation.


❗ Feature Selection Methods

 - 14 feature selection methods based on filter approaches were used.

 - 선정기준: simplicity, computational efficiency, popularity in literature

  • Fisher score
  • Relief
  • T-score
  • Chi-square
  • Wilcoxon
  • Gini index
  • Mutual information maximisation
  • Mutual information feature selection
  • Minimum redundancy maximum relevance
  • Conditional informax feature extraction
  • Joint mutual information
  • Conditional mutual information maximisation
  • Interaction capping
  • Double input symmetric relevance


❗ Classifiers

 - 12 machine learning based classification methods were considered.

 - supervised learning task로 training set, validation set으로 나눔

 - 10 fold cross validation was used

 - predictive performance evaluation: AUC

  • Bagging
  • Bayesian
  • Boosting
  • Decision trees
  • Discriminant analysis
  • Generalised linear models
  • Multiple adaptive regression splines
  • Nearest neighbours
  • Neural networks
  • Partial least square and principle component regression
  • Random forests
  • Support vector machines




Analysis

Predictive Performance of Feature Selection Methods

    - feature의 개수를 (n = 5, 10, 15, 20, ..., 50) 점차 늘려가며 AUC 값들의 중앙값 계산


Results

  a total of 440 radiomic features were extracted from the segmented tumor regions

Predictive performance of feature selection and classification methods

• AUC was used for assessing predictive performance of different feature selection and classification methods.

✅ Classification

   👍 Random Forest showed the highest predictive performance as a classifier.

      (AUC = 0.66 ± 0.03)

   👎 Decision Tree had the lowest predictive performance.

      (AUC = 0.54 ± 0.04)

✅ Feature selection

   👍 Wilcoxon test based methods showed the highest predictive performance

      (AUC = 0.65 ± 0.02)

   👎 Chi-square & Conditional informax feature extraction displayed the lowest predictive performance. (AUC = 0.60 ± 0.03)

Stability of the feature selection and classification methods

✅ Feature selection

   👍 Mutual Information Maximisation was the most stable (stability = 0.94 ± 0.02)

   👍 Relief was the second best (stability = 0.91 ± 0.05)

   👎 GINI(GINI index), JMI(Joint mutual information), CHSQ(Chi-square), DISR(Double input symmetric relevance), CIFE(Conditional informax feature extraction) showed relatively low stability.

✅ Classification

    - RSD(Relative standard deviation) were used for measuring empirical stability.

   👍 Bayesian classifier was the best (RSD = 0.86%)

   👍 Generalised linear models was the second best (RSD = 2.19%)

   👍 Partial least square and principle component regression was the third best (RSD = 2.24%)

   👎 Boosting had the lowest stability among the classification methods.


Stability and Predictive Performance



✅ 👍 Feature selection methods

 Wilcoxon (stability = 0.84 ± 0.05, AUC = 0.65 ± 0.02)

 Mutual information feature selection (stability = 0.8 ± 0.03, AUC = 0.63 ± 0.03)

 Minimum redundancy maximum relevance (stability = 0.74 ± 0.03, AUC = 0.63 ± 0.03)

 Fisher score (stability = 0.78 ± 0.08, AUC = 0.62 ± 0.04)

are preferred as their stability and predictive performance was higher than corresponding median values(stability=0.735, AUC=0.615) across all feature selection methods.

✅ 👍 Classification methods

 RF (RSD = 3.52%, AUC = 0.66 ± 0.03)

 BY (RSD = 0.86%, AUC = 0.64 ± 0.05)

 BAG (RSD = 5.56%, AUC = 0.64 ± 0.03)

 GLM (RSD = 2.19%, AUC = 0.63 ± 0.02)

 PLSR (RSD = 2.24%, AUC = 0.63 ± 0.02)

showed that the stability and predictive performance was higher than the corresponding median values(RSD = 5.93%, AUC = 0.61).


Experimental Factors Affecting the Radiomics Based Survival Prediction

 - 3 experimental factors (feature selection methods, classification methods, and the number of selected features) 의 effect를 quantify하기 위해 AUC score에 대한 ANOVA 실시

 - ANOVA result: all 3 factors and their interactions are significant.

 - Classification method was the most dominant source of variability (34.21%)

 - Feature selection accounted for 6.25%

 - Classification X Feature selection interaction explained 23.03%

 - Size of the selected feature subset only shared 1.65% of the total variance



Discussion

 Feature selection methods는 크게 3 카테고리로 나눌 수 있음

(1) filter methods 

     - This paper only investigated filter methods as these are classifier independent.

    • simple feature ranking methods based on some heuristic scoring criterion

    • computationally efficient

    • high generalisability and scalability

(2) wrapper methods

    • classifier dependent

    ⇨ may produce feature subsets that are overly specific to the classifiers, hence low            generalisability

    • search through the whole feature space and identify a relevant and non-redundant          feature subset.

    • computationally expensive 

(3) embedded methods

    • classifier dependent

    ⇨ lacks in the generalisability

    • incorporate feature selection as a part of training process 

    • computationally efficient as compared to the wrappers. 


Filter Methods

 - J : scoring criterion (relevance index)

 - Y : class labels

 - X : set of all features

 - Xk : the feature to be evaluated

 - S : the set of already selected features


위 내용을 작성할 때 Parmar, C., Grossmann, P., Bussink, J. et al. Machine Learning methods for Quantitative Radiomic Biomarkers. Sci Rep 5, 13087 (2015). 해당 논문을 참고하였음.


댓글

이 블로그의 인기 게시물

민감도와 특이도 - 예측력 측도

민감도와 특이도   Predicted Total Pos Neg Actual Pos a b a+b Neg c d c+d 민감도 sensitivity 는 조건부확률로 실제 양성일 때, 양성으로 판정할 확률이다. $$ Sensitivity = P(Y=1|X=1) , Specificity = P(Y=0|X=0) $$ 이를 위 표에 적용해보면 민감도는 $$ \frac{a}{a+b} $$ 특이도 specificity 는 실제 음성일 때 음성으로 판정할 확률이다. 위 표로부터 구한 특이도는 다음과 같다. $$ \frac{c}{c+d} $$ 따라서 민감도와 특이도가 높을수록 예측력이 좋다고 할 수 있다.

통계 : Dummy Variable Trap

 Regression model을 만들다 보면, 독립변수로 명목형 변수를 사용할 때가 많다. 지역이나 성별, biomaker유/무 등이 대표적으로 많이 사용되는 명목형 변수로, 특히 medical 저널에서는 성별을 covariate으로 취급하여 성별에 따른 종속변수의 차이를 보고자 할 때가 많다. 명목형 변수들은 Dummy Variable로 바꾸어서 regression model을 만드는데, one hot 인코딩 방식으로 더미 변수들을 만든다. 예를 들면,  ❕ male=0, female=1 ❕ Biomarker유=1, Biomarker무=0 이런식으로 만든다. 이와 같은 binary 데이터들은 더미 변수로 만드는 것도 쉽고 큰 문제가 되지 않는다. one-hot encoding의 이유는, 0과 1 대신 빨간색, 녹색, 파란색을 1,2,3 으로 코딩해버리면, 적합한 모델은 숫자가 더 큰 3을 빨간색보다 더 중요한 인자로 생각하게 되기 때문이다. 학력처럼 순서형이면 각 숫자에 의미가 있지만 여러 컬러처럼 단순 명목형일 때에는 one-hot 인코딩을 해야한다.  (단, 컬러가 연함 - 진함과 같이 순서형이라면 굳이 one-hot 인코딩을 안해도 된다.) one-hot encoding 의 예시는 아래 세 가지 카테고리로 보면 더 명확하게 이해가 된다. (image출처: https://towardsdatascience.com/encoding-categorical-variables-one-hot-vs-dummy-encoding-6d5b9c46e2db) 그럼 binary 말고 여러 카테고리가 있는 변수의 경우는 어떤 식으로 더미 변수를 만들까? Regression 모델을 만들 때에는 k개의 카테고리가 있다면, k-1개의 더미변수를 생성하게 된다. 위 그림에서 Blue일 때, d1, d2, d3 모두 0으로 코딩해도 Red, Green과 차이가 있으므로, d1과 d2만 만들어도 무방하다는 의미이다. 즉, d1과 d2가 0일 때, Blue가...

GLM과 이탈도 - 모형검정 & 적합도검정

 GLM에서 이탈도 이탈도 Deviance란 무엇일까? GLM에서 이탈도 deviance는 선형회귀모형에서 잔차제곱합에 대응하는 값으로 표본과 적합된 모형 간 편차를 정량화 하는 값이다. GLM에서 이탈도로 할 수 있는 검정이 무엇이 있을까? 1. 모형의 유의성 검정 2. 모형의 적합도평가 GLM 모형을 적합하면 다음과 같은 결과가 나오는데, 이탈도 값은 2개가 나온다. ①Null deviance ②Residual deviance 이 두 개의 이탈도 값이다. ①Null deviance는 절편만 있는 모형에서의 이탈도 ②Residual deviance는 β  가 포함된 모형에서의 이탈도를 의미한다. 먼저 ① 모형의 유의성 검정 을 살펴보겠다. 유의성 검정에서의 귀무가설과 대립가설은 다음과 같다. 귀무가설:  β1  =  β2  = 0 (Null model) 대립가설: 적어도 하나는 0이 아니다. 이를 검정하기 위한 검정통계량 LR(=가능도비 검정통계량) 은 다음과 같다. LR = [Null deviance] - [Residual deviance] = 2(현재모형 하에서 로그가능도 함수의 최댓값 - 귀무가설 하에서 로그가능도 함수의 최댓값) 따라서 이를 위 적합 결과에 적용하면 $$ LR = 225.76 - 189.12 = 36.64 $$ $$ 1 - pchisq(36.64, df=2) = 1.1(10)^{-8} $$ 결론: 유의확률이 매우 작으므로 귀무가설을 기각하고 대립가설을 선택할 수 있다. 모형의 유의성 검정에 대해 살펴보았으니 이번에는 두 번째인 ② 모형의 적합도 검정 을 살펴보자. 적합도 검정에서의 귀무가설과 대립가설은 다음과 같다. 귀무가설: 현재의 단순한 모형 M (=모형M에 포함되지 않은 모수들은 모두 0이다.) 대립가설: 포화모형 saturated model 이를 검정하기 위한 검정통계량은 이탈도 통계량으로 다음과 같이 정의한다. 이탈도 통계량 = 2[포화모형 로그가능도함수의 최댓값 ...