Évolution des appuis à la souveraineté
Problématique
La page qui suivante reproduit le code complet permettant d'estimer l'appui à la souveraineté du Québec de 1968 à 2015 en fonction de la génération, du revenu et de la scolarité des répondants. Les données sont tirées des Études électorales canadiennes. Attention, le code n'a pas été parfaitement harmonisé.

Nous avons choisi de travailler avec ces données d’abord parce qu’elles sont réputées fiables, mais aussi parce qu’elles sont accessibles et ouvertes pour la période 1968-2015. Chaque fois, nous avons retenu uniquement les répondants francophones du Québec. Nous avons pu ainsi faire des nouveaux traitements qui sont disponibles à l’adresse suivante. Après l’élimination de la non-réponse sur la question de la souveraineté (y compris les répondants discrets), les nombres de répondants sont ceux-ci : 493 en 1968; 264 en 1980; 491 en 1988; 358 en 1993; 779 en 1997; 431 en 2008, 918 en 2011 et 2012 en 2015. Sauf pour 1980, la marge va de plus ou moins 5 pourcent à plus ou moins 2 pourcent, 19 fois sur 20. Nous ne prétendons pas que l'image de l'évolution des appuis soit exacte. Malgré, les incontestables limites, on peut néanmoions suivre les tendances générales
Pré-requis pour l'ensemble des années


#Années de passage d'une génération à l'autre



#Par exemple les gens du baby-boom sont nés pendant la période 1946 -1965, exclusivement.
#Les pré avant 1946
#Les x après de 1966 à 1981
#Les milléniaux de 1982 à maintenant.

#Pour les babyboomers, la division retenu et celle de Statistique Canada: http://www12.statcan.gc.ca/census-recensement/2011/as-sa/98-311-x/98-311-x2011003_2-fra.cfm

Ensemble des années

library(questionr)

#Année 1968
#Récupération du fichier du CES 1968

genDate1=1946
genDate2=1966
genDate3=1982

base1968=read.csv("http://dimension.usherbrooke.ca/voute/CES-E-1968_F1.csv", header=TRUE, encoding="latin1")

#Vérication du n, et du nom des variables
#dim(base1968)
#names(base1968)
#table(base1968$age)


#Libellé des variables
#var001 #Québec 4
#var357: 2 pour lfrancais
#var306 # en faveur de la souveraineté du québec 1 et 2
#var404 #income 0 10000 14999 468 16.9%
#1 UNDER 1000 79 2.9%
#2 1000 1999 216 7.8%
#3 2000 2999 214 7.7%
#4 3000 3999 264 9.5%
#5 4000 4999 297 10.7%
#6 5000 5999 344 12.4%
#7 6000 6999 317 11.5%
#8 7000 7999 268 9.7%
#9 8000 9999

#Sélection pour ne retenir que les francophones du Québec
base1968Qc=subset(base1968, base1968$var001==4) #Sélection du Québec seulement
dim(base1968Qc)
base1968Qc=subset(base1968Qc, base1968Qc$var357==2) #Sélection des francophones seulement
dim(base1968Qc)
base1968Qc=subset(base1968Qc, base1968Qc$age<100)
table(base1968Qc$age)

#Construction de la variable génération
base1968Qc$ageRepondant=1968-base1968Qc$age

base1968Qc$generationRepondantGr[base1968Qc$ageRepondant < genDate1]="1-PRE"
base1968Qc$generationRepondantGr[base1968Qc$ageRepondant >=genDate1 & base1968Qc$ageRepondant < genDate2]="2-BB"
base1968Qc$generationRepondantGr[base1968Qc$ageRepondant >=genDate2 & base1968Qc$ageRepondant < genDate3]="3-X"
base1968Qc$generationRepondantGr[base1968Qc$ageRepondant >=genDate3]="4-MIL"


table(base1968Qc$generationRepondantGr)


#Construction de la variable revenu
base1968Qc$revenu[base1968Qc$var404==1]="1-DÉFAVORISEE"
base1968Qc$revenu[base1968Qc$var404==2]="1-DÉFAVORISEE"
base1968Qc$revenu[base1968Qc$var404==3]="1-DÉFAVORISEE"
base1968Qc$revenu[base1968Qc$var404==4]="1-DÉFAVORISEE"
base1968Qc$revenu[base1968Qc$var404==5]="2-MOYENNE"
base1968Qc$revenu[base1968Qc$var404==6]="2-MOYENNE"
base1968Qc$revenu[base1968Qc$var404==7]="3-FAVORISEE"
base1968Qc$revenu[base1968Qc$var404==8]="3-FAVORISEE"
base1968Qc$revenu[base1968Qc$var404==9]="3-FAVORISEE"
table(base1968Qc$revenu)

#Construction de la variable souveraineté 1: oui, 2 discretion, 3 non.
base1968Qc$souverainete[base1968Qc$var306==1]="1-POUR"
base1968Qc$souverainete[base1968Qc$var306==2]="1-POUR"
base1968Qc$souverainete[base1968Qc$var306==4]="3-CONTRE"
base1968Qc$souverainete[base1968Qc$var306==5]="3-CONTRE"

base1968Qc$scolarite[base1968Qc$var334<=12]="NON-UNIVERSITAIRES"
base1968Qc$scolarite[base1968Qc$var334>=13]="UNIVERSITAIRES"
table(base1968Qc$scolarite)


