利用SPSS的正交设计

正交实验设计Orthogonalexperimentaldesign（一）概念1.定义是利用一套规格化的“正交表”，将实验因素、各水平之间的组合均匀搭配、合理安排，用较少的、有代表性的处理组合数，提供充分有用信息的一种高效、快速的多因素实验设计方法。2.因素和水平因素factor：实验结果或观察指标的影响因素。水平level：各因素的不同状态。3.主效应和交互作用主效应maineffect—每个因素对观察指标的作用。交互作用interaction—一些因素取不同水平时对另一些因素的作用。两个因素间的交互作用，称为一级交互作用，记为A×B。三个因素间的交互作用，称为二级交互作用，记为A×B×C。三个因素以上的交互作用统称为高级交互作用。经验证明，二级交互作用大多可以忽略，高级交互作用则全部可以忽略。Ln(mk)（二）正交表的表示方法正交表实验次数水平数因素个数L4(23)是指最多可以安排3个两水平的因素作4次实验的正交表。实验次数因素及其水平ABC第一次第二次第三次第四次A1A1A2A2B1B2B1B2C1C2C2C1L8(27)正交表实验次数因素及其水平ABC第一次第二次第三次第四次第五次第六次第七次第八次A1A1A1A1A2A2A2A2B1B1B2B2B1B1B2B2C1C2C1C2C1C2C1C2正交表中的每一列可有三种安排：（1）安排某个因素（2）安排某种交互作用（3）空列（做误差用）交互作用表：每张正交表均有对应的交互作用表。L4(23)正交表列号列号12312321L8(27)正交表的交互作用表列号列号123456712321567476174523654321（三）表头设计根据分析需要，选用合适的正交表，将各个因素安排在正交表的各列的过程，称为表头设计。表头设计原则：根据因素个数、水平数合理安排主效应和因素之间的交互作用。必须考虑对正交试验数据作方差分析时误差的计算途径。例：研究A、B、C三个因素的主效应，三个因素均有2个水平。同时要研究交互作用A×B和A×C。根据研究目的确定主要因素。至少正交表要有5列。故选L8(27)列号1234567因素ABA×BCA×C6、7列是空列，作误差估计用。1.根据研究目的确定主要因素。2.确定每个因素的水平数。3.确定各因素间的交互作用。4.确定实验次数。5.综合上述情况作表头设计。6.按表头设计要求进行试验，收集数据。7.对数据作方差分析。（四）正交实验及分析的基本步骤交互作用的确定凭经验来确定。借助于表头设计及结果计算：不清楚交互作用时，尽量空出有交互作用的列号，以防发生意外混杂。举例有研究者以大鼠作实验，观察指标是细胞色素b5（nmol/mg），主要想了解正氟醚的作用，同时要看一下用生理盐水和用戊巴比妥作为诱导药物对正氟醚的毒性作用有何影响，以及不同诱导剂对不同性别大鼠作用有何不同。试作正交试验设计及分析。1.确定因素和水平数因素和水平数因素水平12ABC生理盐水不用正氟醚雄性大鼠戊巴比妥用正氟醚雌性大鼠3.表头设计：查1937页附表18-1，用L8(27)。2.交互作用：A×B、A×C列号1234567因素ABA×BCA×C试验号A1B2A×B3C4A×C567Cytb12345678111122221122112211222211121212121212212112211221122121120.660.230.110.130.640.950.150.12L8(27)正交表及试验结果4.统计分析：方差分析用SPSS16.0计算。步骤：AnalyzeGeneralLinealModelunivariate练习：研究高频呼吸机A、B、C、D四个参数对通气量的影响，每个参数有高低两个水平，试进行正交设计，并列出试验安排及结果记录表。要求研究A、B、C、D、A×B、B×D列号1234567因素ABA×BCDB×D如果考虑交互作用B×C、B×D列号1234567因素BCB×CDB×DA对正交实验结果的解释例1：有研究者以大鼠作实验，观察指标是细胞色素b5（nmol/mg），主要想了解正氟醚的作用，同时要看一下用生理盐水和用戊巴比妥作为诱导药物对正氟醚的毒性作用有何影响，以及不同诱导剂对不同性别大鼠作用有何不同。试作正交试验设计及分析。试验号A1B2A×B3C4A×C567Cytb12345678111122221122112211222211121212121212212112211221122121120.660.230.110.130.640.950.150.12L8(27)正交表及试验结果试验号A1B2A×B3C4A×C5Cytb1234567811112222112211221122221112121212121221210.660.230.110.130.640.950.150.1211.342.631.101.472.1221.950.662.191.821.17方差分析结果P值均大于0.05，各因素对大鼠细胞色素b5均无明显的影响。有重复两水平正交实验例：在正氟醚对大鼠肝重影响的研究中，A、B、C的意义同前，但观察指标为肝重和体重之比，需要分析交互作用A×B、A×C、B×C，不考虑二级交互作用。每重试验重复3次。试做正交设计及方差分析。试验号A1B2A×B3C4A×C5B×C612345678111122221122112211222211121212121212212112211221L8（27）正交表三次重复实验的肝重/体重（%）实验次数X1X2X3∑X123456785.005.266.305.426.135.876.024.655.525.687.025.606.465.505.704.605.385.835.905.705.216.205.485.4415.9016.7719.2216.7217.8017.5717.2014.68水平A1B2A×B3C4A×C5B×C6ⅠⅡ68.6167.2568.0467.8264.5571.3170.1265.7467.3768.4965.1070.76每一因素各水平的肝重/体重（%）结果分析用SPSS11.0进行方差分析。TestsofBetween-SubjectsEffectsDependentVariable:肝体比4.170a6.6954.322.008769.0811769.0814782.841.0007.707E-0217.707E-02.479.4982.017E-0312.017E-03.013.9121.90411.90411.841.003.7991.7994.971.0405.227E-0215.227E-02.325.5761.33511.3358.301.0102.73417.161775.984246.90323SourceCorrectedModelInterceptABABCACBCErrorTotalCorrectedTotalTypeIIISumofSquaresdfMeanSquareFSig.RSquared=.604(AdjustedRSquared=.464)a.三水平正交设计最简单的三水平正交设计：L9（34）任意两列间的交互作用出现于其他两列。试验号列号1234123456789111222333123123123123231312123312231L9（34）列号12345678(1)3424(2)142313(3)1267811913(4)1012579121011813(5)17L27（313）两列间的交互作用三水平正交设计例：在胃蛋白酶生产过程中，产品质量指标是残留蛋白，影响因素是水解温度、水解时间、加烟酸量、烘室温度，每个因素均取三水平。可忽略交互作用。因素：水解温度℃：A1=43A2=46A3=49水解时间h：B1=3.5B2=4.0B3=4.5加盐酸量%：C1=2.0C2=2.6C3=3.2烘室温度℃：D1=55D2=60D3=65试验号A1B2C3D4残留蛋白x1x2x3∑x1234567891112223331231231231232313121233122311.50.70.40.81.10.90.70.33.21.30.50.60.61.20.90.70.33.01.40.20.30.71.00.80.80.43.14.21.41.22.13.32.62.21.09.3L9（34）水平A1B2C3D4ⅠⅡⅢ6.88.012.58.55.713.17.812.86.716.86.24.3每一因素各水平的残留蛋白（mol）TestsofBetween-SubjectsEffectsDependentVariable:PROTEIN17.540a82.193159.993.00027.603127.6032014.297.0002.00721.00373.216.0003.10221.551113.189.0002.34921.17485.703.00010.08225.041367.865.000.24718.01445.3902717.78726SourceCorrectedModelInterceptABCDErrorTotalCorrectedTotalTypeIIISumofSquaresdfMeanSquareFSig.RSquared=.986(AdjustedRSquared=.980)a.方差分析结果：4个主效应都有高度统计学意义最佳配方的确定：把其作用有统计意义的因素用合适的方法组合起来。上例的最佳配方：对于胃蛋白酶的质量，应使残留蛋白越少越好，所以最佳配方是A1B2C3D3如果有些因素并无统计学意义，那么在组成最佳配方时就可以不考虑它们，或者任意取什么水平都可以。在两水平正交表上安排四水平的因素在四水平的正交表中只能安排四水平的因素，如果在研究工作中有些因素是两水平的，有些因素是四水平的，则在同一张正交表同时安排两水平和四水平的因素，其方法是把两水平正交表改为可安排四水平的正交表。如在L8（27）正交表中安排一个四水平因素和3个2水平因素。把L8（27）改为L8（4×24）。方法如下：把1，2列的水平数看成有序的对子，（1,1、1,2、2,1、2,2）。把1，2列合并成1列，并按下列规则改写水平数：1，1→11，2→22，1→32，2→4同时把1，2列的交互作用列3划掉。L8(27)正交表试验号12345671234567811112222112211221122221112121212121221211221122112212112L8(4×24)正交表试验号1,2,34567123456781122334412121212121221211221122112212112例：在蛇毒抑瘤实验中，有三个试验因素：蛇毒成分A、瘤株B、剂量C。A有四个水平，B、C各有两个水平：A1=Ⅻ峰、A2=Ⅰ峰、A3=Ⅷ峰、A4=原毒，B1=小鼠肉瘤(S180)、B2=小鼠艾氏腹水癌,C1=剂量0（对照）、C2=剂量0.075mg/kg。不考虑交互作用。观察指标为瘤重。试做正交试验，找出最优组合条件。每次试验用5头小鼠做重复试验。1.设计L8(4×24)正交表试验号1A2B3C瘤重（g）12345计123456781122334412121212121221211.50.61.70.70.61.00.41.02.00.60.90.90.71.10.60.91.60.81.40.80.50.60.91.12.00.51.61.00.70.90.61.81.70.51.70.71.00.70.51.58.83.07.34.13.54.33.06.3A1A2A3A4B1B2C1C21.181.140.780.931.130.881.340.68每一因素各水平的平均瘤重方差分析结果，A间、B间、C间均为P＜0.01蛇毒抑瘤作用以成分Ⅷ、瘤株艾氏腹水瘤、剂量0.075mg/kg的组合条件为最好。练习题1：长梗胶股蓝总皂甙提取工艺研究。主要研