This works for me:
data(engel)
plot(engel$x, engel$y, xlab="household income", ylab="food expenditure",
cex=.5)
taus <- c(.05,.1,.25,.75,.9,.95)
lapply(taus, function(tau) abline(coef(rq(y~x, data=engel, tau=tau))))
HTH,
Andy
> -----Original Message-----
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of
> [EMAIL PROTECTED]
> Sent: Monday, July 19, 2004 9:53 AM
> To: [EMAIL PROTECTED]
> Subject: [R] why won't rq draw lines?
>
>
> I've been trying to draw quantile linear regression lines
> across a scatterplot of my data using
>
> attach(forrq)
> plot(PREGNANT,DAY8,xlab="pregnant EPDS",ylab="postnatal
> EPDS",cex=.5)
> taus <- c(.05,.1,.25,.75,.9,.95)
> xx <- seq(min(PREGNANT),max(PREGNANT),100)
> for(tau in taus){
> f <- coef(rq(DAY8~PREGNANT,tau=tau))
> yy <- (f[1]+f[2]*xx)
> lines(xx,yy)
> }
> which is simply the method from the help file with my dataset
> attached, and the variable names substituted where
> appropriate. I get the scatterplot, but no lines. Any ideas
> about what's going on? or wrong?
>
>
> ______________________________________________________________
> _________
> Most doctors use http://www.Doctors.net.uk e-mail.
> Move to a free professional address with spam and virus protection.
>
> ______________________________________________
> [EMAIL PROTECTED] mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
>
______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