#Croisement de la souveraineté par le revenu, avec test
cprop(table(base1968Qc$souverainete, base1968Qc$revenu))
T1968RCSQ=chisq.test(table(base1968Qc$souverainete, base1968Qc$revenu))
T1968RA=addmargins(table(base1968Qc$souverainete, base1968Qc$revenu))
T1968R=round(cprop(table(base1968Qc$souverainete, base1968Qc$revenu)))


#Croisement de la souveraineté par la génération, avec test
table(base1968Qc$generationRepondantGr)
T1968GA=addmargins(table(base1968Qc$souverainete, base1968Qc$generationRepondantGr))
cprop(table(base1968Qc$souverainete, base1968Qc$generationRepondantGr))
T1968GCSQ=chisq.test(table(base1968Qc$souverainete, base1968Qc$generationRepondantGr))
#Construction du tableau finale
T1968G=round(cprop(table(base1968Qc$souverainete, base1968Qc$generationRepondantGr)))

#Croisement de la souveraineté par la scolarité, avec test
cprop(table(base1968Qc$souverainete, base1968Qc$scolarite))
T1968SCSQ=chisq.test(table(base1968Qc$souverainete, base1968Qc$scolarite))
T1968SA=addmargins(table(base1968Qc$souverainete, base1968Qc$scolarite))
T1968S=round(cprop(table(base1968Qc$souverainete, base1968Qc$scolarite)))
#Fin du traitement pour l'année 1968




#Année 1980, pour l'année 1974 trop peu de répondants sur la question de la souveraineté


base19741980=read.csv("http://dimension.usherbrooke.ca/voute/CES-E-1974-1980_F1.csv", header=TRUE, encoding="latin1")
#names(base19741980)
dim(base19741980)
base1980Qc=subset(base19741980,base19741980$v1510==2)
table(base1980Qc$v3007)#Vérification qu'ils sont bien de la période référendaire de 1980

table(base1980$v3011, base1980$v1510)

#Construction de la variable souveraineté
base1980Qc$souverainete[base1980Qc$v3011==1]="1-POUR"
base1980Qc$souverainete[base1980Qc$v3011==2]="1-POUR"
base1980Qc$souverainete[base1980Qc$v3011==3]="3-CONTRE"
base1980Qc$souverainete[base1980Qc$v3011==4]="3-CONTRE"

table(base1980Qc$souverainete)

#Construction de la variable age
table(base1980Qc$v3039)
base1980Qc=subset(base1980Qc,base1980Qc$v3039!=0)
base1980Qc$ageRepondant=1980-base1980Qc$v3039
base1980Qc$generationRepondantGr[base1980Qc$ageRepondant < genDate1]="1-PRE"
base1980Qc$generationRepondantGr[base1980Qc$ageRepondant >= genDate1 & base1980Qc$ageRepondant < genDate2]="2-BB"
base1980Qc$generationRepondantGr[base1980Qc$ageRepondant >= genDate2 & base1980Qc$ageRepondant < genDate3]="3-X"
base1980Qc$generationRepondantGr[base1980Qc$ageRepondant >= genDate3]="4-MIL"
table(base1980Qc$generationRepondantGr)

#Construction de la variable revenu
table(base1980Qc$v1516)
base1980Qc$revenu[base1980Qc$v1516==1]="1-DEFAVORISES"
base1980Qc$revenu[base1980Qc$v1516==2]="1-DEFAVORISES"
base1980Qc$revenu[base1980Qc$v1516==3]="1-DEFAVORISES"
base1980Qc$revenu[base1980Qc$v1516==4]="2-MOYENNE"
base1980Qc$revenu[base1980Qc$v1516==5]="2-MOYENNE"
base1980Qc$revenu[base1980Qc$v1516==6]="2-MOYENNE"
base1980Qc$revenu[base1980Qc$v1516==7]="3-FAVORISES"
base1980Qc$revenu[base1980Qc$v1516==8]="3-FAVORISES"
table(base1980Qc$revenu)

table(base1980Qc$v1502)

#Construction de la variable scolarité
base1980Qc$scolarite[base1980Qc$v1502<=21]="NON-UNIVERSITAIRES"
base1980Qc$scolarite[base1980Qc$v1502==22]="UNIVERSITAIRES"
base1980Qc$scolarite[base1980Qc$v1502==23]="UNIVERSITAIRES"
base1980Qc$scolarite[base1980Qc$v1502==24]="UNIVERSITAIRES"
table(base1980Qc$scolarite)

#Croisement des variables indépendantes par la variable dépendante
#Revenu avec souveraineté
table(base1980Qc$souverainete,base1980Qc$revenu)
T1980RA=addmargins(table(base1980Qc$souverainete,base1980Qc$revenu))
cprop(table(base1980Qc$souverainete,base1980Qc$revenu))
T1980RCSQ=chisq.test(table(base1980Qc$souverainete,base1980Qc$revenu))
T1980R=round(cprop(table(base1980Qc$souverainete,base1980Qc$revenu)))

#Revenu avec souveraineté
cprop(table(base1980Qc$souverainete,base1980Qc$generationRepondantGr))
table(base1980Qc$souverainete,base1980Qc$generationRepondantGr)
T1980GCSQ=chisq.test(table(base1980Qc$souverainete,base1980Qc$generationRepondantGr))
T1980GA=addmargins(table(base1980Qc$souverainete,base1980Qc$generationRepondantGr))
T1980G=round(cprop(table(base1980Qc$souverainete,base1980Qc$generationRepondantGr)))

