Alternative sensory eval - FP
Quoc Cuong Nguyen
2022-11-14
Clean
rm(list = ls()); graphics.off()Set R packages path and load packages
.libPaths("C:/Rpack")
library(FactoMineR)
library(cluster)
library(SensoMineR)
library(openxlsx)Import data
FCP <- read.xlsx("FPdata.xlsx", sheet = 1, rowNames = T)
head(FCP)## Wrapping_CM Green_CM Floral_CM Citrus_CM Vanilla_CM
## Angel 7.45 0.000000e+00 1.60 8.881784e-16 6.00
## Aromatics Elixir 5.30 1.110223e-15 3.25 0.000000e+00 0.60
## Chanel N5 9.20 8.500000e-01 6.20 -6.661338e-16 0.75
## Cinema 9.10 4.500000e-01 6.05 5.000000e-01 6.75
## Coco Mademoiselle 8.90 8.500000e-01 8.80 3.050000e+00 0.40
## J'adore EP 7.80 4.950000e+00 9.15 2.050000e+00 0.75
## Woody_CM Fruity_CR Greedy_CR Oriental_CR Citrus_CR
## Angel 4.500000e-01 6.600000e+00 7.10 0.15 0.90
## Aromatics Elixir 2.100000e+00 2.500000e-01 0.15 7.90 0.65
## Chanel N5 3.885781e-16 -1.776357e-15 0.10 8.55 0.15
## Cinema -1.110223e-16 3.450000e+00 4.30 0.15 0.90
## Coco Mademoiselle 3.500000e-01 1.900000e+00 1.25 9.35 0.15
## J'adore EP 0.000000e+00 8.850000e+00 1.50 4.20 0.20
## Floral_CR Vanilla_CR Wrapping_CR Citrus_GV Greedy_GV Woody_GV
## Angel 4.10 7.05 4.00 2.70 4.90 3.30
## Aromatics Elixir 4.75 0.65 7.95 0.95 0.10 4.35
## Chanel N5 7.45 0.45 7.95 0.65 0.30 1.35
## Cinema 7.20 4.00 1.75 3.85 7.00 0.15
## Coco Mademoiselle 9.55 1.10 3.20 2.25 5.65 0.25
## J'adore EP 9.25 0.60 4.60 6.80 2.00 0.20
## Vanilla_GV Fruity_GV Green_GV Wrapping_GV Oriental_GV
## Angel 3.95 1.55 5.551115e-17 4.85 4.70
## Aromatics Elixir 0.10 0.40 4.163336e-17 5.45 5.10
## Chanel N5 0.15 0.50 2.000000e-01 6.60 2.75
## Cinema 0.15 7.90 -3.469447e-17 4.30 5.55
## Coco Mademoiselle 0.10 7.15 5.000000e-02 5.85 7.75
## J'adore EP 0.20 5.45 2.775558e-17 3.95 3.10
## Marine_GV Green_MLD Marine_MLD Oriental_MLD
## Angel 5.000000e-02 6.661338e-16 1.110223e-15 4.25
## Aromatics Elixir 1.500000e-01 9.500000e-01 -4.440892e-16 4.60
## Chanel N5 1.000000e+00 7.500000e-01 1.000000e-01 1.00
## Cinema -5.551115e-17 4.000000e-01 7.771561e-16 6.00
## Coco Mademoiselle 5.551115e-17 4.000000e-01 5.500000e-01 0.50
## J'adore EP 2.220446e-16 1.350000e+00 3.950000e+00 0.10
## Woody_MLD Floral_MLD Vanilla_MLD Wrapping_MLD
## Angel 5.000000e-01 0.10 9.200000e+00 5.00
## Aromatics Elixir 5.000000e-01 1.50 -1.332268e-15 9.60
## Chanel N5 1.000000e-01 5.30 5.000000e-02 8.40
## Cinema -3.885781e-16 3.80 8.750000e+00 8.50
## Coco Mademoiselle 1.942890e-16 9.85 1.332268e-15 0.30
## J'adore EP 5.000000e-01 9.90 1.100000e+00 1.65
## Citrus_MLD Fruity_MLD Greedy_MLD Marine_NMA
## Angel 1.776357e-15 1.500000e-01 9.850000e+00 0.10
## Aromatics Elixir 1.000000e-01 1.776357e-15 0.000000e+00 0.20
## Chanel N5 1.000000e-01 1.776357e-15 -2.664535e-15 0.10
## Cinema 3.500000e-01 8.500000e+00 5.350000e+00 3.55
## Coco Mademoiselle 3.500000e-01 9.950000e+00 6.500000e-01 4.50
## J'adore EP 5.250000e+00 9.800000e+00 1.350000e+00 3.70
## Citrus_NMA Woody_NMA Oriental_NMA Greedy_NMA Floral_NMA
## Angel -1.734723e-17 4.30 9.50 4.20 5.40
## Aromatics Elixir 2.000000e-01 3.00 6.70 0.05 9.25
## Chanel N5 1.040834e-17 0.25 8.80 3.75 6.30
## Cinema 1.387779e-17 3.40 4.90 3.40 8.35
## Coco Mademoiselle 6.938894e-18 3.30 4.10 7.55 9.65
## J'adore EP 4.510281e-17 4.60 3.95 0.15 10.00
## Fruity_NMA Wrapping_NMA Oriental_PR Wrapping_PR Floral_PR
## Angel 3.500000e+00 9.30 7.65 8.65 2.45
## Aromatics Elixir 4.440892e-16 8.65 8.20 9.25 3.05
## Chanel N5 -4.440892e-16 8.90 3.30 9.10 8.20
## Cinema 3.200000e+00 7.75 4.00 6.00 6.75
## Coco Mademoiselle 7.500000e+00 5.45 2.60 4.10 6.55
## J'adore EP 3.850000e+00 7.65 0.65 3.15 7.05
## Woody_PR Greedy_PR Marine_PR Fruity_PR Citrus_PR
## Angel 0.50 9.65 1.500000e-01 0.95 0.15
## Aromatics Elixir 3.25 1.15 -1.665335e-16 4.00 0.20
## Chanel N5 0.70 0.35 2.000000e-01 2.20 0.25
## Cinema 0.75 0.75 2.500000e-01 0.25 0.10
## Coco Mademoiselle 0.70 3.20 2.000000e-01 0.95 0.25
## J'adore EP 1.05 0.65 7.500000e-01 4.90 0.45
## Greedy_RL Fruity_RL Marine_RL Green_SD
## Angel 1.000000e+01 3.000000e-01 5.000000e-02 0.10
## Aromatics Elixir 5.000000e-02 -2.220446e-15 1.040834e-17 0.15
## Chanel N5 -4.440892e-16 -1.776357e-15 2.081668e-17 0.10
## Cinema 2.500000e+00 6.500000e+00 1.000000e-01 0.10
## Coco Mademoiselle 5.050000e+00 5.350000e+00 -1.734723e-17 0.15
## J'adore EP 4.440892e-16 5.000000e+00 2.428613e-17 6.15
## Wrapping_SD Marine_SD Citrus_SD Oriental_SM Floral_SM
## Angel 9.75 0.10 0.05 9.65 2.30
## Aromatics Elixir 7.45 0.15 4.20 9.40 3.70
## Chanel N5 5.40 0.10 6.55 0.35 7.55
## Cinema 8.90 0.10 0.10 2.85 4.30
## Coco Mademoiselle 8.05 0.05 0.05 0.05 8.45
## J'adore EP 2.60 0.15 0.10 0.10 7.70
## Fruity_SM Woody_SM Wrapping_SM Marine_SM Green_SM Citrus_SM
## Angel 2.05 1.20 9.65 0.05 0.05 0.05
## Aromatics Elixir 0.50 1.45 8.75 0.05 0.10 0.00
## Chanel N5 0.65 0.10 8.80 0.05 0.10 2.00
## Cinema 3.40 0.05 6.75 0.05 0.05 0.05
## Coco Mademoiselle 2.25 0.10 4.10 0.05 1.10 3.75
## J'adore EP 4.25 0.05 5.55 2.85 2.00 1.30
## Wrapping_SO Floral_SO Greedy_SO Vanilla_SO Green_SO Fruity_SO
## Angel 8.55 3.35 9.85 3.00 0.40 3.05
## Aromatics Elixir 9.20 4.65 0.10 1.60 1.30 0.55
## Chanel N5 9.65 6.45 0.95 0.10 0.15 2.05
## Cinema 3.95 8.05 2.35 1.25 0.75 6.70
## Coco Mademoiselle 6.15 9.00 0.50 2.05 0.55 7.85
## J'adore EP 2.55 7.60 0.45 0.60 0.40 9.20
## Marine_SQ Floral_SQ Wrapping_SQ Greedy_SQ Green_SQ
## Angel 0.60 5.10 3.45 5.75 0.45
## Aromatics Elixir 0.35 2.75 2.95 0.30 0.20
## Chanel N5 0.10 3.85 2.35 1.05 0.95
## Cinema 1.60 2.70 2.35 1.35 0.15
## Coco Mademoiselle 0.65 4.20 2.60 3.10 0.15
## J'adore EP 0.45 2.95 2.90 2.30 0.65
## Vanilla_ST Citrus_ST Wrapping_ST Woody_ST
## Angel 9.800000e+00 6.938894e-17 1.000000e+01 2.220446e-16
## Aromatics Elixir 1.850000e+00 0.000000e+00 8.500000e+00 4.700000e+00
## Chanel N5 4.300000e+00 8.326673e-17 8.600000e+00 -8.881784e-16
## Cinema 2.950000e+00 1.387779e-17 4.440892e-16 1.250000e+00
## Coco Mademoiselle 1.450000e+00 8.326673e-17 8.500000e-01 1.221245e-15
## J'adore EP 4.440892e-16 -4.163336e-17 -2.664535e-15 1.110223e-15
## Fruity_ST Oriental_ST
## Angel 1.776357e-15 4.000000e-01
## Aromatics Elixir 1.100000e+00 9.100000e+00
## Chanel N5 4.350000e+00 8.200000e+00
## Cinema 9.000000e+00 8.881784e-16
## Coco Mademoiselle 9.100000e+00 2.400000e+00
## J'adore EP 9.000000e+00 -1.332268e-15Run MFA
MFA
res.mfa <- MFA(FCP, group = c(6,7,9,10,8,8,3,4,8,6,5,6), type = rep("s",12), name.group = c("CM","CR","GV","MLD","NMA","PR","RL","SD","SM","SO","SQ","ST"), graph = F)
round(res.mfa$eig,3)## eigenvalue percentage of variance cumulative percentage of variance
## comp 1 7.901 31.088 31.088
## comp 2 5.589 21.989 53.076
## comp 3 3.184 12.525 65.602
## comp 4 1.830 7.200 72.802
## comp 5 1.606 6.320 79.121
## comp 6 1.241 4.884 84.006
## comp 7 1.088 4.281 88.286
## comp 8 0.946 3.722 92.008
## comp 9 0.878 3.454 95.461
## comp 10 0.644 2.533 97.995
## comp 11 0.510 2.005 100.000Plot MFA results
windows(10,10)
plot.MFA(res.mfa, choix = "ind", graph.type = "classic")
windows(10,10)
plot.MFA(res.mfa, choix = "var", graph.type = "classic")
Excess nb of variables
windows(10,10)
plot.MFA(res.mfa, choix = "var", select = "cos2 0.8")
plot.MFA(res.mfa, choix = "var", select = "contrib 10")
dimdesc(res.mfa)## $Dim.1
##
## Link between the variable and the continuous variables (R-square)
## =================================================================================
## correlation p.value
## Floral_SO 0.9144486 3.122001e-05
## Floral_MLD 0.9135088 3.291967e-05
## Fruity_ST 0.8622927 3.080797e-04
## Fruity_MLD 0.8218378 1.039482e-03
## Green_SD 0.7905379 2.207113e-03
## Fruity_RL 0.7820988 2.648216e-03
## Fruity_SO 0.7808247 2.720246e-03
## Floral_SM 0.7804696 2.740587e-03
## Floral_CM 0.7786151 2.848695e-03
## Green_SM 0.7531895 4.683229e-03
## Floral_PR 0.7495261 5.007423e-03
## Floral_CR 0.7483746 5.112747e-03
## Green_CM 0.7282594 7.236187e-03
## Marine_MLD 0.7156150 8.872644e-03
## Floral_NMA 0.6794342 1.508604e-02
## Marine_PR 0.6790919 1.515693e-02
## Fruity_SM 0.6427840 2.416990e-02
## Citrus_MLD 0.6329578 2.716306e-02
## Green_MLD 0.5930702 4.209904e-02
## Greedy_SO -0.5906471 4.316157e-02
## Greedy_PR -0.5953423 4.111974e-02
## Woody_GV -0.6114113 3.464998e-02
## Oriental_NMA -0.6188564 3.191422e-02
## Wrapping_SO -0.6593233 1.968315e-02
## Wrapping_MLD -0.6691841 1.731705e-02
## Vanilla_ST -0.7607128 4.067401e-03
## Woody_SM -0.7629432 3.897220e-03
## Wrapping_SM -0.7730777 3.190917e-03
## Wrapping_SD -0.7856073 2.457364e-03
## Wrapping_PR -0.8157756 1.215528e-03
## Wrapping_ST -0.8374420 6.758509e-04
## Oriental_SM -0.8819701 1.475046e-04
## Oriental_PR -0.9117125 3.637026e-05
##
## $Dim.2
##
## Link between the variable and the continuous variables (R-square)
## =================================================================================
## correlation p.value
## Greedy_GV 0.8060096 0.0015463882
## Greedy_SQ 0.8032239 0.0016523509
## Vanilla_MLD 0.7877336 0.0023469241
## Greedy_CR 0.7786013 0.0028495156
## Greedy_RL 0.7666430 0.0036268930
## Greedy_MLD 0.7495639 0.0050039928
## Vanilla_CR 0.7440603 0.0055222838
## Fruity_GV 0.7425696 0.0056693819
## Greedy_SO 0.7144352 0.0090383163
## Vanilla_SO 0.6697183 0.0171951601
## Greedy_PR 0.6104338 0.0350212111
## Wrapping_SQ 0.6098463 0.0352457009
## Fruity_NMA 0.6007757 0.0388435480
## Fruity_CR 0.5962001 0.0407542560
## Woody_CM -0.6491457 0.0223644921
## Woody_ST -0.6529025 0.0213455921
## Oriental_CR -0.6713523 0.0168261591
## Woody_PR -0.6725649 0.0165561219
## Wrapping_CR -0.6819819 0.0145660629
## Citrus_SD -0.8128357 0.0013087455
## Oriental_ST -0.8381417 0.0006622529
##
## $Dim.3
##
## Link between the variable and the continuous variables (R-square)
## =================================================================================
## correlation p.value
## Citrus_CR 0.7650391 0.003742272
## Marine_GV 0.6632631 0.018711148
## Green_SQ 0.6616582 0.019102741
## Citrus_ST 0.6582644 0.019950560
## Wrapping_NMA 0.6033982 0.037777645
## Green_CM 0.5923072 0.042431626
## Fruity_CR 0.5872962 0.044662030
## Wrapping_CM -0.6416764 0.024494697
## Marine_NMA -0.6545676 0.020905004
## Oriental_GV -0.6576580 0.020104897res.hcpc <- HCPC(res.mfa) # cluster analysis## [1] "Click on the graph to cut the tree"
Group structure
Information shared between groups
windows(10,10)
plot.MFA(res.mfa, choix = "group", graph.type = "classic")
plot.MFA(res.mfa, choix = "ind", partial = "all", habillage = "group")
Relation between each panelist in final product space
res.mfa$group$correlation## Dim.1 Dim.2 Dim.3 Dim.4 Dim.5
## CM 0.8420411 0.8076817 0.7147114 0.5142463 0.2279936
## CR 0.9386000 0.8238809 0.6996418 0.5302432 0.4402471
## GV 0.8191236 0.8752828 0.8207007 0.4298697 0.7613124
## MLD 0.9041807 0.8610971 0.5506487 0.3983491 0.7618902
## NMA 0.7205395 0.8164912 0.6208808 0.4069858 0.7666184
## PR 0.9676199 0.8536618 0.6201820 0.5699056 0.7013625
## RL 0.8302574 0.8705330 0.3269824 0.7330304 0.2283961
## SD 0.8707191 0.9005601 0.5925741 0.4393523 0.2383805
## SM 0.9358547 0.6504157 0.7123539 0.4556074 0.6241345
## SO 0.9430347 0.9545336 0.5949590 0.6398699 0.5407440
## SQ 0.6718014 0.8696481 0.7288451 0.6976033 0.4543576
## ST 0.9147933 0.8202435 0.8463279 0.5811781 0.5139779