[问答] 自定义函数问题：单独令a为0.1，能运行正确结果，但function(0.1)错误。 [推广有奖]

0关注
0粉丝

高中生

77%

还不是VIP/贵宾

威望: 0 级
论坛币: 2 个
通用积分: 1.0000
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 229 点
帖子: 10
精华: 0
在线时间: 45 小时
注册时间: 2018-7-27
最后登录: 2023-4-22

楼主

EmmaSwan13 发表于 2019-3-28 15:56:09 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

错误提示：
> 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)
}