求局部多项式回归的程序！ [推广有奖]

0关注
0粉丝

高中生

45%

还不是VIP/贵宾

威望: 0 级
论坛币: 89 个
通用积分: 0
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 538 点
帖子: 24
精华: 0
在线时间: 12 小时
注册时间: 2008-12-28
最后登录: 2010-8-2

楼主

guluzhe 发表于 2010-4-24 08:59:59 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求局部多项式回归的程序！谢谢！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：局部多项式回归局部多项式多项式程序局部多项式

相关帖子

沙发

oddgod 发表于 2010-5-11 22:02:32

lls<-function(x, y, h, grid, method=1,iker = 1)
{ # SIMPLE BUT NAIVE IMPLEMENTATION of
  # local constant or local linear with kernel symmetric beta kernel
  # (1-u^2)^iker and bandwidth h, evaluated at "grid".
  #method = 0, local linear fit;
  #    = 1, local constant fit
  #iker is an integer.
  #result: grid, est
  ngrid<-length(grid)
  res <- rep(0,ngrid)

  for (j in 1:ngrid)
  { w <- locw(x,grid[j],h)          #Epanechnikov kernel
x0 <- x[ w > 0]                #local data
y0 <- y[ w > 0]
w <- w[w > 0]
w <- w^iker
coef <- lsfit(x0-grid[j],y0,w)$coef
if(method==1)
   res[j] <- lsfit(x0-grid[j],y0,w)$coef[1]
else
   res[j] <- lsfit(x0-grid[j],y0,w,intercept=FALSE)$coef[1]
}
  list(grid=grid,est=res)
  }

ep=function(t)
(3/4)*(1-t^2)*(1-t^2>=0)
locw=function(T,t,h){
  y=T-t;Kh=ep(y/h)/h;
  s1=t(Kh)%*%y;s2=t(Kh)%*%(y^2);
  return(Kh*(s2-s1*y)/sum(Kh*(s2-s1*y)));
  }