similar to: wald test: compare quantile regression estimators from different samples

Displaying 20 results from an estimated 6000 matches similar to: "wald test: compare quantile regression estimators from different samples"

2012 Jul 28
4
quantreg Wald-Test
Dear all, I know that my question is somewhat special but I tried several times to solve the problems on my own but I am unfortunately not able to compute the following test statistic using the quantreg package. Well, here we go, I appreciate every little comment or help as I really do not know how to tell R what I want it to do^^ My situation is as follows: I have a data set containing a
2006 Dec 02
2
nonlinear quantile regression
Hello, I?m with a problem in using nonlinear quantile regression, the function nlrq. I want to do a quantile regression o nonlinear function in the form a*log(x)-b, the coefficients ?a? and ?b? is my objective. I try to use the command: funx <- function(x,a,b){ res <- a*log(x)-b res } Dat.nlrq <- nlrq(y ~ funx(x, a, b), data=Dat, tau=0.25, trace=TRUE) But a can?t solve de problem,
2010 Oct 07
3
quantile regression
Dear all, I am a new user in r and I am facing some problems with the quantile regression specification. I have two matrix (mresultb and mresultx) with nrow=1000 and ncol=nsim, where I specify (let's say) nsim=10. Hence, the columns in my matrix represents each simulation of a determined variable. I need to regress each column of mresultb on mresultx. My codes are the following:
2011 Oct 16
1
nlrq {quantreg}
Dear all, I sent an email on Friday asking about nlrq {quantreg}, but I haven't received any answer. I need to estimate the quantile regression estimators of a model as: y = exp(b0+x'b1+u). The model is nonlinear in parameters, although I can linearise it by using log.When I write: fitnl <- nlrq(y ~ exp(x), tau=0.5) I have the following error: Error in match.call(func, call = cll) :
2005 Apr 15
2
negetative AIC values: How to compare models with negative AIC's
Dear, When fitting the following model knots <- 5 lrm.NDWI <- lrm(m.arson ~ rcs(NDWI,knots) I obtain the following result: Logistic Regression Model lrm(formula = m.arson ~ rcs(NDWI, knots)) Frequencies of Responses 0 1 666 35 Obs Max Deriv Model L.R. d.f. P C Dxy Gamma Tau-a R2 Brier 701 5e-07 34.49
2012 May 28
2
R quantreg anova: How to change summary se-type
He folks=) I want to check whether a coefficient has an impact on a quantile regression (by applying the sup-wald test for a given quantile range [0.05,0.95]. Therefore I am doing the following calculations: a=0; for (i in 5:95/100){ fitrestricted=rq(Y~X1+X2,tau=i) tifunrestrited=rq(Y~X1+X2+X3,tau=i) a[i]=anova(fitrestricted,fitunrestricted)$table$Tn) #gives the Test-Value } supW=max(a) As anova
2010 Oct 13
1
(no subject)
Dear all, I have just sent an email with my problem, but I think no one can see the red part, beacuse it is black. So, i am writing again the codes: rm(list=ls()) #remove almost everything in the memory set.seed(180185) nsim <- 10 mresultx <- matrix(-99, nrow=1000, ncol=nsim) mresultb <- matrix(-99, nrow=1000, ncol=nsim) N <- 200 I <- 5 taus <- c(0.480:0.520) h <-
2009 Aug 17
3
Help understanding lrm function of Design library
Hi, I'm developing an experiment with logistic regression. I've come across the lrm function in the Design library. While I understand and can use the basic functionality, there are a ton of options that go beyond my knowledge. I've carefully read the help page for lrm, but don't understand many of the arguments and optional return values. (penalty, penalty.matrix,
2009 Jun 09
1
Non-linear regression/Quantile regression
Hi, I'm relatively new to R and need to do a quantile regression. Linear quantile regression works, but for my data I need some quadratic function. So I guess, I have to use a nonlinear quantile regression. I tried the example on the help page for nlrq with my data and it worked. But the example there was with a SSlogis model. Trying to write dat.nlrq <- nlrq(BM ~ I(Regen100^2),
2006 Oct 25
1
Quantile Regression
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 --
2011 Dec 05
1
extract cov matrix in summary.rq and use as a matrix.
Dear all, I need to extract the covariance matrix of my quantile regression estimation to use in a test. My regression is: qf2_1 <- summary(rq(wb2 ~ apv2 + vol2, tau = phi2[1]), cov = TRUE) I can extract the covaraince matrix by using: qf2_1 [3]. However, if I try to use it in the test, it does not work. I only need to transform qf2_1[3] in a matrix 3x3. I have already tried:
2011 Nov 05
2
linear against nonlinear alternatives - quantile regression
Dear all, I would like to know whether any specification test for linear against nonlinear model hypothesis has been implemented in R using the quantreg package. I could read papers concerning this issue, but they haven't been implemented at R. As far as I know, we only have two specification tests in this line: anova.rq and Khmaladze.test. The first one test equality and significance of
2006 Jul 23
1
Warning Messages using rq -quantile regressions
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, ...)
2004 Sep 30
1
polr (MASS) and lrm (Design) differences in tests of statistical signifcance
Greetings: I'm running R-1.9.1 on Fedora Core 2 Linux. I tested a proportional odds logistic regression with MASS's polr and Design's lrm. Parameter estimates between the 2 are consistent, but the standard errors are quite different, and the conclusions from the t and Wald tests are dramatically different. I cranked the "abstol" argument up quite a bit in the polr
2008 Jan 01
2
Non-Linear Quantile Regression
Please, I have a problem with nonlinear quantile regression. My data shows a large variability and the quantile regression seemed perfect to relate two given variables. I got to run the linear quantile regression analysis and to build the graph in the R (with quantreg package). However, the up part of my data dispersion seems a positive exponential curve, while the down part seems a negative
2006 Oct 02
1
a question regarding 'lrm'
Hi List, I don't understand why 'lrm' doesn't recognize the '~.' formula. I'm pretty sure it was working before. Please see below: I'm using R2.3.0, WinXP, Design 2.0-12 thanks, ...Tao > dat <- data.frame(y=factor(rep(1:2,each=50)), x1=rnorm(100), x2=rnorm(100), x3=rnorm(100)) > lrm(y~., data=dat, x=T, y=T) Error in terms.formula(formula, specials =
2011 Aug 05
1
Goodness of fit of binary logistic model
Dear All, I have just estimated this model: ----------------------------------------------------------- Logistic Regression Model lrm(formula = Y ~ X16, x = T, y = T) Model Likelihood Discrimination Rank Discrim. Ratio Test Indexes Indexes Obs 82 LR chi2 5.58 R2 0.088 C 0.607 0
2011 Aug 23
3
Change Variable Labels in Quantile Plot
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
2009 Aug 21
1
Possible bug with lrm.fit in Design Library
Hi, I've come across a strange error when using the lrm.fit function and the subsequent predict function. The model is created very quickly and can be verified by printing it on the console. Everything looks good. (In fact, the performance measures are rather nice.) Then, I want to use the model to predict some values. I get the following error: "fit was not created by a Design
2004 Feb 04
1
Fitting nonlinear (quantile) models to linear data.
Hello. I am trying to fit an asymptotic relationship (nonlinear) to some ecological data, and am having problems. I am interested in the upper bound on the data (i.e. if there is an upper limit to 'y' across a range of 'x'). As such, I am using the nonlinear quantile regression package (nlrq) to fit a michaelis mention type model. The errors I get (which are dependant on