我想在R中使用分位数估计H-D的非线性关系,选择了3个分位数(0.1,0.5,0.9),可为什么估计出来的参数只有一组呢?不应该是每个分位数都会有一组参数吗?请各位帮忙!以下附上运行命令及输出结果。
library(quantreg)
rq01<-nlrq(H~1.3+a*((1-exp(-b*D))^c), tau=c(0.1,0.5,0.9), data=wholeplot,start=c(a= 27.49283027, b= 0.03288259, c= 1.09875097) , trace=TRUE)
[size=13.3333330154419px]print(rq01)
[size=13.3333330154419px]Nonlinear quantile regression
[size=13.3333330154419px] model: H ~ 1.3 + a * ((1 - exp(-b * D))^c)
[size=13.3333330154419px] data: wholeplot
[size=13.3333330154419px] tau: 0.1 0.5 0.9
[size=13.3333330154419px]deviance: 4343.479
[size=13.3333330154419px] a b c
[size=13.3333330154419px]27.70278258 0.03259709 1.09796018
AIC(rq01)
[1] 35689.77 24075.63 35689.77
attr(,"edf")
[1] 3