#Scolarité avec souveraineté
cprop(table(base1980Qc$souverainete,base1980Qc$scolarite))
table(base1980Qc$souverainete,base1980Qc$scolarite)
T1980SCSQ=chisq.test(table(base1980Qc$souverainete,base1980Qc$scolarite))
T1980SA=addmargins(table(base1980Qc$souverainete,base1980Qc$scolarite))
T1980S=round(cprop(table(base1980Qc$souverainete,base1980Qc$scolarite)))
#Fin du traitement pour l'année 1980




#Récupération du fichier du CES 1988
base1988=read.csv2("http://dimension.usherbrooke.ca/voute/ces1988.csv", header=TRUE, encoding="latin1")

#Vérication du n, et du nom des variables
dim(base1988)
#names(base1988)

#Sélection pour ne retenir que les francophones du Québec
base1988Qc=subset(base1988, base1988$province==4) #Sélection du Québec seulement
dim(base1988Qc)
base1988Qc=subset(base1988Qc, base1988Qc$n17==1) #Sélection des francophones seulement
dim(base1988Qc)
base1988Qc=subset(base1988Qc, base1988Qc$n1<9000)

#Construction de la variable génération
base1988Qc$generationRepondantGr[base1988Qc$n1 < genDate1]="1-PRE"
base1988Qc$generationRepondantGr[base1988Qc$n1 >= genDate1 & base1988Qc$n1 < genDate2]="2-BB"
base1988Qc$generationRepondantGr[base1988Qc$n1 >= genDate2 & base1988Qc$n1 < genDate3]="3-X"
base1988Qc$generationRepondantGr[base1988Qc$n1 >= genDate3]="4-MIL"


#Construction de la variable revenu
#var404 #income 0 10000 14999 468 16.9%
#1 UNDER 1000 79 2.9%
#2 1000 1999 216 7.8%
#3 2000 2999 214 7.7%
#4 3000 3999 264 9.5%
#5 4000 4999 297 10.7%
#6 5000 5999 344 12.4%
#7 6000 6999 317 11.5%
#8 7000 7999 268 9.7%
#9 8000 9999
base1988Qc$revenu[base1988Qc$n19==1]="1-DÉFAVORISEE"
base1988Qc$revenu[base1988Qc$n19==2]="1-DÉFAVORISEE"
base1988Qc$revenu[base1988Qc$n19==3]="2-MOYENNE"
base1988Qc$revenu[base1988Qc$n19==4]="2-MOYENNE"
base1988Qc$revenu[base1988Qc$n19==5]="2-MOYENNE"
base1988Qc$revenu[base1988Qc$n19==6]="3-FAVORISEE"
base1988Qc$revenu[base1988Qc$n19==7]="3-FAVORISEE"
base1988Qc$revenu[base1988Qc$n19==8]="3-FAVORISEE"
base1988Qc$revenu[base1988Qc$n19==9]="3-FAVORISEE"

table(base1988Qc$revenu)

#Construction de la variable souveraineté 1: oui, 2 discretion, 3 non.
#b10
#What is your opinion on Quebec independence? Do you support it completely, support it somewhat, oppose it somewhat, oppose it completely or don't you really have an opinion on it?

#1 Support Completely 72 8,7%
#2 Support Somewhat 162 19,6%
#3 Oppose Somewhat 139 16,8%
#4 Oppose Completely 254 30,8%
#8 No Opinion 199 24,1%
#9 Refused
base1988Qc$souverainete[base1988Qc$b10==1]="1-POUR"
base1988Qc$souverainete[base1988Qc$b10==2]="1-POUR"
base1988Qc$souverainete[base1988Qc$b10==3]="3-CONTRE"
base1988Qc$souverainete[base1988Qc$b10==4]="3-CONTRE"
table(base1988Qc$souverainete)


#Construction de la variable scolarité
table(base1988Qc$n3)
base1988Qc$scolarite[base1988Qc$n3<=7]="NON-UNIVERSITAIRES"
base1988Qc$scolarite[base1988Qc$n3>=8 & base1988Qc$n3<12]="UNIVERSITAIRES"

#Croisement de la souveraineté par le revenu, avec test
cprop(table(base1988Qc$souverainete, base1988Qc$revenu))
T1988RCSQ=chisq.test(table(base1988Qc$souverainete, base1988Qc$revenu))
T1988RA=addmargins(table(base1988Qc$souverainete, base1988Qc$revenu))
T1988R=round(cprop(table(base1988Qc$souverainete, base1988Qc$revenu)))

#Croisement de la souveraineté par la génération, avec test
table(base1988Qc$generationRepondantGr)
T1988GA=addmargins(table(base1988Qc$souverainete, base1988Qc$generationRepondantGr))
cprop(table(base1988Qc$souverainete, base1988Qc$generationRepondantGr))
T1988GCSQ=chisq.test(table(base1988Qc$souverainete, base1988Qc$generationRepondantGr))

#Construction du tableau finale
T1988G=round(cprop(table(base1988Qc$souverainete, base1988Qc$generationRepondantGr)))

#Croisement de la souveraineté par la scolarité, avec test
cprop(table(base1988Qc$souverainete, base1988Qc$scolarite))
T1988SCSQ=chisq.test(table(base1988Qc$souverainete, base1988Qc$scolarite))
T1988SA=addmargins(table(base1988Qc$souverainete, base1988Qc$scolarite))
T1988S=round(cprop(table(base1988Qc$souverainete, base1988Qc$scolarite)))

#Fin du traitement pour l'année 1988








## Année 1993
base1993=read.csv2("http://dimension.usherbrooke.ca/voute/ces1993.csv", header=TRUE, encoding="latin1")

