如何用R语言实现RCT分层区组随机化,很多新手对此不是很清楚,为了帮助大家解决这个难题,下面小编将为大家详细讲解,有这方面需求的人可以来学习下,希望你能有所收获。
划重点:分层区组随机化,一般多中心的RCT常用,以中心分层后,各层中再区组从而保证各组受试者数量均衡。
调用的package是blockrand,包的内容,就只有这一个功能,就是做分层区组随机化,就喜欢这种直白的包,没有眼花缭乱的东西。
简单说明一下,为了演示的直观,编了一个虚拟场景。假设有个方案的Sample-size是32,Two-arm Designs,两个中心(center A & B),区组长度设为4,随机化代码如下:
输入:
library(blockrand)C_A <- blockrand(n=20, # 该层患者数 id.prefix='A', # 层的编号 block.prefix='A', # 区组号前缀 stratum='center_A', # 层的label block.sizes = 2, # 固定区组长度 num.levels = 2 # 实验分组数量 )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
以上展示的是一个最基本的例子,关于随机化其实根据方案和实际情况的不同,状况也是多种多样的,实际应用靠的还是举一反三。
看完上述内容是否对您有帮助呢?如果还想对相关知识有进一步的了解或阅读更多相关文章,请关注亿速云行业资讯频道,感谢您对亿速云的支持。
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。