PHÂN TÍCH hàm LƯỢNG cà PHÊ và CAFFEINE TRONG cà PHÊ RANG XAY BẰNG PHƯƠNG PHÁP UV VIS kết hợp PHÂN TÍCH dữ LIỆU đa BIẾN

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA HÓA HỌC BM HĨA PHÂN TÍCH PHẠM VÕ PHƯƠNG DUN PHÂN TÍCH HÀM LƯỢNG CÀ PHÊ VÀ CAFFEINE TRONG CÀ PHÊ RANG XAY BẰNG PHƯƠNG PHÁP UV-VIS KẾT HỢP PHÂN TÍCH DỮ LIỆU ĐA BIẾN KHĨA LUẬN TỐT NGHIỆP THÀNH PHỐ HỒ CHÍ MINH – 2019 ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA HĨA HỌC BM HĨA PHÂN TÍCH PHẠM VÕ PHƯƠNG DUN PHÂN TÍCH HÀM LƯỢNG CÀ PHÊ VÀ CAFFEINE TRONG CÀ PHÊ RANG XAY BẰNG PHƯƠNG PHÁP UV-VIS KẾT HỢP PHÂN TÍCH DỮ LIỆU ĐA BIẾN KHĨA LUẬN TỐT NGHIỆP GVHD: PGS.TS Nguyễn Ánh Mai HVCH Nguyễn Phúc Thịnh CN Lê Văn Duy THÀNH PHỐ HỒ CHÍ MINH – 2019 LỜI CẢM ƠN Qua năm học tập rèn luyện trường Đại học Khoa học Tự nhiên TPHCM, em gửi lời cảm ơn đến tất thầy cô khoa Hóa học nói chung mơn Hóa phân tích nói riêng kiến thức kinh nghiệm mà thầy cô truyền đạt Em xin gửi lời cảm ơn chân thành đến cô PGS.TS Nguyễn Ánh Mai, nhà giáo tận tụy, ln hết lịng sinh viên Gửi lời cảm ơn đặc biệt đến chị Nguyễn Phúc Thịnh, anh Lê Văn Duy người trực tiếp hướng dẫn, giúp em giải thắc mắc khó khăn q trình làm đề tài viết báo cáo Cảm ơn thầy Nguyễn Hoàng Phú hỗ trợ em mặt dụng cụ, hóa chất Cảm ơn chị Nguyễn Thị Kim Hạnh người sáng lập công ty The Yellow Chair Specialty Coffee tư vấn cho em kiến thức cà phê cung cấp số mẫu cà phê Cảm ơn cô Hồ Thị Phước nhận lời phản biện đưa góp ý chân thành khóa luận em Cảm ơn bạn sinh viên chun ngành hóa phân tích, hỗ trợ trình thực đề tài Cuối cùng, vơ biết ơn gia đình quan tâm, động viên, chỗ dựa tinh thần cững để có điều kiện học tập cách tốt Cảm ơn tất người! Sinh viên Phạm Võ Phương Duyên GIỚI THIỆU CHUNG Cà phê thức uống tiêu dùng phổ biến giới Việt Nam quốc gia xuất cà phê lớn thứ giới sau Brazil, thành viên Hiệp Hội Cà Phê Thế Giới (ICO) Cà phê mặt hàng nông sản chủ lực Việt Nam, đứng thứ hai (sau mặt hàng gạo) kim ngạch xuất Chính ngành cà phê có vai trị lớn kinh tế nước nhà Theo lịch sử hình thành, cà phê chia thành nhiều loại Ngày nay, tính giá trị thương mại điều kiện ni trồng, hai loại cà phê trồng nhiều cà phê Arabica cà phê Robusta Từ cà phê đến hạt cà phê rang trải qua nhiều giai đoạn ảnh hưởng lớn đến chất lượng thành phần hóa học cà phê Thị trường cà phê Việt Nam thường Robusta đa phần người Việt Nam thích vị đắng chúng, bên cạnh đó, giá thành Robusta thấp Arabica (thường thấp gấp hai lần) Tuy nhiên, thông thường để tăng hương vị đậm đà cho cà phê, hầu hết nhà sản xuất cà phê rang thường cho thêm số chất phụ gia vào trình rang Việc thêm vào cà phê chất phụ gia xuất phát từ nhu cầu thị trường nhiên số biến tướng việc xảy người sản xuất muốn hạ giá thành sản phẩm để tăng lợi nhuận Việc làm xuất tình trạng cà phê bẩn xuất thị trường Do nhu cầu tìm kiếm giải pháp để phân biệt cà phê với thành phần làm giả đời Trong đề tài này, phương pháp phân tích trắc quang kết hợp với phân tích liệu đa biến thực để xác định hàm lượng cà phê hỗn hợp trộn lẫn cà phê, bắp đậu nành Đồng thời phương pháp sử dụng để xác định hàm lượng caffeine mẫu cà phê, hướng đến việc xác định nhanh hàm lượng caffeine mẫu thay cho phương pháp sắc ký lỏng truyền thống MỤC LỤC DANH MỤC HÌNH ẢNH A DANH MỤC BẢNG BIỂU B DANH MỤC VIẾT TẮT C TỔNG QUAN 1 Thành phần hóa học cà phê .1 1.1 Các hợp chất tạo vị đắng – Caffeine Trigonelline 1.2 Các hợp chất tạo vị béo – Lipid .2 1.3 Các hợp chất tạo vị – Carbohydrate phản ứng Maillard 1.4 Các hợp chất tạo vị chua – Các acid hữu Phân loại cà phê 2.1 Phân biệt theo cảm quan 2.2 Phân biệt theo thành phần hóa học Các phương pháp phân tích cà phê 3.1 Phương pháp sắc kí lỏng hiệu cao HPLC-UV (theo TCVN 9723:2013)12 3.2 Phương pháp sắc ký khí đầu dị ion hóa lửa GC-FID13 3.3 Phương pháp NIR kết hợp với chemometrics11 Giới thiệu phương pháp phân tích liệu đa biến 4.1.Phương pháp thành phần PCA (Principal Component Analysis) 4.2 Phương pháp bình phương tối thiểu phần PLS (Projection to Latent Structures - Partial Least Square) 11 4.3 Các cơng cụ mơ tả mơ hình liệu đa biến 12 Hiện trạng cà phê Việt Nam 15 THỰC NGHIỆM 17 Dụng cụ hóa chất 17 1.1 Thiết bị, dụng cụ phần mềm .17 i 1.2 Hóa chất 17 1.3 Lấy mẫu, xử lí sơ bảo quản mẫu 17 Quy trình chiết mẫu 18 2.1.Quy trình chiết mẫu để đo trắc quang 18 2.2.Quy trình xử lý mẫu cho việc xác định caffeine phương pháp HPLC – UV 19 Quy trình thiết lập mơ hình PCA PLS 20 KẾT QUẢ VÀ BIỆN LUẬN 22 Mở rộng mơ hình xác định hàm lượng cà phê phương pháp trắc quang kết hợp với liệu đa biến 22 1.1.Mơ hình phân loại mẫu đơn - mơ hình thành phần PCA 22 1.2.Phương pháp định lượng - Mơ hình bình phương tối thiểu phần PLS 24 Thiết lập mơ hình xác định hàm lượng caffeine phương pháp trắc quang kết hợp với liệu đa biến 32 2.1.So sánh phương pháp chiết quy trình xác định hàm lượng caffeine phương pháp HPLC-UV 32 2.2.Xác định hàm lượng caffeine phương pháp HPLC-UV 33 2.3.Xây dựng mơ hình PLS 34 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 40 Kết luận 40 Hướng phát triển 40 PHỤ LỤC D TÀI LIỆU THAM KHẢO E ii DANH MỤC HÌNH ẢNH Hình Công thức cấu tạo caffeine Hình Cơng thức cấu tạo trigonelline (trái) and nicotinic acid (phải) Hình Cơng thức cấu tạo triacylglycerol Hình Sự hình thành phân tử Melanoidin trình rang8 Hình Phản ứng phân hủy CGA trình rang Hình Phản ứng hình thành Chlorogenic lactone từ CGA Hình Hình dạng hạt Arabica Robusta Hình Quy trình phân tích caffeine phương pháp HPLC-UV Hình Ba loại vấn đề phân tích liệu khác giải MVDA Hình 10 Mơ phương pháp chiếu mơ hình PCA 10 Hình 11 Mô cân trung điểm 10 Hình 12 Hai giai đoạn mơ hình PLS 12 Hình 13 Sự đánh đổi lực giải thích R2 lực dự đốn Q2 12 Hình 14 Mơ tượng Overfitting 12 Hình 15 Quy trình thiết kế mơ hình MVDA 15 Hình 16 Quy trình xử lý mẫu phương pháp trắc quang 18 Hình 17 Chương trình gradient dung mơi 20 Hình 18 Quy trình xử lý mẫu theo TCVN 9723:2013 20 Hình 19 Biểu đồ “Score” “Loading” mơ hình PCA chồng lên 22 Hình 20 Phân bố phổ hấp thu mẫu cà phê 23 Hình 21 Phổ hấp thu mẫu bắp đậu nành 24 Hình 22 Sự phân bố phổ hấp thu mẫu cà phê chọn 25 Hình 23 Đồ thị Hotelling’s T2Range 26 Hình 24 Đồ thị DModX 26 Hình 25 Biểu đồ “Overview” mơ hình 27 Hình 26 Đồ thị phân phối chuẩn (µ,σ)15 29 Hình 27 Sắc ký đồ mẫu cà phê chiết theo phương pháp 32 Hình 28 Đường chuẩn caffeine xây dựng từ phương pháp HPLC – UV 34 Hình 29 Giá trị VIP theo số sóng 34 Hình 30 Đồ thị Hotelling’s T2Range 35 Hình 31 Phổ hấp thu điểm mẫu tính chất với "19R" "31R" 36 Hình 32 Biều đồ “Overview” mơ hình xác định hàm lượng caffeine 36 A DANH MỤC BẢNG BIỂU Bảng Thành phần hóa học hạt cà phê Arabica Robusta xanh2 Bảng Thành phần caffeine chlorogenic acid hạt cà phê rang10 Bảng Một số cơng cụ mơ tả cho mơ hình MVDA 13 Bảng Thông số thiết bị UV-VIS hai chùm tia 18 Bảng Các điều kiện tối ưu quy trình chiết mẫu cà phê 19 Bảng So sánh phân phối xác xuất tổng phổ cà phê phân phối chuẩn 23 Bảng Bảng giá trị RMSEE RMSEP mơ hình 27 Bảng So sánh giá trị RMSEE RMSEP hai mơ hình 28 Bảng Kết dự đoán hàm lượng cà phê kiểm tra 29 Bảng 10 Kết mơ hình tính hàm lượng cà phê mẫu thị trường 31 Bảng 11 Hàm lượng phần trăm caffeine theo phương pháp chiết 33 Bảng 12 Khả dự đốn mơ hình 35 Bảng 13 R2Y Q2 mơ hình 36 Bảng 14 Bộ kiểm tra mơ hình xác định hàm lượng caffeine 37 Bảng 15 Kết kiểm tra mẫu thi trường mơ hình xác định hàm lượng caffeine 38 B DANH MỤC VIẾT TẮT Kí hiệu Tên đầy đủ Ý nghĩa Chlorogenic acid, Caffeine Các hợp chất có cà phê Distance to the Model in X Khoảng cách đến mơ hình space khơng gian tọa độ biến X MVDA Multivariate Data Analysis Phân tích liệu đa biến NIR Near Infrared Spectrometer Phổ hồng ngoại trường gần PC Principal Component Biến ngầm (Latent Variable) Principal Component Analysis Phân tích thành phần CGA, CAF DmodX PCA PLS Pred RMSEE RSD RMSEP Var VIP Projections to Latent Structures Phân tích hồi quy tuyến tính (Partial Least Square) phần Prediction Giá trị dự đoán Root Mean Square Error of Căn bình phương trung bình sai số Estimation ước tính Relative Standard Deviation Độ lệch chuẩn tương đối Root Mean Square Error of Căn bình phương trung bình sai số Prediction dự đoán Variable Giá trị thiết kế Variable Influence of Thông số biến ảnh hưởng đến mô Projections hình C 50-80 % >80 % 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 30.2 30.8 31.9 34.1 34.9 35.0 35.3 41.4 48.5 48.6 49.0 49.5 49.6 49.8 49.9 50.0 50.9 56.9 65.1 66.9 67.0 68.3 69.9 70.0 70.5 72.7 76.4 77.9 78.3 79.0 79.9 80.6 80.8 89.2 89.2 89.3 90.1 90.4 90.9 100.0 100.0 100.0 28.7 29.4 29.3 32.2 35.0 37.9 35.9 37.1 50.4 46.7 51.3 48.2 47.3 55.5 48.7 47.3 45.4 53.1 68.8 67.1 60.7 67.5 68.0 80.2 66.2 76.3 80.7 66.1 83.4 78.2 80.4 78.0 72.6 70.4 88.5 86.8 96.5 83.6 84.4 95.1 97.9 92.3 5.1 4.5 8.2 5.6 0.4 8.2 1.7 10.3 3.9 4.0 4.6 2.6 4.8 11.5 2.5 5.4 10.9 6.7 5.7 0.4 9.4 1.2 2.8 14.5 6.0 5.0 5.6 15.2 6.5 0.9 0.6 3.2 10.2 21.0 0.9 2.9 7.1 7.6 7.1 4.9 2.1 7.8 30 61 62 63 64 65 100.0 100.0 100.0 100.0 100.0 99.3 97.6 109.4 98.8 106.4 0.8 2.4 9.4 1.2 6.4 Kết xác định hàm lượng cà phê mua thị trường 1.2.2.2 Một số mẫu cà phê có nhãn hiệu thu thập thị trường kiểm tra hàm lượng cà phê mơ hình cho kết Bảng 10 Bảng 10 Kết mơ hình tính hàm lượng cà phê mẫu thị trường STT Tên mẫu Hãng sản xuất Hàm lượng cà phê dự đốn từ mơ hình PLS (%) Thunder No 1864 Café 105.0 ± 2.2 Thunder Culi Robusta 1864 Café 93.9 ± 1.7 Sơn Tùng Robusta Sơn Tùng 73.4 ± 2.5 Highlands Di sản Highlands Coffee 106.9 ± 7.3 Chinh phục S Trung Nguyên 49.4 ± 3.8 Sáng tạo Culi Robusta Trung Nguyên 38.1 ± 3.3 Nes Việt nguyên chất NESCAFÉ 117.9 ± 3.6 Nes Việt nguyên chất NESCAFÉ 117.1 ± 2.1 Đức Thành Đức Thành 93.9 ± 3.1 10 Cà phê Chocolate Phúc Long 111.8 ± 2.4 11 Cà phê Buôn Mê Thuột Đ Phương Vy 40.1 ± 1.5 12 Phương Vy Robusta Phương Vy 34.61 ± 0.92 13 Bùi Văn Thọ Robusta Bùi Văn Ngọ 103.3 ± 2.0 14 Bùi Văn Thọ Culi Bùi Văn Ngọ 107.1 ± 2.2 15 Robusta Thượng Hạng Light Coffee 106.8 ± 2.4 16 Robusta Đặc Biệt Light Coffee 98.0 ± 3.8 31 Thiết lập mơ hình xác định hàm lượng caffeine phương pháp trắc quang kết hợp với liệu đa biến 2.1 So sánh phương pháp chiết quy trình xác định hàm lượng caffeine phương pháp HPLC-UV Trong khuôn khổ đề tài phương pháp chiết caffeine từ mẫu cà phê sử dụng dung môi EtOH:H2O sử dụng để xây dựng chuẩn, để có sở khẳng định phương pháp hiệu cần tiến hành so sánh hiệu suất chiết với phương pháp xử lý mẫu trọng tài phương pháp xử lý mẫu theo TCVN 9723:2013 Quy trình xử lý mẫu trình bày mục 6.2.2 Theo TCVN 9723:2013 quy trình chiết mẫu sử dụng MgO để loại hợp chất tạp tannin có mẫu chiết MgO trước sử dụng nung lò gia nhiệt lên 600oC giữ Hai phương pháp chiết mẫu thực dựa mẫu có hàm lượng caffeine thấp, trung bình cao cho kết sau: Hình 27 Sắc ký đồ mẫu cà phê chiết theo phương pháp Hình 27 thể so sánh mẫu cà phê xử lý phương pháp chiết với EtOH: H2O (màu hồng) xử lý theo TCVN 9723:2013 (màu đen) Phương pháp xử 32 lý mẫu theo TCVN loại bỏ mẫu tốt hơn, phương pháp sử dụng nước nóng chiết 90 oC nên khó kiểm sốt q trình Bảng 11 Hàm lượng phần trăm caffeine theo phương pháp chiết Chiết với EtOH:H2O Chiết theo TCVN Sai lệch tương đối Hàm lượng CAF RSD Hàm lượng CAF RSD (%) (%) (%) (%) 7R 1.945 ± 0.026 1.3 1.750 ± 0.077 4.4 11.1 31R 1.433 ± 0.071 4.9 1.270 ± 0.059 4.7 12.8 R1 2.193 ± 0.056 2.6 2.000 ± 0.039 1.9 9.7 Mẫu (%) Kết hàm lượng caffeine từ phương pháp chiết theo TCVN cho giá trị RSD lớn (~5 %) phương pháp chiết với EtOH:H2O (1-5 %) Hàm lượng caffeine chiết theo TCVN thấp so với phương pháp chiết hỗn hợp ethanol nước (thấp ~10 %), điều cho thấy phương pháp chiết EtOH:H2O cho hiệu suất chiết tốt Vì mục tiêu nghiên cứu sử dụng dịch chiết EtOH:H2O để xác định caffeine nên khảo sát dùng cho mục đích chuyển đổi kết phương pháp Nói cách khác muốn ước lượng hàm lượng caffeine xác định TCVN lấy hàm lượng caffeine xác định phương pháp cộng với 10 % giá trị Việc xây dựng mơ hình dự đoán hàm lượng caffeine thực dịch chiết EtOH:H2O, loại dịch chiết để dự đoán thành phần cà phê cà phê rang xay 2.2 Xác định hàm lượng caffeine phương pháp HPLC-UV Đường chuẩn caffeine dùng phương pháp HPLC-UV có độ tuyến tính tốt với R2 = 0.9993 Giá trị giới hạn phát hiên, LOD, giới hạn định lượng, LOQ, phương pháp 0.016 µg.mL-1 0.048 µg.mL-1 33 1.2E+06 Diện tích peak 1.0E+06 y = 101338x + 7134.7 R² = 0.9993 8.0E+05 6.0E+05 4.0E+05 2.0E+05 0.0E+00 Nồng độ (mg.L-1) 10 12 Hình 28 Đường chuẩn caffeine xây dựng từ phương pháp HPLC – UV 2.3 Xây dựng mơ hình PLS Tương tự q trình xây dựng mơ hình định lượng phần trăm cà phê, mơ hình định lượng hàm lượng caffeine hàm lượng caffeine biến đáp ứng Các mẫu cà phê phân tích hàm lượng caffeine phương pháp HPLC – UV, sau mẫu hỗn hợp trộn cà phê, bắp đậu nành tính tốn hàm lượng caffeine để đưa vào mơ hình Tổng cộng mơ hình chứa 193 điểm mẫu chuẩn 2.3.1 Chọn khoảng bước sóng thích hợp Để loại bỏ biến khơng tực đóng góp đáng kể vào việc xây dựng mơ hình ta tiến hành lựa chọn số sóng theo thông số VIP Từ liệu ban đầu, ta tiến hành chọn lọc số sóng theo thơng số VIP Những biến có giá trị VIP > thường ảnh hưởng quan trọng đến mơ hình Do ta tiến hành loại bỏ thơng số sóng có giá trị < 1, từ lập mơ hình thống kê Hình 29 Giá trị VIP theo số sóng 34 Các số sóng có VIP > nằm chủ yếu vùng 250 – 280 nm (vùng phổ hấp thu caffeine) vùng 300-350 nm (vùng phổ hấp thu acid hữu cơ) Các vùng có thơng số sóng VIP > vùng hấp thu đặc trưng hợp chất cà phê nên việc xây dựng mơ hình theo thơng số VIP hợp lý (Hình 29) Bảng 12 Khả dự đốn mơ hình Cách chọn số sóng R2Y Q2 RMSEE (%) RMSEP (%) 250-450 nm 0.971 0.968 5.88 3.07 250-350 nm 0.992 0.991 2.10 0.50 Từ kết mơ hình ta nhận thấy cách chọn theo thơng số VIP có giá trị RMSEE RMSEP thấp chứng tỏ khả dự đốn mơ hình tốt Vì ta chọn cách lấy bước sóng theo thơng số VIP để xây dựng mơ hình 2.3.2 Loại bỏ điểm lệch thô Hotelling’s T2Range DModX Tiến hành xem xét xác điểm lệch thô gây sai lệch mơ hình đồ thị Hotelling’s T2Range ta thấy có điểm có khả gây sai lệch cho mơ hình “19R”, “31R”, “B1”, “B5” “19R80B10N10” (Hình 30) Hình 30 Đồ thị Hotelling’s T2Range Các điểm “19R”, “31R”, “B1”, “B5” điểm gốc đại diện cho mẫu cà phê (19R, 31R) bắp (B1, B5) có mẫu 100 % nguyên chất loại thành phần trộn nên xem điểm gốc mơ hình Nếu loại điểm mơ hình tính định hướng đường chuẩn nên ta không loại bỏ điểm Để giải thích cho lệch thơ điểm ta xem xét lại phổ hấp thu điểm mẫu có tính chất (Hình 31) Phổ hấp thu điểm mẫu “19R” “31R” có giá trị thấp so với điểm lại, nên kết mơ hình trả hợp lý 35 Tương tự phổ hấp thu hai điểm “B1” “B5” có giá trị độ hấp thu thấp điểm mẫu bắp lại Điểm “19R80B10N10” sau thử loại bỏ khỏi mơ hình khả dự đốn mơ hình khơng thay đổi nên định giữ lại điểm Bảng 13 R2Y Q2 mơ hình Mơ hình R2 Y Q2 Mơ hình ban đầu 0.992 0.991 Mơ hình sau loại điểm 0.992 0.991 2.0 1.8 Độ hấp thu (Abs) 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 250 270 290 310 330 350 Bước sóng (nm) 15CU 17CU R2 31R 28R 25R 19R 4R Hình 31 Phổ hấp thu điểm mẫu tính chất với "19R" "31R" Sau trình thiết lập loại điểm lệch thơ, kết thu mơ hình tốt có PC với R2= 0.992 Q2= 0.991, RMSEE= 2.10 % RMSEP= 0.50 % Hình 32 Biều đồ “Overview” mơ hình xác định hàm lượng caffeine 36 Đánh giả khả dự đốn mơ hình thơng qua mẫu kiểm tra 2.3.3 Bộ mẫu kiểm tra gồm 30 mẫu người thực tự trộn để có mẫu hỗn hợp thành phần cà phê, bắp đậu nành xác định hàm lượng caffeine phương pháp HPLC-UV Đồng thời dịch chiết mẫu đem quét phổ máy trắc quang để lấy phổ hấp thu đưa vào mô hình Kết thu thể qua Bảng 14 Trong 30 mẫu kiểm tra có 21 mẫu có kết dự đoán hàm lượng caffeine tốt (RSD < %), mẫu có kết dự đốn khơng tốt (RSD ~ 6-7 %) mẫu có kết dự đoán (RSD ~ 10 %) Bảng 14 Bộ kiểm tra mơ hình xác định hàm lượng caffeine CAF (%) CAF (%) xác định HPLC-UV theo mô hình PLS 0.455 0.442 2.8 0.750 0.757 0.9 1.285 1.318 2.6 0.201 0.198 1.8 0.636 0.638 0.4 1.116 1.268 13.6 1.332 1.331 0.1 0.841 0.865 2.8 1.593 1.528 4.1 10 1.612 1.619 0.4 11 1.676 1.523 9.2 12 0.160 0.170 6.6 13 0.454 0.435 4.3 14 1.441 1.288 10.6 15 1.362 1.278 6.2 16 0.529 0.548 3.5 17 0.858 0.885 3.2 18 1.373 1.377 0.3 19 1.522 1.612 6.0 STT % RSD 37 20 0.338 0.337 0.2 21 0.838 0.832 0.7 22 0.644 0.555 13.8 23 1.092 1.043 4.5 24 2.040 2.138 4.8 25 0.569 0.609 7.1 26 0.954 0.957 0.2 27 1.516 1.538 1.5 28 1.750 1.885 7.7 29 1.350 1.352 0.1 30 1.494 1.538 2.9 Như vậy, kết dự đốn mơ hình có tốt ta nhận thấy sử dụng mơ hình để xác định hàm lượng caffeine mẫu cà phê thị trường Các mẫu thị trường (16 mẫu) kiểm tra hàm lượng caffeine phương pháp HPLC- UV sau liệu caffeine sử dụng để so sánh với kết mơ hình Bảng 15 Kết kiểm tra mẫu thi trường mơ hình xác định hàm lượng caffeine STT Tên mẫu Hãng sản xuất CAF (%) CAF (%) (**) 6.6 % RSD Thunder No 1864 Café 1.920 (***) 1.793 Thunder Culi Robusta 1864 Café 1.679 1.600 4.7 Sơn Tùng Robusta Sơn Tùng 1.288 1.187 7.8 Highlands Di sản Highlands Coffee 1.727 1.565 9.4 Chinh phục S Trung Nguyên 1.056 1.042 1.4 Sáng tạo Culi Robusta Trung Nguyên 0.799 0.871 9.0 Nes Việt nguyên chất NESCAFÉ 1.931 1.920 0.6 Nes Việt nguyên chất NESCAFÉ 2.098 2.046 2.5 Đức Thành Đức Thành 1.657 1.761 6.3 10 Cà phê Chocolate Phúc Long 1.963 1.940 1.2 38 11 Cà phê Buôn Mê Thuột Đ Phương Vy 0.835 0.784 6.1 12 Phương Vy Robusta Phương Vy 0.613 0.637 3.9 13 Bùi Văn Thọ Robusta Bùi Văn Ngọ 1.844 1.795 2.7 14 Bùi Văn Thọ Culi Bùi Văn Ngọ 1.864 1.812 2.8 15 Robusta Thượng Hạng Light Coffee 1.863 1.975 6.0 16 Robusta Đặc Biệt Light Coffee 1.744 1.874 7.4 (**)Xác định phương pháp HPLC-UV (***)Xác định từ mơ hình PLS Với mẫu cà phê thị trường cho kết dự đoán tốt (sai lệch