base1993Qc=subset(base1993, base1993$cpsprov==4) #Sélection des répondants du Québec seulement
base1993Qc=subset(base1993Qc, base1993Qc$cpso15==1) #Sélection des répondant francophones
dim(base1993Qc)

#Construction de la variable age
table(base1993Qc$cpsage)
base1993Qc=subset(base1993Qc,base1993Qc$cpsage<9000) #élimination des valeurs manquantes ici 9999
base1993Qc$generationRepondantGr[base1993Qc$cpsage < genDate1]="1-PRE"
base1993Qc$generationRepondantGr[base1993Qc$cpsage >= genDate1 & base1993Qc$cpsage < genDate2]="2-BB"
base1993Qc$generationRepondantGr[base1993Qc$cpsage >= genDate2 & base1993Qc$cpsage < genDate3]="3-X"
base1993Qc$generationRepondantGr[base1993Qc$cpsage >= genDate3]="4-MIL"


#Construction de la variable génération
#Variable CPSG11: Opinion on Quebec Sovereignty
#Question littérale
#What is your opinion on Quebec sovereignty, that is, Quebec is no longer a part of Canada? Are you very favourable, somewhat favourable, somewhat opposed, or very opposed.
#0 Neither 25 2,6%
#1 Very favourable 160 16,7%
#3 Somewhat favour 234 24,4%
#5 Somewhat opposed 202 21,1%
#7 Very opposed 338 35,2%

base1993Qc$souverainete[base1993Qc$cpsg11==1]="1-POUR"
base1993Qc$souverainete[base1993Qc$cpsg11==3]="1-POUR"
base1993Qc$souverainete[base1993Qc$cpsg11==5]="3-CONTRE"
base1993Qc$souverainete[base1993Qc$cpsg11==7]="3-CONTRE"
table(base1993Qc$souverainete)

#Construction de la variable revenu
#Variable CPSO18A: Total Household Income> Category
#Question littérale
#We don't need the exact amount; could tell me which of these broad categories it falls into...
#1 <$20,000 225 33,4%
#2 $20,000-$29,000 150 22,3%
#3 $30,000-$39,000 85 12,6%
#4 $40,000-$49,000 69 10,2%
#5 $50,000-$59,000 46 6,8%
#6 $60,000-$69,000 22 3,3%
#7 $70,000-$79,000 14 2,1%
#8 $80,000-$89,000 14 2,1%
#9 $90,000-100,000 12 1,8%
#10 $100,000+ 37 5,5%


#À partir des catégories
base1993Qc$revenu[base1993Qc$cpso18a==1]="1-DEFAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18a==2]="1-DEFAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18a==3]="2-MOYENNE"
base1993Qc$revenu[base1993Qc$cpso18a==4]="2-MOYENNE"
base1993Qc$revenu[base1993Qc$cpso18a==5]="3-FAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18a==6]="3-FAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18a==7]="3-FAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18a==8]="3-FAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18a==9]="3-FAVORISEE"
#À partir des valeurs rapportées par les répondants
base1993Qc$revenu[base1993Qc$cpso18<30]="1-DEFAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18>=30 & base1993Qc$cpso18<50]="2-MOYENNE"
base1993Qc$revenu[base1993Qc$cpso18>=50 & base1993Qc$cpso18<80]="3-FAVORISEE"
base1993Qc$revenu[base1993Qc$cpso18>=80 & base1993Qc$cpso18<400 ]="3-FAVORISEE"
table(base1993Qc$revenu)



#Constrution de la variable scolarité
table(base1993Qc$cpso3)
base1993Qc$scolarite[base1993Qc$cpso3<=7]="NON-UNIVERSITAIRES"
base1993Qc$scolarite[base1993Qc$cpso3>=8 & base1993Qc$cpso3<12]="UNIVERSITAIRES"


#Croisement de la variable souveraineté par revenu et test
table(base1993Qc$souverainete,base1993Qc$revenu)
T1993RA=addmargins(table(base1993Qc$souverainete,base1993Qc$revenu))
cprop(table(base1993Qc$souverainete,base1993Qc$revenu))
T1993RCSQ=chisq.test(table(base1993Qc$souverainete,base1993Qc$revenu))
T1993R=round(cprop(table(base1993Qc$souverainete,base1993Qc$revenu)))


#Croisement de la variable souveraineté par génération et test
table(base1993Qc$generationRepondantGr)
table(base1993Qc$souverainete,base1993Qc$generationRepondantGr)
T1993GA=addmargins(table(base1993Qc$souverainete,base1993Qc$generationRepondantGr))
cprop(table(base1993Qc$souverainete,base1993Qc$generationRepondantGr))
T1993GCSQ=chisq.test(table(base1993Qc$souverainete,base1993Qc$generationRepondantGr))
T1993G=round(cprop(table(base1993Qc$souverainete,base1993Qc$generationRepondantGr)))

#Croisement de la variable souveraineté par génération et test
table(base1993Qc$scolarite)
table(base1993Qc$souverainete,base1993Qc$scolarite)
T1993SA=addmargins(table(base1993Qc$souverainete,base1993Qc$scolarite))
cprop(table(base1993Qc$souverainete,base1993Qc$scolarite))
T1993SCSQ=chisq.test(table(base1993Qc$souverainete,base1993Qc$scolarite))
T1993S=round(cprop(table(base1993Qc$souverainete,base1993Qc$scolarite)))

#Fin du traitement pour 1993







