418
Le logiciel R - Effectuer des analyses statistiques
> densite <- function(x,n=1000,eps=0.01) {
+ lx <- length(x)
+ res <- vector("integer",lx)
+ obs <- replicate(n,X()) # On génère des observations venant de
# X.
+ for (i in 1:lx) { res[i] <+
length(which((x[i]-eps) <= obs & obs <= (x[i]+eps)))/n
+ }
+ return(res)
+ }
Si l’on trace la courbe qui, ` a chaque portion du support de X, associe
la valeur y de la concentration relative d’observations dans cette portion
(c’est-` a-dire la fr´ equence de valeurs appartenant ` a cette portion), on obtient la
courbe ci-dessous.
> curve(densite,xlim=c(-1,10))
0
2
4
6
8
1 0
0.000
0.004
0.008
x
densite(x)
Si l’on augmente ` a la fois le nombre n d’observations g´ en´ er´ ees par X tandis
que l’on diminue , cette courbe deviendra de plus en plus lisse.
> curve(densite(x,n=1000000,eps=0.001),xlim=c(-1,10))