[发明专利]一种基于行限制覆盖阵列的药物临床试验设计方法

权利要求书

1.一种基于行限制覆盖阵列的药物临床试验设计方法，其特征在于，包含以下步骤：

步骤1、将临床试验方案转化为行限制覆盖阵列设计问题，若共有k种药物需要进行组合测试，每种药物有v种剂量，实验者只能同时服用w种药物，且满足则将问题转化为构造带有行限制的覆盖阵列CARL(t，k，v：w)，其中，t为覆盖强度，取值为正整数；

步骤2、构造满足参数要求的带有行限制的覆盖阵列，生成行限制覆盖阵列，具体步骤如下：

步骤2.1、生成初始覆盖阵列A＝CA(t，w，v)，w列，N₁行；

步骤2.2、在初始覆盖阵列的基础上，使用矩阵的填充方法扩充阵列，构造满足参数要求的带有行限制的覆盖阵列B＝CARL(t，k，v：w)；

步骤2.3、进一步降阶得到新的带有行限制的覆盖阵列C＝CARL(t-1′k-1，v：w-1)；

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

所述步骤2.2的具体步骤如下：

步骤2.2.1、基于初始覆盖阵列A，生成(t+1)N₁×k的空矩阵B，从第一行开始，每N₁行记为一组，共t+1组；

步骤2.2.2、记s＝k-w，依次将如下位置记为空值“-”：第1组的1～s列，第2组的s+1～2s列，…，第t+1组的ts+1～(t+1)s列，由于

有k≥(t+1)s，保证如上操作的可行性；

步骤2.2.3、填充矩阵B，在B的每一组非空值“-”的位置依次有序将覆盖阵列A的各列循环填入，得到的B为CARL(t，k，v：w)，N₂行，N₂＝(t+1)N₁。

2.根据权利要求1所述的基于行限制覆盖阵列的药物临床试验设计方法，其特征在于，步骤2.1中，所述生成初始覆盖阵列的方法包括组合设计数学构造、贪心法和进化算法。

3.根据权利要求1所述的基于行限制覆盖阵列的药物临床试验设计方法，其特征在于，所述步骤2.3的具体步骤如下：

步骤2.3.1、在矩阵B中任取一列，在其中任取非空值“-”的元素，设为x；

步骤2.3.2、提取B中所有在此列上取值为x的行，并删除该列，生成新矩阵C，C是一个降阶的带有行限制的覆盖阵列CARL(t-1，k-1，v：w-1)，N₃行。

4.根据权利要求3所述的基于行限制覆盖阵列的药物临床试验设计方法，其特征在于，步骤2.3.2中，选取适当的元素x，确保所述矩阵C的行数满足N₃≤tN₁/v。