Programmer fonction répartition empirique
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:
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
