Displaying 1 result from an estimated 1 matches for "supw".
Did you mean:
sunw
2012 May 28
2
R quantreg anova: How to change summary se-type
...egression
(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 is using the summary.rq function I want to change the Standard
error method used (default: se="nid" leads to mistakes, I prefer se="ker").
Do you know how to handle this information in the anova syntax?
Thank you very much
Stefan
--
View this message in context...