歷史文章
筆者前面介紹了簡(jiǎn)單隨機(jī)化分組方法在臨床試驗(yàn)隨機(jī)化分組中的應(yīng)用，同時(shí)提出了“隨機(jī)抽樣法”來(lái)解決簡(jiǎn)單隨機(jī)分組中可能存在的組間例數(shù)不平衡問(wèn)起，一并給出了具體的SAS實(shí)現(xiàn)代碼。接下來(lái)，筆者將介紹區(qū)組隨機(jī)化分組在臨床試驗(yàn)分組中的應(yīng)用及其SAS實(shí)現(xiàn)。
針對(duì)簡(jiǎn)單隨機(jī)化分組可能出現(xiàn)的組間例數(shù)不平衡問(wèn)題，統(tǒng)計(jì)學(xué)家對(duì)簡(jiǎn)單隨機(jī)化分組做了下“改良”--既然簡(jiǎn)單隨機(jī)化可能得到有些組受試者例數(shù)較多而有些組受試者例數(shù)較少，那么是不是可以人為指定每個(gè)組有相同的受試者例數(shù)，只需要保證在分配這個(gè)過(guò)程保持隨機(jī)就好。比方說(shuō)，現(xiàn)在有30個(gè)受試者需要分到A、B、C三組中，我們先按照受試者的年齡從小到大排序并分成5組，針對(duì)每組中的6位受試者，使用簡(jiǎn)單隨機(jī)方法，向A、B、C三組分別分配2例即可，這樣就保證了每個(gè)分組的例數(shù)都為10。
區(qū)組隨機(jī)化
定義：將隨機(jī)加以約束，使各處理組的分配更加平衡，滿足研究要求。在一個(gè)區(qū)間內(nèi)包含一個(gè)預(yù)定的處理分組數(shù)目和比例。
區(qū)組：由若干特征相似的試驗(yàn)對(duì)象組成。如同一窩的動(dòng)物、批號(hào)相同的試劑、體重相近的受試者等。臨床試驗(yàn)中區(qū)組隨機(jī)化常將受試者進(jìn)入研究的時(shí)間序列作為劃分區(qū)組的依據(jù)。
區(qū)組的長(zhǎng)度：區(qū)組中對(duì)象的數(shù)目。
優(yōu)點(diǎn)：區(qū)組隨機(jī)化分組，避免簡(jiǎn)單隨機(jī)化分組可能產(chǎn)生的不平衡現(xiàn)象，不僅提高統(tǒng)計(jì)學(xué)效率，而且保證分配率不存在時(shí)間趨勢(shì)。
使用條件：區(qū)組的長(zhǎng)度不宜太小，太小則形成不隨機(jī)。一般區(qū)組的長(zhǎng)度至少要求為組數(shù)的2倍以上。區(qū)組的長(zhǎng)度也不宜太大，太大易使分段內(nèi)不均衡，如果只有兩個(gè)組別（試驗(yàn)組和對(duì)照組），區(qū)組的長(zhǎng)度一般可取4-8，如果有4個(gè)組別則區(qū)組的長(zhǎng)度至少為8。
需要說(shuō)明的是，區(qū)組長(zhǎng)度還與試驗(yàn)的療程長(zhǎng)短有關(guān)--對(duì)于療程較短的疾病，患者入組快，結(jié)束快，區(qū)組長(zhǎng)短影響不大，而對(duì)于療程比較長(zhǎng)的疾病，區(qū)組長(zhǎng)度不宜過(guò)大。
一個(gè)栗子
一個(gè)臨床試驗(yàn)有160名受試者，現(xiàn)在需要將其使用區(qū)組隨機(jī)化方法分配到A、B兩個(gè)試驗(yàn)組。先按照受試者進(jìn)入研究的時(shí)間排序，同時(shí)設(shè)置區(qū)組長(zhǎng)度為8（即每個(gè)區(qū)組內(nèi)包含8名受試者），160例受試者劃分成20個(gè)區(qū)組。對(duì)劃分得到的每個(gè)區(qū)組，利用簡(jiǎn)單隨機(jī)化方法即可得到分組結(jié)果--先給每個(gè)區(qū)組內(nèi)的受試者編號(hào)（1-8號(hào)），然后在從隨機(jī)數(shù)字表中隨機(jī)地從某個(gè)位置開(kāi)始連續(xù)取8個(gè)隨機(jī)數(shù)給每個(gè)區(qū)組內(nèi)的每一個(gè)受試者，并根據(jù)大小將每個(gè)隨機(jī)數(shù)編上序號(hào)，規(guī)定每個(gè)區(qū)組內(nèi)隨機(jī)數(shù)序號(hào)在1-4范圍內(nèi)的受試者進(jìn)入A組（實(shí)驗(yàn)組），隨機(jī)數(shù)序號(hào)在5-8范圍內(nèi)的受試者進(jìn)入B組（對(duì)照組）。比如，針對(duì)第一個(gè)區(qū)組，我們得到其分組結(jié)果如下表：
如此這般，重復(fù)以上過(guò)程，得到區(qū)組2至區(qū)組20的剩余19個(gè)區(qū)組的分組結(jié)果。
區(qū)組隨機(jī)化分組的SAS實(shí)現(xiàn)
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 區(qū)組隨機(jī)分組
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-;
*--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ū)組長(zhǎng)度，即每個(gè)區(qū)組內(nèi)的受試者例數(shù)，通常為分組數(shù)目的2-4倍;
%let block_seed=20170709; *隨機(jī)種子，通常取值為當(dāng)前日期點(diǎn);
proc datasets lib=work kill noprint;run;quit;
*--step0-劃分區(qū)組，并得到受試者的隨機(jī)數(shù)序號(hào)--;
proc plan seed=&block_seed.;
? ? factors block=&block_split_no.
? ? length=&block_length. ;
? ? output out=block_export00;
quit;
*--step1-對(duì)各區(qū)組，根據(jù)隨機(jī)數(shù)序號(hào)，將受試者分配到目標(biāo)分組中--;
data block_export01;
? ? id=_N_;
? ? set block_export00;
? ? group=ceil(length / (&block_length./&block_num_groups.));
run;
proc print noobs;run;
