Discussion :

[A. D.]

Bonjour,

Déjà, il serait bon de jeter un oeil aux messages d'erreur fournis par R.

Code :

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

1
2
3
4
5
6
> hist(x[x!=0],border="darkred",xlims=hlims,ylim=vlims,nclass=20,col=histcol[1])
Messages d'avis :
1: In title(main = main, sub = sub, xlab = xlab, ylab = ylab, ...) :
  "xlims" n'est pas un paramètre graphique
2: In axis(1, ...) : "xlims" n'est pas un paramètre graphique
3: In axis(2, ...) : "xlims" n'est pas un paramètre graphique

Ici, tout est dit : vous avez ajouté un "s" en trop à la fin du paramètre "xlim".

Pour ce qui est du problème des "breaks", je ne sais pas exactement à quoi correspond ce paramètre ni ce que vous cherchez à faire, mais quand je tape la commande que vous indiquez

Code :

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

seq(from=min(hlims),by=diff(bp1$breaks[1:2]))

, j'obtiens uniquement la valeur zéro.


Bonne continuation


Cordialement,

A.D.

[pitipoisson]

Bonjour,

Pour ce qui est du problème des "breaks", je ne sais pas exactement à quoi correspond ce paramètre ni ce que vous cherchez à faire, mais quand je tape la commande que vous indiquez

Code :

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

seq(from=min(hlims),by=diff(bp1$breaks[1:2]))

, j'obtiens uniquement la valeur zéro.

C'est tiré du dernier exemple que j'avais donné, mais il manque effectivement le

Code :

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

..., to=max(hlims),...

sans lequel la fonction seq ne sait pas jusqu'à quelle valeur doit aller la séquence (et fait donc par défaut une séquence de longueur 1) !

[A. D.]

Ah oui, je n'avais pas rejeté un oeil à votre code.

Bref, en fait donc il n'y a pas de problème, mis à part les petites étourderies qu'il a commise, le code fonctionne


Cordialement,

A.D.