根据R中的特定分布将范围划分到面元

# 生成正态分布的随机数
set.seed(123)
data <- rnorm(1000, mean = 50, sd = 10)

# 定义面元边界
bins <- seq(min(data), max(data), by = 5)

# 使用cut函数将数据分到面元中
data_bins <- cut(data, breaks = bins, include.lowest = TRUE)

# 查看每个面元中的数据点数量
bin_counts <- table(data_bins)
print(bin_counts)

# 绘制直方图
hist(data, breaks = bins, main = "Normal Distribution Histogram", xlab = "Value")

# 使用等频面元
library(Hmisc)
data_quantiles <- cut2(data, g = 20)  # 将数据分为20个等频面元
quantile_counts <- table(data_quantiles)
print(quantile_counts)

# 绘制直方图
hist(data, breaks = data_quantiles, main = "Equal Frequency Histogram", xlab = "Value")