您好,登錄后才能下訂單哦!
如何用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ì)億速云的支持。
免責(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)容。