通信专业学生跨学科请教S-PLUS一个小程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

本人专业是通信与信息系统，偶然发现时间序列分析方法对信号分析非常有用，下面的程序来自范剑青老师的《非线性时间序列-建模、预报及应用》6.2.6局部线性平滑
我的问题是，grid代表什么？实验数据是一个时间序列，我最开始将grid理解为估计值的时刻，method=1也就是local constan fit时，数据靠谱；method=0时，local linear fit，程序算出的数据灰常不靠谱，，，也可能是我理解的错的离谱。

还有一个问题，对平滑后的实验数据进行预测，请问有没有多元局部线性回归预测的S-PLUS实现

感觉是个很基础的问题，但是姐有点急了，希望各位统计专业的好同志能够伸出援手，谢谢了
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 <- 1 - (x-grid[j])^2/h^2; #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,interp=F)$coef[1]
}
list(grid=grid,est=res)
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：PLUS 通信专业 跨学科 Plu 小程序 请教 程序 专业 学生 跨学科

汗，这个问题太外行，大家不屑于回答么

局部线性平滑方法在R中有一个专门的包，此处的grid可以理解为时刻，或称之为时间格点。