# année 1997
#Importation de la base de données du CES 1997
base1997=read.csv2("http://dimension.usherbrooke.ca/voute/ces1997.csv", header=TRUE, encoding="latin1")
dim(base1997)
#names(base1997)
table(base1997$province)

#Sélection des Québécois francophones
base1997Qc=subset(base1997, base1997$province==24) #Sélection du Québec seulement
base1997Qc=subset(base1997Qc, base1997Qc$cpsm14==1) #Sélection des francophones seulement
dim(base1997Qc)


#Construction de la variable génération
table(base1997Qc$cpsage)
base1997Qc=subset(base1997Qc,base1997Qc$cpsage<9000) #élimination des valeurs manquantes ici 9999
base1997Qc$generationRepondantGr[base1997Qc$cpsage < genDate1]="1-PRE"
base1997Qc$generationRepondantGr[base1997Qc$cpsage >= genDate1 & base1997Qc$cpsage < genDate2]="2-BB"
base1997Qc$generationRepondantGr[base1997Qc$cpsage >= genDate2 & base1997Qc$cpsage < genDate3]="3-X"
base1997Qc$generationRepondantGr[base1997Qc$cpsage >= genDate3]="4-MIL"
table(base1997Qc$generationRepondantGr)

#Construction de la variable revenu
base1997Qc$revenu[base1997Qc$cpsm16a==1]="1-DEFAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16a==2]="1-DEFAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16a==3]="2-MOYENNE"
base1997Qc$revenu[base1997Qc$cpsm16a==4]="2-MOYENNE"
base1997Qc$revenu[base1997Qc$cpsm16a==5]="3-FAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16a==6]="3-FAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16a==7]="3-FAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16a==8]="3-FAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16a==9]="3-FAVORISEE"

base1997Qc$revenu[base1997Qc$cpsm16 <30]="1-DEFAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16 >=30 & base1997Qc$cpsm16 < 50]="2-MOYENNE"
base1997Qc$revenu[base1997Qc$cpsm16 >=50 & base1997Qc$cpsm16 < 80]="3-FAVORISEE"
base1997Qc$revenu[base1997Qc$cpsm16 >=80 & base1997Qc$cpsm16 < 400 ]="3-FAVORISEE"
table(base1997Qc$revenu)

#Construction de la variable souveraineté
table(base1997$cpsj3a)
#j3a
#Are you very favourable, somewhat favourable, somewhat opposed
#[n] [r]   or very opposed to Québec sovereignty
#<1> very favourable
#<2> somewhat favourable
#<3> somewhat opposed
#<4> very opposed
#<0> neither favourable or opposed
#<8> don't know   <9> refused

base1997Qc$souverainete[base1997Qc$cpsj3a==1]="1-POUR"
base1997Qc$souverainete[base1997Qc$cpsj3a==2]="1-POUR"
base1997Qc$souverainete[base1997Qc$cpsj3a==3]="3-CONTRE"
base1997Qc$souverainete[base1997Qc$cpsj3a==4]="3-CONTRE"


#Construction de la variable scolarité
table(base1997Qc$cpsm3)
base1997Qc$scolarite[base1997Qc$cpsm3<=7]="NON-UNIVERSITAIRES"
base1997Qc$scolarite[base1997Qc$cpsm3>=8 & base1997Qc$cpsm3<12]="UNIVERSITAIRES"


#Croisement de la variable souveraineté par la variable revenu
table(base1997Qc$souverainete,base1997Qc$revenu)
T1997RA=addmargins(table(base1997Qc$souverainete,base1997Qc$revenu))
cprop(table(base1997Qc$souverainete,base1997Qc$revenu))
T1997RCSQ=chisq.test(base1997Qc$souverainete, base1997Qc$revenu)
T1997R=round(cprop(table(base1997Qc$souverainete,base1997Qc$revenu)))


#Croisement de la variable souveraineté par la variable génération
table(base1997Qc$souverainete,base1997Qc$generationRepondantGr)
T1997GA=addmargins(table(base1997Qc$souverainete,base1997Qc$generationRepondantGr))
cprop(table(base1997Qc$souverainete,base1997Qc$generationRepondantGr))
T1997GCSQ=chisq.test(base1997Qc$souverainete, base1997Qc$generationRepondantGr)
T1997G=round(cprop(table(base1997Qc$souverainete,base1997Qc$generationRepondantGr)))

#Croisement de la variable souveraineté par la variable scolarité
table(base1997Qc$souverainete,base1997Qc$scolarite)
T1997SA=addmargins(table(base1997Qc$souverainete,base1997Qc$scolarite))
cprop(table(base1997Qc$souverainete,base1997Qc$scolarite))
T1997SCSQ=chisq.test(base1997Qc$souverainete, base1997Qc$scolarite)
T1997S=round(cprop(table(base1997Qc$souverainete,base1997Qc$scolarite)))

#Fin du traitement pour 1997






#Élection fédérale de 2008
base2008=read.csv2("http://dimension.usherbrooke.ca/voute/ces2008.csv", header=TRUE, encoding="latin1")


#Sélection des répondants du Québec et francophones
base2008Qc=subset(base2008, base2008$ces08_PROVINCE==24)
base2008Qc=subset(base2008Qc, base2008$ces08_CPS_S17==5)
dim(base2008Qc)

