Chapitre 6. Programmation en R
219
$divers$b
[1] "titi"
Il suffit de quelques lignes de code pour obtenir le r´ esultat, ce qui fait
de R un langage ayant de r´ eelles qualit´ es d’introspection. Et cela n’est rien
par rapport `
a tout ce que le R sait faire dans ce contexte. Attention, il n’est
nullement question ici de vous inciter ` a rentrer tout de suite dans ce type
de d´ eveloppement, mais plutˆ ot de vous faire comprendre les possibilit´ es
offertes par ce langage.
6.2.3 Application ` a la probl´ ematique
Fort des quelques explications th´ eoriques propos´ ees ci-dessus, notre utilisateur d´ ebutant propose alors les codes des fonctions relatives ` a sa probl´ ematique
de r´ egression lin´ eaire simple.
1 resume . reg1 <- function (y ,x) {
2
aEst <- cov (x ,y)/ var (x)
3
bEst <- mean (y)- aEst * mean (x)
4
return ( list ( aEst = aEst , bEst = bEst , cor = cor (x ,y )))
5 }
6
7 affiche . reg1 <- function (y ,x) {
8
aEst <- cov (x ,y)/ var (x)
9
bEst <- mean (y)- aEst * mean (x)
10
plot (x ,y)
11
abline (a= bEst ,b= aEst )
12 }
Remarque
Notons au passage que dans les anciennes versions de R, il ´ etait autoris´ e
d’´ ecrire
return(aEst=aEst, bEst=bEst,cor=cor(x,y))
mais que cet usage est en voie d’ˆ etre abandonn´ e dans les futures versions.
Apr` es avoir charg´ e ces fonctions soit par un copier-coller, soit ` a l’aide de la
commande source(), l’utilisateur fait ses premiers tests sur un exemple sans
r´ eel int´ erˆ et.
> y <- rnorm(10);x <- 1:10
> resume.reg1(y,x)
$aEst
[1] -0.218551