利用R（不能使用内置函数）编写一个函数, 求数据 y <- sample(rnorm(10000), 1000, replace = TRUE) 的均值, 标准差, 偏度与峰度.

# 定义函数，输入为数据 y
my_summary <- function(y) {
# 求均值
mean_val <- sum(y) / length(y)

# 求标准差
sd_val <- sqrt(sum((y - mean_val)^2) / (length(y) - 1))

# 求偏度
skewness_val <- sum((y - mean_val)^3) / (length(y) * sd_val^3)

# 求峰度
kurtosis_val <- sum((y - mean_val)^4) / (length(y) * sd_val^4) - 3

# 返回结果
return(c(mean_val, sd_val, skewness_val, kurtosis_val))
}

# 生成数据
y <- sample(rnorm(10000), 1000, replace = TRUE)

# 调用函数，输出结果
my_summary(y)