similar to: quantiles and dataframe

Hello First, thanks for the help for an earlier question about error handling! I have problem filtering a dataset. I'm trying to filter the data in the y columns based on the values in the x column, e.g.: x y1 y2 yn 1.0 1 NA 3 2.0 1 NA 11 2.0 2 NA NA 3.0

Hi All, I am looking at applications of percentiles to time sequenced data. I had just been using the quantile function to get percentiles over various periods, but am more interested in if there is an accepted (and/or R-implemented) method to apply weighting to the data so as to weigh recent data more heavily. I wrote the following function, but it seems quite inefficient, and not really very

I have a set of values and their corresponding weights. I can use the function weighted.mean to calculate the weighted mean, I would like to be able to similarly calculate the weighted median and quantiles? Is there a function in R that can do this? thanks, Spencer [[alternative HTML version deleted]]

Dear list, I have to look examine hourly time - series and would like to plot variable section of them using plot.zoo. Hourly time series data which looks like this: YYYY MM DD HH P-uk P-kor P-SME EPOT EREA RO R1 R2 RGES S-SNO SI SSM SUZ SLZ 2003 1 1 1 0.385 0.456 0.021 0.000 0.000 0.000 0.013 0.223 0.235 0.01 0.38

Hi All, I am running quantile regression in a "for loop" starting with 1 variable and adding a variable at a time reaching a maximum of 20 variables. I get the following warning messages after my "for" loop runs. Should I be concerned about these messages? I am building predictive models and am not interested in inference. Warning messages: 1: In

The help page and vignette for summary.rq(quantreg) mention that there are three different bootstrap methods available for the se="bootstrap" argument, but I can't figure out how to select a particular method. For example, if I want to use the "xy-pair bootstrap" how do I indicate this in summary.rq? Tom -- View this message in context:

Hej, i need the which() funktion to find the positions of an entry in a matrix. the entries i'm looking for are : seq(begin,end,0.01) and there are no empty spaces i'm searching in the right range. so i was looking for the results R can find and i recieved this answer. for (l in

Hi, how is it possible to retrieve the corresponding tau value for each observed data pair (x(t) y(t), t=1,...,n) when doing a quantile regression like rq.fit <- rq(y~x,tau=-1). Thank you for your help. Jaci --

I am new to R and am interested in using the program to fit quantile regression models to data collected from a multi-stage probability sample of the US population. The quantile regression package, rq, can accommodate person weights. However, it is not clear to me that boot.rq is appropriate for use with multi-stage samples (i.e., is capable of sampling primary sampling units instead of survey

I have spent hours on this ---looked through the quantreg manual and r-help site--- still couldn't figure out the answer. Can someone please help me on this? I plot the result from quantile regression and want to change the variable labels: temp<-rq(dep~inc+age50, data=newdata, tau=1:9/10) temp2<-plot(summary(temp)) dimnames(temp2)[[1]]<-c("Intercept", "Per Capita

Motivated by Deepayan's recent inquiries about the efficiency of the R 'quantile' function: http://tolstoy.newcastle.edu.au/R/devel/05/11/3305.html http://tolstoy.newcastle.edu.au/R/devel/06/03/4358.html I decided to try to revive an old project to implement a version of the Floyd and Rivest (1975) algorithm for finding quantiles with O(n) comparisons. I used

Hello Everyone, I'm currently learning about quantile regressions. I've been using an optimizer to compare with the rq() command for quantile regression. When I run the code, the results show that my coefficients are consistent with rq(), but the intercept term can vary by a lot. I don't think my optimizer code is wrong and suspects it has something to do with the starting

Hi, I am trying to estimate a quantile regression using panel data. I am trying to use the model that is described in Dr. Koenker's article. So I use the code the that is posted in the following link: http://www.econ.uiuc.edu/~roger/research/panel/rq.fit.panel.R While this code run perfectly, it does not work for my data providing a warning message: In rq.fit.sfn(D, y, rhs = a) : tiny

Dear all, I am trying to compare the estimated coefficients of a quantile regression model between two different samples. It is a Wald test, but I cannot find one way to do that in R.The samples are collected conditional on a specific characteristic and I would like to test whether such characteristic indeed affect the estimators. The problem in the test anova.rq is that the response variable

Dear Users, I loaded my dataset as following: presu <- read.table("C:/_Ricardo/Paty/qtdata_f.txt", header=TRUE, sep="\t", na.strings="NA", dec=".", strip.white=TRUE) dep<-presu[,3]; exo<-presu[,4:92]; When I try: rq(dep ~ exo, ...) or mle.stepwise(dep ~ exo, ...) I got the same error: > rq(dep ~ exo) Error in model.frame(formula, rownames,

Hi, how to measure the goodness of fit, when using the rq() function of quantreg? I need something like an R^2 for quantile regression, a single number which tells me if the fit of the whole quantile process (not only for a single quantile) is o.k. or not. Is it possible to compare the (conditional) quantile process with the (unconditional) empirical distribution function? Perhaps with a Chi^2

Dear all, I'm trying to get the variance/covarince matrix after fitting a quantile regression model (either linear or non linear), in order to get the variance of my predictions and be able to calculate the median squared error. The commands working for the lm models (corr=T or vcov=T) do not seem to work for the rq models. Could you advise me a way of getting it? Best regards

Hello, I?m wondering if anyone can offer advice on the out-of-memory error I?m getting. I?m using R2.12.2 on Windows XP, Platform: i386-pc-mingw32/i386 (32-bit). I am using the quantreg package, trying to perform a quantile regression on a dataframe that has 11,254 rows and 5 columns. > object.size(subsetAudit.dat) 450832 bytes > str(subsetAudit.dat) 'data.frame': 11253 obs.

Dear all. I have a question on plotting the coefficients from a series of mutivariate quantile regressions. The following code plots the coefficients for each RHS variable x1 and x2. Is there a way to plot only the coefficients on x1? In the data I am using, I have a large number of fixed effects and do want to plot the coefficients on these fixed effects. quant.plot <-

I am a new to using quantile regressions in R. I have estimated a set of coefficients using the method="br" algorithm with the rq command at various quantiles along the entire distribution. My data set contains approximately 2,500 observations and I have 7 predictor variables. I receive the following warning message: Solution may be nonunique in: rq.fit.br(x, y, tau = tau, ...)