Discussion :

[habasque]

Bonjour,

Je souhaiterais réaliser une régression quantile en utilisant la fonction "rqss" du package quantreg. Le but est de modéliser la distribution "haute" de mes données. Grosso modo je souhaiterais représenter la courbe rouge qui est sur la figure jointe.

Le code est le suivant (en se basant sur l'exemple fournit dans la documentation du package):

Code :

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

1
2
3
4
5
 
x11()
plot(dataset$CHLA,dataset$SST,xlab = "CHLA", ylab = "SST", xlim=c(0,3),ylim=c(9,30))
fit <- rqss(SST ~ qss(CHLA), tau = 0.90, data=dataset)
lines(dataset$CHLA[-1],fit$coef[1]+fit$coef[-1])

Mais l'erreur suivante survient :

Code :

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

xy.coords(x, y) : 'x' and 'y' lengths differ

Quelqu'un aurait une idée d'où cela peut provenir ?

Merci d'avance,