星点设计效应面优化法

1、药学实验设计优化法(3)星点设计-效应面优化法参考书目及文献：1星点设计效应面优化法及其在药学中的应用21中的参考文献3design and analysis of experiments (dc montgomery ed.)4中心多点等距设计法优化醋酸地塞米松聚丙交酯微球的制备工艺.药学学报,1999,34(5):3873915实验设计中多指标的优化：星点设计和总评“归一值”的应用. 中国药学杂志,2000,35(8):530533概述ccd设计表的组成代码值的计算特殊的ccd设计代码值与实际操作物理量的对应关系星点设计概述ccd设计表的组成代码值的计算特殊的ccd设计代码值与实际操作物理

2、量的对应关系星点设计名称辨析：central composite design（ccd）球面设计？中心复合设计？中心组合设计？效应面设计？星点设计！主要特点：多因素五水平，由二水平的析因设计加上中心点和星点组成。五水平：，1，0 x1x2k=2x1k=3x2x3 星点/轴点概述ccd设计表的组成代码值的计算特殊的ccd设计代码值与实际操作物理量的对应关系(1) 析因设计部分：2k或2k1/2（1，1）(2) 星点部分：star point, axial point 为各因素的极值水平(3) 中心点因素数x11+10(,0,0)(0, , ,0)(0,0, ) x1 x2 xk 0 0 0 0

3、0 0 x1 x2 +1 + 1+1 11 + 1 1 1+ 0 0 0 + 0 0 0 析因设计部分f4星点部分次数2k中心点次数5k23455(1/2)f48163216星点数46810101.4141.6822.0002.3782.000中心点数（均一精密）567106中心点数（正交）89121710实验总次数（均一精密）1320315232实验总次数（正交）1623365936概述ccd设计表的组成代码值的计算特殊的ccd设计代码值与实际操作物理量的对应关系x11+10（f）1/4f2k或2k1/2析因设计部分实验次数因素数k=2, =1.414k=3, =1.682k=4, =2.0

4、00 x1 x2 x3 1 +1 + 1 +12 +1 +1 13 +1 1 +14 +1 1 15 1 + 1 +16 1 +1 17 1 1 +18 1 1 1 9 +1.682 0010 1.682 0 011 0 +1.682 012 0 1.682 013 0 0 +1.682 14 0 0 1.68215 0 0 0 20 x1x2k=2x1k=3x2x3概述ccd设计表的组成代码值的计算特殊的ccd设计代码值与实际操作物理量的对应关系等距设计equiradical designx1x2k=2=1.414( 2 )k=3?=1.682, equiradical?=1.732( 3

5、)x2k=2x2k=2x1x1box-behnken面中心立方设计面中心立方设计概述ccd设计表的组成代码值的计算特殊的ccd设计代码值与实际操作物理量的对应关系x11+10k=2 1.414 1 0 +1 +1.414k=3 1.682 1 0 +1 +1.682例.三因素ccd设计1.6821011.682x1 0.0100 0.0290 0.0550 0.0810 0.1000 x2 0.0500 0.0606 0.0750 0.0894 0.1000 x3 0.0500 0.1134 0.2000 0.2866 0.3500中心值(0.1000.010)/2+0.010+1(1.682

6、) 0.010+1.682(1.682) 0.1000.010 x1 x2 x3 x1 x2 x3 loading l e md span p1 t85 od 1 1 1 1 0.0290 0.0606 0.1134 6.95 61.23 26.63 1.37 4.53 132.36 0.48252 1 11 0.0810 0.0606 0.1134 3.06 26.96 12.42 1.74 15.57 25.10 0.20763 1 11 0.0290 0.0894 0.1134 7.64 67.37 78.12 0.97 2.09 156.61 0.44304 +1 11 0.0810

7、0.0894 0.1134 3.72 32.76 16.91 2.19 9.61 50.96 0.26555 1 1 1 0.0290 0.0606 0.2866 20.24 70.61 51.26 1.62 2.72 113.35 0.49386 1 1 1 0.0810 0.0606 0.2866 12.81 44.68 12.75 1.83 6.54 61.58 0.43887 1 1 1 0.0290 0.0894 0.2866 23.41 81.66 76.15 1.17 2.44 164.01 0.51938 1 1 1 0.0810 0.0894 0.2866 14.42 50.

8、30 19.78 1.26 3.27 127.83 0.5118 9 1.732 00 0.0100 0.0750 0.2000 16.90 84.50 99.31 0.94 2.12 287.42 010 1.732 0 0 0.1000 0.0750 0.2000 4.54 22.72 11.40 1.00 12.01 43.06 0 11 0 1.732 0 0.0550 0.0500 0.2000 9.18 45.89 14.37 1.64 10.58 37.48 0.428412 0 1.732 0 0.0550 0.1000 0.2000 11.36 56.80 28.60 1.7

9、0 3.35 145.66 0.476113 0 0 1.732 0.0550 0.0750 0.0500 1.75 34.92 21.55 2.58 7.45 47.98 014 0 0 1.732 0.0550 0.0750 0.3500 23.45 67.01 30.23 2.02 2.94 131.69 0.505415 0 0 0 0.0550 0.0750 0.2000 10.89 54.46 24.66 1.98 23.17 63.29 0.410620均值od值的计算overall desirability 总评“归一值”？步骤：（1）各效应分别求“归一值”（desirability)hassan法：对于越小越好的效应：dmin=(ymaxyi)/(ymaxymin)对于越大越好的效应： dmax=(yiymin)/(ymaxymin)释放度数据p1、t85用harrington方法 p1 = 12 时,d = 0. 75 ;p1 = 4 时,d = 0. 25 t85 = 24 时,d = 0. 8 ;t85 = 60 时,d = 0. 2将d 在转化为无量纲的指标: y= - ln ( -lnd) , y= b0 + b1 y（2）各“归一值”求几何平均值od=(d1,d2,dn)1/nrsm与均匀设计