Université Abdelmalek Essaâdi Année universitaire 2025-2026
(portail Licence SSD ). Faculté des Sciences et Techniques - Tanger
Simulations des variables aléatoires
TP 2 sous R
Exercice 1 (loi normale) Comment utiliser un package ?
1. Taper norm(4) et rnorm(4).
2. Simuler une variable alátoire de loi normale N (0, 1), sa densité, sa fonction de répartition,
et sa fonction quantile (utiliser : dnorm(x, mean = 0, sd = 1, log = FALSE) ; pnorm(q,
mean = 0, sd = 1, [Link] = TRUE, log.p = FALSE) ; qnorm(p, mean = 0, sd = 1,
[Link] = TRUE, log.p = FALSE) ; rnorm(n, mean = 0, sd = 1) ;)
3. Fixer Taper rnorm(100) puis rnorm(100). Taper [Link](20) ; rnorm(100) puis rnorm(100)
puis [Link](20) ; rnorm(100). Que peut-on déduire ?
4. Générer 20 nombres aléatoires suivant une loi normale de paramètres 3 2 (fixer la graine
([Link]) puis utlisrer rnorm(nb, mean,sd)).
5. Représenter la densité gaussiennes sur l’intervalle [−3, 3].
6. Taper ?mvrnorm. Que peut-on déduire ?
7. Taper library(MASS) ; ?mvrnorm.
8. Générer quatre vecteurs
selon
une loi multinormale de moyenne (0,1) et de matrice de
1 0
variance-covariance Σ = .
0 2
9. Calculer les statistiques de base des trois varaiable du jeu de données ethanol(diponible sous
dans le package lattice).
10. Taper ?apply puis calculer les quartiles de chacune des trois variables (utiliser apply avec
quantile).
11. Calculer tous les déciles de chacune des trois variables (préciser les proba dans le précédent
code).
Exercice 2 Exportation de graphiques
1. Tracer, en rouge, la fonction sinus sur [−3π, 3π] en donnant au graphe obtenu le titre
suivant : ”Graphe de la fonction sinus” (voir title ou l’option main de plot).
2. Personnalisation d’un graphique
(a) Tracer la densité de la loi normal N (0, 1) sur [−5, 5] .
(b) Ajouter l’axe des abscisses (voir abline).
(c) Colorier en bleu l’aire sous la courbe à droite correspondant à la probabilité de 5%
voir(polygon). On pourra se servir des codes suivants
X= x >=qnorm(0.95)
Z=c(x[X],rev(x[X]))
W=c(rep(0,sum(X)),rev(dnorm(x[X])))
polygon(Z,W,col=”blue”)
(d) Ajouter une flèche désignant l’aire coloriée (arrows).
(e) Indiquer au bout de la flèche α = 5% (voir text et expression).
(f ) Exporter, sur votre bureau, le graphique obtenu en format jpeg (voir jpeg).
3. Comparaison de distributions
(a) Tracer la courbe de la loi normale entre -4 et 4.
(b) Tracer sur le même graphe les lois de Student à 5 et 30 degrés de liberté. Utiliser la
fonction curve et une couleur pour chaque courbe.
(c) Ajouter une légende en haut à gauche (ou à droite) pour différencier chaque distrubition
(voir legend).
(d) Exporter les courbes, sur votre bureau, en format pdf (voir pdf ).
Exercice 3 (Importation des données sous R)
1. Créer un [Link] contenant les données suivantes
individu ;sexe ;age ;taille
Karim ;M ;18 ;1,60
Fatima ;F ;16 ;1,62
Omar ;M ;20 ;1,70
Yassine ;M ;19 ;1,60
Sanae ;F ;21 ;1,80
Ayoub ;M ;17 ;1,60
2. Taper ?[Link] puis importer les données précédentes sous R.
3. Créer le tableau K concaténant le tableau précédent et la ligne ”Karim M 19 NA” (utiliser
rbind).
4. Créer le tableau Z concaténant le tableau précédent et la colonne ”moyenne 12 13 11 16 13
14 NA” (utiliser cbind).
5. Remplacer les valeurs manquantes par les moyennes correspondantes et visualiser le tableau
obtenu.
Exercice 4 (Intervalle de confiance)
1. Importer le tableau de données hassan sous R.
2. Calculer la moyenne de la variable B.
3. Construire une intervalle de confiance de cette moyenne (voir [Link])
4. Construire à la main cet intervalle en utilisant le fait que
√ √
P x − t499 (0.975) × σ
b/ n, x + t499 (0.975) × σ
b/ n ≈ 0.95
Exercice 5 (Graphique en 3D)
1. Tracer la courbe de la fonction f définie par
p
cos x2 + y 2
f (x, y) = π p
x2 + y 2
on pourra se servir des codes suivants
x = seq(-10, 10, length= 30)
y = x f = f unction(x, y) {r = sqrt(x2 + y 2 ); π ∗ sin(r)/r}
z = outer(x, y, f )
z[[Link](z)] = 1
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = ”lightblue”)
ou
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = ”lightblue”,ltheta = 120, shade =
0.75, ticktype = ”detailed”,xlab = ”X”, ylab = ”Y”, zlab = ”Sinc( r )”)