search for: pcinc

...nse's quantile prediction as a function of the quantile index? (If I understand correctly). So, if I use the following code the predict command seems to work fine dyear<-dummy(ekc$year)[,-1] dstate<-dummy(ekc$state)[,-1] dekc<-cbind(ekc, dyear, dstate) z.nox<-rq(nox~dyear+dstate+pcinc+I(pcinc^2)+I(pcinc^3), tau=-1, data=dekc) zp.nox <- predict(z.nox,newdata=list(pcinc=ekc$pcinc, dyear=dummy(ekc$year)[,-1], dstate=dummy(ekc$state)[,-1]), type="stepfun") but when I am going to do the plot plot(zp.nox,do.points = FALSE, xlab = expression(tau), ylab = expression(Q...

...there any way to do it ? I have post the part of the code Thanks Dimitris ####################################################### quant<-c(0.25, 0.5, 0.75) qrnox<-c("qrnox1", "qrnox2", "qrnox3") for (i in 1:3){ qrnox[i]<-rq(nox~factor(year)+factor(state)+pcinc+I(pcinc^2)+I(pcinc^3), tau=quant[i], data=exmp) } summary(qrnox1) -- View this message in context: http://www.nabble.com/Problem-to-male-an-Index-in-looping-tp19706138p19706138.html Sent from the R help mailing list archive at Nabble.com.

...is happening and I use both predict and predict.lm but there is no difference. Explicitly the code is: fe.nox <- lm(nox~ state.1 + state.2 + state.3 + state.4 + state.5 + state.6 + state.7 + state.8 + state.9 + time.1 + time.2 + time.3 + time.4 + time.5 + time.6 + time.7 + pcinc + I(pcinc^2) + I(pcinc^3), data=ekc) p.fe.nox<-predict.lm(fe.nox, new, interval = "prediction", level=0.95, incl.non.slopes=FALSE) Any Help would be highly appreciated Thanks Dimitris -- View this message in context: http://www.nabble.com/incl.non.slopes%3DFALSE-does-not-work-at...