Sử dụng phương pháp quang phổ hấp phụ cận hồng ngoại (NIRS) để xác định thành phần hóa học của khô dầu đỗ tương

ĐINH VĂN MƯỜI – Sử dụng phưỡng pháp quang phổ hấp phụ cận hồng ngoại .... 1 SỬ DỤNG PHƯƠNG PHÁP QUANG PHỔ HẤP PHỤ CẬN HỒNG NGOẠI (NIRS) ĐỂ XÁC ĐỊNH THÀNH PHẦN HÓA HỌC CỦA KHÔ DẦU ĐỖ TƯƠNG Đinh Văn Mười 1,Vũ Chí Cương2*, Nguyễn Đức Chuyên2, Nguyễn Sức Mạnh2 và Bùi Thị Thu Hiền2. 1Viên Chăn Nuôi 2Trung tâm nghiên cứu và phát triển chăn nuôi miền Núi *Tác giả liên hệ: Vũ Chí Cương - Viện Chăn nuôi - Từ Liêm - Hà Nội Tel: (04) 38.386.127/ 0912.121.506; Fax: (04) 38.389.775; Email: vuchicuong@gmail.com ABSTRACT The chemical composition of soybean cake predicted by Near infrared reflectance spectroscopy Near infrared reflectance spectroscopy (NIRS) was used to predict the chemical composition of soybean cake (dry matter - DM, crude protein -CP, crude fibre - CF and Fat). Samples were analysed by reference methods and spectra were collected using a NIRS machine in a range of wave lengths of 1100–2500 nm. A linear regression model was used to establish the relashionship of two sets of data (one from the chemical composition and one from NIRS). It was revealed that the NIRS based calibration equations can be accurately predicted chemical compositions of soybean cake. The following equations can be used for determination of DM (CK-PT), CP (Pr-PT); Fat (Mỡ- PT), CF (Xơ-PTL) of soybean cake with the accuracy of more than 95%: CK-PT = 2,627 + 0,9710 CK-NIRS; Pr-PT= 0,770 + 0,9825 Pr-NIRS; Mỡ-PT = 0,0717 + 0,9582 Mỡ-NIRS; Xơ-PTL = - 0,0301 + 1,003 Xơ-NIRS Key words: NIRS, calibration, prediction, equation. ĐẶT VẤN ĐỀ Để xác định thành phần hoá học của thức ăn, giúp xây dựng được khẩu phần ăn hợp lí đáp ứng đủ nhu cầu dinh dưỡng cho gia súc, gia cầm, người ta thường sử dụng phương pháp phân tích hóa học trong phòng thí nghiệm. Phương pháp này tốn nhiều thời gian và kinh phí. Hiện nay một phương pháp thay thế đang được sử dụng để xác định thành phần hoá học của thức ăn là phương pháp sử dụng quang phổ hấp phụ cận hồng ngoại (Near infrared reflectance spectroscopy – NIRS). Dùng NIRS để xác định thành phần hoá học của thức ăn hiện đã được công nhận là phương pháp phòng thí nghiệm có độ chính xác cao (Boval và cs, 2004), được AOAC chính thức công nhận để ước tính protein thô, ADF (AOAC 989.03), ẩm độ (AOAC 991.01; Barton và Windham, 1988), tinh bột, đường polysaccharides không phải tinh bột, mỡ, dầu, năng lượng trao đổi, tồn dư thuốc bảo vệ thực vật, độc tố trong ngũ cốc (Wrigley,1999) chất khô ở các loại cỏ làm thức ăn gia súc (Murray, 1993). NIRS còn được dùng để kiểm tra các loại thực phẩm (De Boever và cs, 1987), protein bị nhiệt làm biến tính, tồn dư nấm mốc và các chất phụ gia trộn trong nguyên liệu (Givens và Deaville, 1999). NIRS có những ưu điểm nổi bật hơn hẳn so với các phương pháp truyền thống khác: nhanh (45 giây/1 mẫu), nhiều mẫu đồng thời, chuẩn bị mẫu đơn giản (sấy khô và nghiền nhỏ, hoặc tươi), tránh được ô nhiễm môi trường do hoá chất, tránh được độc hại, không cần nhiều lao động, giá thành phân tích rẻ. VIỆN CHĂN NUÔI - Tạp chí Khoa học Công nghệ Chăn nuôi - Số 17-Tháng 4-2009 2 Xuất phát từ các ứng dụng trên của NIRS chúng tôi tiến hành đề tài: "Sử dụng kỹ thuật quang phổ hấp phụ cận hồng ngoại (NIRS) để xác định thành phần hoá học của khô dầu đỗ tương." với mục tiêu: Xây dựng các phương trình chẩn đoán thành phần hoá học của của khô dầu đỗ tương từ các số liệu phổ hấp phụ trên máy NIRS và số liệu phân tích thông thường. VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU Vật liệu nghiên cứu Khô dầu đỗ tương đã xác định thành phần hóa học tại Phòng phân tích Viện Chăn nuôi Địa điểm nghiên cứu Bộ môn dinh dưỡng thức ăn chăn nuôi và đồng cỏ, Trung tâm Thực nghiệm và Bảo tồn vật nuôi Thụy Phương, Phòng phân tích Viện Chăn nuôi. Thời gian nghiên cứu Đề tài được tiến hành từ tháng 5/2005 đến tháng 12/2007 Nội dung nghiên cứu Xác định phổ hấp phụ cận hồng ngoại của khô dầu đỗ tương đã xác định thành phần hóa học. Xây dựng đường hồi qui hiệu chỉnh chẩn đoán thành phần hoá học của khô dầu đỗ tương từ các số liệu về phổ hấp phụ cận hồng ngoại và thành phần hóa học. Áp dụng phương trình hồi qui tìm được cho các mẫu khô dầu đậu tương lấy ngẫu nhiên, chạy phổ trên máy NIRS nhưng không dùng để xây dựng phương trình hồi qui hiệu chỉnh để kiểm tra độ tin cậy của phương trình. Phương pháp nghiên cứu Phổ hấp phụ cận hồng ngoại của thức ăn được xác định với máy NIRS loại NIR Systems 5000 Monochromator của hãng Foss, USA với bước sóng từ 1100 - 2500 nm. Phổ và thành phần hóa học sau đó được xử lý bằng phần mềm WinISI. Từ kết quả chạy trên máy NIRS và nhờ phần mềm Winisi chúng ta có được hai bộ số liệu về thành phần hóa học phân tích trong phòng thí nghiệm và thành phần hóa học chẩn đoán trên máy NIRS. Từ hai bộ số liệu này, sử dụng kỹ thuật hồi qui đa chiều bậc một để xây dựng được phương trình hồi qui hiệu chỉnh chẩn đoán thành phần hóa học của thức ăn. Đây sẽ là phương trình hồi qui dùng để kiểm tra các mẫu thức ăn khác để biết độ chính xác của phương trình vừa tạo ra. Sở dĩ phải dùng phương trình này vì các phương trình có sẵn trong máy là phương trình được tạo ra cho các thức ăn ôn đới nên độ chính xác không cao. Mô hình toán học để xây dựng phương trình chẩn đoán hiệu chỉnh với NIRS trên cơ sở thành phần hóa học của thức ăn phân tích trong phòng thí nghiệm và thành phần hóa học của thức ăn chẩn đoán bằng phương trình có sẵn trong máy NIRS là mô hình hồi quy tuyến tính đa chiều bậc một: Y = bo + b1x1 + b2x2 ++ bixi + bpxp Trong đó: Y là kết quả phân tích phòng thí nghiệm của chất A nào đó x là kết quả chẩn đoán với máy NIRS b là hệ số Các phương trình chẩn đoán hiệu chỉnh với NIRS được áp dụng cho các TĂ mà số liệu NIRS ĐINH VĂN MƯỜI – Sử dụng phưỡng pháp quang phổ hấp phụ cận hồng ngoại .... 3 của nó không dùng để lập các phương trình chẩn đoán hiệu chỉnh với NIRS để kiểm tra độ tin cậy. Hai bộ số liệu tính từ phương trình chẩn đoán hiệu chỉnh với NIRS và bộ số liệu phân tích phòng thí nghiệm được so sánh bằng T-student (Paired test). Chỉ các phương trình cho hai kết quả không sai khác về thống kê (P > 0,05) mới được chấp nhận để dùng. Xử lý số liệu Số liệu được tính toán xác định giá trị trung bình và phân tích phương sai ANOVA trên phần mềm Genstat, phiên bản Discovery Edition (VSN International, 2007). KẾT QUẢ VÀ THẢO LUẬN Khô dầu đỗ tương Kết quả chẩn đoán chất khô (CK), protein thô (Pr), xơ, mỡ của khô dầu đỗ tương Kết quả chẩn đoán CK, Pr, xơ, mỡ của khô dầu đỗ tương từ NIRS được trình bày ở Bảng 1, Bảng 2. Bảng 1. Kết quả chẩn đoán CK, Pr, mỡ, xơ của khô dầu đỗ tương trên máy NIRS Chỉ tiêu n. Mean SD Est. Min Est. Max SEC SECV 1-VR CK 83 89,4921 0,7968 87,1017 91,8824 0,3484 0,3969 0,7496 Pr 82 46,3801 1,4383 42,0653 50,6949 0,6911 0,7426 0,7384 Mỡ 83 1,2233 0,3047 0,3091 2,1374 0,1202 0,1286 0,8258 Xơ 85 6,3659 0,5086 4,8403 7,8916 0,2277 0,2489 0,7631 Kết quả này cho thấy, chẩn đoán CK, Pr, mỡ, xơ của khô dầu đỗ tương bằng NIRS có SEC (Standard error of calibration - sai số của phép đo) chấp nhận được (0,1202 đến 0,6911%), và SECV (sai số của phép đo sau khi đã hiệu chỉnh -the standard error for cross validation) cũng đủ nhỏ để không gây sai số quá lớn của phép đo (0,1286 đến 0,7426%). Bảng 2. So sánh kết quả phân tích CK, Pr, mỡ, xơ của khô dầu đỗ tương với kết quả của NIRS Chỉ tiêu n. Mean SE StDev Min Max PT* 88 89,359 0,179 1,681 75,732 91,52 NIRS 88 89,318 0,173 1,627 75,732 91,831 Ck Sai khác (%) 0,3764 PT 83 46,324 0,167 1,518 41,737 50,39 NIRS 83 46,315 0,153 1,393 40,976 50,133 Pr Sai khác (%) 1,1978 PT 85 1,2373 0,0358 0,3297 0,65 2,46 NIRS 85 1,2257 0,0303 0,2793 0,754 2,009 Mỡ Sai khác (%) 8,233 PT 85 6,3659 0,0552 0,5086 4,62 7,32 NIRS 85 6,3659 0,0497 0,4581 4,893 7,262 Xơ Sai khác (%) 2,727 Chú thích: PT*: kết quả phân tích Với SEC và SECV không lớn nên không có sự sai khác lớn về các giá trị: Mean ± SD, Min và Max của CK, Pr, mỡ, xơ của khô dầu đỗ tương phân tích và chạy trên máy NIRS. Sai số tính bằng % giữa kết quả phân tich phòng thí nghiệm (Lab hoặc PT) và kết quả trên NIRS không VIỆN CHĂN NUÔI - Tạp chí Khoa học Công nghệ Chăn nuôi - Số 17-Tháng 4-2009 4 lớn và tương ứng là: 0,3764; 1,1978; 8,233; 2,727% cho CK, Pr, mỡ, xơ của khô dầu đỗ tương. *Xây dựng phương trình chẩn đoán CK, Pr, mỡ, xơ của khô dầu đỗ tương Sử dụng thuật toán hồi qui trên MINITAB với hai bộ số liệu là kết quả vừa có trên máy NIRS và kết quả phân tích chúng tôi xây dựng được phương trình để chẩn đoán CK, Pr, mỡ, xơ của khô dầu đỗ tương trong Bảng 3. Bảng 3: Phương trình hồi quy chẩn đoán DM, CP, Fat, CF của khô dầu đỗ tương với NIRS Kết quả xây dựng phương trình chẩn đoán để chẩn đoán CK, Pr, mỡ, xơ của khô dầu đỗ tương từ NIRS cho thấy dù số mẫu chưa nhiều nhưng phương trình chẩn đoán CK, Pr, mỡ, xơ của khô dầu đỗ tương bằng NIRS với số mẫu > 80 có R2 khá cao, tương ứng là (92,0; 79,8; 68,3 và 80,9) cho CK, Pr, mỡ, xơ. Phương trình chẩn đoán CK, Pr, mỡ, xơ của khô dầu đỗ tương bằng NIRS với số mẫu = 65 có R2 khá, tương ứng là (81,0; 76,8; 62,3 và 80,6) cho CK, Pr, mỡ, xơ. Áp dụng phương trình hồi qui cho các mẫu khô dầu đỗ tương tương khác để kiểm tra độ tin cậy của phương trình chẩn đoán của khô dầu đỗ tương *Đối với chất khô Sau khi áp dụng phương trình có số thứ tự 5 trên 23 mẫu ngẫu nhiên chúng tôi có kết quả ở Bảng 4 và Đồ thị 1. Bảng 4 So sánh kết quả phân tích CK khô dầu đỗ tương trong phòng thí nghiệm và sử dụng phương trình hiệu chỉnh với NIRS Chỉ tiêu n Mean SE StDev Min Max R2 (%) r Lab 88,654 0,605 2,9 75,732 91,09 NIRS 88,639 0,57 2,734 76,163 89,824 Phương trình số 5 Sai khác (%) 23 0,551 93,8 0,97 Kết quả cho thấy, với 23 mẫu và chỉ với phương trình có dung lương mẫu là 65, kết quả chẩn đoán bằng phương trình hiệu chỉnh với NIRS cho CK của khô dầu đỗ tương đã rất tốt: R2 (%) là 93.8%, sai khác % giữa hai giá trị Mean tính được và thực của 23 mẫu là đủ nhỏ: 0,551 %. Hai bộ số liệu tính từ phương trình và số liệu phân tích 23 mẫu ngẫu nhiên có hồi qui tuyến tính chặt chẽ và đều có phân bố chuẩn và gần như trùng khít nhau. Như vậy, phương trình số 5 có thể dùng để ước tính CK của khô dầu đỗ tương. TT n Phương trình R2 (%) R P 1 88 CK-PT = 0,827 + 0,9912 CK-NIRS 92,0 0,959 < 0,01 2 83 Pr- PT = 1,191 + 0,9745 Pr –NIRS 79,8 0,895 < 0,01 3 85 Mỡ- PT b = 0,03841 + 0,9782 Mỡ-NIRS 68,3 0,829 < 0,01 4 85 Xơ- PT = - 0,0001 + 1,000 Xơ-NIRS 80,9 0,901 < 0,01 5 65 CK- PT = 2,627 + 0,9710 CK-NIRS 81,0 0,901 < 0,01 6 65 Pr- PT = 0,770 + 0,9825 Pr-NIRS 76,8 0,878 < 0,01 7 65 Mỡ- PT = 0,0717 + 0,9582 Mỡ-NIRS 62,3 0,793 < 0,01 8 65 Xơ- PT = - 0,0301 + 1,003 Xơ- NIRS 80,6 0,899 < 0,01 ĐINH VĂN MƯỜI – Sử dụng phưỡng pháp quang phổ hấp phụ cận hồng ngoại .... 5 DM_NIRS D M _L ab 90.087.585.082.580.077.575.0 92.5 90.0 87.5 85.0 82.5 80.0 77.5 75.0 S 0.723816 R-Sq 94.1% R-Sq(adj) 93.8% DM_Lab = - 2.553 + 1.029 DM_NIRS CP_NIRS CP _ La b 50494847464544434241 50 48 46 44 42 40 S 0.778508 R-Sq 86.0% R-Sq(adj) 85.1% CP_Lab = 0.913 + 0.9852 CP_NIRS Data Fr eq ue nc y 9590858075 9 8 7 6 5 4 3 2 1 0 Mean StDev N 88.65 2.900 23 88.64 2.734 23 Variable DM_lab DM_NIRS Histogram of DM_lab, DM_NIRS Normal Data Fr eq ue nc y 5048464442 4 3 2 1 0 Mean StDev N 46.42 2.020 18 46.19 1.902 18 Variable CP_Lab CP_NIRS Histogram of CP_Lab, CP_NIRS Normal Đồ thị 1: Hồi quy giữa giá trị chẩn đoán (NIRS) với giá trị phân tích phòng thí nghiệm của khô đỗ (% CK) và phân bố tần số của hai bộ số liệu *Đối với protein thô Sau khi áp dụng phương trình có số thứ tự 6 trên 18 mẫu ngẫu nhiên chúng tôi có kết quả ở Bảng 5 và Đồ thị 2. Bảng 5. So sánh kết quả phân tích Pr của khô dầu đỗ tương trong phòng thí nghiệm và sử dụng phương trình hiệu chỉnh chẩn đoán của NIRS Chỉ tiêu n Mean SE StDev Min Max R2 (%) r PT 46,424 0,476 2,02 41,737 50,39 NIRS 46,192 0,448 1,902 41,029 49,246 Phương trình số 6 Sai khác (%) 18 1,39 0,213 0,903 0,0302 2,77 85,1 0,927 Kết quả cho thấy, với 18 mẫu và chỉ với phương trình có dung lượng mẫu là 65, kết quả chẩn doán bằng phương trình hiệu chỉnh với NIRS cho Pr của khô dầu đỗ tương đã rất tốt: R2 (%) là 81,5 %, sai khác % giữa hai giá trị Mean tính được và thực của 18 mẫu là 1,39% là đủ nhỏ, nhỏ hơn 5% nên chấp nhận được để dùng. VIỆN CHĂN NUÔI - Tạp chí Khoa học Công nghệ Chăn nuôi - Số 17-Tháng 4-2009 6 Fat_NIRS Fa t_ La b 2.001.751.501.251.00 2.25 2.00 1.75 1.50 1.25 1.00 0.75 0.50 S 0.107790 R-Sq 88.5% R-Sq(adj) 87.8% Fat_Lab = - 0.3399 + 1.217 Fat_NIRS CF_NIRS CF _L ab 7.06.56.05.55.0 7.5 7.0 6.5 6.0 5.5 5.0 4.5 S 0.262890 R-Sq 82.2% R-Sq(adj) 81.2% CF_Lab = - 0.0304 + 1.014 CF_NIRS Data Fr eq ue nc y 2.22.01.81.61.41.21.00.8 7 6 5 4 3 2 1 0 Mean StDev N 1.355 0.3087 20 1.392 0.2385 20 Variable Fat_Lab Fat_NIRS Histogram of Fat_Lab, Fat_NIRS Normal Data Fr e qu en cy 7.57.06.56.05.55.04.5 8 7 6 5 4 3 2 1 0 Mean StDev N 6.303 0.6064 20 6.248 0.5423 20 Variable CF_Lab CF_NIRS Histogram of CF_Lab, CF_NIRS Normal Đồ thị 2: Hồi quy giữa giá trị chẩn đoán (NIRS) với giá trị phân tích phòng thí nghiệm của khô dầu đỗ tương (% Pr) và phân bố tần số của hai bộ số liệu Hai bộ số liệu tính từ phương trình và số liệu phân tích 18 mẫu ngẫu nhiên có hồi qui tuyến tính chặt chẽ và đều có phân bố chuẩn và gần như trùng khít nhau. Như vậy, phương trình số 6 có thể dùng để ước tính Protein của khô dầu đỗ tương. *Đối với mỡ và xơ Sau khi áp dụng hai phương trình có số thứ tự 7, 8 trên 20 mẫu ngẫu nhiên chúng tôi có kết quả ở Bảng 6 và các Đồ thị 3. Bảng 6. So sánh kết quả phân tích mỡ, xơ của khô dầu đỗ tương trong phòng thí nghiệm và sử dụng phương trình chẩn đoán của NIRS Chỉ tiêu n. Mean SE StDev Min Max R2 (%) r Mỡ-PT 20 1,3551 0,069 0,3087 0,77 2,24 Mỡ-NIRS 20 1,3924 0,0533 0,2385 0,8622 1,9967 Phương trình số 7 Sai khác (%) 20 7,28 1,32 5,91 0,407 18,43 87,8 0,94 Xơ-PT 20 6,303 0,136 0,606 4,62 7,32 Xơ-NIRS 20 6,248 0,121 0,542 4,961 7,119 Phương trình số 8 Sai khác (%) 20 3,38 0,56 2,504 0,0849 7,874 81,2 0,907 Kết quả cho thấy, với 20 mẫu và chỉ với phương trình có dung lương mẫu là 65, kết quả chẩn đoán bằng phương trình hiệu chỉnh với NIRS cho mỡ và xơ của khô dầu đỗ tương đã rất tốt: R2 (%) là 87,8 và 81,2. Tuy nhiên sai khác % giữa hai giá trị Mean tính được và thực của 20 mẫu là 3.38 % cho xơ là chấp nhận được, nhưng sai khác này với mỡ là lớn hơn 5 % (7.28%). ĐINH VĂN MƯỜI – Sử dụng phưỡng pháp quang phổ hấp phụ cận hồng ngoại .... 7 Đồ thị 3. Hồi qui giữa giá trị chẩn đoán (NIRS) với giá trị phân tích phòng thí nghiệm (% mỡ, xơ) và phân bố tần số của bốn bộ số liệu Hai bộ số liệu tính từ phương trình và số liệu phân tích 20 mẫu ngẫu nhiên cho xơ có hồi qui tuyến tính chặt chẽ và đều có phân bố chuẩn và gần như trùng khít nhau. Hai bộ số liệu tính từ phương trình và số liệu phân tích 20 mẫu ngẫu nhiên cho mỡ có hồi qui tuyến tính chặt chẽ và đều có phân bố chuẩn và gần như trùng khít nhau. Như vậy, phương trình số 8 có thể dùng để ước tính xơ của khô dầu đỗ tương, phương trình 7 có thể tạm dùng để ước tính mỡ của khô dầu đỗ tương. Kiểm tra sự sai khác của các số trung bình Để kiểm tra lại cả bốn phương trình chúng tôi đã dùng T-student (Paired test) để so sánh từng cặp số liệu về CK, Pr, mỡ, xơ của khô dầu đỗ tương. Bộ số liệu thứ nhất là bộ số liệu tính được từ phương trình chẩn đoán, bộ số liệu kia là kết quả phân tích. Kết quả so sánh T-student (Paired test) Chỉ tiêu n Mean SE StDev Minimum Maximum P Lab 23 88,654 0,605 2,9 75,732 91,09 CK NIRS 23 88,639 0,57 2,734 76,163 89,824 0,985 Lab 18 46,424 0,476 2,02 41,737 50,39 Pr NIRS 18 46,192 0,448 1,902 41,029 49,246 0,725 Lab 20 1,3551 0,069 0,3087 0,77 2,24 Mỡ NIRS 20 1,3924 0,0533 0,2385 0,8622 1,9967 0,671 Lab 20 6,303 0,136 0,606 4,62 7,32 Xơ NIRS 20 6,248 0,121 0,542 4,961 7,119 0,763 Kết quả của phép thử T-student cho thấy: các giá trị (P>0,05) rất nhiều, có nghĩa là các giá trị CK, Pr, mỡ và xơ của khô dầu đỗ tương tính theo các phương trình và giá trị phân tích cho CK, Pr, mỡ và xơ là không khác nhau có ý nghĩa thống kê hay đúng hơn chúng tương tự nhau. Như vậy, có thể dùng các phương trình trên NIRS để chẩn đoán CK, Pr, mỡ và xơ của khô dầu đỗ tương với độ chính xác > 95%. KẾT LUẬN VÀ ĐỀ NGHỊ Kết luận Có thể dùng các phương trình sau để xác định chất khô, Protein, mỡ và xơ của khô dầu đỗ tương với độ chính xác xấp xỉ 95%. CK-PT = 2,627 + 0,9710 CK-NIRS; Pr - PT= 0,770 + 0,9825 Pr-NIRS; Mỡ-PT = 0,0717 + 0,9582 Mỡ-NIRS; Xơ - PTL = - 0,0301 + 1,003 Xơ-NIRS Đề nghị Cho áp dụng kết quả nghiên cứu để xác định thành phần hóa học của thức ăn để giảm chi phí phân tích và làm thí nghiệm in vivo. Tiếp tục nghiên cứu theo hướng này cho thức ăn khác và tăng độ chính xác của các phương trình hiện có, đặc biệt là các phương trình cho mỡ, chất khô và khoáng. VIỆN CHĂN NUÔI - Tạp chí Khoa học Công nghệ Chăn nuôi - Số 17-Tháng 4-2009 8 TÀI LIỆU THAM KHẢO Barton. F.E., II and Windham. W.R, (1988). Determination of acid detergent fibre and crude protein in forages by near infrared reflectance spectroscopy: collaborative study. Journal of the Association of official Analytical Chemists, 71: pp. 1162-1167. Boval.M., Coates. D.B., Lacomte.P., Decruyenaere.V and Archimede.H, (2004). Faecal Near-Infrared Reflectance Spectroscopy (NIRS) to access chemical composition, in vivo digestibility and intake of tropical grass by Creolo cattle. Animal Feed Science and Technology, Vol 114, Issues 1-4, 3 May 2004, Pp. 19-29. De Boever. J. L,.Cottyn. B. G; De Brander. D. L and Buysse. F. S, (1987). Traitement de la paille.1. Effet de l’ammoniac sur la compostion, la digestibilité et la valeur alimentaire. Revue de l'Agric, 40: pp.347. Given. D. I and Deaville, E. R, (1999). The current and future role of near infrared reflectance spectroscopy in animal nutrition: a riview. Aust. J. Agric. Res, 1999, 50. pp.1131-1145. Murray. I, (1993). Forage analysis by near infrared reflectance spectroscopy. In A. Davies, R.D. Baker, S.A. Grant & A.S. Laidlaw, eds. Sward management handbook. British Grassland Society, Reading, UK. pp. 285-312 Wrigley.C.W, (1999). Potential methodologies and strategies for the rapid assessment of feed-grain quality. Australian Journal of Agricultural Research, 50: pp.789-805. *Người phản biện: TS. Đoàn Thị Khang ; TS. Mai Văn Sánh