similar to: variance covariance matrix of parameter estimate using nlrq

I have been quantreg library for a number of projects but have just hit a snag. I am using nlrq to examine an asymptotic relationship between 2 variables at the 99th percentile. It performs as expected, however when I try to extract the coefficients along with se and significance I am running into problems. The problem is that for the nlrq regression Dat.nlrq, summary(Dat.nlrq) reports a different

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

Dear List, Apologies if this is a known problem ... I wasn't able to find it on the bug list, but it is a problem that does not seem to occur with a MAC build of R 2.0, so perhaps this problem has already been addressed for the future. I am getting *sporadic* errors when refitting the same model to the same data set, using nlrq() in the nlrq package. The algorithm is not stochastic, so I

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,

Hi there, Can anyone please tell me if it is possible to limit parameters in nlrq() to 'upper' and 'lower' bounds as per nls()? If so how?? Many thanks in advance

Dear r-help members I posted this message already yesterday, but don't know whether it reached you since I joined the group only yesterday. I would like to estimate the boxcox transformed model (y^t - 1)/t ~ b0 + b1 * x. Unfortunately, R returns with an error message when I try to perform this with the call nls( I((y^t - 1)/t) ~ I(b0 + b1*x), start = c(t=1,b0=0,b1=0), data = mydataframe)

Hi, I'm attempting to calculate the 0.25 and 0.97 quantiles for tree height (0-50 meters) against tree age (0-300 years) and I am running into some difficulty with the plotted grafic. I've run the examples in the quantreg help and can get those to work properly and by plugging in my data I can also get the lines plotted on my dataset. Unfortunately I'm running into a problem with the

Johannes, You can minimize an model expression by just putting the ~ on the left and everything else on the righthand side, but I don't think that this is really what you want. In the NLS expression this would ignore the jacobian of the transformation from errors to response, and in nlrq there is the same problem, however you can adjust for the jacobian by rescaling by the geometric mean of

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

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

All, This worked: mBW <- function( ... ) ... # matrix-valued function BaconWatts <- function(formula, mmf=mBW, # model matrix function(x, bp, g) data, plot=T, tau=0.5 ) { ... m.nl <- nlrq(y ~ b0 + mBW(x,bp,g) %*% c(b1,b2), tau=tau, start=par0, trace=T )$m ... } For some reason the following reports a failure to find the

Hi All, Could any one tells me if R or S has the capacity to fit nonlinear regression with Huber's M estimation? Any suggestion is appreciated. I was aware of 'rlm' in MASS library for robust linear regression and 'nls' for nonlinear least squares regression, but did not seem to be able to find robust non-linear regression function. Thanks and regards, Ray Liu

I used the nlrq function in the package "quantreg". There is a lambda in the output when I set trace=TRUE. With different start point, the results are converged, but the last lambda is different. I want to know the meaning "lambda=1" and "lambda=0". Many Thanks! Examples of output 1. Where the last lambda=1: 108.6581 : 0.3 8.0 iter 0 value 108.658087 final

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

Dear R Community! We analyse the impact of climbing activity on cliff vegetation. During our fieldwork, we recorded 90 Transects in 3 climbing sites. The aim is to see, if the plant cover (response: Cover) is influenced only by crevice availability (predictor: Cracs), or, additional, by the distance to the climbing route (predictor: Distance). Six plots are nested within one Transect

Hi, I am running a simple linear model with (say) 5 independent variables. Is there a simple way of getting the variance-covariance matrix of the coeffcient estimates? None of the values of the lm() seem to provide this. Thanks in advance, Ritwik Sinha rsinha@darwin.cwru.edu Grad Student Case Western Reserve University [[alternative HTML version deleted]]

From: "Michael S. Tsirkin" <mst at redhat.com> Date: Sun, 10 Aug 2014 21:45:59 +0200 > On Sun, Aug 10, 2014 at 11:30:35AM +0300, Razya Ladelsky wrote: ... > And, did your tests actually produce 100% load on both host CPUs? ... Michael, please do not quote an entire patch just to ask a one line question. I truly, truly, wish it was simpler in modern email clients to

From: "Michael S. Tsirkin" <mst at redhat.com> Date: Sun, 10 Aug 2014 21:45:59 +0200 > On Sun, Aug 10, 2014 at 11:30:35AM +0300, Razya Ladelsky wrote: ... > And, did your tests actually produce 100% load on both host CPUs? ... Michael, please do not quote an entire patch just to ask a one line question. I truly, truly, wish it was simpler in modern email clients to

I need to install Puppy linux as a guest. Usually I boot from the Puppy CD and have so a life system which I can than install to hard disk. How do I do that as a Xen guest? bye Ronald _______________________________________________ Xen-users mailing list Xen-users@lists.xensource.com http://lists.xensource.com/xen-users

I would appreciate confirmation that the function vcov(model.nlme) gives the var-cov matrix of the fixed effects in an nlme model. Presumably the random-effects var-cov matrix is given by cov(ranef (model.nlme)? Rob Forsyth