请求将地震数据处理的R代码改写为Base R实现

1. 将深度数据划分为因子分组 #

用Base R内置的cut()函数对地震深度进行区间划分，并直接转换为因子类型。示例中按0-300km、300-600km、600km以上分为三组，你可根据需求调整breaks和labels参数：

# 划分深度区间并转为因子
quakes$depth_group <- cut(quakes$depth, 
                          breaks = c(0, 300, 600, Inf),
                          labels = c("浅源地震", "中源地震", "深源地震"))

2. 绘制分组箱线图 #

使用Base R原生的boxplot()函数，通过公式语法mag ~ depth_group指定按深度分组绘制震级的箱线图，还可自定义标题、坐标轴标签和配色：

# 绘制不同深度分组的震级箱线图
boxplot(mag ~ depth_group, 
        data = quakes,
        main = "不同深度分组的地震震级分布",
        xlab = "深度分组",
        ylab = "震级",
        col = c("#66c2a5", "#fc8d62", "#8da0cb"))

3. 执行Kruskal-Wallis检验 #

直接调用Base R的kruskal.test()函数，同样用公式语法传入数据，即可完成非参数检验：

# 执行Kruskal-Wallis检验
kw_result <- kruskal.test(mag ~ depth_group, data = quakes)
# 输出检验结果
print(kw_result)

以上代码完全依赖Base R内置函数，替代了dplyr的分组处理、ggplot2的绘图功能，输出结果与tidyverse版本完全一致，满足作业要求。

内容的提问来源于stack exchange，提问作者tetra4892