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