生物检定统计法(二)
http://www.100md.com
大众医药网
| 标题 | 生物检定统计法(二) |
| 附录序号 | 附录ⅩⅣ |
| 内容全文 | 生物检定统计法 三、量反应平行线测定法 药物对生物体所引起的反应随着药物剂量的增加产生的量变可以测量者,称量反应。量反应检定用平行线测定法,要求在一定剂量范围内,S和T的对数剂量x和反应或反应的特定函数y呈直线关系,当S和T的活性组分基本相同时,两直线平行。 本药典量反应检定主要用(2.2)法、(3.3)法或(2.2.2)法、(3.3.3)法,即S、T(或U)各用2个剂量组或3个剂量组,统称(k.k)法或(k·k·k)法;如果S和T的剂量组数不相等,则称(k. k')法;前面的k代表S的剂量组数,后面的k或k'代表T的剂量组数。一般都是按(k·k)法实验设计,当S或T的端剂量所致的反应未达阈值,或趋于极限,去除此端剂量后,对数剂量和反应的直线关系成立,这就形成了(k·k')法。例如(3.3)法设计就可能形成(2.3)或(3.2)法等。因此,(k·k')法中的k只可能比k'多一组或少一组剂量。(k·k')法的计算结果可供重复试验时调节剂量或调整供试品估计效价时参考。无论是(k·k) 法、(k.k')法或(k.k.k)法,都以K代表S和T的剂量组数之和,故K=k+k、K=k+k'或K=k+k+k。 本药典平行线测定法的计算都用简算法,因此对各种(k·k)法要求: (1) S和T相邻高低剂量组的比值(r)要相等,一般r用1:0.8~1:0.5,logr=I; (2) 各剂量组的反应个数(m)应相等。 1.平行线测定的实验设计类型 根据不同的检定方法可加以限制的因级数采用不同的实验设计类型。本药典主要用下面三种实验设计类型。 (1) 随机设计 剂量组内不加因级限制,有关因子的各级随机分配到各剂量组。本设计类型的实验结果只能分离不同剂量(剂间)所致变异。如绒促性素的生物检定。 (2) 随机区组设计 将实验动物或实验对象分成区组,一个区组可以是一窝动物、一只双碟、 或一次实验.在剂量组内的各行间加以区组间(如窝间、碟间、实验次序间)的因级限制。随机区组设计要求每一区组的容量(如每一窝动物的受试动物只数、每一只双碟能容纳的小杯数等)必须和剂量组数相同,这样可以使每一窝动物或每一只双碟都能接受到各个不同的剂量。因此随机区组设计除了从总变异中分离剂间变异之外,还可以分离区组间变异,减小实验误差。例如抗生素杯碟法效价测定。 (3) 交叉设计 同一动物可以分两次进行实验者适合用交叉设计。交叉设计是将动物分组,每组可以是一只动物,也可以是几只动物,但各组的动物只数应相等。标准品(S)和供试品(T)对比时,一组动物在第一次试验时接受(S)的一个剂量,第二次试验时则接受(T)的一个剂量,如此调换交叉进行,可以在同一动物身上进行不同试品、不同剂量的比较,以去除动物间差异对实验误差的影响,提高实验精确度,节约实验动物。 (2.2)法S和T各两组剂量,用双交叉设计,将动物分成四组;对各组中的每一只动物都标上识别号。每一只动物都按给药次序表进行两次实验。 双交叉设计两次实验的给药次序表
───────┬─────┬─────┬─────┬──────
│ 第一组 │ 第二组 │ 第三组 │ 第四组
───────┼─────┼─────┼─────┼──────
第一次实验 │ds<[1]> │ds<[2]> │dT<[1]> │ dT<[2]>
第二次实验 │dT<[2] │dT<[1]> │ds<[2]> │ ds<[1]>
───────┴─────┴─────┴─────┴──────
2.平行线测定法的方差分析和可靠性测验随机设计和随机区组设计的方差分析和可靠性测验 (1) 将反应值或其规定的函数(y)按S和T的剂量分组列成方阵表 见表二。
表二 剂量分组方阵表
──────┬───────────────────────────┬────
│ (S)和(T)的剂量组 │ 总 和
├─────┬─────┬─────┬───┬─────┤
│ (1) │ (2) │ (3) │ … │ (k) │∑y<[m]>
──────┼─────┼─────┼─────┼───┼─────┼────
行 1 │ y<[1]>(1)│ y<[1]>(2)│ y<[1]>(3)│ … │y<[1]>(k) │∑y<[1]>
间 2 │ y<[2]>(1)│ y<[2]>(2)│ y<[2]>(3)│ … │y<[2]>(k) │∑y<[2]>
组 3 │ y<[3]>(1)│ y<[3]>(2)│ y<[3]>(3)│ … │y<[3]>(k) │∑y<[3]>
内 │ │ │ │ │ │ m
│ y<[m]>(1)│ y<[m]>(2)│ y<[m]>(3)│ … │y<[m]>(k) │∑y<[m]>
──────┼─────┼─────┼─────┼───┼─────┼────
总和∑y(k)│ ∑y(1) │ ∑y(2) │ ∑y(3) │ … │∑y(k) │ ∑y
──────┴─────┴─────┴─────┴───┴─────┴────方阵中,K为S和T的剂量组数和,m为各剂量组内y的个数,如为随机区组设计,m为行间或组内所加的因级限制;n为反应的总个数,n=mK。 (2) 特异反应剔除和缺项补足 特异反应剔除 在同一剂量组内的各个反应中,如出现个别特大或特小的反应,应按下法判断其是否可以剔除。 设y<[a]>示特异反应值(或其规定的函数),y<[m]>为与y<[a]>相对的另一极端的反应值,y<[2]>、y<[3]>为与y<[a]>最接近的两个反应值,y<[m-1]>、y<[m-2]>为与y<[m]>最接近的两个反应值,m是该剂量组内的反应个数,将各数值按大小次序排列如下: y<[a]>、y<[2]>、y<[3]>…y<[m-2]>、y<[m-1]>、y<[m]> 如y<[a]>为特大值,则依次递减,y<[m]>最小;如y<[a]>为特小值则依次递升,y<[m]>最大。按(10)~(12)式计算J值。 当m=3~7时, y<[2]>-y<[a]> J<[1]>=──────── (10) y<[m]>-y<[a]> 当m=8~13时, y<[3]>-y<[a]> J<[2]>=──────── (11) y<[m-1]>-y<[a]> 当m=14~20时, y<[3]>-y<[a]> J<[3]>=──────── (12) y<[m-2]>-y<[a]> 如J的计算值大于J值表(表三)中规定的相应数值时,y<[a]>即可剔除。 表三 剔除特异反应的J值表
┏━━━┯━━━┯━━━┯━━━┯━━━┯━━━┓
┃ m │ 3 │ 4 │ 5 │ 6 │ 7 ┃
┠───┼───┼───┼───┼───┼───┨
┃J<[1]>│ 0.98 │ 0.85 │ 0.73 │ 0.64 │ 0.59 ┃
┠───┼───┼───┼───┼───┼───╄━━━┓
┃m │ 8 │ 9 │ 10 │ 11 │ 12 │ 13 ┃
┠───┼───┼───┼───┼───┼───┼───┨
┃J<[2]>│ 0.78 │ 0.73 │ 0.68 │ 0.64 │ 0.61 │ 0.58 ┃
┠───┼───┼───┼───┼───┼───┼───╄━━━┓
┃m │ 14 │ 15 │ 16 │ 17 │ 18 │ 19 │ 20 ┃
┠───┼───┼───┼───┼───┼───┼───┼───┨
┃J<[3]>│ 0.60 │ 0.58 │ 0.56 │ 0.54 │ 0.53 │ 0.51 │ 0.50 ┃
┗━━━┷━━━┷━━━┷━━━┷━━━┷━━━┷━━━┷━━━┛缺项补足 因反应值被剔除或因故反应值缺失造成缺项,致m不等时,根据实验设计类型做缺项补足,使各剂量组的反应个数m相等。 随机设计 对缺失数据的剂量组,以该组的反应均值补入,缺1个反应补一个均值,缺2个反应补2个均值。 随机区组设计 按(13)式计算,补足缺项。 KC+mR-G 缺项y=────────── (13) (k-2)(m-1) 式中 C为缺项所在剂量组内的反应值总和; R为缺项所在行的反应值总和; G为全部反应值总和。 如果缺1项以上,可以分别以y<[1]>、y<[2]>、y<[3]>……等代表各缺项, 然后在计算其中之一时,把其他缺项y直接用符号y<[1]>、y<[2]>…当作未缺项代入(13)式,这样可得与缺项数相同的方程组,解方程组即得。 随机区组设计当剂量组内安排的区组数较多时,也可将缺项所在的整个区组除去。 随机设计的实验结果中,如在个别剂量组多出1~2个反应值,可按严格的随机原则去除,使各剂量组的反应个数m相等。 不论哪种实验设计,每补足一个缺项,就需把S2的自由度减去1,缺项不得超过反应总个数的5%。 (3) 方差分析 方阵表(表二)的实验结果,按(14)~(21)式计算各项变异的差方和、自由度(f)及误差项的方差(s<2>)。 随机设计 按(14)、(15)式计算差方和(总)、差方和(剂间)。按(20)式计算差方和(误差)。按(18)或(21)式计算s<2>。 随机区组设计 按(14)~(17)式计算差方和(总)、差方和(剂间)、差方和(IX组间)、差方和(误差)。按(18)或(19)式计算s<2>。
(∑y)<2>
差方和(总)=∑y<2> - ───── (14)
mK
f(总)=mK-1
∑[∑y(k)]<2> (∑y)<2>
差方和(剂间)=────────-────── (15)
m mK
f(剂间)=K-1
∑[∑y<[m]>]<2> (∑y)<2>
差方和(区组间)=───────-─────── (16)
K mK
f(区组间)=m-1
差方和(误差)=差方和(总)-差方和(剂间)-差方和(区组间) (17)
f(误差)=f(总)-f(剂间)-f(区组间)=(K-1)(m-1)
各变异项差方和
各变异项方差=───────── (18)
各变异项自由度
差方和(误差)
误差项方差(s<2>)=────────
f(误差)
Km∑y<2>-k·∑[∑y(k)]<2>-m·∑(∑y<[m]>)<2>+(∑y)<2>
或s<2>=─────────────────────────────── (19)
Km(K-1)(m-1)
f=(K-1)(m-1)
差方和(误差)=差方和(总)-差方和(剂间) (20)
f(误差)=f(总)-f(剂间)=K(m-1)
m∑y<2>-∑[∑y(k)]<2>
s<2>=────────────── (21)
Km(m-1)
f=K(m-1)(4) 可靠性测验 通过对剂间变异的分析,以测验S和T的对数剂量和反应的关系是否显著偏离平行直线。(2.2)法和(2.2.2)法的剂间变异分析为试品间、回归、偏离平行三项,其他(k·k)法还需再分析二次曲线、反向二次曲线等。 可靠性测验的剂间变异分析 (k·k)法、(k·k')法。按表四计算各变异项的m·∑C<[i]><2>及∑(C〈i〉·∑y(k)),按(22)式计算各项变异的差方和。
[∑(C<[i]>·∑y(k))]<2>
各项变异的差方和=────────────── (22)
m·∑C<[i]><2>
f=1
(k·k.k)法按(23)、(24)式计算试品间差方和。
(2.2.2)法
(S<[2]>+S<[1]>)<2>+(T<[2]>+T<[1]>)<2>+(U<[2]>+U<[1]>)<2>
差方和(试品间)=──────────────────────────────
(∑y)<2>
3m - ───── (23)
mK
f=2
(3.3.3)法
(S<[1]>+S<[2]>+S<[3]>)<2>+(T<[1]>+T<[2]>+T<[3]>)<2>+(U<[1]>+U<[2]>+U<[3]>)<2>
差方和(试品间)=─────────────────────────────────────
(∑y)2
3m -─── (24)
mk
f=2按表五计算回归、二次曲线、反向二次曲线各项变异的m·∑C<[i]><2>及∑[C<[i]>·∑y(k)];按(22)式计算差方和(回归)、差方和(二次曲线)。 按(25)式计算差方和(偏离平行)及差方和(反向二次曲线)。
2∑[∑(C<[i]>·∑y(k))]<2>
差方和(偏离平行)、差方和(反向二次曲线)=──────────── (25)
∑(m·∑C<[i]><2>)f=2 按(18)式计算各项变异的方差。
表四 (K·K)法、(K·K')法可靠性测验正交多项系数表
━━━┯━━━━━━┯━━━━━━━━━━━━━━━━━┯━━━━━┯━━━━━━━━━━━━━━
│ │ ∑y(k)的正交多项系数(Ci) │ │
方法 │变 异 来 源 ├────────┬────────┤ m.∑Ci<2>│ ∑[Ci·∑y(k)]
│ │S1 S2 S3 S4 │ T1 T2 T3 T4 │ │
━━━┿━━━━━━┿━━━━━━━━┿━━━━━━━━┿━━━━━┿━━━━━━━━━━━━━━
(2.2)│试品间 │ -1 -1 │ 1 1 │ 4m │ T2+T1-S2-S1
│回归 │ -1 1 │ -1 1 │ 4m │ T2-T1+S2-S1
│偏离平行 │ 1 -1 │ -1 1 │ 4m │ T2-T1-S2+S1
───┼──────┼────────┼────────┼─────┼──────────────
(3.3)│试品间 │ -1 -1 -1 │ 1 1 1 │ 6m │ T3+T2+T1-S3-S2-S1
│回归 │ -1 0 1 │ -1 0 1 │ 4m │ T3-T1+S3-S1
│偏离平行 │ 1 0 -1 │ -1 0 1 │ 4m │ T3-T1-S3+S1
│二次曲线 │ 1 -2 1 │ 1 -2 1 │ 12m │ T3-2T2+T1+S3-2S2+S1
│反向二次曲线│ -1 2 -1 │ 1 -2 1 │ 12m │ T3-2T2+T1-S3+2S2-S1
───┼──────┼────────┼────────┼─────┼──────────────
(4.4)│试品间 │ -1 -1 -1 -1 │ 1 1 1 1 │ 8m │ T4+T3+T2+T1-S4-S3-S2-S1
│回归 │ -3 -1 1 3 │ -3 -1 1 3 │ 40m │3T4+T3-T2-3T1+3S4+S3-S2-3S1
│偏离平行 │ 3 1 -1 -3 │ -3 -1 1 3 │ 40m │3T4+T3-T2-3T1-3S4-S3+S2+3S1
│二次曲线 │ 1 -1 -1 1 │ 1 -1 -1 1 │ 8m │ T4-T3-T2+T1+S4-S3-S2+S1
│反向二次曲线│ -1 1 1 -1 │ 1 -1 -1 1 │ 8m │ T4-T3-T2+T1-S4+S3+S2-S1
───┼──────┼────────┼────────┼─────┼──────────────
(3.2)│试品间 │ -2 -2 -2 │ 3 3 │ 30m │ 3(T2+T1)-2(S3+S2+S1)
│回归 │ -2 0 2 │ -1 1 │ 10m │ T2-T1+2(S3-S1)
│偏离平行 │ 1 0 -1 │ -2 2 │ 10m │ 2(T2-T1)-S3+S1
│二次曲线 │ 1 -2 1 │ 0 0 │ 6m │ S3-2S2+S1
───┼──────┼────────┼────────┼─────┼──────────────
(4.3)│试品间 │ -3 -3 -3 -3 │ 4 4 4 │ 84m │ 4(T3+T2+T1)-3(S4+S3+S2+S1)
│回归 │ -3 -1 1 3 │ -2 0 2 │ 28m │ 2(T3-T1)+3(S4-S1)-S2+S3
│偏离平行 │ 3 1 -1 -3 │ -5 0 5 │ 70m │ 5(T3-T1)-3(S4-S1)-S3+S2
│二次曲线 │ 3 -3 -3 3 │ 2 -4 2 │ 60m │2(T3+T1)-4T2+3(S4-S3-S2+S1)
│反向二次曲线│ -1 1 1 -1 │ 1 -2 1 │ 10m │ T3-2T2+T1-S4+S3+S2-S1
━━━┷━━━━━━┷━━━━━━━━┷━━━━━━━━┷━━━━━┷━━━━━━━━━━━━━━注:表中字母T、S后面的数字1、2、3、4均表示其的下标 用(2.3)法及(3.4)法时,分别将(3.2)法及(4.3)法中S和T的正交多项系数互换即得。 表中S<[1]>、S<[2]>…T<[1]>、T<[2]>…在量反应分别为标准品和供试品每一剂量组内的反应值或它们规定函数的总和(相当于表二的∑y(k)各项)。所有足序1、2、3… 都是顺次由小剂量到大剂量,C<[i]>是与之相应的正交多项系数。m·∑C<[i]><2>是该项变异各正交多项系数的平方之和与m的乘积, ∑[C<[i]>·∑y(k)]为S<[1]>、S<[2]>…T<[1]>、T<[2]>…分别与该项正交多项系数乘积之和。 将方差分析结果列表进行可靠性测验。例如随机区组设计(3.3) 法可靠性测验结果列表,见表六。 表六中概率P是以该变异项的自由度为分子,误差项(S<2>)的自由度为分母,查F值表(表七),将查表所得F值与表六F项下的计算值比较而得。 当F计算值大于P=0.05或P=0.01的查表值时,则P<0.05或P<0.01,即为在此概率水平下该项变异有显著意义。
表五 (K·K·K)法可靠性测验正交多项系数表
━━━━┯━━━━┯━━━━━━━━━━━━━━━━━━━━━━┯━━━━━┯━━━━━━━━━━━━━━━
│ 变 异 │ ∑y(K)的正交多项系数(Ci) │ │
方 法 │ ├───────┬───────┬──────┤m·∑Ci<2>│ ∑[Ci. ∑y(k)]
│ 来 源 │S1 S2 S3 │T1 T2 T3 │ U1 U2 U3 │ │
────┼────┼───────┼───────┼──────┼─────┼───────────────
(2.2.2)│ 回 归 │-1 1 │-1 1 │-1 1 │ 6m │ S2-S1+T2-T1+U2-U1
├────┼───────┼───────┼──────┼─────┼───────────────
│ 偏 离 │ 1 -1 │-1 1 │ │ 4m │ T2-T1-S2+S1
│ 平 行 │ 1 -1 │ │-1 1 │ 4m │ U2-U1-S2+S1
│ │ │ 1 -1 │-1 1 │ 4m │ U2-U1-T2+T1
────┼────┼───────┼───────┼──────┼─────┼───────────────
(3.3.3)│ 回 归 │ -1 0 1 │ -1 0 1 │ -1 0 1 │ 6m │ U3-U1+T3-T1+S3-S1
├────┼───────┼───────┼──────┼─────┼───────────────
│ 偏 离 │ 1 0 -1 │ -1 0 1 │ │ 4m │ T3-T1-S3+S1
│ 平 行 │ 1 0 -1 │ │ -1 0 1 │ 4m │ U3-U1-S3+S1
│ │ │ 1 0 -1 │ -1 0 1 │ 4m │ U3-U1-T3+T1
├────┼───────┼───────┼──────┼─────┼───────────────
│二次曲线│ 1 -2 1 │ 1 -2 1 │ 1 -2 1 │ 18m │U3-2U2+U1+T3-2T2+T1+S3-2S2+S1
├────┼───────┼───────┼──────┼─────┼───────────────
│ 反向二 │ -1 2 -1 │ 1 -2 1 │ │ 12m │T3-2T2+T1-S3+2S2-S1
│ 次曲线 │ -1 2 -1 │ │ 1 -2 1 │ 12m │U3-2U2+U1-S3+2S2-S1
│ │ │ -1 2 -1 │ 1 -2 1 │ 12m │U3-2U2+U1-T3+2T2-T1
━━━━┷━━━━┷━━━━━━━┷━━━━━━━┷━━━━━━┷━━━━━┷━━━━━━━━━━━━━━━注:表中字母S、T、U后面的数字1、2、3均表示其的下标
表六 随机区组设计(3.3)法可靠性测验结果
──────┬─────┬────┬───────┬─────┬──
变异来源 │ f │ 差方和 │ 方 差 │ F │ P
──────┼─────┼────┼───────┼─────┼──
试品间 │ 1 │公式(22)│ 差方和/f │ 方差/s<2>│
回归 │ 1 │公式(22)│ 差方和/f │ 方差/s<2>│
偏离平行 │ 1 │公式(22)│ 差方和/f │ 方差/s<2>│
二次曲线 │ 1 │公式(22)│ 差方和/f │ 方差/s<2>│
反向二次曲线│ 1 │公式(22)│ 差方和/f │ 方差/s<2>│
──────┼─────┼────┼───────┼─────┼──
剂间 │ k-1 │公式(15)│ 差方和/f │ 方差/s<2>│
区组间 │ m-1 │公式(16)│ 差方和/f │ 方差/s<2>│
误差 │(k-1)(m-1)│公式(17)│差方和/f(s<2>)│ │
──────┼─────┼────┼───────┼─────┼──
总 │ mK-1 │公式(14)│ │ │
──────┴─────┴────┴───────┴─────┴──随机设计没有区组间变异项。 可靠性测验结果判断 可靠性测验结果,回归项应非常显著(p〈0.01)。 (2.2)法、(2.2.2)法偏离平行应不显著(p>0.05)。 其他(k.k)法、(k·k·k)法偏离平行、二次曲线、反向二次曲线各项均应不显著(p>0.05)。 试品间一项不作为可靠性测验的判断标准,试品间变异非常显著者,重复试验时,应参考所得结果重新估计T的效价或重新调整剂量试验。 双交叉设计的方差分析和可靠性测验 (1) 双交叉设计实验结果的方阵表 将动物按体重随机分成四组,各组的动物数(m)相等,四组的动物总数为4m。对四组中的每一只动物都加以识别标记,按双交叉设计给药次序表进行实验,各组的每一只动物都给药两次,共得2×4m个反应值。将S、T各两个剂量组两次实验所得反应值排列成表, 见表八。
表七 F值表
━━━━━┯━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
│ f1 (分子的自由度)
├────┬────┬────┬────┬────┬────┬────┬────┬────
│ 1 │ 2 │ 3 │ 4 │ 6 │ 12 │ 20 │ 40 │ ∞
──┬──┼────┼────┼────┼────┼────┼────┼────┼────┼────
│ 1 │ 161 │ 200 │ 216 │ 225 │ 234 │ 244 │ 248 │ 251 │ 254
│ │4052 │4999 │5403 │5625 │5859 │6106 │6208 │6286 │6366
│ 2 │ 18.51 │ 19.00 │ 19.16 │ 19.25 │ 19.33 │ 19.41 │ 19.44 │ 19.47 │ 19.50
│ │ 98.49 │ 99.00 │ 99.17 │ 90.25 │ 99.33 │ 99.42 │ 99.45 │ 99.48 │ 99.50
│ 3 │ 10.13 │ 9.55 │ 9.28 │ 9.12 │ 8.94 │ 8.74 │ 8.66 │ 8.60 │ 8.53
│ │ 34.12 │ 30.82 │ 29.46 │ 28.71 │ 27.91 │ 27.05 │ 26.69 │ 26.41 │ 26.12
│ 4 │ 7.71 │ 6.94 │ 6.59 │ 6.39 │ 6.16 │ 5.91 │ 5.80 │ 5.71 │ 5.63
│ │ 21.20 │ 18.00 │ 16.69 │ 15.98 │ 15.21 │ 14.37 │ 14.02 │ 13.74 │ 13.46
│ 5 │ 6.61 │ 5.79 │ 5.41 │ 5.19 │ 4.95 │ 4.68 │ 4.56 │ 4.46 │ 4.36
│ │ 16.26 │ 13.27 │ 12.06 │ 11.39 │ 10.67 │ 9.89 │ 9.55 │ 9.29 │ 9.02
f2 │ 6 │ 5.99 │ 5.14 │ 4.76 │ 4.53 │ 4.28 │ 4.00 │ 3.87 │ 3.77 │ 3.67
│ │ 13.74 │ 10.92 │ 9.78 │ 9.15 │ 8.47 │ 7.72 │ 7.39 │ 7.14 │ 6.88
分│ 7 │ 5.59 │ 4.74 │ 4.35 │ 4.12 │ 3.87 │ 3.57 │ 3.44 │ 3.34 │ 3.23
母│ │ 12.25 │ 9.55 │ 8.45 │ 7.85 │ 7.19 │ 6.47 │ 6.15 │ 5.90 │ 5.65
的│ 8 │ 5.32 │ 4.46 │ 4.07 │ 3.84 │ 3.58 │ 3.28 │ 3.15 │ 3.05 │ 2.93
自│ │ 11.26 │ 8.65 │ 7.59 │ 7.01 │ 6.37 │ 5.67 │ 5.36 │ 5.11 │ 4.86
由│ 9 │ 5.12 │ 4.26 │ 3.86 │ 3.63 │ 3.37 │ 3.07 │ 2.93 │ 2.82 │ 2.71
度│ │ 10.56 │ 8.02 │ 6.99 │ 6.42 │ 5.80 │ 5.11 │ 4.80 │ 4.56 │ 4.31
│ 10 │ 4.96 │ 4.10 │ 3.71 │ 3.48 │ 3.22 │ 2.91 │ 2.77 │ 2.67 │ 2.54
│ │ 10.04 │ 7.56 │ 6.55 │ 5.99 │ 5.39 │ 4.71 │ 4.41 │ 4.17 │ 3.91
│ 15 │ 4.54 │ 3.68 │ 3.29 │ 3.06 │ 2.79 │ 2.48 │ 2.33 │ 2.21 │ 2.07
│ │ 8.68 │ 6.36 │ 5.42 │ 4.89 │ 4.32 │ 3.67 │ 3.36 │ 3.12 │ 2.87
│ 20 │ 4.35 │ 3.49 │ 3.10 │ 2.87 │ 2.60 │ 2.28 │ 2.12 │ 1.99 │ 1.84
│ │ 8.10 │ 5.85 │ 4.94 │ 4.43 │ 3.87 │ 3.23 │ 2.94 │ 2.69 │ 2.42
│ 30 │ 4.17 │ 3.32 │ 2.92 │ 2.69 │ 2.42 │ 2.09 │ 1.93 │ 1.79 │ 1.62
│ │ 7.56 │ 5.39 │ 4.51 │ 4.02 │ 3.47 │ 2.84 │ 2.55 │ 2.29 │ 2.01
│ 40 │ 4.08 │ 3.23 │ 2.84 │ 2.61 │ 2.34 │ 2.00 │ 1.84 │ 1.69 │ 1.51
│ │ 7.31 │ 5.18 │ 4.31 │ 3.83 │ 3.29 │ 2.66 │ 2.37 │ 2.11 │ 1.81
│ 60 │ 4.00 │ 3.15 │ 2.76 │ 2.52 │ 2.25 │ 1.92 │ 1.75 │ 1.59 │ 1.39
│ │ 7.08 │ 4.98 │ 4.13 │ 3.65 │ 3.12 │ 2.50 │ 2.20 │ 1.93 │ 1.60
│ ∞ │ 3.84 │ 2.99 │ 2.60 │ 2.37 │ 2.09 │ 1.75 │ 1.57 │ 1.40 │ 1.00
│ │ 6.64 │ 4.60 │ 3.78 │ 3.32 │ 2.80 │ 2.18 │ 1.87 │ 1.59 │ 1.00
━━┷━━┷━━━━┷━━━━┷━━━━┷━━━━┷━━━━┷━━━━┷━━━━┷━━━━┷━━━━
上行,P=0.05;下行,P=0.01
表八 双交叉实验结果
━┯━━━━━━━━━━━━━━━━━━┯━━━━━━━━━━━━━━━━━━┯━━━━━━━━━━━━━━━━━━┯━━━━━━━━━━━━━━━━━━┯━━━━━
│ 第 一 组 │ 第 二 组 │ 第 三 组 │ 第 四 组 │
├─────┬─────┬──────┼─────┬─────┬──────┼─────┬─────┬──────┼─────┬─────┬──────┤
│ 第(1)次 │ 第(2)次 │ 两次 │ 第(1)次 │ 第(2)次 │ 两次 │ 第(1)次 │ 第(2)次 │ 两次 │ 第(1)次 │ 第(2)次 │ 两次 │
├─────┼─────┤ ├─────┼─────┤ ├─────┼─────┤ ├─────┼─────┤ │
│ ds<[1]> │ dT<[2]> │ 反应和 │ ds<[2]> │ dT[1] │ 反应和 │ dT<[1]> │ dS<[2]> │ 反应和 │ dT<[2]> │ dS<[1]> │ 反应和 │
─┼─────┼─────┼──────┼─────┼─────┼──────┼─────┼─────┼──────┼─────┼─────┼──────┼─────
│ys<[1]>(1)│yT<[2]>(2)│ y(1)+y(2) │ys<[2]>(1)│yT<[1]>(2)│ y(1)+y(2) │yT<[1]>(1)│ys<[2]>(2)│ y(1)+y(2) │yT<[2]>(1)│ys<[1]>(2)│ y(1)+y(2) │
y│ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │
│ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │
│ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │ . │总和
─┼─────┼─────┼──────┼─────┼─────┼──────┼─────┼─────┼──────┼─────┼─────┼──────┼─────
│S<[1]>(1) │ │ │ │ │ │ │ │ │ │S<[1]>(2) │ │ S<[1]>
∑│ │ │ │ S<[2]>(1)│ │ │ │S<[2]>(2) │ │ │ │ │ S<[2]>
│ │ │ │ │ T<[1]>(2)│ │T<[1]>(1) │ │ │ │ │ │ T<[1]>
│ │ T<[2]>(2)│ │ │ │ │ │ │ │ T<[2]>(1)│ │ │ T<[2]>
━┷━━━━━┷━━━━━┷━━━━━━┷━━━━━┷━━━━━┷━━━━━━┷━━━━━┷━━━━━┷━━━━━━┷━━━━━┷━━━━━┷━━━━━━┷━━━━━
(2) 缺项补足 表八中如有个别组的1个反应值因故缺失,均作该只动物缺失处理,在组内形成两个缺项。此时,可分别用两次实验中该组动物其余各反应值的均值补入;也可在其余三组内用严格随机的方法各去除一只动物,使各组的动物数相等。每补足一个缺项,误差(Ⅰ)和误差(Ⅱ)的方差S<2>和S<2>Ⅱ的自由度都要减去1。缺项不得超过反应总个数的5%。同一组内缺失的动物不得超过1只。 (3) 方差分析 双交叉设计的总变异中,包含有动物间变异和动物内变异。对表八的2×4m个反应值进行方差分析时,总变异的差方和(总)按(26)式计算。 (∑y)<2> 差方和(总)=∑y<2>-────── (26) 2×4m f总=2×4m-1 动物间变异是每一只动物两次实验所得反应值的和(表八每组动物的第三列)之间的变异,其差方和按(27)式计算。 ∑[y(1)+y(2)]<2> (∑y)<2> 差方和(动物间)=────────── - ────── (27) 2 2×4m f(动物间)=4m-1 总变异中分除动物间变异,余下为动物内变异。 动物间变异和动物内变异的分析 将表八中S和T各剂量组第(1)次实验所得反应值之和S<[1]>(1)、S<[2]>(1)、T<[1]>(1)、T<[2]>(1)及第(2)次实验反应值之和S<[1]>(2)、S<[2]>(2)、T<[1]>(2)、T<[2]>(2)按表九双交叉设计正交系数表计算各项变异的m·∑C<[i]><2>及∑(C<[i]>.y),按公式(22)计算各项变异的差方和。 总变异的差方和减去动物间变异的差方和,再减去动物内各项变异的差方和,余项为误差(Ⅰ)的差方和,按(28)式计算。 差方和(误差Ⅰ)=差方和(总)-差方和(动物间)-差方和(试品间)-差方和(回归)-差方和(次间)-差方和(次间×偏离平行) (28) f(误差Ⅰ)=f(总)-f(动物间)-f(试品间)-f(回归)-f(次间) -f(次间×偏离平行)=4(m-1) 误差(Ⅰ)的方差S<2>,用以计算实验误差S<[M]>、FL,及进行动物内各项变异(* 标记者)的F测验。 误差(Ⅱ)的差方和为动物间变异的差方和减去表九中其余三项变异(无* 标记者)的差方和,按(29)式计算。 差方和(误差Ⅱ)=差方和(动物间)-差方和(偏离平行)-差方和(次间×试品间) -差方和(次间×回归) (29) f(误差Ⅱ)=f(动物间)-f(偏离平行)-f(次间×试品间)-f(次间×回归)=4(m-1) 误差(Ⅱ)的方差S<2>Ⅱ用以进行上述三项变异的F测验。 (4) 可靠性测验 将方差分析及F测验的结果列表,如表十。 表十中的概率P,计算同表六,但表的上半部分是以S<2>Ⅱ的自由度为分母,表的下半部分以S<[2]>的自由度为分母,查F值表(表七),将查表所得的F值与表+F项下的计算值比较而得。 可靠性测验结果判断 回归、偏离平行、试品间三项的判断标准同(2.2)法。 次间×试品间、次间×回归、次间×偏离平行三项中,如有F测验非常显著者,说明该项变异在第一次和第二次实验的结果有非常显著的差别,对出现这种情况的检定结果,下结论时应慎重,最好复试。
表九 双交叉设计正交系数表
━━━━━━━━┯━━━━━━━━━━┯━━━━━━━━━━┯━━━━━┯━━━━━━━━━━━━━━━━━━━━━━━━
│ 第(1)次实验 │ 第(2)次实验 │ │
├──────────┼──────────┤ │ 变 异 来 源
│S1(1)S2(1)T1(1)T2(1)│S1(2)S2(2)T1(2)T2(2)│m·∑Ci<2>│ ∑(Ci·∑y)
├──────────┴──────────┤ │
│ 正 交 多 项 系 数 Ci │ │
────────┼──────────┬──────────┼─────┼────────────────────────
试 品 间* │ -1 -1 1 1 │ -1 -1 1 1 │ 8m │T2(1)+T1(1)-S2(1)-S1(1)+T2(2)+T1(2)-S2(2)-S1(2)
回 归* │ -1 1 -1 1 │ -1 1 -1 1 │ 8m │T2(1)-T1(1)+S2(1)-S1(1)+T2(2)-T1(2)+S2(2)-S1(2)
偏离平行 │ 1 -1 -1 1 │ 1 -1 -1 1 │ 8m │T2(1)-T1(1)-S2(1)+S1(1)+T2(2)-T1(2)-S2(2)+S1(2)
次 间* │ -1 -1 -1 -1 │ 1 1 1 1 │ 8m │T2(2)+T1(2)+S2(2)+S1(2)-T2(1)-T1(1)-S2(1)-S1(1)
次间×试品间 │ 1 1 -1 -1 │ -1 -1 1 1 │ 8m │T2(2)+T1(2)-S2(2)-S1(2)-T2(1)-T1(1)+S2(1)+S1(1)
次间×回归 │ 1 -1 1 -1 │ -1 1 -1 1 │ 8m │T2(2)-T1(2)+S2(2)-S1(2)-T2(1)+T1(1)-S2(1)+S1(1)
次间×偏离平行* │ -1 1 1 -1 │ 1 -1 -1 1 │ 8m │T2(2)-T1(2)-S2(2)+S1(2)-T2(1)+T1(1)+S2(1)-S1(1)
━━━━━━━━┷━━━━━━━━━━┷━━━━━━━━━━┷━━━━━┷━━━━━━━━━━━━━━━━━━━━━━━━
注:表中字母S、T后面的数字1(1)、2(2)、1(2)、2(1)均表示其的下标。①各项变异的自由度均为1。有*号标记的四项为动物内变异,其余三项为动物间变异。
表十 双交叉设计可靠性测验结果
────────┬────┬─────┬────────┬──────┬───
变异来源 │ f │ 差方和 │ 方 差 │ F │ P
────────┼────┼─────┼────────┼──────┼───
偏离平行 │ 1 │ 公式(22) │ 差方和/f │方差/S<2>Ⅱ │
次间×试品间 │ 1 │ 公式(22) │ 差方和/f │方差/S<2>Ⅱ │
次间×回归 │ 1 │ 公式(22) │ 差方和/f │方差/S<2>Ⅱ │
误差 (Ⅱ) │ 4(m-1) │ 公式(29) │差方和/f(S<2>Ⅱ)│ │
────────┼────┼─────┼────────┼──────┼───
动物间 │ 4m-1 │ 公式(27) │ 差方和/f │ 方差/S<2>│
试品间 │ 1 │ 公式(22) │ 差方和/f │ 方差/S<2>│
回归 │ 1 │ 公式(22) │ 差方和/f │ 方差/S<2>│
次间 │ 1 │ 公式(22) │ 差方和/f │ 方差/S<2>│
次间×偏离平行 │ 1 │ 公式(22) │ 差方和/f │ 方差/S<2>│
误差 (Ⅰ) │ 4(m-1) │ 公式(28) │差方和/f(S<2>) │ │
────────┼────┼─────┼────────┼──────┼───
总 │2×4m-1 │ 公式(26) │ │ │
────────┴────┴─────┴────────┴──────┴────
3.效价(PT)及可信限(FL)计算各种(k·K)法都按表十一计算V、W、D、A、B、g等数值,代入公式(30)~(33)及(3)、(8)计算R、P<[T]>、S<[M]>以及R、P<[T]>的FL和FL%等。 Ⅳ R=D·antilog──── (30) W I ┌───────────── S<[M]>=─────── √ms<2>[(1-g)AW<2>+BV<2>] (31) W<2>(1-g) logR R的FL=antilog[─────± t·S<[M]> ] (32) 1 - g logR P<[T]>的FL=A<[T]>·antilog [ ────± t·S<[M]> ] (33) 1 - g (2.2) 法双交叉设计 计算方法同上述(2.2)法。 双交叉设计各剂量组都进行两次试验,S和T每一剂量组的反应值个数为组内动物数的两倍(2m)。 (1)双交叉设计用S和T各组剂量两次试验所得各反应值之和(表八中的S<[1]>、S<[2]>、T<[1]>、T<[2]>)按表十一(2.2)法公式计算V、W、D、g等数值。
表十一 量反应平行线检定法的计算公式①
━━━┯━━━━━━━━┯━━━━━━━━┯━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┯━━━━━━━
│ │ │ 效 价 计 算 用 数 值 │Sm计算用数值
方法│ S │ T ├────────────────┬──────────────────┬──────┼─┬─┬───
│ │ │ V │ W │ D │A │B │ g
───┼────────┼────────┼────────────────┼──────────────────┼──────┼─┼─┼───
│ │ │ 1 │ 1 │ ds2 │ │ │t<2>s<2>
m 2.2 │ds1 ds2 │dT1 dT2 │──(T1+T2-S1-S2) │──(T2-T1+S2-S1) │─── │1 │1 │───
│ │ │ 2 │ 2 │ dT2 │ │ │ W<2>
───┼────────┼────────┼────────────────┼──────────────────┼──────┼─┼─┼───
│ │ │ 1 │ 1 │ ds3 │2 │1 │t<2>s<2>
m 3.3 │ds1 ds2 ds3 │dT1 dT2 dT3 │──(T1+T2+T3-S1-S2-S3) │──(T3-T1+S3-S1) │─── │─│─│───
│ │ │ 3 │ 4 │ dT3 │3 │4 │ 4W<2>
───┼────────┼────────┼────────────────┼──────────────────┼──────┼─┼─┼───
│ │ │ 1 │ 1 │ ds4 │1 │1 │t<2>s<2>
m 4.4 │ds1 ds2 ds3 ds4 │dT1 dT2 dT3 dT4 │──(T1+T2+T3+T4-S1-S2-S3-S4) │──[(T3-T2+S3-S2)+3(T4-T1+S4-S1)] │─── │─│─│───
│ │ │ 4 │ 20 │ dT4 │2 │10│ 10W<2>
───┼────────┼────────┼────────────────┼──────────────────┼──────┼─┼─┼───
│ │ │ 1 1 │ 1 │ ds3 1 │ │ │
m 3.2 │ds1 ds2 ds3 │dT1 dT2 │──(T2+T1)-──(S1+S2+S3) │──[(T2-T1)+2(S3-S1)] │──·───│ │ │
│ │ │ 2 3 │ 5 │ dT2 ── │ │ │
│ │ │ │ │ r │5 │2 │2t<2>s<2>
───┼────────┼────────┼────────────────┼──────────────────┼──────┤─│─│───
│ │ │ 1 1 │ 1 │ ds2 ── │6 │5 │ 5W<2>
m 2.3 │ds1 ds2 │dT1 dT2 dT3 │──(T1+T2+T3)-──(S1+S2) │──[2(T3-T1)+(S2-S1)] │──· r │ │ │
│ │ │ 3 2 │ 5 │ dT3 │ │ │
───┼────────┼────────┼────────────────┼──────────────────┼──────┼─┼─┼───
│ │ │ 1 1 │ 1 │ ds4 1 │ │ │
m 4.3 │ds1 ds2 ds3 ds4 │dT1 dT2 dT3 │──(T1+T2+T3)-──(S1+S2+S3+S4)│──[2(T3-T1)+(S3-S2)+3(S4-S1)] │── ·── │ │ │
│ │ │ 3 4 │ 14 │ ─ │ │ │
│ │ │ │ │ dT3 r │7 │1 │ t<2>s<2>
───┼────────┼────────┼────────────────┼──────────────────┼──────┤─│─│───
│ │ │ 1 1 │ 1 │ ds3 ──│12│7 │ 7W<2>
m 3.4 │ds1 ds2 ds3 │dT1 dT2 dT3 dT4 │──(T1+T2+T3+T4)-──(S1+S2+S3)│──[2(S3-S1)+(T3-T2)+3(T4-T1)] │───·r │ │ │
│ │ │ 4 3 │ 14 │ dT4 │ │ │
───┼────────┼────────┼────────────────┼──────────────────┼──────┼─┼─┼───
│ │ │ 1 │ 1 │ ds2 │ │ │
m2.2.2│ds1 ds2 │dT1 dT2 │ ──(T1+T2-S1-S2) │──(T2-T1+U2-U1+S2-S1) │ ─── │ │ │
│ │ │ 2 │ 3 │ dT2 │ │2 │2t<2>s<2>
│ │ │ 1 │ │ ds2 │1 │─│───
│ │dU1dU2 │──(U1+U2-S1-S2) │ │ ─── │ │3 │ 3W<2>
│ │ │ 2 │ │ dU2 │ │ │
───┼────────┼────────┼────────────────┼──────────────────┼──────┼─┼─┼───
│ │ │ 1 │ 1 │ ds3 │ │ │
m3.3.3│ds1ds2ds3 │dT1 dT2 dT3 │──(T1+T2+T3-S1-S2-S3) │ ──(T3-T1+U3-U1+S3-S1) │─── │ │ │
│ │ │ 3 │ 6 │ dT3 │2 │1 │t<2>s<2>
│ │ │ │ │ │─│─│───
│ │ │ 1 │ │ ds3 │3 │6 │ 6W<2>
│ │dU1 dU2 dU3 │ ──(U1+U2+U3-S1-S2-S3) │ │ ─── │ │ │
│ │ │ 3 │ │ dU3 │ │ │
━━━┷━━━━━━━━┷━━━━━━━━┷━━━━━━━━━━━━━━━━┷━━━━━━━━━━━━━━━━━━┷━━━━━━┷━┷━┷━━━
注:表中字母S、T、U后面的字数1、2、3均表示其的下标。① 表中ds、dT分别为S和T的剂量,下角1、2、3是顺次由小剂到大剂。(2) 参照(31)式计算S<[M]>,因每只动物进行两次实验,式中m用2m代替,(2.2)法A=1,B=1,S<[M]>的公式为 Ⅰ ┌──────────── S<[M]>=───────√2ms<2>[(1-g)W<2>+V<2>] (34) W<2> (1-g) 式中 S<2>为表十中误差(Ⅰ)的方差 s<2>·t<2>·2m g=─────── W<2> 例2 量反应平行线测定随机设计(3.3.3)法 绒促性素(HCG)效价测定──小鼠子宫增重法 S为绒促性素标准品 ds<[1]>∶0.135u/鼠 ds<[2]>∶0.225u/鼠 ds<[3]>∶0.375u/鼠 T为绒促性素 估计效价A<[T]>∶2500u/mg dT<[1]>∶0.135u/鼠 dT<[2]>∶0.225u/鼠 dT<[3]>∶0.375u/鼠 U为绒促性素粉针,标示量AU∶500u/安瓿 dU<[1]>∶0.144u/鼠 dU<[2]>∶0.240u/鼠 dU<[3]>∶0.400U/鼠 r=1:0.6 I=0.222 反应(y):10g体重的子宫重(mg) 测定结果见表2-1。
表2-1 HCG效价测定结果
━━━━━━┯━━━━━━━━━━━┯━━━━━━━━━━━┯━━━━━━━━━━━┯━━━━━
剂 量 │ds<[1]> ds<[2]>ds<[3]>│dT<[1]> dT<[2]>dT<[3]>│dU<[1]> dU<[2]>dU<[3]>│
u/ 鼠 │ 0.135 0.225 0.375 │ 0.135 0.225 0.375│ 0.144 0.240 0.400 │
──────┼───────────┼───────────┼───────────┼─────
│ 9.31 33.70 15.10 │20.80 25.70 35.60 │ 26.20 10.00 55.00 │
│17.50 56.80 47.20 │16.40 6.37 48.40 │ 10.00 40.20 41.70 │
│21.90 44.60 51.80 │ 5.66 38.30 41.90 │ 19.22 22.30 15.40 │
│14.60 32.30 47.30 │ 9.50 46.80 44.70 │ 22.00 40.50 53.60 │
│ 8.20 16.70 49.90 │ 9.27 43.40 29.80 │ 20.70 50.90 53.70 │
│11.00 6.17 47.20 │ 7.56 27.80 38.80 │ 23.20 23.50 33.00 │
│24.40 41.50 47.10 │15.40 26.00 37.40 │ 18.70 19.60 44.30 │
y │16.80 36.20 45.10 │20.30 27.20 33.70 │ 12.60 27.20 44.70 │
│29.90 9.83 46.40 │11.50 27.30 35.40 │ 20.90 30.30 23.00 │
│ 8.95 20.00 52.90 │22.20 11.90 47.90 │ 19.10 58.80 31.60 │
│17.80 22.00 32.50 │20.60 33.40 14.60 │ 19.40 55.30 49.20 │
│18.00 60.60 56.40 │13.90 29.00 49.80 │ 14.50 40.70 55.30 │
│13.70 6.43 39.50 │12.60 6.43 14.50 │ 11.40 35.40 23.80 │
│ 8.82 26.00 8.08 │ 7.25 27.80 42.00 │ 16.20 15.20 21.80 │
│17.80 34.80 37.10 │15.80 17.70 11.50 │ 20.80 28.70 36.00 │
──────┼───────────┼───────────┼───────────┼─────
│238.68 447.63 623.58│208.74 395.10 526.00│274.92 498.60 582.10│ ∑y ∑y(k)
│ │ │ │
│S<[1]> S<[2]> S<[3]> │ T<[1]> T<[2]> T<[3]>│U<[1]> U<[2]> U<[3]>│ 3795.35
━━━━━━┷━━━━━━━━━━━┷━━━━━━━━━━━┷━━━━━━━━━━━┷━━━━━(3.3.3)法,K=9;每组15只小鼠,m=15 (1) 按(14)、(15)、(20)式计算各项的差方和
(3795.35)<2>
差方和(总)=9.31<2>+17.50<2>+…+23.80<2>+21.80<2>+36.00<2> - ───── =29 868.26
9×15
f(总)=9×15-1=134
238.68<2>+477.63<2>+…+582.10<2> (3795.35)<2>
差方和(剂间)=──────────────── - ──────=12 336.55
15 9×15
f(剂间)=9-1=8差方和(误差)=29 868.26-12 336.55=17 531.71 f(误差)=134-8=126 (2) 剂间变异分析及可靠性测验 按(24)式及表五(3.3.3)法分析。
(238.68+447.63+623.58)<2>+(208.74+395.10+526.00)<2>
差方和(试品间)=─────────────────────────
3×15
(274.92+498.60+582.10)<2> (3795.35)<2>
+ ───────────── - ──────
3×15 9×15
=633.23f(试品间)=2 各项分析结果见表2-2、表2-3。 结论:回归非常显著,偏离平行、二次曲线、反向二次曲线均不显著,实验结果成立。 (3) 效价(P<[T]>)(P<[U]>)及可信限(FL)计算 按表十一(3.3.3)法及(30)~(33)、(3)、(8)式计算。 r=1:0.6 I=0.222s<2>=139.14 f=126 t=1.98 P<[T]>及其FL计算 1 V=──(208.74+395.10+526.00-238.68-447.63-623.58)=-60.017 3 1 W=──(526.00-208.74+623.58-238.68+582.10-274.92)=168.223 6
表2-2 HCG(3.3.3)法剂间变异分析
━━━━━━┯━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┯━━━━━┯━━━━━━━━━┯━━━━━━━━━━━━━━━━━━━━━━━━━━
│ ∑V(k) │ │ │ 差 方 和
├───┬───┬───┬───┬───┬───┬───┬───┬───┤ │ ├───────────┬──────────────
│S<[1]>│S<[2]>│S<[3]>│T<[1]>│T<[2]>│T<[3]>│U<[1]>│U<[2]>│U<[3]>│ 分 母 │ │ │
│ │ │ │ │ │ │ │ │ │ │┌────────│ │
变异来源 │ │ │ │ │ │ │ │ │ │ │∑[C<[i]>. ∑y(k)]│[∑(C<[i]>.∑y(k))]<2>│ 2∑[∑(C<[i]>·∑y(k))]<2>
│238.68│447.63│623.58│208.74│395.10│526.00│274.92│498.60│582.10│ │ │───────────│──────────────
├───┴───┴───┴───┴───┴───┴───┴───┴───┤m·∑Ci<2>│ │ m·∑C<[i]><2> │ ∑(m·∑C<[i]><2>)
│ 正交多项系数C<[i]> │ │ │ │
──────┼───┬───┬───┬───┬───┬───┬───┬───┬───┼─────┼─────────┼───────────┼──────────────
回归 │ -1 │ 0 │ 1 │ -1 │ 0 │ 1 │ -1 │ 0 │ 1 │ 15×6 │ 1009.34 │ 11319.64 │
──────┼───┼───┼───┼───┼───┼───┼───┼───┼───┼─────┼─────────┼───────────┼──────────────
│ │ │ │ │ │ │ │ │ │ │ │ │
│ 1 │ 0 │ -1 │ -1 │ 0 │ 1 │ │ │ │ 15×4 │ -67.64 │ │ 119.08
偏离平行 │ 1 │ 0 │ -1 │ │ │ │ -1 │ 0 │ 1 │ 15×4 │ -77.72 │ │
│ │ │ │ 1 │ 0 │ -1 │ -1 │ 0 │ 1 │ 15×4 │ -10.08 │ │
│ │ │ │ │ │ │ │ │ │ │ │ │
──────┼───┼───┼───┼───┼───┼───┼───┼───┼───┼─────┼─────────┼───────────┼──────────────
二次曲线 │ 1 │ -2 │ 1 │ 1 │ -2 │ 1 │ 1 │ -2 │ 1 │ 15×18 │ -228.64 │ 193.62 │
──────┼───┼───┼───┼───┼───┼───┼───┼───┼───┼─────┼─────────┼───────────┼──────────────
│ │ │ │ │ │ │ │ │ │ │ │ │
│ -1 │ 2 │ -1 │ 1 │ -2 │ 1 │ │ │ │ 15×12 │ -22.46 │ │ 71.0
反向二次曲线│ -1 │ 2 │ -1 │ │ │ │ 1 │ -2 │ 1 │ 15×12 │ -107.18 │ │
│ │ │ │ -1 │ 2 │ -1 │ 1 │ -2 │ 1 │ 15×12 │ -87.72 │ │
│ │ │ │ │ │ │ │ │ │ │ │ │
━━━━━━┷━━━┷━━━┷━━━┷━━━┷━━━┷━━━┷━━━┷━━━┷━━━┷━━━━━┷━━━━━━━━━┷━━━━━━━━━━━┷━━━━━━━━━━━━━━
表2-3 HCG效价测定(3.3.3)法可靠性测验结果
───────┬────┬─────┬──────┬────┬───
变异来源 │ f │ 差方和 │ 方 差 │ F │ P
───────┼────┼─────┼──────┼────┼───
试品间 │ 2 │ 633.2│ 316.6 │ 2.28│ >0.05
回归 │ 1 │ 11 319.64│11 319.64 │ 81.35│ <0.01
偏离平行 │ 2 │ 119.08│ 59.54 │ <1 │ >0.05
二次曲线 │ 1 │ 193.62│ 193.62 │ 1.39│ >0.05
反向二次曲线 │ 2 │ 71.00│ 35.50 │ <1 │ >0.05
───────┼────┼─────┼──────┼────┼───
剂间 │ 8 │ 12 336.55│1542.07 │ 11.08 │ <0.01
误差 │ 126 │ 17 531.71│139.14(s<2>)│ │
───────┼────┼─────┼──────┼────┼───
总 │ 134 │ 29 868.26│ │ │
───────┴────┴─────┴──────┴────┴───
139.14×1.98<2>×15
g=───────── =0.048
6×168.223<2>
0.375 -60.017
R<[T]>=────·antilog(─────×0.222)=0.833
0.375 168.223
P<[T]>=2500×0.833=2082.5u/mg
0.222
SM<[T]>= ──────────×168.223<2>(1-0.048)
┌─────────────────────────
│ 2 1
│15×139.14[(1-0.048)─×168.223<2>+ ─(-60.017)2]=0.051 29
√ 3 6
log0.833
R<[T]>的FL=antilog[─────±1.98×0.051 29]=0.653~1.043
1-0.048
P<[T]>的FL=2500(0.653~1.043)=1632.5~2607.5u/mg
2607.5-1632.5
P<[T]>的FL%=[───────×100]%=23.4%
2×2082.5
P<[U]>及其FL计算
1
V=──(274.92+498.60+582.10-238.68-447.63-623.58)=15.243
3
W=168.223
g=0.048
0.375 15.243
R<[U]>=────. antilog(─────×0.222)=0.982
0.400 168.223
P<[U]>=500×0.982=491.0u/安瓿
0.222
S<[MU]>=──────────×168.223<2>(1-0.048)
┌─────────────────────────
│ 2 1
│15×139.14[(1-0.048)─×168.223<2>+ ─×15.243<2>]=0.050
│ 3 6
√
log0.982
R<[U]>的FL=antilog[─────±1.98×0.050 51]=0.779~1.235
(1-0.048)
P<[U]>的FL=500(0.779~1.235)=389.5~617.5u/安瓿
617.5-389.5
P<[U]>的FL%=[───────×100]%=23.2%
2×491.0
按(21)式计算S<2>
15(9.31<2>+17.50<2>+…+21.80<2>+36.00<2>)-(238.68<2>+447.63<2>+… +582.10<2>)
S<2>=─────────────────────────────────────── =139.14
9×15(15-1)
和表2-3结果相同。(注):由于计算机中无平方表示符,故用“数字<2>”表示平方。 |