評(píng)定飼料營(yíng)養(yǎng)物質(zhì)相對(duì)生物學(xué)效價(jià)的SAS程序

韓進(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

1 多項(xiàng)式比較

用GLM中Contrast方法確定不同維生素D水平與肉雞股骨灰分重量的線性關(guān)系.

1.1 VD3水平與肉雞股骨灰分重量的關(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線性影響肉雞股骨灰分重量(P<0.0001)，二次曲線影響差異不顯著(P=0.2911).數(shù)據(jù)表明：1～21 d肉雞VD3食入量與股骨灰分重量之間存在顯著線性關(guān)系.結(jié)果列于表2.

1.2 25-OH-D3水平與肉雞股骨灰分重量的關(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.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線性影響肉雞股骨灰分重量(P<0.0001))，二次曲線影響差異不顯著(P=0.4489).數(shù)據(jù)表明：1～21 d肉雞25-OH-D3食入量與股骨灰分重量之間存在顯著線性關(guān)系.結(jié)果列于表2.

表2 VD3和25-OH-D3與21日齡肉雞股骨灰分重量的關(guān)系

續(xù)表2

2 求解斜率及其標(biāo)準(zhǔn)誤

(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)誤.

3 求解相對(duì)生物學(xué)效價(jià)及其標(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倍.