%==========================================================================

% 函數名：pillar_cii
% 函數介紹： 畫基因分類信息指數_基因個數柱圖。d是一個行向量。
% 輸入參數： d是基因的分類信息指數。
%           r1是第一組樣本的個數。
%           r2是第二組樣本的個數。
% 輸出參數： alfa是柱圖縱軸的突跳點。
%           y是柱圖中每根柱的高，是基因的分布情況＝對應每個橫軸區間的基因個數。

%==========================================================================
function [ alfa, y ] = pillar_cii ( d, r1, r2 )
d = transpose ( d );
[m, n] = size ( d );
A = min ( d );
B = max ( d );
fenshu = round ( 1.87 * (( r1 + r2 - 1 ) ^ 0.4 )) + 1;
t = ( B - A ) / fenshu;
f = A + t / 2 : t : B - t / 2;
x = A : t : B;
%==============================

% 求直方圖的高。
[y] = hight ( x, d, fenshu, m );

%==============================
% 畫柱圖。
bar ( f, y, 'k' );%黑色柱圖
ylabel ('基因個數');
xlabel ('基因分類信息指數');
%==============================

% 確定基因個數突跳點。
[alfa] = sharp_jump ( f, y, fenshu, t );
%==============================
% end