阿雷西博天文台的科学家发表了关于近地小行星的主要研究

最近芬兰的一项研究报告称，满月期间的自杀率升高，但仅限于绝经前妇女和冬季。这无法在奥地利的一项研究中复制，并引发了关于芬兰的发现是否为假阳性或芬兰和奥地利之间是否存在未解释的调节变量不同的讨论。本研究的目的是提供另一个来自瑞典的数据，瑞典在地理上与芬兰的可比性高于奥地利。我们还研究了讨论的调节变量纬度和夜间人工亮度。有 48,537 起自杀事件可供分析。冬季满月季的自杀比例与绝经前女性(23.3%)和完整样本(24.7%)中预期的 25% 没有显着差异。泊松回归模型中满月季的发病风险比对于绝经前女性为 0.96(95% CI：0.90-1.02)，对于完整样本为 1.01(95% CI：0.99-1.04)。根据贝叶斯因子分析，证据支持零假设(无关联)而不是替代假设(某些关联)。当我们分别按纬度和人工夜间亮度分割数据时，我们发现了类似的结果。与奥地利的研究一致，瑞典的绝经前妇女在冬季满月期间的自杀率并没有增加。芬兰研究的结果很可能是假阳性，

介绍

尽管缺乏令人信服的证据，但人们对月相(即月相周期)对自杀行为的推定影响反复出现兴趣。事实上，有大量研究表明月球对自杀行为没有影响，而且这个问题被认为很久以前就已经解决了，1985在纬度、农村或性别方面没有差异 [ 1 ]。1992 年的一篇评论同样得出结论，“没有足够的证据表明会合月球周期与自杀或自杀未遂之间存在关系”[ 2]。从那时起，月相与自杀行为或精神​​病入院之间的关联被重新检查了几次，而没有改变被广泛接受的无效结论 [ 3 , 4 , 5 ]。

然而，最近，Meyer-Rochow 等人在芬兰北部奥卢区进行了一项研究。(2021 年)(以下简称 MR)发现满月期间自杀率升高，但仅在绝经前妇女中且仅在冬季 [ 6 ]。这一发现是月相与自杀之间的统计学显着关联，这是在涉及性别、年龄和季节的亚组分析中发现的。这种数据分析似乎是探索性的，可能是对亚组的事后定义，从而大大增加了假阳性结果的风险 [ 7 , 8 , 9 , 10 ]。因此，Plöderl 和 Hengartner [ 11](以下简称 PH)认为，这一惊人且相当令人惊讶的发现可能是一个统计伪影，即由于抽样变异性而导致的假阳性机会发现。为了支持他们的论点，PH 试图用来自奥地利的更高功率数据复制 MR 的发现，但发现在冬季满月阶段绝经前妇女的自杀率没有增加。匿名审阅者和 Meyer-Rochow 等人对这次失败的复制做出了回应。[ 12] 认为，奥地利和芬兰北部结果之间的差异可能是由无法解释的调节变量来解释的，例如纬度、昼夜长度、夜间人工亮度和维生素 D 缺乏。原则上，他们的主要论点是 MR 报告的主要发现不能被 PH 与奥地利的数据复制，因为芬兰北部比奥地利更北，因此，冬季满月对绝经前妇女的影响仅在芬兰北部很明显。

我们目前研究的主要目标是用瑞典的数据复制 MR 的主要发现，瑞典在地理上与芬兰北部相比比奥地利更具可比性。此外，我们还研究了一些建议调节变量的影响，以解释奥地利和芬兰研究结果之间的差异，即纬度和夜间人工亮度。因此，本研究的结果还应该增加新的见解，这些新见解与假设涉及假定的会合月球周期对心理健康的影响的因素有关。

自杀数据

自杀数据提取自瑞典国家卫生和福利委员会于 2022 年提供的匿名瑞典死因登记处。根据瑞典国家健康与福利委员会的故意自残相关代码(X60–X84 或 E950–E959)确定自杀。国际疾病分类 (ICD)，(第 8-10 版)。我们只包括对自杀意图毫无疑问的案例。登记处包括了 1980 年至 2021 年间的所有此类病例。主要分析中排除了有关确切死亡日期和年龄的数据缺失的病例，有关纬度和夜间亮度的亚组分析排除了有关城市的数据缺失的病例。瑞典共有290个市镇，市镇数据指死者死亡时的居住地，根据瑞典税务局的民事登记。这个变量为自杀提供了最细粒度的地理信息。在协议中，我们打算使用发生自杀事件的市镇，而不是居住市镇。这偏离了协议，因为我们可以获得仅在 2015-2021 年发生自杀事件的市政当局的信息。然而，这不应该影响结果，因为大多数自杀发生在居住地。例如，在瑞典最大的城市斯德哥尔摩居住的自杀死亡者中，80% 的人死于斯德哥尔摩，13% 的人死于斯德哥尔摩附近的城市，7% 的人死于更远的城市。

