Quantile regression à 2 paramètres

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:

---

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:

---

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

---

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

Merci d'avance,

Citation:

---

Envoyé par habasque

Bonjour,

Mais l'erreur suivante survient :

Code:

---

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

---

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

---

Salut,

Code:

---

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:

---

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

---

Pourtant voici l'exemple fournit par le package :

Code:

---

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)])

---

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:

---

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...

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

Code:

---

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.