music$AGE <- recode(music$AGE,'15:18="15-18"; 19:25="19-25"; 26:35="26-35"; 36:200="36 and more"; ', as.factor.result=TRUE)
colnames(music)[6]="Musical.knowledge"
colnames(music)[7]="Budget"
colnames(music)[8]="Use.of.internet"
colnames(music)[9]="Play.an.instrument"
colnames(music)[10]="Listen.to.radio"
colnames(music)[11]="Access.to.music"
colnames(music)[16]="New.bands.from.friends"
music$New.bands.from.friends <- factor(music$New.bands.from.friends,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[17]="New.bands.from.websites"
music$New.bands.from.websites <- factor(music$New.bands.from.websites,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[18]="New.bands.from.radio"
music$New.bands.from.radio <- factor(music$New.bands.from.radio,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[19]="New.bands.from.tv.shows"
music$New.bands.from.tv.shows <- factor(music$New.bands.from.tv.shows,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[20]="New.bands.from.magazines"
music$New.bands.from.magazines <- factor(music$New.bands.from.magazines,labels=c('very few','few','average','quite a lot'))
colnames(music)[21]="New.bands.from.parents"
music$New.bands.from.parents <- factor(music$New.bands.from.parents,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[22]="New.bands.from.posters"
music$New.bands.from.posters <- factor(music$New.bands.from.posters,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[23]="New.bands.from.movies"
music$New.bands.from.movies <- factor(music$New.bands.from.movies,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[24]="New.bands.from.concerts"
music$New.bands.from.concerts <- factor(music$New.bands.from.concerts,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[25]="Go.to.enough.concerts"
music$go.to.concert.of.known.bands <- factor(music$go.to.concert.of.known.bands,labels=c('very few','few','average','quite a lot','a lot'))
colnames(music)[27]="Influenced.by.radio"

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

summary(music)

#MCA
music.MCA<-music[, c("Kind.of.music", "Musical.knowledge", "Kind.of.material", 
  "How.many.new.bands.by.month", "New.bands.from.friends", "New.bands.from.websites", "New.bands.from.radio", 
  "New.bands.from.tv.shows", "New.bands.from.magazines", "New.bands.from.parents", 
  "New.bands.from.posters", "New.bands.from.movies", "New.bands.from.concerts", 
  "go.to.concert.of.known.bands", "Budget", "listen.to.whole.album", "Why", 
  "Play.an.instrument", "Go.to.enough.concerts", "AGE", "SEX", "SPC")]

res<-MCA(music.MCA, ncp=5, quali.sup=16:22, level.ventil=0.05, graph = FALSE)
plot.MCA(res, axes=c(1, 2), col.ind="black", col.ind.sup="blue", 
  col.var="darkred", col.quali.sup="darkgreen", label=c("ind", "ind.sup", 
  "quali.sup", "var", "quanti.sup"), invisible=c(""))
plot.MCA(res, axes=c(1, 2), choix="var", col.var="darkred", 
  col.quali.sup="darkgreen", label=c("var", "quali.sup"), invisible=c(""))

#dimensions description
dimdesc(res,axes=c(1,2))

#Classification
res.hcpc = HCPC(res)
res.hcpc$desc.var

#dynGraph
dynGraph(res)

#EnQuireR (all the variables are active here since the algorithm for missing values 
#does not work with supplementary individuals yet)
res.enmca = ENMCA(music.MCA, level.ventil=0.05, report=TRUE)