%==========================================================================
%
% 函數名：choose
% 函數介紹：根據設定的顯著性水平，選出特征基因。
% 輸入參數：data是基因表達譜數據的原始表達量。數據的第一行是樣本標記行。0表示一類樣本，1表示另一類樣本。
%          p_value是基因的犯錯概率向量，是單邊檢驗值。1*表達譜中的基因個數。
% 輸出參數：
%          feature_gene是特征基因的原始表達量。
%          feature_order是特征基因的順序號。
% 閾值設定：顯著性水平alfa=0.05
%==========================================================================
function [ feature_gene, feature_order ] = choose ( data, p_value )

[ row, column ] = size ( data );
%去掉樣本標記行
gene_data = data ( 2 : row, :);
m = row - 1;

alfa = input ('請輸入顯著性水平值：') ;%顯著性水平

p = 1;
for i = 1: m
    if p_value ( i ) <= alfa             
        feature_gene ( p, :) = gene_data( i, : );
        feature_order ( p ) = i;
         p = p+1;
    end
end

if p == 1%說明顯著性水平alfa太小，沒有特征基因。
    feature_gene = nan;
    feature_order = nan;
end

%end