月相

每天的月亮照度是使用 R 的 lunar Package [ 13 ] 计算的，它为每个日期和月相类别提供连续的照度测量。四类措施将月亮的照度分成相等大小的四分之一(满月、残月、新月、上盈月);八类措施将月亮的照明分成八类等长。与 MR 类似，我们将四级月相变量分解为具有三类(满月、新月、亏月/打蜡月组合)的变量。对于敏感性分析，我们使用了八级月相。

其他变量

季节的定义方式与 MR 相同：冬季(11 月、12 月、1 月)、春季(2 月、3 月、4 月)、夏季(5 月、6 月、7 月)和秋季(8 月、9 月、10 月)。年龄组也以与 MR 相同的方式定义，即 45 岁以下的女性被归类为绝经前，而 45 岁或以上的女性被归类为绝经后。对于自杀的纬度，我们使用 R 的 swemaps 包 [ 14 ] 来提取瑞典城市的地理信息，即每个城市对应的多边形。然后我们计算了每个城市的地理中心，并用它来计算每个自杀事件的纬度(见 补充图)。瑞典每个城市的夜间亮度是通过将城市的地理信息与世界人工夜空亮度图集联系起来计算得出的 [ 15 ]。我们使用了为每个城市提取的数据点的平均值(见 补充图)。

统计分析

我们使用与 MR 相似的统计方法： χ 2-测试用于调查自杀和月相的性别差异。多项检验用于比较三个月相中观察到的自杀与预期的自杀，也就是说，我们预计在亚组(男性、女性，绝经前妇女)在冬天。我们使用与 MR 类似的泊松回归模型来调查每日自杀是否因月相而异，以亏月/上盈月为参考类别，根据季节进行调整，以冬季为参考类别。除了MR使用的统计方法，我们计算了贝叶斯因子 (BF)，以量化有多少证据支持零假设(与其他月相相比，满月期间的自杀发生率为零)而不是替代假设(满月期间或多或少的自杀)。对于这些二项式检验，我们使用了 BayesFactor 包 [16 ]。我们根据 Wagenmakers [ 17 ]描述的 Raftery (1995) 解释了 BF 。由于大量分析，预计统计显着的发现仅是偶然发生的(对于p < 0.05的阈值，大约 20 分之一 )，因此我们计划使用顺序 Bonferroni alpha 调整。我们使用 R [ 18 ] 进行所有分析。R 代码可通过开放科学框架 (OSF) https://osf.io/cfq3k/获得。

二次分析

对于二次分析，我们通过运行亚组分析研究了两个讨论过的潜在调节变量(纬度和夜间亮度)的影响。为此，我们使用中值对夜间亮度变量进行了二分法。对于纬度，我们使用瑞典的平均纬度。这与协议有偏差，我们在协议中声明使用城市纬度的中值，但瑞典南部的中值太远了，因为大多数城市都位于南部。因此，这将阻止估计瑞典北部和南部之间的差异。对整个样本和绝经前妇女分别进行亚组分析。

敏感性分析

我们重复了主要分析，仅限于纬度与奥卢相当甚至更北的城市(以便与 MR 具有可比性)。我们决定将纬度的界限也包括 Skellefteå 市。这个城市位于奥卢以南约 50 公里处，但它是较大的城市之一，因此将增加统计能力(见 补充图)。在另一项敏感性分析中，我们使用了 1988 年 3 月至 2011 年 6 月的数据，类似于 MR。最后，我们使用八个月相级别重复分析，正如 PH 所做的那样。

预注册和数据可用性

研究方案于 2022 年 4 月 14 日在 OSF https://osf.io/cfq3k/上上传。本研究使用的数据最初由瑞典国家卫生和福利委员会于 2022 年 7 月 11 日提供，数据分析。出于道德原因，未提供来自瑞典死因登记处的原始数据，只有隶属于瑞典国家自杀研究和预防中心 (JW, SH, GH) 的作者可以访问。上传研究协议后，R 代码由 MP 在 1000 名自杀者的虚拟数据集上开发，并由 JW 首次应用于最终数据集 2022-07-13。支持独立研究人员复制的方法学信息可从通讯作者处获得。

样品特征

1980 年 1 月 1 日至 2021 年 12 月 31 日期间发生了 54,942 起自杀事件。其中，已知死亡日期为 48,537 人(88.3%)，其中 14,326 人(29.5%)为女性自杀，34,211 人(70.5%)为男性自杀。在 10 例中缺少年龄数据，在 41 例中缺少市镇数据。有 5580 名绝经前妇女的自杀事件可用于主要分析。