Discussion :

[Lhaneslo]

Bonjour,
Je me donne un vecteur X de points sur [0,1] et un nombre D d'intervalles.

Voici mon programme; je note N_(j+1) le nombre de X_i qui tombent dans l'intervalle I_j=[(j-1)/D,j/D]. Quelqu'un peut-il me dire où est l'erreur?
N est le vecteur qui contient le nombre de points que contient chaque intervalle.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
comptage<- function(X,D){
Y<-sort(X)
 
N<-numeric(D+1)
N[1]<-0
s<-0
k<-1
for(j in 1:D){
	while((j-1)/D <= Y[k] & Y[k]<j/D){k<-k+1}
	N[j+1]<-k-1-s
	s<-s+N[j+1]
	}
}

L'erreur suivante s'affiche alors:

Code :

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

1
2
Erreur dans while ((j - 1)/D <= Y[k] & Y[k] < j/D) { : 
  valeur manquante là où TRUE / FALSE est requis

[Gakusei]

Bonjour,

l'indice k dépasse la taille de la liste Y, donc Y[k] renvoie NA donc la comparaison devient aussi NA et while renvoie ton erreur.
pour s"en assurer il suffit d'afficher k dans la boucle.

[Lhaneslo]

Bien vu!
Du coup je vais faire tourner ma boucle de j=1 à D-1. Ensuite je peux écrire N[D+1]=length(Y)-s. Et j'ai mon vecteur N!

Merci beaucoup!

Bonjour,

l'indice k dépasse la taille de la liste Y, donc Y[k] renvoie NA donc la comparaison devient aussi NA et while renvoie ton erreur.
pour s"en assurer il suffit d'afficher k dans la boucle.

[Gakusei]

Je pense qu'il faut revoir ton script:

Code :

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

1
2
3
4
5
6
 
N[j+1]<-k-1-s
s<-s+N[j+1]
<=>
N[j+1]<-k-1-s
s<-k-1

Il faut que tu définisses bien le rôle de chaque variable.

J'ai simplifié ton code
avec interval les bornes supérieures de ton intervalle
et cmpt l'indice des intervalles:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
comptage<- function(X,D){
  Y<-sort(X)
  interval=1:D/D
  k=0
  s=1
  cmpt=1
  N=interval*0
  for(j in interval){
    print(paste("intervalle:",j))
    while(Y[k+s]<j & (k+s)< length(Y)){k<-k+1; }
    print(paste("nmbre de valeurs:",k))
    N[cmpt]=k
    s<-s+k
    k=0
    cmpt=cmpt+1
  }
  return(N)
}