10-rbase.Rmd

# Aller plus loin avec les objets et la programmation fonctionnelle

Ce qui a été présenté dans ce module repose sur les fonctions du package `tidyverse`. Cette approche tend à se généraliser depuis quelques années, mais quand on cherche la réponse à un problème sur Internet, on trouve d'autres façons de programmer en R, qui font appel aux fonctions du package `base` et non du `tidyverse`
$\Rightarrow$ Cette partie donne quelques clés de compréhension.

## Les objets dans R, plus de détails

Rappel : en informatique, un objet est défini par : ses *attributs* et ses *méthodes* (fonctions). Dans l'exemple du jeu d'échec, chaque pièce peut être vue comme un objet :

- sa position sur le plateau constitue ses attributs
- sa façon de se déplacer peut être vue comme une fonction qui ne s'applique qu'à ce type de pièce, donc une méthode

R est un langage orienté objet ; ces objets permettent de structurer les données selon leurs caractéristiques $\Rightarrow$ on retrouve les données dans les attributs. Les méthodes sont en général transparentes pour l'utilisateur (cf. utilisation des fonctions `summary`, `plot`...). Les objets les plus courants sont les suivants :

- **Vecteurs** : suite unidimensionnelle de **valeurs** ayant le même type.
- **Facteurs** : vecteur qui prend un nombre limité de modalités (exemple : sexe). Il est défini par les niveaux (`levels`) et les libellés associés (`labels`).
- **Matrice** et **arrays** : suites multidimensionnelles de **valeurs** (matrices=dimension 2 ; array=dimension n).
A la différence d'une dataframe, les valeurs d'une matrice sont toutes du même type.
Les arrays peuvent être très puissants pour gérer des millésimes.
- **Liste** : ensemble d'**objets** différents. On peut stocker un vecteur alphanumérique + une matrice numérique dans une liste.
- **Tableaux** (`data.frame`) : Objet qui ressemble le plus aux tables Excel, SAS ou SPSS... : description d'individus statistiques (observations, en ligne) par des caractéristiques (variables, en colonnes).
* **Fonctions** : Objets particuliers qui donnent un *résultat* à partir de paramètres en entrée.
- **Autres objets** : Il existe un très grand nombre d'objets *ad hoc* dans R. Par exemple 
    + *ts* (time serie) pour les séries temporelles, 
    + *lm* (linear model) qui contient tous les résultats d'une régression linéraire... 
    + des graphiques
    + On peut même en définir de nouveaux soi-même !


## Créer une nouvelle fonction en R

La fonction est un objet comme les autres, qu'on crée avec l'opérateur d'affectation. Elle est définie par des paramètres et elle se termine par la fonction `return()`. On reprend l'exemple du calcul de l'IMC

```{r}

calcul_IMC <- function (poids, taille)
{
  ## La taille est exprimée en mètres
  imc <- poids / taille ^ 2
  return (imc)
}
calcul_IMC (poids = 80, taille = 1.89)
calcul_IMC (poids = 60, taille = 1.55)
```

## Les boucles conditionnelles

Les commandes `if` et `else` sont bien entendues utilisables. Le "then" n'existe pas : il est implicite après les accolades.

```{r}
diag_IMC <- function(poids,taille)
{
  imc <- poids / taille ^ 2
  if (imc < 18.5) {diag <- "maigre"}
  else if (imc < 25) {diag <- "normal"}
       else {diag <- "surpoids"}
  return (diag)
}
diag_IMC (poids=60,taille=1.89)
diag_IMC (poids=80,taille=1.89)
diag_IMC (poids=80,taille=1.55)
```


## Les boucles

On peut utiliser les boucles classiques : `repeat`, `while`, `for` :

```{r}
for (pp in seq(from = 50, to = 100, by = 5))
{
  print(paste ("Taille = 1,70m, poids =", pp, "Diagnotic :",
               diag_IMC (poids = pp, taille = 1.70)))
}
```


## Pour aller plus loin 

### matrices et arrays

Les matrices et les arrays permettent des calculs rapides et efficaces, et peuvent être très pratiques et optimisent le stockage des données. Ils demandent cependant plus de réflexion en amont quant à leur utilisation. On accède aux éléments avec les [].

Un hypercube de trois dimensions peut être représenté comme suit :
![](images/grille.jpg){width=50%}
On peut par exemple créer une matrice à 10 lignes et 10 colonnes remplie avec un tirage aléatoire selon une loi normale centrée réduit.
De la même façon on peut créer un hypercube avec la fonction avec 10 lignes, 5 colonnes et de profondeur 3, toujours avec un tirage aléatoire selon une loi normale

```{r}
mat <- matrix(rnorm(50), ncol = 5, nrow = 10)
arr <- array(rnorm(150),dim = c(10,5,3))
mat
arr
```

*Pourquoi s'embêter avec ça ?* Parce qu'on peut appliquer des fonctions facilement sur les lignes, colonnes et autres dimensions grâce à la fonction  ```apply()```.
Exemple : résultats de validations croisées par bloc, simulations de loi selon différents paramètres. Et on calcule facilement des statistiques "marginales".

Par, exemple, sur une matrice, on peut calculer des statistiques par lignes :

```{r}
apply(mat, MARGIN = 1, FUN=mean)
```

Ou par colonnes :

```{r}
apply(mat, MARGIN = 2, FUN=mean)
```

Sur notre hypercube de type ```array```, on peut aussi calculer des stats sur ses différentes dimensions :

```{r}
apply (arr, MARGIN = 3, FUN=mean)
apply (arr, MARGIN = c(2,3), FUN = mean)
```

### Inspection d'un objet : la régression

```{r,echo=FALSE}
library ("ggplot2")
reg <- lm (data = iris, formula = Sepal.Length ~ Petal.Length)
ggplot (data = iris, aes(x = Petal.Length, y = Sepal.Length)) +
  geom_point ()  +
  geom_smooth (method = "lm", se = FALSE) +
  ggtitle ("Illustration de la régression linéaire simple")
```


La régression linéaire consiste à exprimer une variable Y en fonction d'une variable X dans une fonction linéaire. C'est à dire qu'on cherche a et b tels que : $$ Y = a \cdot X + b + \epsilon$$ 
où $\epsilon$ est le résidu de la régression. On utilise dans cet exemple la table des iris de Fisher, existant dans R base qu'il suffit d'appeler avec `data(iris)` (il existe d'autres dataframe inclus dans les packages et qui sont utilisés en exemple dans l'aide).

```{r}
data ("iris")
str (iris)
```

Faire la régression de la Sepal.Length sur Petal.length à l'aide de la fonction `lm()`

```{r}
lm (data = iris, formula = Sepal.Length ~ Petal.Length)
```
On a les paramètres a et b mais on aimerait en savoir plus... Au moins la qualité d'ajustement (le $R^2$ par exemple), et un graphique des résidus pour détecter une éventuelle structure. Pour cela, stocker le résultat dans un nouvel objet, et explorez-le avec les fonctions `str()`, `summary()` et `plot()`


```{r}
reg <- lm(data = iris, formula = Sepal.Length ~ Petal.Length)
str (reg)
summary (reg)
plot (reg)
```

Les **méthodes** summary, print et plot sont implémentées pour tous les objets en R, et on peut les utiliser pour avoir un premier aperçu de ce que l'on obtient avec la fonction.