错误报告 为什么？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

                10^at
            else at)[-c(1, ny + 1)]
        }
        abline(h = at, col = col, lty = lty, lwd = lwd)
    }
}
<bytecode: 0x04e5a070>
<environment: namespace:graphics>
$est
[1] 0
>   }
错误: 意外的'}'在"  }"里
>
>
> 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)));     
+

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：错误报告 function abline RETURN line 错误报告

帖子不完整，请编译一下。

> 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)));     
+   }
>

友情帮顶

伊树书华 发表于 2012-5-14 17:34
> lls

试了一下，没有问题。

继续帮顶