#Construction de la variable génération
table(base2008Qc$YEARofBIRTH)
base2008Qc=subset(base2008Qc,base2008Qc$YEARofBIRTH<9000) #élimination des valeurs manquantes ici 9999
base2008Qc$generationRepondantGr[base2008Qc$YEARofBIRTH < genDate1]="1-PRE"
base2008Qc$generationRepondantGr[base2008Qc$YEARofBIRTH >= genDate1 & base2008Qc$YEARofBIRTH < genDate2]="2-BB"
base2008Qc$generationRepondantGr[base2008Qc$YEARofBIRTH >= genDate2 & base2008Qc$YEARofBIRTH < genDate3]="3-X"
base2008Qc$generationRepondantGr[base2008Qc$YEARofBIRTH >= genDate3]="4-MIL"

#Construction de la variable revenu
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==1]="1-DEFAVORISEE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==2]="1-DEFAVORISEE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==3]="1-DEFAVORISEE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==4]="2-MOYENNE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==5]="2-MOYENNE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==6]="2-MOYENNE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==7]="2-MOYENNE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==8]="3-FAVORISEE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18B==9]="3-FAVORISEE"

base2008Qc$revenu[base2008Qc$ces08_CPS_S18A<40]="1-DEFAVORISEE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18A >=40 & base2008Qc$ces08_CPS_S18A < 60]="2-MOYENNE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18A >=60 & base2008Qc$ces08_CPS_S18A < 80]="2-MOYENNE"
base2008Qc$revenu[base2008Qc$ces08_CPS_S18A >=80 & base2008Qc$ces08_CPS_S18A < 400 ]="3-FAVORISEE"
table(base2008Qc$revenu)

#Variable ces08_CPS_S18B_ISR: ISR Income Categories
#Valeurs
#Catégories
#N
#1 Less than $20,000 70 14,6%
#2 $20,000-$29,999 65 13,6%
#3 $30,000-$39,999 56 11,7%
#4 $40,000-$49,999 35 7,3%
#5 $50,000-$59,999 57 11,9%
#6 $60,000-$69,999 36 7,5%
##7 $70,000-$79,999 27 5,6%
#8 $80,000-$89,999 35 7,3%
#9 $90,000-$99,999 15 3,1%
#10 $100,000-$109,999 18 3,8%
#11 $110,000-119,999 9 1,9%
#12 More than $120,000 56 11,7%
#98 don't know 130
#99 refused 320

#Construction de la variable scolarité
base2008Qc$souverainete[base2008Qc$ces08_CPS_Q9==1]="1-POUR"
base2008Qc$souverainete[base2008Qc$ces08_CPS_Q9==3]="1-POUR"
#base2008Qc$souverainete[base2008Qc$ces08_CPS_Q9==8]="2-DISCRET"
#base2008Qc$souverainete[base2008Qc$ces08_CPS_Q9==9]="2-DISCRET"
base2008Qc$souverainete[base2008Qc$ces08_CPS_Q9==5]="3-CONTRE"
base2008Qc$souverainete[base2008Qc$ces08_CPS_Q9==7]="3-CONTRE"
table(base2008Qc$ces08_CPS_S3)

base2008Qc$scolarite[base2008Qc$ces08_CPS_S3<=7]="NON-UNIVERSITAIRES"
base2008Qc$scolarite[base2008Qc$ces08_CPS_S3>=8 & base2008Qc$ces08_CPS_S3<12]="UNIVERSITAIRES"
table(base2008Qc$scolarite)


#Croisement de la variable souveraineté par revenu
cprop(table(base2008Qc$souverainete,base2008Qc$revenu))
T2008RA=addmargins(table(base2008Qc$souverainete,base2008Qc$revenu))
T2008RCSQ=chisq.test(table(base2008Qc$souverainete,base2008Qc$revenu))
T2008R=round(cprop(table(base2008Qc$souverainete,base2008Qc$revenu)))


#Croisement de la variable souveraineté par âge
table(base2008Qc$souverainete,base2008Qc$generationRepondantGr)
T2008GA=addmargins(table(base2008Qc$souverainete,base2008Qc$generationRepondantGr))
cprop(table(base2008Qc$souverainete,base2008Qc$generationRepondantGr))
T2008GCSQ=chisq.test(table(base2008Qc$souverainete,base2008Qc$generationRepondantGr))
T2008G=round(cprop(table(base2008Qc$souverainete,base2008Qc$generationRepondantGr)))

#Croisement de la variable souveraineté par scolarité
table(base2008Qc$souverainete,base2008Qc$scolarite)
T2008SA=addmargins(table(base2008Qc$souverainete,base2008Qc$scolarite))
cprop(table(base2008Qc$souverainete,base2008Qc$scolarite))
T2008SCSQ=chisq.test(table(base2008Qc$souverainete,base2008Qc$scolarite))
T2008S=round(cprop(table(base2008Qc$souverainete,base2008Qc$scolarite)))
#Fin du traitement pour 2008




#Année 2011
base2011=read.csv2("http://dimension.usherbrooke.ca/voute/ces2011guay.csv", header=TRUE, encoding="latin1")

#names(base2011)

base2011Qc=subset(base2011, base2011$PROVINCE==24)
base2011Qc=subset(base2011Qc, base2011Qc$CPS11_90==5)
dim(base2011Qc)


table(base2011Qc$CPS11_78)
base2011Qc=subset(base2011Qc, base2011Qc$CPS11_78<9000) #élimination des gens nés en l'an 10000 (certainement la  non-réponse)

