有偿求助！！请修改一下我写的程序。。

R编写的局部多项式回归程序，刚接触R，不知道如何修改。。
麻烦大家帮我修改一下。。谢谢！！！！！

LPR=function(n,c,p)
{
    e=rnorm(n,0,0.25)
    x<-seq(0,1,length=n)     
    y=sin(2*pi*x)+e
    z=sin(2*pi*x)
    plot(x,y)
    lines(x,z,lwd=2)
    h=c*n^(-1/5)
    mx=rep(0,n)
    for(i in 1:n)
     {
         mx=sum(((y-qx(x,x,p)%*%(solve(t(X(x,x,n,p))%*%W(x,x,h,c)%*%X(x,x,n,p))%*%t(X(x,x,n,p))%*%W(x,x,h,c)%*%y))^2)*K(x,x,h))         
     }
    lines(x,mx)
}
##########################################################
qx=function(t.val,t,p)
{
    qx=rep(1,p+1)  
           for(i in 2:p+1)
             qx=qx*((t.val-t)^(i-1))
         
}
##########################################################
X=function(t.val,t,n,p)
{
   X=matrix(1,nrow=n,ncol=p+1)
            for(i in 1:n)
              for(j in 2:p+1)
                 X[i,j]=X[i,j]*((t.val-t)^(j-1))
         
}
##########################################################
W=function(t.val,t,h,c)
{
    h.d=c*n^(-1/5)
           h.s=1/h.d
           diag=matrix(0,n,n)
           for(i in 1:n)
              diag[i,i]=h.s*K(t.val,t,h)
           
           
}
############################################
K=function(t.val,t,h)
{
    return(kernel((t.val-t)/h))
}
############################################
kernel<-function(x)
{
    return ((15/16)*(1-x^2)*(abs(x)<=1))
}

先关注一个，对此不是很了解，等有时间了，好好研究一下。

好多错误，例如第一个循环没有i 呢？如果你需要局部多项式回归的R  code ，我可以
帮助你。

这个我毕业论文刚做过。

你的程序问题确实挺多的。很多细节语法没注意

比如，你的循环for(i in 2:p+1)应该改成for(2  in 2:(p+1))
另外你的t.val是向量还是单个数值？如果是向量，那么你的编写的核函数肯定运行出错。分段函数不允许整个向量直接计算。

6# 楚韵荆风
谢谢了，我还是个初学者，请问如果想深入学习R软件应该看哪本书?

错误于solve.default(t(X) %*% W %*% X) :
Lapack例行程序dgesv: 系统正好是奇异的
怎么会这样？？

6# 楚韵荆风

你的这个程序有错！！

他的肯定有问题，他说你的循环for(i in 2:p+1)应该改成for(2  in 2:(p+1))
后者明显就是错的，你的是对的，我已经给你发了这个CODE，你是想做哪类模型？
一般非参数？变系数？还是单指标模型？或是偏线性模型？