投稿：地震越来越频繁了吗？ 黄湘云

[XiangyunHuang]

非常感谢您的PR, 如果您是在为主站投稿, 请将PR的标题改为"投稿:标题+作者的形式",如: "投稿: 数据通灵术 杜亚磊" 并保留下面的内容.

• 投稿者请注意
  • [x] 主编审核确认接收
  • [x] 主编安排审稿人
  • [x] 审稿、修改
  • [x] 已发其他公众号的文章，请在投稿前注明相关公众号信息，并建议给统计之都公众号开白名单

至此，投稿部分的工作结束

• 编辑部分工作
  • [x] 文字编辑做校对工作。需校对文章分类，目前的文章分类包括：COS访谈，R会议，R语言，推荐文章，新闻动态，机器学习，统计之都，统计图形，统计应用，统计模型，统计计算，统计软件，职业事业。不得擅自增加分类。
  • [x] 论坛帖子负责人发帖，粘贴链接，提供 forum_id 加入文章
  • [x] 合并 PR，发布文章，转交微信端编辑（项目管理中标注PR编号跟进）
  • [ ] 微信编辑、发布，阅读原文附主站链接，文末附作者介绍

投稿指南在这里，有任何问题，可以直接在PR留言，其他问题请联系: editor@cos.name。

[XiangyunHuang]

@hongtaoh 鸿涛，你有兴趣来看看这篇文章吗？提一些修改建议等。

[hongtaoh]

@hongtaoh 鸿涛，你有兴趣来看看这篇文章吗？提一些修改建议等。

可以啊，但我貌似看不了，不知道怎么看。

[XiangyunHuang]

@hongtaoh 我们是用 Github Pull Requests Review 来审阅的，文章在主站的预览效果，见这 https://deploy-preview-1066--physicist-ralph-51373.netlify.app/。这篇文章是 R Markdown 文档生成的，对 R Markdown 文档提修改意见就行。

[hongtaoh]

@XiangyunHuang 好的！我看完了。整体而言我觉得非常好，图很精美，代码也很简洁。我但不知道如何直接在 rmarkdown 上提意见，就把我的意见写在这里了：

• 3.2 节讲完之后，从图中看出次数确实增加。但一个解释原因是检测能力变强，你也提到了。但你下面讲的，我不知道你为什么突然要开始分析震级的频率分布。你的逻辑可能是，有可能，1990 年之后虽然次数增加，但有不少是小震级地震，而 1990 年之前很多是大震级地震。但你要把这个逻辑在分析震级之前说明白，这样读者读起来就容易一些。

• 你下面的那几个图来看，好像并不容易判断是否 1990 年之后小震级地震变多。根据图 3.2，我在想，你在看震级按年份的分布时，是不是可以只关注 6-7 级这个区间，这样更容易看到，是不是1990 年之后，小震级的出现次数越来越多，这样才能佐证是否 1990 年之后，随着科技发展，小震级被发现的频率增加。

• 分析中国地震变化时，在原始数据无法读取的情况下你做的各种尝试，我觉得非常好。

• 你的图 4.3 显示，如果只关注中国监测到的 6 级以上地震，貌似频率没有增加。但图 3.1 显示世界范围内，6 级以上地震是增加了。这个差别你也许可以探讨一下。

[XiangyunHuang]

@hongtaoh 谢谢，修改建议都很好，预计这两周我会找时间改下。

[XiangyunHuang]

文章的组织也会再调整下，开篇提出问题，回答问题，然后再讲分析思路和探索过程，最后，简单总结一下。最近，手头事情比较多，文章修改的进度预计会比较慢。

[XiangyunHuang]

找到了美国地震局数据的下载入口 https://earthquake.usgs.gov/earthquakes/search/ 到时候把美国的数据更新一下，这样中美相互对比和校验，从数据中获得的一些地震活动规律，更加确认了。

[XiangyunHuang]

找到了美国地震局数据的下载入口 https://earthquake.usgs.gov/earthquakes/search/ 到时候把美国的数据更新一下，这样中美相互对比和校验，从数据中获得的一些地震活动规律，更加确认了。

原想和中国地震台网那样，有多少数据都下载下来，实际上，美国地震局数据下载限制为 20000 行/次。因此，更新美国地震数据，仍然仅考虑 6 级及以上全球范围 1976 年至今的数据。

[XiangyunHuang]

因此，更新美国地震数据，仍然仅考虑 6 级及以上全球范围 1976 年至今的数据。

在中国地震台网，6 级及以上的地震关注度高，会标记红色。可以对比看看中美 6 级及以上的全球范围的数据。

[XiangyunHuang]

中美对比暂时放弃了。检查数据质量是个无底洞呀，最终取决于政府选择性公开的质量。比如美国会公开一些中国及其他地区疑似的核爆炸或大爆炸导致的地震，反正就是没有美国人自己的，当然，也有人会说这正好把脏数据过滤了呀，嗯，是选择性过滤掉了。毕竟一个国家的地质构造信息还是挺机密的。中国政府发布的地震位置的经纬度精度很有限。所以，本文数一数地震的次数及其分布规律还是有一些意义的，当然还有一些可视化的分析方法，其他的别太当真。无论怎么说，比起新闻报道的内容还是强多了，客观得多了。

[XiangyunHuang]

@hongtaoh 你提的意见非常好，我在开头新加了一节，说明后续的分析思路、一些分析结论和未来还可以继续往哪些方面分析。

1. 补充了本文结构，交代分析思路和结论。
2. 只关注 6-7 级，似乎不够，地震监测能力要结合震级、和震深一起来看，找到能力变化的时间点changepoint，然后分析这个点前后的变化情况。这个我没有做下来，中国的数据不够。美国的数据，因每次下载的数据量有限制，我又没有找到一个批量爬数据的办法，最后，我只下载了 6 级及以上的。
3. 谢谢。
4. 解释这个差别，要先找到检测能力的拐点，并把前后的数据找到。原因见 2，我没有能力把数据下载下来。

[XiangyunHuang]

找 changepoint 分离效应是数据分析里很重要的事，后续，如果可以把数据弄下来，再就这个话题深入分析一波，更加聚焦一点，本文太发散了，似乎什么都写了一点，这篇就当是踩坑积累经验教训了。changepoint 是否有在官网发布呢？目前，我也没找到。

[XiangyunHuang]

批量下载数据的秘方藏在链接里，改时间段就可以了。

https://earthquake.usgs.gov/fdsnws/event/1/query.csv?starttime=2022-12-01%2000:00:00&endtime=2022-12-31%2023:59:59&minmagnitude=0&orderby=time

[XiangyunHuang]

@hongtaoh 我应该改完了。麻烦有空的时候，再来看看哈 😄

[hongtaoh]

@hongtaoh 我应该改完了。麻烦有空的时候，再来看看哈 😄

刚才看了。这次没有上次看得仔细，我只是大体上看了一下全篇。抛开内容不说，整个挖掘数据、寻找数据的过程就很感人，看得出来花了很多时间在上面，恭喜最后把它做完！这个探索精神是很有感染力的。

作图很精美，估计大家看到后会更有学 R 的动力。

[XiangyunHuang]

@hongtaoh 谢谢鼓励😄。那文章就进入编辑环节了。