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,

[Sengar]

Bonjour,

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 ?

Salut,

Code :

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

lines(dataset$CHLA[-1],fit$coef[1]+fit$coef[-1])

Que veux tu plotter avec ca?
Ca me parait bizarre de plotter la variable d'intérêt en fonction des coefficients (les deux n'ayant à priori pas la même taille, d'où l'erreur renvoyée...).

Après je connais mal le package donc peut être fit$coef correspond à autre chose?



Edit: ok j'ai compris tu voulais juste plotter la regression quantile fittée.

Du coup ce serait plutôt quelque chose dans ce genre la:

Code :

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

lines(dataset$CHLA,fit$coef[1]+fit$coef[-1] * dataset$CHLA)

[habasque]

Pourtant voici l'exemple fournit par le package :

Code :

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

1
2
3
4
5
6
7
8
 
n <- 200
x <- sort(rchisq(n,4))
z <- x + rnorm(n)
y <- log(x)+ .1*(log(x))^2 + log(x)*rnorm(n)/4 + z
plot(x, y-z)
f.N  <- rqss(y ~ qss(x, constraint= "N") + z)
lines(x[-1], f.N $coef[1] + f.N $coef[-(1:2)])

[habasque]

Merci, effectivement je veux plotter la courbe fittée.
Par contre, le code proposé ne fonctionne pas.
J'ai le message d'avis suivant :

Code :

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

1
2
3
 
In fit$coef[-1] * dataset_cjm_model_sans_na_presence$CHLA :
  la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court

Et j'obtiens le graph joint...

[habasque]

En fait il manquait unique sur le 1er terme de la commande lines

Code :

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

1
2
3
4
5
6
x11()
plot(dataset$CHLA,dataset$SST, xlab = "CHLA", ylab = "SST", xlim=c(0,5),ylim=c(9,30), main="Quantile regression tx = 95%")
fit <- rqss(SST ~ qss(CHLA, constraint="N"), tau = 0.95, data=dataset)
lines(unique(dataset$CHLA[-1]), fit$coef[1] + fit$coef[-1], col="blue")
fit2 <- rqss(SST ~ qss(CHLA, constraint="N"), tau = 0.01, data=dataset)
lines(unique(dataset$CHLA[-1]), fit2$coef[1] + fit2$coef[-1]), col="black")

Le résultat est sur la figure jointe.