base2011Qc$generationRepondantGr[base2011Qc$CPS11_78 < genDate1]="1-PRE"
base2011Qc$generationRepondantGr[base2011Qc$CPS11_78 >= genDate1 & base2011Qc$CPS11_78 < genDate2]="2-BB"
base2011Qc$generationRepondantGr[base2011Qc$CPS11_78 >= genDate2 & base2011Qc$CPS11_78 < genDate3]="3-X"
base2011Qc$generationRepondantGr[base2011Qc$CPS11_78 >= genDate3]="4-MIL"
table(base2011Qc$generationRepondantGr)


#Construction de la variable souveraineté
#Êtes-vous très favorable, plutôt favorable, plutôt opposé(e)  ou très  #opposé(e) a la souveraineté du Québec, c'est-a-dire que le Québec NE FAIT #PLUS PARTIE DU CANADA?

#1 Très favorable
#3 Plutôt favorable
#5 Plutôt  opposé
#7  Très  opposé
#98 Ne  sais  pas
#99 Refus


table(base2011$CPS11_75)
base2011Qc$souverainete[base2011Qc$CPS11_75==1]="1-POUR"
base2011Qc$souverainete[base2011Qc$CPS11_75==3]="1-POUR"
#base2011Qc$souverainete[base2011Qc$CPS11_75==98]="2-DISCRET"
#base2011Qc$souverainete[base2011Qc$CPS11_75==99]="2-DISCRET"
base2011Qc$souverainete[base2011Qc$CPS11_75==5]="3-CONTRE"
base2011Qc$souverainete[base2011Qc$CPS11_75==7]="3-CONTRE"
table(base2011Qc$souverainete)

#Construction de la variable revenu
base2011Qc$revenu[base2011Qc$income.gr==1]="1-DEFAVORISES"
base2011Qc$revenu[base2011Qc$income.gr==2]="2-MOYENNE"
base2011Qc$revenu[base2011Qc$income.gr==3]="3-FAVORISES"
base2011Qc$revenu[base2011Qc$income.gr==4]="3-FAVORISES"
base2011Qc$revenu[base2011Qc$income.gr==5]="3-FAVORISES"

table(base2011Qc$CPS11_92)
CPS11_93

#Construction de la variable revenu
base2011Qc$revenu[base2011Qc$CPS11_93==1]="1-DEFAVORISEE"
base2011Qc$revenu[base2011Qc$CPS11_93==2]="2-MOYENNE"
base2011Qc$revenu[base2011Qc$CPS11_93==3]="2-MOYENNE"
base2011Qc$revenu[base2011Qc$CPS11_93==4]="2-MOYENNE"
base2011Qc$revenu[base2011Qc$CPS11_93==5]="3-FAVORISEE"


base2011Qc$revenu[base2011Qc$CPS11_92<30]="1-DEFAVORISEE"
base2011Qc$revenu[base2011Qc$CPS11_92>=30 & base2011Qc$CPS11_92<90]="2-MOYENNE"
base2011Qc$revenu[base2011Qc$CPS11_92>=110 & base2011Qc$CPS11_92<400 ]="3-FAVORISEE"

table(base2011Qc$revenu)

#Construction de la variable scolarité
table(base2011Qc$CPS11_79)
base2011Qc$scolarite[base2011Qc$CPS11_79<=7]="NON-UNIVERSITAIRES"
base2011Qc$scolarite[base2011Qc$CPS11_79>=8 & base2011Qc$CPS11_79<12]="UNIVERSITAIRES"


#Croisement de la variable souveraineté par génération
table(base2011Qc$souverainete,base2011Qc$generationRepondantGr)
T2011GCSQ=chisq.test(table(base2011Qc$souverainete,base2011Qc$generationRepondantGr))
T2011GA=addmargins(table(base2011Qc$souverainete,base2011Qc$generationRepondantGr))
T2011G=round(cprop(table(base2011Qc$souverainete,base2011Qc$generationRepondantGr)))

#Croisement de la variable souveraineté par revenu
table(base2011Qc$souverainete,base2011Qc$revenu)
T2011RA=addmargins(table(base2011Qc$souverainete,base2011Qc$revenu))
cprop(table(base2011Qc$souverainete,base2011Qc$revenu))
T2011RCSQ=chisq.test(table(base2011Qc$souverainete,base2011Qc$revenu))
T2011R=round(cprop(table(base2011Qc$souverainete,base2011Qc$revenu)))

#Croisement de la variable souveraineté par scolarité
table(base2011Qc$souverainete,base2011Qc$scolarite)
T2011SA=addmargins(table(base2011Qc$souverainete,base2011Qc$scolarite))
cprop(table(base2011Qc$souverainete,base2011Qc$scolarite))
T2011SCSQ=chisq.test(table(base2011Qc$souverainete,base2011Qc$scolarite))
T2011S=round(cprop(table(base2011Qc$souverainete,base2011Qc$scolarite)))

#Fin du traitement pour l'année 2011




#Année 2015
base2015=read.csv2("http://dimension.usherbrooke.ca/voute/ces2015guay.csv", header=TRUE, encoding="latin1")
dim(base2015)
#names(base2015)
table(base2015$province)
table(base2015$age)

#Sélection des Québécois francophones
base2015Qc=subset(base2015, base2015$province==11)
base2015Qc=subset(base2015Qc, base2015Qc$first.lang==2)
dim(base2015Qc)

