筆者前面介紹了簡單隨機化分組方法在臨床試驗隨機化分組中的應(yīng)用，同時提出了“隨機抽樣法”來解決簡單隨機分組中可能存在的組間例數(shù)不平衡問起，一并給出了具體的SAS實現(xiàn)代碼。接下來，筆者將介紹區(qū)組隨機化分組在臨床試驗分組中的應(yīng)用及其SAS實現(xiàn)。

針對簡單隨機化分組可能出現(xiàn)的組間例數(shù)不平衡問題，統(tǒng)計學(xué)家對簡單隨機化分組做了下“改良”--既然簡單隨機化可能得到有些組受試者例數(shù)較多而有些組受試者例數(shù)較少，那么是不是可以人為指定每個組有相同的受試者例數(shù)，只需要保證在分配這個過程保持隨機就好。比方說，現(xiàn)在有30個受試者需要分到A、B、C三組中，我們先按照受試者的年齡從小到大排序并分成5組，針對每組中的6位受試者，使用簡單隨機方法，向A、B、C三組分別分配2例即可，這樣就保證了每個分組的例數(shù)都為10。

區(qū)組隨機化

定義：將隨機加以約束，使各處理組的分配更加平衡，滿足研究要求。在一個區(qū)間內(nèi)包含一個預(yù)定的處理分組數(shù)目和比例。

區(qū)組：由若干特征相似的試驗對象組成。如同一窩的動物、批號相同的試劑、體重相近的受試者等。臨床試驗中區(qū)組隨機化常將受試者進(jìn)入研究的時間序列作為劃分區(qū)組的依據(jù)。

區(qū)組的長度：區(qū)組中對象的數(shù)目。

優(yōu)點：區(qū)組隨機化分組，避免簡單隨機化分組可能產(chǎn)生的不平衡現(xiàn)象，不僅提高統(tǒng)計學(xué)效率，而且保證分配率不存在時間趨勢。

使用條件：區(qū)組的長度不宜太小，太小則形成不隨機。一般區(qū)組的長度至少要求為組數(shù)的2倍以上。區(qū)組的長度也不宜太大，太大易使分段內(nèi)不均衡，如果只有兩個組別（試驗組和對照組），區(qū)組的長度一般可取4-8，如果有4個組別則區(qū)組的長度至少為8。

需要說明的是，區(qū)組長度還與試驗的療程長短有關(guān)--對于療程較短的疾病，患者入組快，結(jié)束快，區(qū)組長短影響不大，而對于療程比較長的疾病，區(qū)組長度不宜過大。

一個栗子

一個臨床試驗有160名受試者，現(xiàn)在需要將其使用區(qū)組隨機化方法分配到A、B兩個試驗組。先按照受試者進(jìn)入研究的時間排序，同時設(shè)置區(qū)組長度為8（即每個區(qū)組內(nèi)包含8名受試者），160例受試者劃分成20個區(qū)組。對劃分得到的每個區(qū)組，利用簡單隨機化方法即可得到分組結(jié)果--先給每個區(qū)組內(nèi)的受試者編號（1-8號），然后在從隨機數(shù)字表中隨機地從某個位置開始連續(xù)取8個隨機數(shù)給每個區(qū)組內(nèi)的每一個受試者，并根據(jù)大小將每個隨機數(shù)編上序號，規(guī)定每個區(qū)組內(nèi)隨機數(shù)序號在1-4范圍內(nèi)的受試者進(jìn)入A組（實驗組），隨機數(shù)序號在5-8范圍內(nèi)的受試者進(jìn)入B組（對照組）。比如，針對第一個區(qū)組，我們得到其分組結(jié)果如下表：

區(qū)組1 的隨機化分組結(jié)果

如此這般，重復(fù)以上過程，得到區(qū)組2至區(qū)組20的剩余19個區(qū)組的分組結(jié)果。

區(qū)組隨機化分組的SAS實現(xiàn)

*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-

                                                                      區(qū)組隨機分組

*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-;

*--init--;

%let block_num_subjects=160; *受試者的總例數(shù);

%let block_num_groups=2; *需要分組的組數(shù);

%let block_split_no=20; *劃分的區(qū)組數(shù)目;

%let block_length=8; *區(qū)組長度，即每個區(qū)組內(nèi)的受試者例數(shù)，通常為分組數(shù)目的2-4倍;

%let block_seed=20170709; *隨機種子，通常取值為當(dāng)前日期點;

proc datasets lib=work kill noprint;run;quit;

*--step0-劃分區(qū)組，并得到受試者的隨機數(shù)序號--;

proc plan seed=&block_seed.;

    factors block=&block_split_no.

    length=&block_length. ;

    output out=block_export00;

quit;

*--step1-對各區(qū)組，根據(jù)隨機數(shù)序號，將受試者分配到目標(biāo)分組中--;

data block_export01;

    id=_N_;

    set block_export00;

    group=ceil(length / (&block_length./&block_num_groups.));

run;

proc print noobs;run;