全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 经管代码库
832 1
2019-03-28
错误提示:
> function(0.1)
Error: unexpected numeric constant in "function(0.1"
代码:
n <- 100
X <- list()
#生成X(1),...,X(n)
samples <- sample(1:(10*n), n)
order <- order(samples, decreasing=F)
for(i in 1:n){
  order
  X[[i]] <- samples[[order[[i]]]]
}

#s^2(a)
s2 <- function(a){
  #m(α)
  sumX <- 0
  for(i in (floor(a*n)+1):(n-floor(a*n))){
    sumX <- sumX+X[[i]]
  }
  m <- 1/(n-2*floor(a*n))*sumX
  #Var(a)
  VarX <- 0
  for(i in (floor(a*n)+1):(n-floor(a*n))){
    VarX <- VarX+(X[[i]]-m)^2
  }
  VarX <- VarX/n
  result <- 1/(1-2*a)^2*(VarX+a*(X[[floor(a*n)+1]]-m)^2+a*(X[[n-floor(a*n)]]-m)^2)
  return(result)
}

二维码

扫码加我 拉你入群

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

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

全部回复
2019-3-28 16:54:23
您好,如果您的求助没有解决,请到项目交易发布需求,会有更快更专业的用户帮助您 https://bbs.pinggu.org/prj/
二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群