Discussion :

[Atacama]

Bonjour à tous,

Je sollicite votre aide car je suis débutant sur R et je ne parviens pas à réaliser ce que je souhaiterais.

Mon objectif est d'avoir une ligne de programme qui permet de calculer la fonction de réparition empirique d'un ensemble de données.
Sur R, j'ai pu voir qu'il existe une fonction "ecdf()" qui permet de l'obtenir graphiquement. J'ai également vu que la commande "knots" renvoie aux valeurs en x de la fonction.

De mon côté, je voudrais obtenir les valeurs en Fn(x) qui sont comprises entre 0 et 1.

Pour ça, j'ai tenté d'écrire une ligne de programme qui malheureusement ne fonctionne pas. Elle prend la forme suivante:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
### Valeurs initiales ###
 
a=c(14,15,12,9,7,3,2,18,6,17)
 
### Mise en ordre croissant des valeurs du vecteur a ###
 
b=sort(a)
 
### Calcul de la fonction de répartition empirique ###
 
FRE=function(x) function(t){
t=1:10
sum(x %in% t)/length(x)}
epdf(b)

Je souhaiterais avoir quelques indications pour faire fonctionner ce code.

Je vous remercie.


Atacama

[hoccha]

Bonjour,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
> ttt <- ecdf(c(14,15,12,9,7,3,2,18,6,17))
> knots(ttt)
 [1]  2  3  6  7  9 12 14 15 17 18
> ttt(knots(ttt))
 [1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

Dans ce code, la ligne ttt(knots(ttt)) donne "les valeurs en Fn(x) qui sont comprises entre 0 et 1".
Pour conprendre, il faut savoir que ecdf() est une fonction qui retourne une fonction.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
> class(ttt)
[1] "ecdf"     "stepfun"  "function"

C'est bien ce que vous vouliez faire?

A++

[Atacama]

Merci pour votre réponse rapide. Ca me donne bien les valeurs que je voulais.