Intégration d'un estimateur histogramme

Bonjour à tous,

J'ai écrit la fonction suivante, qui est un histogramme régulier à D intervalles sur [0,1]. Disons X<-runif(100) et D<-10. La fonction est bien écrite et ne me pose pas de problème. Par contre, je ne parviens pas à l'intégrer. Quand j'écris integrate(f,-5,5);, il me répond "evaluation of function gave a result of wrong length".

Quelqu'un pourrait m'aider à intégrer ce type de fonction?
Merci d'avance !

Citation:

---

f<-function(x){
#on compte d'abord le nombre de points qui tombent dans l'intervalle I_j pour j=1..D

#but:créer un vecteur qui donne le nbre de points X_i dans chaque intervalle de longueur 1/D
n<-length(X)
Y<-sort(X)
N<-numeric(D+1)
N[1]<-0
s<-0
k<-1
for(j in 1:D-1){
while((j-1)/D <= Y[k] && Y[k]<j/D){k<-k+1}
N[j+1]<-k-1-s
s<-s+N[j+1]
}
N[D+1]=length(Y)-s
N<-N[-1]

if(x<0 || x>1) return(0)

for(j in 1:D){
if((j-1)/D<=x && x <j/D) return(N[j]*D/n)
}

}

---

Bonjour,

Dans mon cas ça fonctionne si je mets au début de la fonction (car les calculs de la suite sont inutiles dans ce cas) :

Code:

---

if(x<=0 || x>=1) return(x*0)

---

Cette modification résout l'erreur:
"evaluation of function gave a result of wrong length"


Votre fonction a des erreurs du coup elles se répercutent sur la fonction integrate:
"Erreur dans while ((j - 1)/D <= Y[k] && Y[k] < j/D) { :
valeur manquante là où TRUE / FALSE est requis"


Une impression de déjà vue :
http://www.developpez.net/forums/d12...tion-comptage/

Avez-vous regardé le dernier post à propos des variables utilisées?

ma fonction marche sans erreur chez moi;
et si je mets ce que tu me dis au début (la condition if) alors il intègre bien, mais la réponse est mauvaise, il donne 0 comme résultat (alors que la réponse est 1).

Tu as une idée?

Le problème est que votre fonction n'est compatible qu'avec des scalaires et non des vecteurs.
Si vous faîtes:

Code:

---

1 2 3

x=seq(0,1,0.01) f(x)

---

Vous obtenez une liste nul.

Il faut donc la modifier pour des listes 1D.

C'est vrai! ma fonction n'est compatible qu'avec des scalaires.
Mais pour quelle(s) raison(s)? et en quoi cela empêche de pouvoir l'intégrer?
Comment la modifier pour des listes 1D?


Merci encore pour votre aide!

Ce n'est pas compliqué, il suffit de rajouter une boucle for(i in 1:length(x)) au début et remplacer x par x[i].