I. Giới thiệu về số liệu và định dạng của số liệu (hedochoc, sensochoc, sensopanels)…………………………….....................Trang 4 II. Phân tích mô tả số liệu…………………………………Trang 6 III. Đối với số liệu hedochoc: 1. Biểu đồ cột về điểm thị hiếu của người tiêu dùng đối với các sản phẩm chocolates…………………………………...............Trang 9 2. Phân tích phương sai về sự sai khác mức độ ưa thích của người tiêu dùng đối với các sản phẩm. nếu cần thiết sử dụng Tukey để đánh giá kiểm định hậu nghiệm………….................................................Trang 10 IV. Đối với số liệu sensochoc: 1. phân tích phương sai sự khác nhau trên từng tiêu chuẩn đối với hội đồng mô tả. vẽ biểu đồ boxplot để hỗ trợ quá trình phân tích……………………………………………..............Trang 11 2. Đánh giá sự đồng thuận của thành viên hội đồng….......Trang 22 3. Vẽ biểu đồ radar mô tả tiêu chuẩn của sản phẩm…....…Trang 23
TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP TP.HỒ CHÍ MINH … … VIỆN CÔNG NGHIỆ SINH HỌC VÀ THỰC PHẨM BÀI TẬP NHÓM MÔN Thiết kế thí nghiệm xử lý số liệu thực nghiệm ngành công nghệ thực phẩm Đề tài : Thống kê mô tả loại Chocolates Pháp GVHD: TS.Nguyễn Bá Thanh Lớp: ĐHTP10C TP,HCM ngày 22 tháng 10 năm 2015 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Danh sách nhóm Họ tên Mã số sinh viên Diệp Trung Hiếu 14126071 Lê Thị Dung 14132021 Phạm Võ Tân Thuận 14130201 Lê Ngọc Châu 14131841 Nguyễn Thị Tuyết Mai 14129161 Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Mục Lục I Giới thiệu số liệu định dạng số liệu (hedochoc, sensochoc, sensopanels)…………………………… Trang II Phân tích mô tả số liệu…………………………………Trang III Đối với số liệu hedochoc: Biểu đồ cột điểm thị hiếu người tiêu dùng sản phẩm chocolates………………………………… .Trang Phân tích phương sai sai khác mức độ ưa thích người tiêu dùng sản phẩm cần thiết sử dụng Tukey để đánh giá kiểm định hậu nghiệm………… .Trang 10 IV Đối với số liệu sensochoc: phân tích phương sai khác tiêu chuẩn hội đồng mô tả vẽ biểu đồ boxplot để hỗ trợ trình phân tích…………………………………………… Trang 11 Đánh giá đồng thuận thành viên hội đồng… .Trang 22 Vẽ biểu đồ radar mô tả tiêu chuẩn sản phẩm… …Trang 23 Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH I Giới thiệu số liệu định dạng số liệu (hedochoc, sensochoc, sensopanels) Khảo sát tính chất loại chocolates Pháp Mô tả tính chất cảm quan sản phẩm: loại chocolate đánh giá 29 bảng danh sách xếp theo 14 tính chất cảm quan Dữ liệu mức độ ưa thích: loại chocolate đánh giá theo cấu trúc từ (thích) đến 10 (không thích) 222 người tiêu dùng Hedochoc: Dữ liệu gồm dòng 222 cột, dòng tương ứng với loại chocolate cột điểm số ưa thích cho 222 người tiêu dùng tham gia nghiên cứu Từ số liệu ta thấy row.names giá trị định tính, giá trị khác giá trị định lượng Bảng liệu kết khảo sát mức độ thích 222 người tiêu dùng sản phẩm chocolate, từ ta kiểm tra người tiêu dùng có khác hay không điều mong đợi khác người thử; kiểm tra mẫu thử có khác hay không, việc kiểm tra mong đợi phải có khác xem sản phẩm khác nhiều Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Sensochoc: liệu gồm 348 dòng 19 cột, biến chất lượng (Panelist, Session, Rank, Product) 14 tính chất cảm quan Từ bảng số liệu ta thấy giá trị đầu (panelist, session, rank, product) giá trị định tính giá trị lại(cocaA, milkA, caramel, vanilla, )là giá trị định lượng Từ liệu sensochoc ta nêu lên xem có khác biệt tính chất chocolate cocaA, milkA, caramel, vanilla, với panelist hay session, rank, product; so sánh thành phần chocolates với nhóm yếu tố từ ta biết với tính chất tập trung có sản phẩm nào; hay sản phẩm có tính chất trội sản phẩm khác xem Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH II Phân tích mô tả liệu Hedochoc Vars n mean sd se Choc1 222 5.33 2.55 0.17 Choc2 222 5.69 2.22 0.15 Choc3 222 5.77 2.66 0.18 Choc4 222 5.77 2.30 0.15 Choc5 222 5.63 2.05 0.14 Choc6 222 5.91 2.00 0.13 Sensochoc Vars n mean sd se CocoaA 58 7.09 1.82 0.24 MilkA 58 3.59 2.34 0.31 CocoaF 58 8.07 1.57 0.21 MilkF 58 1.57 1.45 0.19 Caramel 58 1.67 2.05 0.27 Vanilla 58 1.10 1.36 0.18 58 3.14 2.14 0.28 Acidity 58 4.66 2.93 0.38 Bitterness 10 58 7.07 1.88 0.25 Astringency 11 58 4.76 2.82 0.37 Crunchy 12 58 5.97 2.52 0.33 Melting 13 58 4.74 2.60 0.34 Sticky 14 58 3.76 2.37 0.31 Granular 15 58 3.45 2.80 0.37 Choc Sweetness Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Vars n mean sd se CocoaA 58 6.55 2.18 0.29 MilkA 58 4.00 2.23 0.29 CocoaF 58 6.91 1.59 0.21 MilkF 58 2.38 1.55 0.20 Caramel 58 2.78 2.43 0.32 Vanilla 58 1.81 1.89 0.25 58 4.62 2.22 0.29 Acidity 58 3.14 2.54 0.33 Bitterness 10 58 4.95 2.31 0.30 Astringency 11 58 3.16 2.48 0.33 Crunchy 12 58 7.71 1.64 0.22 Melting 13 58 4.33 2.30 0.30 Sticky 14 58 3.83 2.22 0.29 Granular 15 58 3.16 2.77 0.36 Vars n mean sd se CocoaA 58 4.67 2.13 0.28 MilkA 58 6.05 2.65 0.35 CocoaF 58 3.38 1.95 0.26 MilkF 58 7.71 2.03 0.27 Caramel 58 6.33 2.52 0.33 Vanilla 58 3.67 2.69 0.35 58 7.60 1.61 0.21 Acidity 58 1.57 1.67 0.22 Bitterness 10 58 1.40 1.45 0.19 Astringency 11 58 1.21 1.39 0.18 Crunchy 12 58 2.98 2.11 0.28 Melting 13 58 7.31 1.73 0.23 Sticky 14 58 5.03 2.63 0.35 Granular 15 58 1.60 1.95 0.26 Choc Sweetness Choc Sweetness Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Vars n mean sd se CocoaA 58 6.26 1.92 0.25 MilkA 58 4.10 2.28 0.30 CocoaF 58 6.69 2.02 0.27 MilkF 58 2.59 1.94 0.25 Caramel 58 2.67 2.34 0.31 Vanilla 58 2.12 2.16 0.28 58 4.29 2.11 0.28 Acidity 58 3.93 2.87 0.38 Bitterness 10 58 5.19 2.50 0.33 Astringency 11 58 3.69 2.50 0.33 Crunchy 12 58 6.10 2.30 0.30 Melting 13 58 4.38 2.47 0.32 Sticky 14 58 4.10 2.52 0.33 Granular 15 58 3.55 2.66 0.35 Vars n mean sd se CocoaA 58 6.79 1.68 0.22 MilkA 58 4.17 2.05 0.27 CocoaF 58 6.79 1.68 0.22 MilkF 58 3.12 2.04 0.27 Caramel 58 3.41 2.66 0.35 Vanilla 58 1.79 2.02 0.26 58 5.22 2.27 0.30 Acidity 58 3.09 2.49 0.33 Bitterness 10 58 4.88 2.19 0.29 Astringency 11 58 3.10 2.29 0.30 Crunchy 12 58 6.64 2.20 0.29 Melting 13 58 4.47 2.46 0.32 Sticky 14 58 3.22 2.76 0.36 Granular 15 58 3.07 2.69 0.35 Choc Sweetness Choc Sweetness Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Vars n mean sd se CocoaA 58 6.36 1.85 0.24 MilkA 58 4.57 2.47 0.32 CocoaF 58 6.22 1.87 0.25 MilkF 58 3.36 2.17 0.28 Caramel 58 3.26 2.42 0.32 Vanilla 58 1.91 2.12 0.28 58 5.62 2.23 0.29 Acidity 58 2.67 2.19 0.29 Bitterness 10 58 4.19 2.16 0.28 Astringency 11 58 2.76 2.57 0.34 Crunchy 12 58 7.33 1.82 0.24 Melting 13 58 4.21 2.32 0.30 Sticky 14 58 3.93 2.58 0.34 Granular 15 58 3.17 2.78 0.36 Choc Sweetness III Đối với liệu hedochoc Vẽ biểu đồ cột điểm thị hiếu người tiêu dùng sản phẩm chocolate: Trang | GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH 9.00 8.00 7.00 Thị hiếu 6.00 5.00 4.00 3.00 2.00 1.00 0.00 choc1 choc2 choc3 choc4 choc5 choc6 Loại chocolate Biểu đồ thị hiếu người tiêu dùng loại chocolate Đánh giá nhận xét : biểu đồ cột điểm thị hiếu nói chung cho loại chocolates Lần lượt từ trái sang phải choc1 đến choc6 Dựa vào giá trị trung bình cho sản phẩm chocolates cho thấy điểm ưa thích sản phẩm điểm khác biệt lớn, tương tự biểu đồ cảm quan ta thấy choc1, choc5 có khác biệt rõ so với sản phẩm khác, sản phẩm ưa thích nhiều.( Ở đánh giá theo thang điểm đến 10, điểm coi ưa thích nhất, điểm 10 không thích) Phân tích phương sai sai khác mức độ ưa thích người tiêu dùng sản phẩm cần thiết sử dụng Tukey để đánh giá kiểm định hậu nghiệm Để phân tích phương sai ta dùng kiểm định ANOVA để đánh giá Analysis of Variance Table Response: thihieu Df Sum Sq Mean Sq F value Pr(>F) loai 42.6 8.5171 1.5976 0.1577 Residuals 1326 7069.2 5.3312 Trang | 10 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Caramel Analysis of Variance Table Response: Caramel Df Sum Sq Mean Sq F value Pr(>F) Panelist 28 818.11 29.218 19.1208 3.152e-08 *** Session 0.03 0.026 0.0169 0.898019 Rank 40.34 8.067 5.2793 0.004162 ** Product 717.56 143.512 93.9161 9.054e-12 *** Panelist:Session 28 171.22 6.115 4.0018 0.002103 ** Panelist:Rank 140 563.07 4.022 2.6320 0.012333 * Panelist:Product 89 274.45 3.084 2.0180 0.051108 Session:Rank 14.86 2.973 1.9453 0.139297 Session:Product 9.61 1.921 1.2573 0.326796 Rank:Product 24 74.31 3.096 2.0263 0.068673 Residuals 17 25.98 1.528 Signif codes: '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' Dựa kết qur phân tích cho thấy có khác biệt Caramel nhóm: Panelist, Rank, Product, Panelis:Session Biểu đồ thể khác biệt: Panelist 10 Biểu đồ khác Caramel Panelist 11 14 17 20 23 26 29 Trang | 19 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Rank 10 Biểu đồ khác Caramel Rank Product 10 Biểu đồ khác Caramel Product choc1 choc2 choc3 choc4 choc5 choc6 Trang | 20 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Panelist:Session 10 Biểu đồ khác Caramel Panelist:Session 1.1 6.1 12.1 18.1 24.1 1.2 6.2 12.2 18.2 24.2 Trang | 21 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Đánh giá đồng thuận thành viên hội đồng Analysis of Variance Table Response: score Df Sum Sq Mean Sq F value Pr(>F) choc 221 2259.1 10.2221 2.3382 < 2.2e-16 *** Residuals 1110 4852.7 4.3718 Signif codes: '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' Tóm lược ước số phân tích: Call: lm(formula = score ~ choc + lchoc) Residuals: Min 1Q Median 3Q Max -6.1667 -1.1667 0.1667 1.3333 5.6667 Coefficients: (221 not defined because of singularities) Nhận xét: Trị số p=2.2e-16 có nghĩa số liệu thống kê có nghĩa Hay có khác biệt đánh giá thành viên hội đồng mức độ yêu thích loại chocolate * Biểu đồ đánh giá đồng thuận thành viên hội đồng loại chocolate Trang | 22 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Vẽ biểu đồ rada mô tả tiêu chuẩn sản phẩm Radar Chart CocoaA MilkA Granular Sticky CocoaF MilkF Melting Caramel Crunchy Vanilla Sweetness Bitterness Acidity choc1 choc2 choc3 Astringency choc4 choc5 choc6 n Đây đồ thị thể Sensory profile loại chocolates: + Choc1 thuộc tính Crunchy, CocoaF, Acidity có nhiều loại chocolates ưa thích nhóm đối tượng nghiên cứu choc1,còn thuộc tính khác mức độ trung bình, không làm giảm giá trị choc1 + Choc2 thuộc tính caramel không ưa thích nhóm sản phẩm này, đặc tính CocoaF chấm điểm thấp độ không chấm điểm cao loai choc2 dây hai thuộc tính ưa thích choc2, thuộc tính khác mức độ trung bình đánh giá + Choc3 thuộc tính CocoaF thuộc tính ưa thích nhóm sản phẩm này, cá thuộc tính lại đánh giá mức độ trung bình, + Choc4 Melting, Granular, CocoaF thuộc tính ưa thích sản phẩm này, Sticky không ưa thích nhóm đối tuogwj nghiên cứu choc4 Trang | 23 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH + Choc5 thuộc tính CocoaF chấm điểm 0, mức độ ưa thích đạt điểm ưa thích nhất, Caramel, Bitterness Granular không ưa thích choc5 + Choc6 có điểm tương quan so với choc1, thuộc tinh CocoaF ưa thích nhóm đối tượng Đánh giá chung: Thuộc tính CocoaF ưa thích nhiều loại sản phẩm nghiên cứu, yếu tố Melting, Granular ưa thích không cao so với CocoaF, thuộc tính lại ảnh hưởng đến đánh giá tham luận viên, người tieu dùng thuộc tính Phần phụ lục (Code lệnh sử dụng) Dữ liệu lấy từ Pakages SensomineR >library(SensoMineR) >library(FactomineR) >library(Rcmdr) > data(chocolates), package="SensoMineR") III Đối với liệu hedochoc >data(hedochoc) >attach(hedochoc) >hedo hedo attach(hedo) Vẽ biểu đồ cột điểm thị hiếu người tiêu dùng sản phẩm chocolate: > a2 a1 a3 a4 a5 a6 thihieu=c(a1,a2,a3,a4,a5,a6) >barplot(thihieu, xlab =”Loại sô-cô-la”, ylab=”Điểm thị hiếu”, col=”blue”, main=”Biểu đồ thị hiếu \n Người tiêu dung loại sô-cô-la”) Phân tích phương sai sai khác mức độ ưa thích người tiêu dùng sản phẩm cần thiết sử dụng Tukey để đánh giá kiểm định hậu nghiệm >loai=c(rep("choc1",222),rep("choc2",222),rep("choc3",222),rep("choc4",222),rep("c hoc5",222),rep("choc6",222)) > data=data.frame(thihieu,loai) Trang | 24 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH >attach(data) > ana=lm(thihieu~loai) > anova(ana) IV Đối với số liệu senschoc > data(sensochoc) >attach(sensochoc) > Panelist Session Rank Product Phantich1=aov(CocoaA~Panelist+Session+Rank+Product+Panelist:Session+Panelist :Rank+Panelist:Product +Session :Rank +Session :Product +Rank:Product) >anova(Phantich1) >boxplot(CocoaA~Panelist, col=”red”, main=”Biểu đồ khác CocoaA Panelist”) >boxplot(CocoaA~Product, col=”yellow”, main=”Biểu đồ khác CocoaA Product”) CocoaF >Phantich2=aov(CocoaF~Panelist+Session+Rank+Product+Panelist:Session+Panelist: Rank+Panelist:Product +Session :Rank +Session :Product +Rank:Product) >anova(Phantich2) >boxplot(CocoaF~Panelist, col=”red”, main=”Biểu đồ khác CocoaF Panelist”) >boxplot(CocoaF~Product, col=”yellow”, main=”Biểu đồ khác CocoaF Product”) MilkA >Phantich3=aov(MilkA~Panelist+Session+Rank+Product+Panelist:Session+Panelist: Rank+Panelist:Product +Session :Rank +Session :Product +Rank:Product) >anova(Phantich3) >boxplot(MilkA~Panelist, col=”red”, main=”Biểu đồ khác MilkA Panelist”) >boxplot(MilkA~Rank, col=”green”, main=”Biểu đồ khác MilkA Rank”) >boxplot(MilkA~Product, col=”yellow”, main=”Biểu đồ khác MilkA Product”) Trang | 25 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH MilkF >Phantich4=aov(MilkF~Panelist+Session+Rank+Product+Panelist:Session+Panelist:R ank+Panelist:Product +Session :Rank +Session :Product +Rank:Product) >anova(Phantich4) >boxplot(MilkF~Panelist, col=”red”, main=”Biểu đồ khác MilkF Panelist”) >boxplot(MilkF~Rank, col=”green”, main=”Biểu đồ khác MilkF Rank”) >boxplot(MilkF~Product, col=”yellow”, main=”Biểu đồ khác MilkF Product”) Caramel >Phantich5=aov(Caramel~Panelist+Session+Rank+Product+Panelist:Session+Panelist :Rank+Panelist:Product +Session :Rank +Session :Product +Rank:Product) >anova(Phantich5) >boxplot(Caramel~Panelist, col=”red”, main=”Biểu đồ khác Caramel Panelist”) >boxplot(Caramel ~Rank, col=”green”, main=”Biểu đồ khác Caramel Rank”) >boxplot(Caramel ~Product, col=”yellow”, main=”Biểu đồ khác Caramel Product”) >boxplot(Caramel ~Panelist:Session, col=”blue”, main=”Biểu đồ khác Caramel Panelist:Session”) Đánh giá đồng thuận thành viên hội đồng *Phân tích phương sai hai chiều ANOVA: >choc=gl(222,6,1332) > lchoc=gl(6,222,1332) > id=1:1332 >score=c(j1,j2,j3,j4,j5,j6,j7,j8,j9,j10,j11,j12,j13,j14,j15,j16,j17,j18,j19,j20,j21, +j22,j23,j24,j25,j26,j27,j28,j29,j30,j31,j32,j33,j34,j35,j36,j37,j38,j39,j40, +j41,j42,j43,j44,j45,j46,j47,j48,j49,j50,j51,j52,j53,j54,j55,j56,j57,j58,j59,j60, +j61,j62,j63,j64,j65,j66,j67,j68,j69,j70,j71,j72,j73,j74,j75,j76,j77,j78,j79, +j80,j81,j82,j83,j84,j85,j86,j87,j88,j89,j90,j91,j92,j93,j94,j95,j96,j97,j98,j99, +j100,j101,j102,j103,j104,j105,j106,j107,j108,j109,j110,j111,j112,j113,j114, +j115,j116,j117,j118,j119,j120,j121,j122,j123,j124,j125,j126,j127,j128,j129,j130,j131,j +132,j133,j134, +j135,j136,j137,j138,j139,j140,j141,j142,j143,j144,j145,j146,j147,j148,j149,j150,j151, +j152,j153,j154,j155,j156,j157,j158,j159,j160,j161,j162,j163,j164,j165,j166,j167, +j168,j169,j170,j171,j172,j173,j174,j175,j176,j177,j178,j179,j180,j181,j182,j183, +j184,j185,j186,j187,j188,j189,j190,j191,j192,j193,j194,j195,j196,j197,j198,j199, +j200,j201,j202,j203,j204,j205,j206,j207,j208,j209,j210,j211,j212,j213,j214,j215, +j216,j217,j218,j219,j220,j221,j222) > data=data.frame(choc,lchoc,score,id) > attach(data) Trang | 26 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH > twoway= lm(score~choc+lchoc) > anova(twoway ) *Các số liệu phân tích khác biệt thành viên hội đồng với nhau: Residuals: Min 1Q Median 3Q Max -6.1667 -1.1667 0.1667 1.3333 5.6667 Coefficients: (221 not defined because of singularities) Estimate Std Error t value Pr(>|t|) (Intercept) 5.167e+00 8.536e-01 6.053 1.94e-09 *** choc[T.2] 5.000e-01 1.207e+00 0.414 0.67881 choc[T.3] 6.667e-01 1.207e+00 0.552 0.58088 choc[T.4] -1.167e+00 1.207e+00 -0.966 0.33403 choc[T.5] -8.333e-01 1.207e+00 -0.690 0.49014 choc[T.6] 1.333e+00 1.207e+00 1.105 0.26961 choc[T.7] -1.500e+00 1.207e+00 -1.243 0.21429 choc[T.8] -2.942e-15 1.207e+00 0.000 1.00000 choc[T.9] -2.779e-15 1.207e+00 0.000 1.00000 choc[T.10] -1.500e+00 1.207e+00 -1.243 0.21429 choc[T.11] -1.000e+00 1.207e+00 -0.828 0.40763 choc[T.12] 3.333e-01 1.207e+00 0.276 0.78250 choc[T.13] -1.167e+00 1.207e+00 -0.966 0.33403 choc[T.14] -2.167e+00 1.207e+00 -1.795 0.07295 choc[T.15] -2.500e+00 1.207e+00 -2.071 0.03859 * choc[T.16] 1.667e-01 1.207e+00 0.138 0.89021 choc[T.17] -1.667e-01 1.207e+00 -0.138 0.89021 choc[T.18] -5.000e-01 1.207e+00 -0.414 0.67881 choc[T.19] -8.333e-01 1.207e+00 -0.690 0.49014 choc[T.20] -1.833e+00 1.207e+00 -1.519 0.12912 choc[T.21] -2.167e+00 1.207e+00 -1.795 0.07295 choc[T.22] -5.000e-01 1.207e+00 -0.414 0.67881 choc[T.23] 1.667e-01 1.207e+00 0.138 0.89021 choc[T.24] -1.333e+00 1.207e+00 -1.105 0.26961 choc[T.25] 2.500e+00 1.207e+00 2.071 0.03859 * choc[T.26] 2.000e+00 1.207e+00 1.657 0.09785 choc[T.27] 2.667e+00 1.207e+00 2.209 0.02738 * choc[T.28] 2.167e+00 1.207e+00 1.795 0.07295 choc[T.29] 2.000e+00 1.207e+00 1.657 0.09785 choc[T.30] 1.667e+00 1.207e+00 1.381 0.16767 choc[T.31] 8.333e-01 1.207e+00 0.690 0.49014 choc[T.32] 1.167e+00 1.207e+00 0.966 0.33403 choc[T.33] 8.333e-01 1.207e+00 0.690 0.49014 choc[T.34] 2.000e+00 1.207e+00 1.657 0.09785 choc[T.35] 1.333e+00 1.207e+00 1.105 0.26961 choc[T.36] 2.167e+00 1.207e+00 1.795 0.07295 choc[T.37] 2.167e+00 1.207e+00 1.795 0.07295 choc[T.38] 1.667e+00 1.207e+00 1.381 0.16767 Trang | 27 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH choc[T.39] choc[T.40] choc[T.41] choc[T.42] choc[T.43] choc[T.44] choc[T.45] choc[T.46] choc[T.47] choc[T.48] choc[T.49] choc[T.50] choc[T.51] choc[T.52] choc[T.53] choc[T.47] choc[T.48] choc[T.49] choc[T.50] choc[T.51] choc[T.52] choc[T.53] choc[T.54] choc[T.55] choc[T.56] choc[T.57] choc[T.58] choc[T.59] choc[T.60] choc[T.61] choc[T.62] choc[T.63] choc[T.64] choc[T.65] choc[T.66] choc[T.67] choc[T.68] choc[T.69] choc[T.70] choc[T.71] choc[T.72] choc[T.73] choc[T.74] choc[T.75] choc[T.76] choc[T.77] choc[T.78] choc[T.79] choc[T.80] 1.667e+00 1.207e+00 1.381 0.16767 3.167e+00 1.207e+00 2.623 0.00883 ** -2.840e-15 1.207e+00 0.000 1.00000 2.833e+00 1.207e+00 2.347 0.01910 * 3.000e+00 1.207e+00 2.485 0.01310 * 2.167e+00 1.207e+00 1.795 0.07295 1.500e+00 1.207e+00 1.243 0.21429 3.500e+00 1.207e+00 2.899 0.00381 ** 6.667e-01 1.207e+00 0.552 0.58088 6.667e-01 1.207e+00 0.552 0.58088 -6.667e-01 1.207e+00 -0.552 0.58088 8.333e-01 1.207e+00 0.690 0.49014 5.000e-01 1.207e+00 0.414 0.67881 6.667e-01 1.207e+00 0.552 0.58088 -8.333e-01 1.207e+00 -0.690 0.49014 6.667e-01 1.207e+00 0.552 0.58088 6.667e-01 1.207e+00 0.552 0.58088 -6.667e-01 1.207e+00 -0.552 0.58088 8.333e-01 1.207e+00 0.690 0.49014 5.000e-01 1.207e+00 0.414 0.67881 6.667e-01 1.207e+00 0.552 0.58088 -8.333e-01 1.207e+00 -0.690 0.49014 5.000e-01 1.207e+00 0.414 0.67881 5.000e-01 1.207e+00 0.414 0.67881 1.333e+00 1.207e+00 1.105 0.26961 5.000e-01 1.207e+00 0.414 0.67881 3.333e-01 1.207e+00 0.276 0.78250 5.000e-01 1.207e+00 0.414 0.67881 -6.667e-01 1.207e+00 -0.552 0.58088 1.000e+00 1.207e+00 0.828 0.40763 5.000e-01 1.207e+00 0.414 0.67881 -1.333e+00 1.207e+00 -1.105 0.26961 -2.862e-15 1.207e+00 0.000 1.00000 -6.667e-01 1.207e+00 -0.552 0.58088 6.667e-01 1.207e+00 0.552 0.58088 5.000e-01 1.207e+00 0.414 0.67881 -3.000e+00 1.207e+00 -2.485 0.01310 * -1.333e+00 1.207e+00 -1.105 0.26961 2.500e+00 1.207e+00 2.071 0.03859 * 1.667e+00 1.207e+00 1.381 0.16767 1.500e+00 1.207e+00 1.243 0.21429 1.667e-01 1.207e+00 0.138 0.89021 1.833e+00 1.207e+00 1.519 0.12912 -1.500e+00 1.207e+00 -1.243 0.21429 -3.333e-01 1.207e+00 -0.276 0.78250 -2.667e+00 1.207e+00 -2.209 0.02738 * 2.500e+00 1.207e+00 2.071 0.03859 * -3.333e-01 1.207e+00 -0.276 0.78250 -2.863e-15 1.207e+00 0.000 1.00000 Trang | 28 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH choc[T.81] -3.333e-01 1.207e+00 -0.276 0.78250 choc[T.82] 1.333e+00 1.207e+00 1.105 0.26961 choc[T.83] 1.333e+00 1.207e+00 1.105 0.26961 choc[T.84] -1.667e-01 1.207e+00 -0.138 0.89021 choc[T.85] 2.833e+00 1.207e+00 2.347 0.01910 * choc[T.86] 8.333e-01 1.207e+00 0.690 0.49014 choc[T.87] 3.500e+00 1.207e+00 2.899 0.00381 ** choc[T.88] -5.000e-01 1.207e+00 -0.414 0.67881 choc[T.89] -1.000e+00 1.207e+00 -0.828 0.40763 choc[T.90] -2.167e+00 1.207e+00 -1.795 0.07295 choc[T.91] 1.667e+00 1.207e+00 1.381 0.16767 choc[T.92] 5.000e-01 1.207e+00 0.414 0.67881 choc[T.93] 3.333e-01 1.207e+00 0.276 0.78250 choc[T.94] 1.667e+00 1.207e+00 1.381 0.16767 choc[T.95] -2.000e+00 1.207e+00 -1.657 0.09785 choc[T.96] -1.833e+00 1.207e+00 -1.519 0.12912 choc[T.97] 1.833e+00 1.207e+00 1.519 0.12912 choc[T.98] 1.667e-01 1.207e+00 0.138 0.89021 choc[T.99] -1.167e+00 1.207e+00 -0.966 0.33403 choc[T.100] -5.000e-01 1.207e+00 -0.414 0.67881 choc[T.101] -5.000e-01 1.207e+00 -0.414 0.67881 choc[T.102] -2.550e-15 1.207e+00 0.000 1.00000 choc[T.103] -6.667e-01 1.207e+00 -0.552 0.58088 choc[T.104] -1.167e+00 1.207e+00 -0.966 0.33403 choc[T.105] 5.000e-01 1.207e+00 0.414 0.67881 choc[T.106] -3.333e-01 1.207e+00 -0.276 0.78250 choc[T.107] 6.667e-01 1.207e+00 0.552 0.58088 choc[T.108] -1.000e+00 1.207e+00 -0.828 0.40763 choc[T.109] 6.667e-01 1.207e+00 0.552 0.58088 choc[T.110] 1.000e+00 1.207e+00 0.828 0.40763 choc[T.111] 6.667e-01 1.207e+00 0.552 0.58088 choc[T.112] 5.000e-01 1.207e+00 0.414 0.67881 choc[T.113] -3.333e-01 1.207e+00 -0.276 0.78250 choc[T.114] 2.167e+00 1.207e+00 1.795 0.07295 choc[T.115] 2.000e+00 1.207e+00 1.657 0.09785 choc[T.116] 2.833e+00 1.207e+00 2.347 0.01910 * choc[T.117] 1.500e+00 1.207e+00 1.243 0.21429 choc[T.118] 3.333e-01 1.207e+00 0.276 0.78250 choc[T.119] -3.333e-01 1.207e+00 -0.276 0.78250 choc[T.120] 1.833e+00 1.207e+00 1.519 0.12912 choc[T.121] -5.000e-01 1.207e+00 -0.414 0.67881 choc[T.122] -2.490e-15 1.207e+00 0.000 1.00000 choc[T.123] 2.000e+00 1.207e+00 1.657 0.09785 choc[T.124] -3.497e-15 1.207e+00 0.000 1.00000 choc[T.125] 1.167e+00 1.207e+00 0.966 0.33403 choc[T.126] 3.333e-01 1.207e+00 0.276 0.78250 choc[T.127] 1.000e+00 1.207e+00 0.828 0.40763 choc[T.128] 6.667e-01 1.207e+00 0.552 0.58088 choc[T.129] 1.500e+00 1.207e+00 1.243 0.21429 Trang | 29 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH choc[T.130] choc[T.131] choc[T.132] choc[T.133] choc[T.134] choc[T.135] choc[T.136] choc[T.137] choc[T.138] choc[T.139] choc[T.140] choc[T.141] choc[T.142] choc[T.143] choc[T.144] choc[T.145] choc[T.146] choc[T.147] choc[T.148] choc[T.149] choc[T.150] choc[T.151] choc[T.152] choc[T.153] choc[T.154] choc[T.155] choc[T.156] choc[T.157] choc[T.158] choc[T.159] choc[T.160] choc[T.161] choc[T.162] choc[T.163] choc[T.164] choc[T.165] choc[T.166] choc[T.167] choc[T.168] choc[T.169] choc[T.170] choc[T.171] choc[T.172] choc[T.173] choc[T.174] choc[T.175] choc[T.176] choc[T.177] choc[T.178] 3.333e-01 1.207e+00 1.000e+00 1.207e+00 8.333e-01 1.207e+00 -8.333e-01 1.207e+00 -1.833e+00 1.207e+00 1.500e+00 1.207e+00 1.333e+00 1.207e+00 1.667e-01 1.207e+00 2.167e+00 1.207e+00 3.333e-01 1.207e+00 -8.333e-01 1.207e+00 -1.167e+00 1.207e+00 -1.667e-01 1.207e+00 -2.648e-15 1.207e+00 -5.000e-01 1.207e+00 2.167e+00 1.207e+00 -5.000e-01 1.207e+00 -6.667e-01 1.207e+00 5.000e-01 1.207e+00 2.500e+00 1.207e+00 -3.333e-01 1.207e+00 -1.667e-01 1.207e+00 2.333e+00 1.207e+00 1.667e-01 1.207e+00 -1.667e-01 1.207e+00 5.000e-01 1.207e+00 1.667e+00 1.207e+00 3.333e-01 1.207e+00 3.333e-01 1.207e+00 -1.000e+00 1.207e+00 1.667e-01 1.207e+00 2.000e+00 1.207e+00 1.000e+00 1.207e+00 2.333e+00 1.207e+00 -4.124e-15 1.207e+00 6.667e-01 1.207e+00 2.667e+00 1.207e+00 1.667e-01 1.207e+00 2.167e+00 1.207e+00 5.000e-01 1.207e+00 -8.333e-01 1.207e+00 -3.333e-01 1.207e+00 1.833e+00 1.207e+00 -3.333e-01 1.207e+00 2.167e+00 1.207e+00 -3.333e-01 1.207e+00 1.500e+00 1.207e+00 -2.695e-15 1.207e+00 -6.667e-01 1.207e+00 0.276 0.78250 0.828 0.40763 0.690 0.49014 -0.690 0.49014 -1.519 0.12912 1.243 0.21429 1.105 0.26961 0.138 0.89021 1.795 0.07295 0.276 0.78250 -0.690 0.49014 -0.966 0.33403 -0.138 0.89021 0.000 1.00000 -0.414 0.67881 1.795 0.07295 -0.414 0.67881 -0.552 0.58088 0.414 0.67881 2.071 0.03859 * -0.276 0.78250 -0.138 0.89021 1.933 0.05350 0.138 0.89021 -0.138 0.89021 0.414 0.67881 1.381 0.16767 0.276 0.78250 0.276 0.78250 -0.828 0.40763 0.138 0.89021 1.657 0.09785 0.828 0.40763 1.933 0.05350 0.000 1.00000 0.552 0.58088 2.209 0.02738 * 0.138 0.89021 1.795 0.07295 0.414 0.67881 -0.690 0.49014 -0.276 0.78250 1.519 0.12912 -0.276 0.78250 1.795 0.07295 -0.276 0.78250 1.243 0.21429 0.000 1.00000 -0.552 0.58088 Trang | 30 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH choc[T.179] choc[T.180] choc[T.181] choc[T.182] choc[T.183] choc[T.184] choc[T.185] choc[T.186] choc[T.187] choc[T.188] choc[T.189] choc[T.190] choc[T.191] choc[T.192] choc[T.193] choc[T.194] choc[T.195] choc[T.196] choc[T.197] choc[T.198] choc[T.199] choc[T.200] choc[T.201] choc[T.202] choc[T.203] choc[T.204] choc[T.205] choc[T.206] choc[T.207] choc[T.208] choc[T.209] choc[T.210] choc[T.211] choc[T.212] choc[T.213] choc[T.214] choc[T.215] choc[T.216] choc[T.217] choc[T.218] choc[T.219] choc[T.220] choc[T.221] choc[T.222] 1.333e+00 1.207e+00 -1.500e+00 1.207e+00 -1.667e-01 1.207e+00 1.667e-01 1.207e+00 -1.333e+00 1.207e+00 1.000e+00 1.207e+00 6.667e-01 1.207e+00 -6.667e-01 1.207e+00 2.167e+00 1.207e+00 1.333e+00 1.207e+00 -8.333e-01 1.207e+00 5.000e-01 1.207e+00 1.333e+00 1.207e+00 1.667e+00 1.207e+00 1.000e+00 1.207e+00 2.667e+00 1.207e+00 -5.000e-01 1.207e+00 1.167e+00 1.207e+00 1.833e+00 1.207e+00 1.667e+00 1.207e+00 -3.434e-15 1.207e+00 2.667e+00 1.207e+00 2.167e+00 1.207e+00 1.167e+00 1.207e+00 6.667e-01 1.207e+00 2.333e+00 1.207e+00 1.333e+00 1.207e+00 6.667e-01 1.207e+00 1.167e+00 1.207e+00 2.333e+00 1.207e+00 2.000e+00 1.207e+00 2.000e+00 1.207e+00 -3.333e-01 1.207e+00 2.500e+00 1.207e+00 5.000e-01 1.207e+00 -6.667e-01 1.207e+00 1.333e+00 1.207e+00 -3.263e-15 1.207e+00 -1.167e+00 1.207e+00 6.667e-01 1.207e+00 8.333e-01 1.207e+00 -1.667e-01 1.207e+00 -2.667e+00 1.207e+00 1.667e-01 1.207e+00 1.105 0.26961 -1.243 0.21429 -0.138 0.89021 0.138 0.89021 -1.105 0.26961 0.828 0.40763 0.552 0.58088 -0.552 0.58088 1.795 0.07295 1.105 0.26961 -0.690 0.49014 0.414 0.67881 1.105 0.26961 1.381 0.16767 0.828 0.40763 2.209 0.02738 * -0.414 0.67881 0.966 0.33403 1.519 0.12912 1.381 0.16767 0.000 1.00000 2.209 0.02738 * 1.795 0.07295 0.966 0.33403 0.552 0.58088 1.933 0.05350 1.105 0.26961 0.552 0.58088 0.966 0.33403 1.933 0.05350 1.657 0.09785 1.657 0.09785 -0.276 0.78250 2.071 0.03859 * 0.414 0.67881 -0.552 0.58088 1.105 0.26961 0.000 1.00000 -0.966 0.33403 0.552 0.58088 0.690 0.49014 -0.138 0.89021 -2.209 0.02738 * 0.138 0.89021 Trang | 31 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Vẽ biểu đồ rada mô tả tiêu chuẩn sản phẩm >library(psych) >webplot = function(data, data.row = NULL, y.cols = NULL, main = NULL, add = F, col = "red", lty = 1, scale = T) { if (!is.matrix(data) & !is.data.frame(data)) stop("Requires matrix or data.frame") if (is.null(y.cols)) y.cols = colnames(data)[sapply(data, is.numeric)] if (sum(!sapply(data[, y.cols], is.numeric)) > 0) { out = paste0("\"", colnames(data)[!sapply(data, is.numeric)], "\"", collapse = ", ") stop(paste0("All y.cols must be numeric\n", out, " are not numeric")) } if (is.null(data.row)) data.row = if (is.character(data.row)) if (data.row %in% rownames(data)) { data.row = which(rownames(data) == data.row) } else { stop("Invalid value for data.row:\nMust be a valid rownames(data) or row-index value") } if (is.null(main)) main = rownames(data)[data.row] if (scale == T) { data = scale(data[, y.cols]) data = apply(data, 2, function(x) x/max(abs(x))) } data = as.data.frame(data) n.y = length(y.cols) min.rad = 360/n.y polar.vals = (90 + seq(0, 360, length.out = n.y + 1)) * pi/180 # if (add == F) { plot(0, xlim = c(-2.2, 2.2), ylim = c(-2.2, 2.2), type = "n", axes = F, xlab = "", ylab = "") title(main) lapply(polar.vals, function(x) lines(c(0, * cos(x)), c(0, * sin(x)))) lapply(1:n.y, function(x) text(2.15 * cos(polar.vals[x]), 2.15 * sin(polar.vals[x]), y.cols[x], cex = 0.8)) lapply(seq(0.5, 2, 0.5), function(x) lines(x * cos(seq(0, * pi, length.out = 100)), x * sin(seq(0, * pi, length.out = 100)), lwd = 0.5, lty = 2, col = "gray60")) lines(cos(seq(0, * pi, length.out = 100)), sin(seq(0, * pi, length.out = 100)), lwd = 1.2, col = "gray50") } r = + data[data.row, y.cols] xs = r * cos(polar.vals) ys = r * sin(polar.vals) xs = c(xs, xs[1]) ys = c(ys, ys[1]) lines(xs, ys, col = col, lwd = 2, lty = lty) } > choc1=subset(sensochoc,Product=="choc1") > summary(choc1) > choc2=subset(sensochoc,Product=="choc2") > summary(choc2) > choc3=subset(sensochoc,Product=="choc3") > summary(choc3) > choc4=subset(sensochoc,Product=="choc4") > summary(choc4) > choc5=subset(sensochoc,Product=="choc5") > summary(choc5) > choc6=subset(sensochoc,Product=="choc6") > summary(choc6) > par(mar = c(1, 1, 6, 1)) > webplot(choc1,main="Radar Chart") > webplot(choc2,add = T, col = "blue", lty = 2) > webplot(choc3,add = T, col = "green") > webplot(choc4,add = T, col = "violet", lty = 2) > webplot(choc5,add = T, col = "yellow", lty = 1) Trang | 32 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH > webplot(choc6,add = T, col = "black", lty = 2) >legend("bottomright",lty=c(1,2,1,2,1,2),col=c("red","blue","green","violet","yellow", "black),c("choc1","choc2","choc3","choc4","choc5","choc6",bty="n")) TÀI LIỆU THAM KHẢO Phân tích số liệu tạo biểu đồ R Nguyễn Văn Tuấn, nhà xuất khoa học kỹ thuật http://statisticstoproveanything.blogspot.com/2013/11/spider-web-plots-in-r.html http://sensominer.free/ http://www.r-bloggers.com/ Trang | 33 [...]... quả phân tích cho thấy bậc từ do ở đây là 1(Df), từ kết quả phân tích phương sai cho thấy P-value=0.1577>0.05, điều này chứng tỏ sự khác nhau của các loai chocolates không có ý nghĩa thống kê, nên ta chọn loại chocolates nào cũng được IV Đối với số liệu senschoc 1 Phân tích phương sai sự khác nhau trên từng tiêu chuẩn đối với hội đồng mô tả Vẽ biểu đồ boxplot để hỗ trợ quá trình phân tích Phân tích. .. 4.3718 Signif codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Tóm lược ước số phân tích: Call: lm(formula = score ~ choc + lchoc) Residuals: Min 1Q Median 3Q Max -6.1667 -1.1667 0.1667 1.3333 5.6667 Coefficients: (221 not defined because of singularities) Nhận xét: Trị số p=2.2e-16 có nghĩa là các số liệu thống kê có nghĩa Hay đúng hơn là có sự khác biệt trong sự đánh giá của các thành viên hội... >legend("bottomright",lty=c(1,2,1,2,1,2),col=c("red","blue","green","violet","yellow", "black),c("choc1","choc2","choc3","choc4","choc5","choc6",bty="n")) TÀI LIỆU THAM KHẢO 1 Phân tích số liệu và tạo biểu đồ bằng R Nguyễn Văn Tuấn, nhà xuất bản khoa học và kỹ thuật 2 http://statisticstoproveanything.blogspot.com/2013/11/spider-web-plots-in-r.html 3 http://sensominer.free/ 4 http://www.r-bloggers.com/ Trang... ana=lm(thihieu~loai) > anova(ana) IV Đối với số liệu senschoc > data(sensochoc) >attach(sensochoc) > Panelist Session Rank Product Phantich1=aov(CocoaA~Panelist+Session+Rank+Product+Panelist:Session+Panelist... thuộc tính Crunchy, CocoaF, và Acidity có nhiều trong loại chocolates này và được ưa thích trên nhóm đối tượng nghiên cứu choc1,còn các thuộc tính khác ở mức độ trung bình, không làm giảm giá trị của choc1 + Choc2 thuộc tính caramel không được ưa thích trên nhóm sản phẩm này, đặc tính CocoaF được chấm điểm thấp nhất và độ ngọt cũng không được chấm điểm cao ở loai choc2 này và dây là hai thuộc tính được... 1.4447 0.2193139 Residuals 17 36.45 2.144 Signif codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Trang | 14 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH Dựa vào kết quả phân tích cho thấy có sự khác biệt về MilkA ở các nhóm: Panelist, Rank, Product và Panelist:Product Ta bỏ qua sự khác biệt ở nhóm Panelist:Product vì nó không quá lớn Biểu đò sự khác biệt: Panelist 0 2 4 6 8 10 Biểu đồ sự khác nhau của... +j200,j201,j202,j203,j204,j205,j206,j207,j208,j209,j210,j211,j212,j213,j214,j215, +j216,j217,j218,j219,j220,j221,j222) > data=data.frame(choc,lchoc,score,id) > attach(data) Trang | 26 GIÁO VIÊN HƯỚNG DẪN: NGUYỄN BÁ THANH > twoway= lm(score~choc+lchoc) > anova(twoway ) *Các số liệu phân tích sự khác biệt giữa các thành viên hội đồng với nhau: Residuals: Min 1Q Median 3Q Max -6.1667 -1.1667 0.1667 1.3333 5.6667 Coefficients: (221 not defined because of singularities) Estimate Std... Rank:Product 24 49.07 2.044 1.1685 0.3761031 Residuals 17 29.74 1.750 Signif codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Dựa vào kết quả phân tích cho thấy có sự khác nhau giữa CocoaF trên 4 nhóm: Panelist, Senssion, Rank, Product Ta bỏ qua sự khác biệt ở Session và Rank Biểu đồ thể hiện sự khác nhau: Panelist 0 2 4 6 8 10 Biểu đồ sự khác nhau của CocoaF trong Panelist 1 3 5 7 9 11 14 17 20... 0.6874402 Session:Product 5 16.36 3.272 1.1544 0.3708950 Rank:Product 24 35.53 1.481 0.5223 0.9293156 Residuals 17 48.19 2.835 Signif codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Dựa vào kết quả phân tích trên, ta chỉ thấy có sự khác biệt về CocoaA trên 3 nhóm: Panelist, Rank, Product Ở nhóm Rank sự khác biệt không quá lơn nên ta có thể bỏ qua Sự khác biệt được biểu hiên qua ba biểu đồ sau... a6 thihieu=c(a1,a2,a3,a4,a5,a6) >barplot(thihieu, xlab =”Loại sô-cô-la”, ylab=”Điểm thị hiếu”, col=”blue”, main=”Biểu đồ thị hiếu của \n Người tiêu dung đối với các loại sô-cô-la”) 2 Phân tích phương sai về sự sai khác mức độ ưa thích của người tiêu dùng đối với các sản phẩm nếu cần thiết sử dụng Tukey để đánh giá kiểm định hậu nghiệm >loai=c(rep("choc1",222),rep("choc2",222),rep("choc3",222),rep("choc4",222),rep("c