[发明专利]递归生成大阶数行限制覆盖阵列的药物临床试验设计方法

权利要求书

1.一种递归生成大阶数行限制覆盖阵列的药物临床试验设计方法，其特征在于，包括以下步骤：

步骤1、将药物临床试验方案转化为行限制覆盖阵列设计问题；

步骤2、构造行数较少的行限制覆盖阵列和一般覆盖阵列，利用递归方法，生成大阶数行限制覆盖阵列；

步骤3、将生成的行限制覆盖阵列转化为实际的药物临床试验方案，

所述步骤1具体包括以下步骤：

步骤1.1、若共有k种药物需要进行组合测试，每种药物有v种剂量，将v分解为v₁、v₂，其中v₁和v₂均为大于等于2的正整数，而实验者只能同时服用w种药物，则将问题转化为构造带有行限制的覆盖阵列CARL(2，k，v：w)，

所述步骤2具体包括以下步骤：

步骤2.1、通过算法生成基础的小阶数行限制覆盖阵列A＝CARL(2，k，v₁：w)，覆盖阵列B＝CA(2，v₁，v₂)，以及行限制覆盖阵列C＝CARL(2，k，v₂：w)，并分别记阵列A，B，C的行数分别为N_A、N_B、N_C；

记矩阵A＝(a_ij)，这里空元素记为“-”，非空元素记矩阵B中的各列分别为β_i，i∈{0，1，…，v₁}，即β_i为N_B维列向量；记矩阵C＝(ci_j)，这里空元素记为“一”，非空元素

步骤2.2、将矩阵A中的每个非空元素ai_j替换为N_B行二元序对空元素“-”替换为N_B行的空元素“-”，记新得的矩阵为D，矩阵D共有N_AN_B行，所有非空元素

步骤2.3、将矩阵C复制v₁份竖放，并将其中第s份中每个非空元素ci_j替换为二元序对(s，ci_j)，空元素“-”替换为N_C行的空元素“-”，这里s＝0，…，v₁-1，记得到的新矩阵为E，共有v₁N_C行；

步骤2.4、将矩阵D，E竖放构成新矩阵，此矩阵是一个行限制覆盖阵列CARL(2，k，v₁v₂：w)，行数为N_AN_B+v₁N_C行。

2.根据权利要求1所述的递归生成大阶数行限制覆盖阵列的药物临床试验设计方法，其特征在于，步骤2.1中，所述算法包括贪心法、退火法和进化算法。

3.根据权利要求1所述的递归生成大阶数行限制覆盖阵列的药物临床试验设计方法，其特征在于，步骤2.2中，所述替换是将A中元素0，1，…，v₁-1及空元素“-”分别代换：

0→(0，β₁)，共N_B行

1→(1，β₂)，共N_B行

…

共N_B行

“-”→N_B行“-”；

进而得到新矩阵D。