Mesure de dénombrement de la population de grues du Canada

Présentation du jeu de données

#---Specifier le répertoire de travai
setwd("C:\\Users\\Pounie\\Documents\\Universite_Laval_21-22\\Gestion\\JeuDeDonnees\\Donnees_30mars\\fichiersR_grues\\Analyse_grues")

#---Importer les données
grues<-read.table("C:/Users/Pounie/Documents/Universite_Laval_21-22/Gestion/JeuDeDonnees/Donnees_30mars/fichiersR_grues/denombrement_grues.txt",
                  header=TRUE)

Analyse

Régression linéaire simple

Transformations

#---Transformation log 10
grues$Nombre_individus.log<-log1p(grues$Nombre_individus)

Régression avec données transformées

#---Analyse de régression
reg.grues.log<-lm(Nombre_individus.log~Annee,data=grues)
#---Résulats de l'analyse
summary(reg.grues.log)
## 
## Call:
## lm(formula = Nombre_individus.log ~ Annee, data = grues)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.18429 -0.50730 -0.06138  0.58657  1.61939 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -0.20006    0.17113  -1.169    0.245    
## Annee        0.69659    0.02758  25.257   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7922 on 98 degrees of freedom
## Multiple R-squared:  0.8668, Adjusted R-squared:  0.8655 
## F-statistic: 637.9 on 1 and 98 DF,  p-value: < 2.2e-16
#---Visulation graphique des condtitions d'application
par(mfrow= c(2,2))
plot(reg.grues.log)

#--- Distance de cook
cookRegGrues<-cooks.distance(reg.grues.log)
plot(cookRegGrues)

summary(cookRegGrues)
##      Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
## 4.750e-06 1.140e-03 5.205e-03 1.005e-02 1.229e-02 7.744e-02
#---Graphique du nombre de grues en fonction du temps 
plot(grues$Nombre_individus.log~grues$Annee,
pch=20, ylab="Nombre de grues du Canada",
xlab="Temps écoulé (années)",
yaxt="n", main="Nombre de grues du Canada en fonction du temps (années)",cex.main=.75)

#---Ajustement des axes
axis(1,at= 1:10)

axis(side = 2,
     ## Rotate the labels.
     las = 2,
     ## Adjust the label position.
     mgp = c(3, 0.75, 0))

#---Ajout de la droite de régression
abline(reg.grues.log, col="red")