全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
708 0
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[] <- samples[[order[]]]
}

#s^2(a)
s2 <- function(a){
  #m(α)
  sumX <- 0
  for(i in (floor(a*n)+1):(n-floor(a*n))){
    sumX <- sumX+X[]
  }
  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[]-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)
}
二维码

扫码加我 拉你入群

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

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

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

分享

扫码加好友,拉您进群