溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊(cè)×
其他方式登錄
點(diǎn)擊 登錄注冊(cè) 即表示同意《億速云用戶服務(wù)條款》

如何用R語言實(shí)現(xiàn)RCT分層區(qū)組隨機(jī)化

發(fā)布時(shí)間:2021-11-24 14:53:22 來源:億速云 閱讀:1325 作者:柒染 欄目:大數(shù)據(jù)

如何用R語言實(shí)現(xiàn)RCT分層區(qū)組隨機(jī)化,很多新手對(duì)此不是很清楚,為了幫助大家解決這個(gè)難題,下面小編將為大家詳細(xì)講解,有這方面需求的人可以來學(xué)習(xí)下,希望你能有所收獲。

劃重點(diǎn):分層區(qū)組隨機(jī)化,一般多中心的RCT常用,以中心分層后,各層中再區(qū)組從而保證各組受試者數(shù)量均衡。

調(diào)用的package是blockrand,包的內(nèi)容,就只有這一個(gè)功能,就是做分層區(qū)組隨機(jī)化,就喜歡這種直白的包,沒有眼花繚亂的東西。

簡(jiǎn)單說明一下,為了演示的直觀,編了一個(gè)虛擬場(chǎng)景。假設(shè)有個(gè)方案的Sample-size是32,Two-arm Designs,兩個(gè)中心(center A & B),區(qū)組長(zhǎng)度設(shè)為4,隨機(jī)化代碼如下:

輸入:

library(blockrand)C_A <- blockrand(n=20, # 該層患者數(shù)                  id.prefix='A', # 層的編號(hào)                  block.prefix='A',  # 區(qū)組號(hào)前綴                  stratum='center_A', # 層的label                  block.sizes = 2, # 固定區(qū)組長(zhǎng)度                   num.levels = 2 # 實(shí)驗(yàn)分組數(shù)量                  )C_B <- blockrand(n=12,                  id.prefix='B',                  block.prefix='B',                   stratum='center_B',                 block.sizes = 2,                  num.levels = 2                 )result <- rbind(C_A,C_B)print(result)

輸出:

    id  stratum block.id block.size treatment1  A01 center_A       A1          2         B2  A02 center_A       A1          2         A3  A03 center_A       A2          2         A4  A04 center_A       A2          2         B5  A05 center_A       A3          6         A6  A06 center_A       A3          6         B7  A07 center_A       A3          6         B8  A08 center_A       A3          6         A9  A09 center_A       A3          6         B10 A10 center_A       A3          6         A11 A11 center_A       A4          2         B12 A12 center_A       A4          2         A13 A13 center_A       A5          2         B14 A14 center_A       A5          2         A15 A15 center_A       A6          6         B16 A16 center_A       A6          6         A17 A17 center_A       A6          6         B18 A18 center_A       A6          6         A19 A19 center_A       A6          6         B20 A20 center_A       A6          6         A21 B01 center_B       B1          4         A22 B02 center_B       B1          4         B23 B03 center_B       B1          4         A24 B04 center_B       B1          4         B25 B05 center_B       B2          4         A26 B06 center_B       B2          4         B27 B07 center_B       B2          4         A28 B08 center_B       B2          4         B29 B09 center_B       B3          4         B30 B10 center_B       B3          4         A31 B11 center_B       B3          4         B32 B12 center_B       B3          4         A

以上展示的是一個(gè)最基本的例子,關(guān)于隨機(jī)化其實(shí)根據(jù)方案和實(shí)際情況的不同,狀況也是多種多樣的,實(shí)際應(yīng)用靠的還是舉一反三。

看完上述內(nèi)容是否對(duì)您有幫助呢?如果還想對(duì)相關(guān)知識(shí)有進(jìn)一步的了解或閱讀更多相關(guān)文章,請(qǐng)關(guān)注億速云行業(yè)資訊頻道,感謝您對(duì)億速云的支持。

向AI問一下細(xì)節(jié)

免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如果涉及侵權(quán)請(qǐng)聯(lián)系站長(zhǎng)郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。

AI