1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 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 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75
| # ************************************
# Read a CSV File (comming from Excel)
# ************************************
ConsoAlim=read.table("d:/user/data/SC/CMP/Acrylamide/BDD&Analyses_R/ConsoGroupesAliments080911.csv",header=TRUE,sep=";")
# ******************************
# Print it in the R Table Editor
# ******************************
fix(ConsoAlim)
# **********************************************
# Print only 4 variables present in the DataFile
# **********************************************
ConsoAlim[,c("CER_J1", "CER_J2", "CER_Hist","CER_HistF")]
# *****************************************
# Counting and comparing previous Variables
# *****************************************
# Détermination du nombre d'enregistrement (N)
NbRec=length(ConsoAlim[,c("SCId")])
# Analyse des profils de consommation J1/J2 et Habituellement
# On crée une variable qui va identifier les différents profils possibles :
# Profil = 0 ==> Aucune consommation
# = 1 ==> Consommation en J1 seulement
# = 2 ==> Consommation en J2 seulement
# = 3 ==> Consommation en J1 et J2, mais pas habituellement ==> Problème !!!!
# = 4 ==> Consommation Habituellement, mais pas en J1/J2
# = 5 ==> Consommation en J1 et Habituellement
# = 6 ==> Consommation en J2 et Habituellement
# = 7 ==> Consommation en J1, J2 et Habituellement
# _______________________________________________________________
# Fonction "ProfilConso" permettant l'Analyse Aliment par Aliment
# _______________________________________________________________
ProfilConso <- function(Pattern,X1,X2,XHist,NbX1,NbX2,NbXH)
{
# Initialisation des variables de type compteur
NbX1=0; NbX2=0; NbXH=0
for(i in 1:NbRec) {
Pattern[i]=0
if (X1[i]=="O") {
Pattern[i]=Pattern[i]+1
NbX1=NbX1+1
}
if (X2[i]=="O") {
Pattern[i]=Pattern[i]+2
NbX2=NbX2+1
}
if (XHist[i]=="O") {
Pattern[i]=Pattern[i]+4
NbXH=NbXH+1
}
}
print(NbRec); print(NbX1); print(NbX2); print(NbXH); print(Pattern); table(Pattern)
}
ProfilConso(ConsoAlim$CER_Pattern, ConsoAlim$CER_J1, ConsoAlim$CER_J2, ConsoAlim$CER_Hist, NbCERJ1, NbCERJ2, NbCERHist)
ProfilConso(ConsoAlim$PAIN_Pattern, ConsoAlim$PAIN_J1, ConsoAlim$PAIN_J2, ConsoAlim$PAIN_Hist, NbPAINJ1, NbPAINJ2, NbPAINHist)
ProfilConso(ConsoAlim$FRI_H_Pattern, ConsoAlim$FRI_H_J1, ConsoAlim$FRI_H_J2, ConsoAlim$FRI_H_Hist, NbFRI_HJ1, NbFRI_HJ2, NbFRI_HHist)
ProfilConso(ConsoAlim$FRI_F_Pattern, ConsoAlim$FRI_F_J1, ConsoAlim$FRI_F_J2, ConsoAlim$FRI_F_Hist, NbFRI_FJ1, NbFRI_FJ2, NbFRI_FHist)
ProfilConso(ConsoAlim$CHI_Pattern, ConsoAlim$CHI_J1, ConsoAlim$CHI_J2, ConsoAlim$CHI_Hist, NbCHIJ1, NbCHIJ2, NbCHIHist)
ProfilConso(ConsoAlim$CHI_AU_Pattern, ConsoAlim$CHI_AU_J1, ConsoAlim$CHI_AU_J2, ConsoAlim$CHI_AU_Hist, NbCHI_AUJ1, NbCHI_AUJ2, NbCHI_AUHist)
ProfilConso(ConsoAlim$BRE_Pattern, ConsoAlim$BRE_J1, ConsoAlim$BRE_J2, ConsoAlim$BRE_Hist, NbBREJ1, NbBREJ2, NbBREHist)
ProfilConso(ConsoAlim$CRAQ_Pattern, ConsoAlim$CRAQ_J1, ConsoAlim$CRAQ_J2, ConsoAlim$CRAQ_Hist, NbCRAQJ1, NbCRAQJ2, NbCRAQHist)
ProfilConso(ConsoAlim$OLI_Pattern, ConsoAlim$OLI_J1, ConsoAlim$OLI_J2, ConsoAlim$OLI_Hist, NbOLIJ1, NbOLIJ2, NbOLIHist)
ProfilConso(ConsoAlim$BISC_PE_Pattern,ConsoAlim$BISC_PE_J1,ConsoAlim$BISC_PE_J2,ConsoAlim$BISC_PE_Hist,NbBISC_PEJ1,NbBISC_PEJ2,NbBISC_PEHist)
ProfilConso(ConsoAlim$BISC_AU_Pattern,ConsoAlim$BISC_AU_J1,ConsoAlim$BISC_AU_J2,ConsoAlim$BISC_AU_Hist,NbBISC_AUJ1,NbBISC_AUJ2,NbBISC_AUHist)
ProfilConso(ConsoAlim$MAIS_Pattern, ConsoAlim$MAIS_J1, ConsoAlim$MAIS_J2, ConsoAlim$MAIS_Hist, NbMAISJ1, NbMAISJ2, NbMAISHist)
ProfilConso(ConsoAlim$CAF_Pattern, ConsoAlim$CAF_J1, ConsoAlim$CAF_J2, ConsoAlim$CAF_Hist, NbCAFJ1, NbCAFJ2, NbCAFHist)
ProfilConso(ConsoAlim$AMAN_Pattern, ConsoAlim$AMAN_J1, ConsoAlim$AMAN_J2, ConsoAlim$AMAN_Hist, NbAMANJ1, NbAMANJ2, NbAMANHist)
ProfilConso(ConsoAlim$PRU_Pattern, ConsoAlim$PRU_J1, ConsoAlim$PRU_J2, ConsoAlim$PRU_Hist, NbPRUJ1, NbPRUJ2, NbPRUHist) |
Partager