使用R语言：请从 words 中每次取 x 个单词，统计辅音结尾的比率，并将其重复 n 次，将其写成 函数。 注意，函数变量有两个，为：“抽取单词个数” x 和 “重复次数” n。请用 x=10, n=5000 测试你的函数并显示它是可运行的。要求运行结果产生一个新表,里面变量 nonvowel_ratio, 同时产生一个直方图，并伴有密度曲线。请问它是什么分布？为什 么？ (15 分)

以下是解答：

首先，我们需要写一个函数来计算每次抽取x个单词时辅音结尾的比率。

consonant_ratio <- function(x) {
  words <- sample(words, x)
  consonant_count <- sum(sapply(words, function(w) length(gsub("[aeiouAEIOU]", "", w))))
  ratio <- consonant_count / (x * nchar(words[1]))
  return(ratio)
}

接着，我们可以将其重复n次，并将每次计算出的比率存入一个向量中。

nonvowel_ratio <- numeric(n)
for (i in 1:n) {
  nonvowel_ratio[i] <- consonant_ratio(x)
}

最后，我们可以将结果写成一个数据框，并绘制直方图和密度曲线。

results <- data.frame(nonvowel_ratio)
hist(results$nonvowel_ratio, freq = FALSE, main = "Distribution of non-vowel ratios", xlab = "Ratio")
lines(density(results$nonvowel_ratio))

根据中心极限定理，当样本数量足够大时，样本均值的分布会趋向于正态分布。因此，这个分布可能是正态分布。