Interfaces utilisateur et rapports statistiques avec R et Sas
Introduction
Chapitre 1. Environnement de travail
Chapitre 2. Le langage R
2.1. Vecteurs
2.1.1. Création
2.1.2. Indexation
2.1.3. Valeurs manquantes
2.1.4. Fonctions prédéfinies
2.2. Matrices
2.3. Listes
2.4. Tableaux de données
2.5. Fonctions personnalisées
2.6. Exercices
Chapitre 3. Rapports statistiques
Chapitre 4. Interfaces utilisateur
Chapitre 5. Bases de données et langage Sql
Annexe A : une applette Java
Annexe B : le logiciel Sas
Page d'accueil Table des matières Niveau supérieur Page précédente Bas de la page Page suivante Démarrage du logiciel Webographie du module Glossaire du module

2.1.4. Fonctions prédéfinies

Il existe des milliers de fonctions prédéfinies. Pour en avoir un premier aperçu, je vous suggère d'examiner la liste de celles qui appartiennent à 3 packages fondamentaux (base, stats et graphics).

library(help="base")  # ou help(base) 
library(help="graphics") 
help(mean)  # ou toute autre fonction...

Dans cette jungle, j'ai sélectionné quelques fonctions pour les mettre en valeur :

# manipulation de vecteurs : length, is.na...
notes = c(15,8,5,NA,10)
length(notes)  # taille du vecteur : 5 ici
is.na(notes)  # valeurs manquantes : c(F,F,F,T,F)
which(!is.na(notes))	# indices des valeurs non manquantes : c(1,2,3,5)
# fonctions de recodage : ifelse, cut...
avis = ifelse(notes<8,"non", ifelse(notes<10,"oral","oui") )
print(avis)
avis = cut(notes, breaks=c(0,8,10,20), labels=c("non","oral","oui"), 
		right=F, include.highest=T)
# fonctions statistiques : mean, sd (ecart-type)...
mean(notes, na.rm=T)  # ne pas oublier d'oter les valeurs manquantes
quantile(notes, na.rm=T)
table(avis)  # les effectifs par avis
# graphiques : boxplot, hist, barplot...
boxplot(notes, main="boite à moustaches")
hist(notes, main = "histogramme")
pie(table(avis), main="camembert")
barplot(table(avis), main="diagramme en barres" )

Il est souvent utile d'examiner l'aide relative aux fonctions car celles-ci peuvent avoir de nombreuses options dont on peut changer les valeurs par défaut.

# quelques options de la fonction quantile : 
help(quantile)	# pour connaitre les options
quantile(notes) # NA car notes a une valeur manquante
quantiles(notes, na.rm=T) # min, max et quartiles
quantiles(notes, na.rm=T, probs=0.25)	# premier quartile
quantiles(notes, na.rm=T, probs=c(0.25,0.75))	# quartiles 1 et 3
# quelques options de la fonction hist : 
hist(notes)
hist(notes, main = "histogramme des notes")	# on change le titre
hist(notes, breaks=c(0,8,10,20)) # avec intervalles fixés  
hist(notes, breaks=c(0,8,10,20), right=FALSE) 
	# les intervalles sont [0,8[ (8  exclu), [8, 10[ (8 inclus)...
hist(notes, breaks=c(0,8,10,20), right=FALSE, 
	labels=c("recale","oral","admis")) # avec legendes pour les barres
hist(notes, breaks=c(0,8,10,20), right=FALSE, 
	labels=c("recale","oral","admis"), col=c("red","orange","green"))
hist(notes, breaks=c(0,8,10,20), right=FALSE, plot=FALSE) 
	# pas de graphique, juste des statistiques	

Pour plus d'informations sur les graphiques, je vous suggère de faire un tour sur le site d'enseignement [PBIL] et d'y chercher "graphiques" avec la zone de texte prévue à cet effet.

Page d'accueil Table des matières Niveau supérieur Page précédente Haut de la page Page suivante Démarrage du logiciel Webographie du module Glossaire du module