#Construction de la variable génération
base2015Qc$ageRepondant=base2015Qc$age
table(base2015Qc$ageRepondant)
base2015Qc=subset(base2015Qc, base2015Qc$ageRepondant>=1900) #élimination des gens nés en l'an 10000 (certainement la  non-réponse)
base2015Qc$generationRepondantGr[base2015Qc$ageRepondant < genDate1]="1-PRE"
base2015Qc$generationRepondantGr[base2015Qc$ageRepondant >= genDate1 & base2015Qc$ageRepondant < genDate2]="2-BB"
base2015Qc$generationRepondantGr[base2015Qc$ageRepondant >= genDate2 & base2015Qc$ageRepondant < genDate3]="3-X"
base2015Qc$generationRepondantGr[base2015Qc$ageRepondant >= genDate3]="4-MIL"
table(base2015Qc$generationRepondantGr)


#Construction de la variable revenu
table(base2015Qc$income.grp)
#?	Moins de $30,000 (1)
#?	$30,000 à $60,000 (2)
#?	$60,000 à $90,000 (3)
#?	$90,000 à $110,000 (4)
#?	Plus de $110,000 (5)
#?	Ne sais pas (6)
#?	Refus (7)


base2015Qc$revenu[base2015Qc$income.gr==1]="1-DEFAVORISES"
base2015Qc$revenu[base2015Qc$income.gr==2]="2-MOYENNE"
base2015Qc$revenu[base2015Qc$income.gr==3]="3-FAVORISES"
base2015Qc$revenu[base2015Qc$income.gr==4]="3-FAVORISES"
base2015Qc$revenu[base2015Qc$income.gr==5]="3-FAVORISES"

#Construction de la variable scolarité
table(base2015Qc$education)
base2015Qc$scolarite[base2015Qc$education<=7]="NON-UNIVERSITAIRES"
base2015Qc$scolarite[base2015Qc$education>=8 & base2015Qc$education<12]="UNIVERSITAIRES"


#sov Êtes-vous très favorable, plutôt favorable, plutôt opposé, ou très opposé à la souveraineté du Québec, c'est-à-dire que le Québec NE FAIT PLUS PARTIE DU CANADA? [AA]
#?	Très favorable à la souveraineté du Québec (1)
#?	Plutôt favorable à la souveraineté du Québec (2)
#?	Plutôt opposé à la souveraineté du Québec (3)
#?	Très opposé à la souveraineté du Québec (4)
#?	Ne sais pas (5)
#?	Refus (6)


base2015Qc$souverainete[base2015Qc$sov==1]="1-POUR"
base2015Qc$souverainete[base2015Qc$sov==2]="1-POUR"
#base2015Qc$souverainete[base2015Qc$sov==5]="2-DISCRET"
#base2015Qc$souverainete[base2015Qc$sov==6]="2-DISCRET"
base2015Qc$souverainete[base2015Qc$sov==3]="3-CONTRE"
base2015Qc$souverainete[base2015Qc$sov==4]="3-CONTRE"


#Croisement de la variable souveraineté par revenu
table(base2015Qc$souverainete,base2015Qc$revenu)
T2015RA=addmargins(table(base2015Qc$souverainete,base2015Qc$revenu))
cprop(table(base2015Qc$souverainete,base2015Qc$revenu))
T2015RCSQ=chisq.test(table(base2015Qc$souverainete,base2015Qc$revenu))
T2015R=round(cprop(table(base2015Qc$souverainete,base2015Qc$revenu)))

#Croisement de la variable souveraineté par génération
table(base2015Qc$souverainete,base2015Qc$generationRepondantGr)
T2015GCSQ=chisq.test(table(base2015Qc$souverainete,base2015Qc$generationRepondantGr))
T2015GA=addmargins(table(base2015Qc$souverainete,base2015Qc$generationRepondantGr))
T2015G=round(cprop(table(base2015Qc$souverainete,base2015Qc$generationRepondantGr)))

#Croisement de la variable souveraineté par scolarité
table(base2015Qc$souverainete,base2015Qc$scolarite)
T2015SCSQ=chisq.test(table(base2015Qc$souverainete,base2015Qc$scolarite))
T2015SA=addmargins(table(base2015Qc$souverainete,base2015Qc$scolarite))
T2015S=round(cprop(table(base2015Qc$souverainete,base2015Qc$scolarite)))
addmargins(table(base2015Qc$souverainete,base2015Qc$scolarite,base2015Qc$generationRepondantGr))



#Liste complète des tableaux générés (G pour génération en %, GA, pour valeur absolues, GCSQ, pour calcul du Chi carré)

T1968G
T1968GA
T1968GCSQ

T1980G
T1980GA
T1980GCSQ

T1988G
T1988GA
T1988GCSQ

T1993G
T1993GA
T1993GCSQ

T1997G
T1997GA
T1997GCSQ

T2008G
T2008GA
T2008GCSQ

T2011G
T2011GA
T2011GCSQ

T2015G
T2015GA
T2015GCSQ

T1968R
T1968RA
T1968RCSQ

T1980R
T1980RA
T1980RCSQ

T1988R
T1988RA
T1988RCSQ

T1993R
T1993RA
T1993RCSQ

T1997R
T1997RA
T1997RCSQ

T2008R
T2008RA
T2008RCSQ

T2011R
T2011RA
T2011RCSQ

T2015R
T2015RA
T2015RCSQ



T1968S
T1968SA
T1968SCSQ

T1980S
T1980SA
T1980SCSQ

T1988S
T1988SA
T1988SCSQ

T1993S
T1993SA
T1993SCSQ

T1997S
T1997SA
T1997SCSQ

T2008S
T2008SA
T2008SCSQ

T2011S
T2011SA
T2011SCSQ

T2015S
T2015SA
T2015SCSQ








Fichier Excel reproduisant les résultats finaux pour les classes et les générations