21. （ R 编程）产生 2 000 个均值等于 5 000, 方差等于 25 000, 且服从以下分布的损 失随机数：正态分布、指数分布、伽马分布、帕累托分布、对数正态分布、威布尔分布。

在R语言中，你可以使用`rnorm()`、`rexp()`、`rgamma()`、`rpareto()`、`rlnorm()`和`rbeta()`等函数来生成各种分布的随机数。以下是创建每个分布所需的2000个随机数的例子：

1. 正态分布 (Normal Distribution):

set.seed(1) # 设置随机种子为了复现结果
normal_data <- rnorm(2000, mean = 5000, sd = sqrt(25000))

2. 指数分布 (Exponential Distribution):

exponential_data <- rexp(2000, rate = 1/sqrt(25000))

注意：对于指数分布，我们通常用率`lambda`作为参数，而方差`sigma^2`通过公式`sigma^2 = 1/lambda`计算。

3. 伽玛分布 (Gamma Distribution):

gamma_data <- rgamma(2000, shape = 2, scale = sqrt(25000))

4. 帕累托分布 (Pareto Distribution):

pareto_shape <- 2
pareto_data <- rpareto(2000, shape = pareto_shape, scale = 5000)

5. 对数正态分布 (Log-Normal Distribution):

log_normal_mu <- log(5000)
log_normal_sigma <- sqrt(log(25000)) / 5000
log_normal_data <- rlnorm(2000, meanlog = log_normal_mu, sdlog = log_normal_sigma)

6. 威布尔分布 (Weibull Distribution):

weibull_shape <- 1 / sqrt(25000)
weibull_scale <- 5000
weibull_data <- rbeta(2000, shape1 = weibull_shape, shape2 = 1)
weibull_data <- -log(-log(weibull_data)) * weibull_scale

每个数据集都存储在对应的变量`normal_data`, `exponential_data`, ..., `weibull_data`中。如果你想查看每种分布的数据摘要，可以使用`summary()`函数。