# association_ex_corr.r

source("mesofas.r")
library(descr)

# Exercice 2.1
load("AccCeinture.rda")
load("AccVitesse.rda")

mesofas.chi2(AccCeinture$Gravité,AccCeinture$Ceinture,table=TRUE)
mesofas.chi2(AccVitesse$Gravité,AccVitesse$Vitesse,table=TRUE)

mesofas.u(AccCeinture$Gravité,AccCeinture$Ceinture)
mesofas.u(AccVitesse$Gravité,AccVitesse$Vitesse)

mesofas.d(AccCeinture$Gravité,AccCeinture$Ceinture)
mesofas.d(AccVitesse$Gravité,AccVitesse$Vitesse)


# Exercice 2.2
load("MatuStatut.rda")

# 1
mesofas.chi2(MatuStatut$Maturité,MatuStatut$Statut)
# 3
mesofas.V(MatuStatut$Maturité,MatuStatut$Statut)
# 4
mesofas.lambda(MatuStatut$Maturité,MatuStatut$Statut)


# Exercice 2.3
load("RevZone.rda")

# 1&2
crosstab(RevZone$Revenu,RevZone$Zone,plot=FALSE,expected=TRUE)

# 3
mesofas.phi(RevZone$Revenu,RevZone$Zone)
mesofas.C(RevZone$Revenu,RevZone$Zone)
mesofas.T(RevZone$Revenu,RevZone$Zone)
mesofas.V(RevZone$Revenu,RevZone$Zone)

# 4
mesofas.lambda(RevZone$Revenu,RevZone$Zone)


# Exercice 2.4
load("Cigognes.rda")
summary(Cigognes)
crosstab(Cigognes$Nids,Cigognes$Enfants,plot=FALSE)

# 1
mesofas.d(Cigognes$Nids,Cigognes$Enfants)

# 2
mesofas.eta(Cigognes$Nids,Cigognes$Enfants)
mesofas.eta(Cigognes$Enfants,Cigognes$Nids)

# 3
mesofas.spearman(Cigognes$Nids,Cigognes$Enfants)
mesofas.pearson(Cigognes$Nids,Cigognes$Enfants)