楼主: efish1
3252 1

[问答] 如何用R来做出local polynomial拟合数据 [推广有奖]

  • 0关注
  • 0粉丝

硕士生

66%

还不是VIP/贵宾

-

威望
0
论坛币
3309 个
通用积分
2.4000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
597 点
帖子
22
精华
0
在线时间
353 小时
注册时间
2009-5-14
最后登录
2024-1-8

50论坛币
用非参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)



最佳答案

hugebear 查看完整内容

locpoly函数会在做拟合时选取一些等间距的grid points(默认值是401个), 然后在这些点上给出局部多项式回归的拟合值。如果你想要的是grid points上的拟合值,用fit$y就可以得到(这些拟合值已经足够细密以反映被估计函数的全貌)。如果你想要的是拟合曲线在解释变量HHH上的拟合值,似乎没有太简便的操作。可以自己根据local polynomial的算法自己写个函数。
关键词:Polynomial nomial Local poly OCA local polynomial R代码 拟合值
沙发
hugebear 发表于 2015-5-7 18:03:47 |只看作者 |坛友微信交流群
locpoly函数会在做拟合时选取一些等间距的grid points(默认值是401个), 然后在这些点上给出局部多项式回归的拟合值。如果你想要的是grid points上的拟合值,用fit$y就可以得到(这些拟合值已经足够细密以反映被估计函数的全貌)。如果你想要的是拟合曲线在解释变量HHH上的拟合值,似乎没有太简便的操作。可以自己根据local polynomial的算法自己写个函数。
已有 1 人评分经验 收起 理由
静水深流 + 20 精彩帖子

总评分: 经验 + 20   查看全部评分

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-3 21:38