如题。请问用R软件做非参数回归时怎样求回归系数?下图是一个实例(支出为因变量收入为自变量)和用到的代码,输入coef(fit)时显示的是Null顿时就方了。还请各位指教Orz
install.packages("np")
library(np)
library(boot)
install.packages("cubature")(似乎是无必要的代码因为一步骤已经包含)
library(cubature)
# 把数据文件放在根目录文件的下面,或指定你的根目录到数据文件所在的文件夹
dat=read.table("income.txt",head=TRUE)
cuh=dat$income
icup=dat$payment
plot(cuh,icup,xlab="城镇居民的可支配收入",ylab="城镇居民的消费支出")
# 使用非参数h核回归方法,通过最小二乘交叉验证方法(cv.ls)选择最优带宽,核函数为高斯核
bw1=npregbw(icup~cuh,regtype="lc",bwmethod="cv.ls",bwscaling=FALSE,ckerorder=2,ckertype="gaussian")$bw # 最优带宽 bw=942.80058
fit=npreg(icup~cuh,bws=bw1,gradients=TRUE,residuals = TRUE,regtype="ll")
summary(fit) # 拟合优度 R^2=0.999914