holiday$age <- as.factor(holiday$age)
holiday$departement <- as.factor(holiday$departement)
colnames(holiday)[6:7]=c("SPC.father", "SPC.mother")
colnames(holiday)[8]="nb.holidays.more.than.three.days"
holiday$nb.holidays.more.than.three.days <- as.factor(holiday$nb.holidays.more.than.three.days)
colnames(holiday)[9:13]=c("holidays.for.four.to.six.days", "holidays.for.one.week", "holidays.for.two.weeks", "holidays.for.three.weeks", "holidays.for.one.month.and.more")
colnames(holiday)[14:17]=c("with.family", "with.friends", "with.boy/girlfriend", "alone")
colnames(holiday)[18:22]=c("destination.area", "destination.border.area", "destination.France", "destination.Europe", "destination.out.of.Europe")
colnames(holiday)[23:27]=c("during.summer", "during.spring", "during.christmas", "during.february", "during.autumn")
colnames(holiday)[42:49]=c("favourite.place.hotel", "favourite.place.bungalow", "favourite.place.rented.flat", "favourite.place.rented.house", "favourite.place.campsite", "favourite.place.family.house", "favourite.place.friends.house", "favourite.place.bed.and.breakfast")
colnames(holiday)[50:52]=c("trip.with.young.people", "trip.with.hotel.club", "trip.without.organization")
colnames(holiday)[53:65]=c("activity.job", "activity.visiting", "activity.strolling", "activity.outdoor.sports", "activity.winter.sports", "activity.beach.sports", "activity.sun.tanning", "activity.books", "activity.tabloid", "activity.tv", "activity.computer", "activity.theme.park", "activity.bathing")
colnames(holiday)[66:78]=c("favourite.activity.job", "favourite.activity.visiting", "favourite.activity.strolling", "favourite.activity.outdoor.sports", "favourite.activity.winter.sports", "favourite.activity.beach.sports", "favourite.activity.sun.tanning", "favourite.activity.books", "favourite.activity.tabloid", "favourite.activity.tv", "favourite.activity.computer", "favourite.activity.theme.park", "favourite.activity.bathing")
colnames(holiday)[79:89]=c("holiday.means.relaxation", "holiday.means.escape", "holiday.means.work", "holiday.means.boredom", "holiday.means.resourcefulness", "holiday.means.routine", "holiday.means.sport", "holiday.means.party", "holiday.means.meetings", "holiday.means.spending", "holiday.means.discoveries")
colnames(holiday)[90]="parents.financial.help"
colnames(holiday)[91]="work.to.earn.money"
colnames(holiday)[92:96]=c("holiday.idea.from.internet", "holiday.idea.from.agency", "holiday.idea.from.family", "holiday.idea.from.friends", "holiday.idea.from.oneself")

for(i in c(4, 9:21, 25:28, 30:33, 36:39, 41, 51, 53, 58, 63:65, 80, 82, 87, 89, 92:95)){
	levels(holiday[,i])=c("yes", "no")
}
for(i in c(5, 22:24, 29, 34:35, 40, 42:50, 52, 54:57, 59:62, 66:79, 81, 83:86, 88, 96)){
	levels(holiday[,i])=c("no", "yes")
}

summary(holiday)

library(FactoMineR)
library(dynGraph)
library(EnQuireR)

holiday.MCA<-holiday[, c("nb.holidays.more.than.three.days", "holidays.for.four.to.six.days", "holidays.for.one.week", "holidays.for.two.weeks",
 "holidays.for.three.weeks", "holidays.for.one.month.and.more", "with.family", "with.friends", "with.boy/girlfriend", "alone",
 "destination.area", "destination.border.area", "destination.France", "destination.Europe", "destination.out.of.Europe", "during.summer",
 "during.spring", "during.christmas", "during.february", "during.autumn", "mountain", "seaside", "countryside",
 "city", "hotel", "bungalow", "rented.flat", "rented.house", "campsite", "family.house",
 "friend.house", "bed.and.breakfast", "youth.hostel", "couchsurfing", "trip.with.young.people", 
"trip.with.hotel.club", "trip.without.organization", "activity.job", "activity.visiting", "activity.strolling", "activity.outdoor.sports", "activity.winter.sports",
 "activity.beach.sports", "activity.sun.tanning", "activity.books", "activity.tabloid", "activity.tv", "activity.computer", "activity.theme.park", "activity.bathing", "parents.financial.help", 
"work.to.earn.money", "holiday.idea.from.internet", "holiday.idea.from.agency", "holiday.idea.from.family", "holiday.idea.from.friends", "holiday.idea.from.oneself", "sex", "age", "departement", "little.siblings", "older.siblings", 
"SPC.father", "SPC.mother", "favourite.place.hotel", "favourite.place.bungalow", "favourite.place.rented.flat", "favourite.place.rented.house", "favourite.place.campsite", "favourite.place.family.house", "favourite.place.friends.house", 
"favourite.place.bed.and.breakfast", "favourite.activity.job", "favourite.activity.visiting", "favourite.activity.strolling", "favourite.activity.outdoor.sports", "favourite.activity.winter.sports", "favourite.activity.beach.sports", 
"favourite.activity.sun.tanning", "favourite.activity.books", "favourite.activity.tabloid", "favourite.activity.tv", "favourite.activity.computer", "favourite.activity.theme.park", "favourite.activity.bathing", "holiday.means.relaxation", 
"holiday.means.escape", "holiday.means.work", "holiday.means.boredom", "holiday.means.resourcefulness", "holiday.means.routine", "holiday.means.sport", "holiday.means.party", "holiday.means.meetings", "holiday.means.spending", "holiday.means.discoveries")]

res<-MCA(holiday.MCA, ncp=5, quali.sup=58:96, graph = FALSE, level.ventil=0.05)
plot.MCA(res, axes=c(1, 2), col.ind="black", col.ind.sup="blue", col.var="darkred", col.quali.sup="darkgreen", label=c("ind.sup", "quali.sup", "var", "quanti.sup"), invisible=c("ind"), title=" Active factors")
plot.MCA(res, axes=c(1, 2), choix="var", col.var="darkred", col.quali.sup="darkgreen", label=c("quali.sup", "var"), invisible=c(""),cex=0.7)

#dimension description
dimdesc(res)

#Hierarchical Clustering on the result of the MCA where we keep the 27 first dimensions (80% total inertia)
res<-MCA(holiday.MCA, ncp=27, quali.sup=58:96, graph = FALSE, level.ventil=0.05)
res.hcpc<-HCPC(res)

#dynGraph
dynGraph(res)

#EnQuireR
ENMCA(holiday.MCA, quali.sup=58:96, level.ventil=0.05, report=TRUE)
graphics.off()