Anova test with quantile regression model

by StephanEr   Last Updated August 10, 2018 13:19 PM

I want to create a anova test with two or more non-linear quantile regression models:

obj.nlrq.025 <- nlrq(Biomass ~ SSlogis(Time,phi1,phi2,phi3), data=dat,tau=0.25)
obj.nlrq.05 <- nlrq(Biomass ~ SSlogis(Time,phi1,phi2,phi3), data=dat,tau=0.5)
obj.nlrq.075 <- nlrq(Biomass ~ SSlogis(Time,phi1,phi2,phi3), data=dat,tau=0.75)

quantreg:anova.rq(obj.nlrq.025$m, obj.nlrq.05$m, obj.nlrq.075$m)

    Error in formula(x)[[2]] : object of type 'closure' is not subsettable

sig.qr <- anova(obj.nlrq.025$m, obj.nlrq.05$m, obj.nlrq.075$m)

    Error in UseMethod("anova") : 
    no applicable method for 'anova' applied to an object of class "nlrqModel"

Both methods anova and anova.rq from quartreg don't work. With normal nls objects I have no problems.

I want to test whether significant difference exists among these quantiles ...(referred from http://dx.doi.org/10.1590/0103-8478cr20170322 PUIATTI, Guilherme Alves et al. Quantile regression of nonlinear models to describe different levels of dry matter accumulation in garlic plants. )



Related Questions



Quantile crossing with Quantile Regression Forests

Updated November 14, 2017 10:19 AM

Linear Properties of the Quantile Function

Updated May 11, 2017 08:19 AM


Quantile regression with spline smoothing

Updated November 22, 2017 13:19 PM