for语句的一个具体应用的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

lar<-function(a,b,x,fun){
w1<-1
s<-0
  if(length(a)!=length(b)) m<-c("有关插值节点的变量输入错误，请检查节点处函数值是否对应");m
      
  else {
         for(i in 1:length(a)){
           w1<- w1*(x-a)
           }
         for(j in 1:length(a)){
            p<-a[j]-a[-j]
            u<-prod(p)
            t<-fun(a[j])
            q<-w1*t/(u*（x-a[j])
            s<-s+q
            j=j+1
            s
           }
        
       }
}
a<-c(0.4,0.5,0.7,0.8,3.2)
b<-c(0.4,1.2,0.6,0.7,0.8)
fun<-function(x) log(x)
lar(a,b,2,fun)

         整个程序是想实现拉格朗日插值算法，红色字体是预计能运算(Xi-X0 )…(Xi-Xi-1 )(Xi-Xi+1 )…(Xi-Xn )
         并且在与Yi,W1做乘积后累加，所以我使用了for循环，但是出现了很大偏差
         Error: object 's' not found
>             j=j+1
>             s
            恳请朋友们帮着检查错误！
　　　　　any suggestion will highly be appreciated. thanks in advance

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：For R语句 LAR 统计 运行 可行性

j在for循环里面已经自动+1了。你再加一就是加二

同意j加1的问题

紫血丹心 发表于 2016-4-3 19:03
同意j加1的问题

确是，这次已经没有报错了