韓進(jìn)誠(chéng),瞿紅俠
(商丘師范學(xué)院 生物與食品學(xué)院,河南 商丘 476000)
在動(dòng)物營(yíng)養(yǎng)與飼料科學(xué)試驗(yàn)中,營(yíng)養(yǎng)物質(zhì)相對(duì)生物學(xué)效價(jià)評(píng)定是一項(xiàng)重要的研究?jī)?nèi)容.以常用的營(yíng)養(yǎng)物質(zhì)為參照物,評(píng)定另一種新型營(yíng)養(yǎng)物質(zhì)的相對(duì)生物學(xué)效價(jià),就可以根據(jù)該生物學(xué)效價(jià)推算新型營(yíng)養(yǎng)物質(zhì)的適宜添加量.在營(yíng)養(yǎng)物質(zhì)相對(duì)生物學(xué)效價(jià)評(píng)定方法中,Littell等(1997)提出的斜率比方法被廣泛應(yīng)用.
本文利用Littell等(1997)提出的斜率比方法,結(jié)合課題組試驗(yàn)結(jié)果(Han等,2016),介紹兩種維生素D的相對(duì)生物學(xué)效價(jià)比較的SAS程序.計(jì)算營(yíng)養(yǎng)物質(zhì)相對(duì)生物學(xué)效價(jià)需要經(jīng)過(guò)多項(xiàng)式比較、求解斜率及其標(biāo)準(zhǔn)誤、求解相對(duì)生物學(xué)效價(jià)及其標(biāo)準(zhǔn)誤等步驟.
以1~21日齡肉雞為試驗(yàn)動(dòng)物,用維生素D3(VD3)為參照物,用肉雞股骨灰分重量作為評(píng)價(jià)指標(biāo),評(píng)估25-羥基維生素D3(25-OH-D3)相對(duì)生物學(xué)效價(jià)(Han等,2016).試驗(yàn)設(shè)計(jì)9個(gè)處理,每個(gè)處理5個(gè)重復(fù).向1~21日齡肉雞基礎(chǔ)飼糧中添加不同水平VD3(0 μg/kg、2.5 μg/kg、5.0 μg/kg、10.0 μg/kg、20.0 μg/kg)和25-OH-D3(1.25 μg/kg、2.5 μg/kg、5.0 μg/kg、10.0 μg/kg),VD3食入量、25-OH-D3食入量、股骨灰分重量的試驗(yàn)結(jié)果列于表1.
表1 不同水平VD3和25-OH-D3對(duì)21日齡肉雞股骨灰分重量的影響
續(xù)表1
用GLM中Contrast方法確定不同維生素D水平與肉雞股骨灰分重量的線(xiàn)性關(guān)系.
(1)SAS程序
Data contrast;
Do t=1 to 5;
Do r=1 to 5;
Input x@@;
Output;
End;
End;
Drop r;
Cards;
0.204 0.131 0.115 0.206 0.170
0.203 0.216 0.217 0.205 0.167
0.267 0.301 0.267 0.300 0.320
0.367 0.416 0.411 0.437 0.353
0.466 0.400 0.487 0.458 0.457
Proc glm;
Class t;
Model x=t;
Means t/tukey;
Contrast′linear′t-2-1 0 1 2;
Contrast′quadratic′t 2-1-2-1 2;
Run;
(2)輸出結(jié)果
Contrast DF Contrast SS Mean Square F Value Pr>F
Linear 1 0.29799200 0.29799200 299.82 <0.0001
Quadratic 1 0.00116851 0.00116851 1.18 0.2911
結(jié)果顯示:VD3線(xiàn)性影響肉雞股骨灰分重量(P<0.0001),二次曲線(xiàn)影響差異不顯著(P=0.2911).數(shù)據(jù)表明:1~21 d肉雞VD3食入量與股骨灰分重量之間存在顯著線(xiàn)性關(guān)系.結(jié)果列于表2.
(1)SAS程序
Data contrast;
Do t=1 to 5;
Do r=1 to 5;
Input x@@;
Output;
End;
End;
Drop r;
Cards;
0.204 0.131 0.115 0.206 0.170
0.225 0.227 0.275 0.201 0.206
0.268 0.284 0.298 0.280 0.318
0.411 0.446 0.371 0.318 0.342
0.443 0.499 0.432 0.463 0.414
Proc glm;
Class t;
Model x=t;
Means t/tukey;
Contrast′linear′t-2-1 0 1 2;
Contrast′quadratic′t 2-1-2-1 2;
Run;
(2)輸出結(jié)果
Contrast DF Contrast SS Mean Square F Value Pr>F
Linear 1 0.25977632 0.25977632 194.80 <0.0001
Quadratic 1 0.00079566 0.00079566 0.60 0.4489
結(jié)果顯示:25-OH-D3線(xiàn)性影響肉雞股骨灰分重量(P<0.0001)),二次曲線(xiàn)影響差異不顯著(P=0.4489).數(shù)據(jù)表明:1~21 d肉雞25-OH-D3食入量與股骨灰分重量之間存在顯著線(xiàn)性關(guān)系.結(jié)果列于表2.
表2 VD3和25-OH-D3與21日齡肉雞股骨灰分重量的關(guān)系
續(xù)表2
(1)SAS程序
Data glm;
Input xs xt y@@;
Output;
Cards;0 0 0.204
0 0 0.131
0 0 0.115
0 0 0.206
0 0 0.170
1.484 0 0.203
1.773 0 0.216
1.553 0 0.217
1.379 0 0.205
1.576 0 0.167
3.393 0 0.267
3.547 0 0.301
3.515 0 0.267
3.990 0 0.300
3.290 0 0.320
9.244 0 0.367
9.510 0 0.416
9.321 0 0.411
9.638 0 0.437
8.890 0 0.353
19.222 0 0.466
18.960 0 0.400
18.720 0 0.487
19.096 0 0.458
19.480 0 0.457
0 0.720 0.225
0 0.687 0.227
0 0.676 0.275
0 0.782 0.201
0 0.791 0.206
0 1.623 0.268
0 1.838 0.284
0 1.920 0.298
0 1.885 0.280
0 1.508 0.318
0 4.170 0.411
0 4.495 0.446
0 4.180 0.371
0 4.345 0.318
0 4.470 0.342
0 8.660 0.443
0 9.000 0.499
0 9.102 0.432
0 9.730 0.463
0 9.840 0.414
Proc glm;
Model y=xs xt;
Run;
(2)輸出結(jié)果
R-Square Coeff Var Root MSE y Mean
0.808955 15.06982 0.047761 0.316933
Parameter Estimate Standard Error t Value Pr>|t|
Intercept 0.2156554039 0.01041173 20.71 <0.0001
xs 0.0138951691 0.00123959 11.21 <0.0001
xt 0.0277155565 0.00256662 10.80 <0.0001
結(jié)果顯示:VD3(即xs)斜率=0.0139,標(biāo)準(zhǔn)誤SE=0.0012;25-OH-D3(即xt)斜率=0.0277,標(biāo)準(zhǔn)誤SE=0.0026;決定系數(shù)R2=0.81,概率值P<0.0001.結(jié)果列于表3.
數(shù)據(jù)表明:以VD3為參照物,用21 d肉雞股骨灰分重量作為指標(biāo),25-OH-D3相對(duì)生物學(xué)效價(jià)RBV=0.0277/0.0139=1.99,即25-OH-D3相對(duì)生物學(xué)效價(jià)約為VD3的199%.
該程序可以計(jì)算VD3和25-OH-D3的斜率及其相應(yīng)的標(biāo)準(zhǔn)誤.但是不能計(jì)算相對(duì)生物學(xué)效價(jià)RBV(1.99)對(duì)應(yīng)的標(biāo)準(zhǔn)誤.
(1)SAS程序
Data nlin;
Input xs xt x0 y@@;
Output;
Cards;
0 0 1 0.204
0 0 1 0.131
0 0 1 0.115
0 0 1 0.206
0 0 1 0.170
1.484 0 0 0.203
1.773 0 0 0.216
1.553 0 0 0.217
1.379 0 0 0.205
1.576 0 0 0.167
3.393 0 0 0.267
3.547 0 0 0.301
3.515 0 0 0.267
3.990 0 0 0.300
3.290 0 0 0.320
9.244 0 0 0.367
9.510 0 0 0.416
9.321 0 0 0.411
9.638 0 0 0.437
8.890 0 0 0.353
19.222 0 0 0.466
18.960 0 0 0.400
18.720 0 0 0.487
19.096 0 0 0.458
19.480 0 0 0.457
0 0.720 0 0.225
0 0.687 0 0.227
0 0.676 0 0.275
0 0.782 0 0.201
0 0.791 0 0.206
0 1.623 0 0.268
0 1.838 0 0.284
0 1.920 0 0.298
0 1.885 0 0.280
0 1.508 0 0.318
0 4.170 0 0.411
0 4.495 0 0.446
0 4.180 0 0.371
0 4.345 0 0.318
0 4.470 0 0.342
0 8.660 0 0.443
0 9.000 0 0.499
0 9.102 0 0.432
0 9.730 0 0.463
0 9.840 0 0.414
Proc nlin;
Parms a=1 a0=1 bs=1 bts=1;
Der.a=1;
Der.a0=x0;
Der.bs=bts*xt+xs;
Der.bts=bs*xt;
Model y=a+a0*x0+bs*(bts*xt+xs);
Run;
(2)輸出結(jié)果
Parameter Estimate Approx Std Error Approximate 95% Confidence Limits
a 0.2314 0.0110 0.2092 0.2535
a0 -0.0662 0.0225 -0.1116 -0.0207
bs 0.0128 0.0012 0.0103 0.0152
bts 1.9884 0.2047 1.5750 2.4019
結(jié)果顯示:以VD3為參照物,用21 d肉雞股骨灰分重量評(píng)定25-OH-D3相對(duì)生物學(xué)效價(jià)RBV=bts=1.9884=1.99,標(biāo)準(zhǔn)誤SE=0.2047=0.20.數(shù)據(jù)表明:通過(guò)兩種維生素D斜率計(jì)算所得效價(jià)(RBV=0.0277/0.0139=1.99)與直接求解所得效價(jià)(RBV=bts=1.99)結(jié)果相同(表3).
表3 以股骨灰分重量評(píng)估1~21日齡肉雞飼糧中25-OH-D3相對(duì)VD3的生物學(xué)效價(jià)
經(jīng)過(guò)多項(xiàng)式比較、求解斜率及其標(biāo)準(zhǔn)誤、求解相對(duì)生物學(xué)效價(jià)及其標(biāo)準(zhǔn)誤等步驟,計(jì)算25-OH-D3相對(duì)于VD3的生物學(xué)效價(jià),可以得出結(jié)論:在1~21 d肉雞飼糧中,用股骨灰分重量作為評(píng)價(jià)指標(biāo),25-OH-D3生物學(xué)效價(jià)約為VD3的1.99倍.