Thêm Thư viện cần thiết
library(readxl)
library(dplyr)
library(FactoMineR)
library(SensoMineR)
library(tidyr)
library(ggplot2)
library(factoextra)
Thêm đường dẫn cho dữ liệu cảm quan
combined_df <- read_excel("D:/Project/SENSORY_CAMQUANTIEU_BKHN/Huan_luyen_hoi_dong-23-7-2023/HoaDescriptors.xlsx")
combined_df <- as.data.frame(combined_df)
Tổng quan về dữ liệu
## Panelist Sample Rep Thuoc.bac
## Length:80 Length:80 Length:80 Min. :0.000
## Class :character Class :character Class :character 1st Qu.:0.000
## Mode :character Mode :character Mode :character Median :0.000
## Mean :0.575
## 3rd Qu.:1.000
## Max. :4.000
## Bac.ha Cam.chanh Chanh.xanh O.mai Beo
## Min. :0.0 Min. :0.00 Min. :0.000 Min. :0.000 Min. :0.000
## 1st Qu.:0.0 1st Qu.:0.00 1st Qu.:0.750 1st Qu.:0.000 1st Qu.:0.000
## Median :1.0 Median :0.00 Median :2.000 Median :2.000 Median :1.000
## Mean :1.3 Mean :0.85 Mean :2.062 Mean :2.125 Mean :1.325
## 3rd Qu.:2.0 3rd Qu.:1.25 3rd Qu.:3.000 3rd Qu.:4.000 3rd Qu.:2.000
## Max. :4.0 Max. :5.00 Max. :5.000 Max. :5.000 Max. :5.000
## Beo.oi Co Thao.moc Huong.hoa
## Min. :0.0000 Min. :0.000 Min. :0.00 Min. :0.000
## 1st Qu.:0.0000 1st Qu.:0.000 1st Qu.:0.00 1st Qu.:0.000
## Median :0.0000 Median :2.000 Median :0.00 Median :1.000
## Mean :0.3125 Mean :2.163 Mean :0.85 Mean :1.575
## 3rd Qu.:0.0000 3rd Qu.:4.000 3rd Qu.:1.25 3rd Qu.:3.000
## Max. :3.0000 Max. :5.000 Max. :4.00 Max. :5.000
## Moc Tieu Thuoc Thong
## Min. :0.0000 Min. :0.000 Min. :0.000 Min. :0.00
## 1st Qu.:0.0000 1st Qu.:0.000 1st Qu.:0.000 1st Qu.:0.00
## Median :0.0000 Median :1.000 Median :0.000 Median :0.00
## Mean :0.4625 Mean :1.288 Mean :0.175 Mean :0.25
## 3rd Qu.:0.0000 3rd Qu.:2.000 3rd Qu.:0.000 3rd Qu.:0.00
## Max. :5.0000 Max. :4.000 Max. :4.000 Max. :2.00
## Len.men Man Vo.cay Go
## Min. :0.0000 Min. :0.00 Min. :0.0000 Min. :0.000
## 1st Qu.:0.0000 1st Qu.:0.00 1st Qu.:0.0000 1st Qu.:0.000
## Median :0.0000 Median :0.00 Median :0.0000 Median :0.000
## Mean :0.3125 Mean :0.40 Mean :0.3625 Mean :0.575
## 3rd Qu.:0.0000 3rd Qu.:0.25 3rd Qu.:0.0000 3rd Qu.:1.000
## Max. :5.0000 Max. :4.00 Max. :5.0000 Max. :5.000
## Dat
## Min. :0.000
## 1st Qu.:0.000
## Median :0.000
## Mean :0.175
## 3rd Qu.:0.000
## Max. :4.000
Đánh giá năng lực hội đồng tiêu
### Test nang luc hoi dong############
res.panelperf <- panelperf(combined_df,firstvar=4,formul="~Sample+Panelist+
+ Rep+Sample:Panelist+Sample:Rep+Panelist:Rep")
coltable(res.panelperf$p.value[order(res.panelperf$p.value[,1]),], col.lower="mistyrose",
main.title= "Năng lực hội đồng Tiêu")
