用非参local polynomial方法。看到很多论文上有非参数回归拟合值与原始值对比的表
如:
时间 原始值 拟合值
1990 12 11.2
1991 13.2 13.4
1992 12.5 12
..........
,但是我只会做非参数回归拟合图。
请教大家编程命令得到拟合数据,谢谢!
附上我的命令
require(KernSmooth) data <- read.table('C:/data.csv', header = TRUE, sep = ',')
h <- dpill(data$LLL, data$HHH) fit <- locpoly(data$LLL, data$HHH, bandwidth = h) plot(data$LLL, data$HHH) line(fit)
|