The need for improved models of zooplankton abundance

技術的要約：動物プランクトン個体数モデルの改善の必要性

問題提起
動物プランクトンの個体数データは、海洋生態系および栄養段階の解析において基本的であるが、特有の統計的課題を提示している。これらのデータセットは、異なるサンプリング・プロトコル（具体的には、異なる垂直サンプリング深度区間、すなわちビン幅）を用いる複数の調査から集計されることが多い。広いビンはより大きな水塊を統合するため、精度の異なる個体数推定値をもたらす。これらを考慮に入れない場合、不均一性が導入され、統計的推論にバイアスが生じる可能性がある。さらに、動物プランクトンの体積あたりの個体数は厳密に正の値であり、強い右裾の長い分布（右偏分布）を示す。この分野における一般的な慣行は、応答変数を対数変換し、正規性を近似するためにガウス線形モデルを適合させることである。しかし、この変換に基づくアプローチは、推論にバイアスを与え、分布の感度を不明瞭にし、元のスケールにおける係数の解釈を困難にする可能性がある。加えて、個体数推定の精度がサンプリング深度区間の関数として明示的にモデル化されることは稀であり、これが過度に自信に満ちた結論を導く可能性がある。

手法
本研究では、メソペラジック（中層）のメソゾプランクトンおよびマイクロネクトン・データベース（MMMD）からの動物プランクトン・サンプリング記録を使用し、1987年の南極海における5種の豊富なカラヌス目カイアシ類（Rhincalanus gigas, Metridia lucens, Calanus simillimus, Pleuromamma robusta, Ctenocalanus vanus）に限定した。応答変数は体積あたりの個体数（個体 $\cdot$ m$^{-3}$）であり、共変量には温度、塩分、溶存酸素、および平均サンプリング深度が含まれた。

著者らは4つの統計的枠組みを比較した：

1. ガウス線形モデル： 対数変換された個体数に対して適合（一般的な慣行）。
2. 一般化線形モデル (GLMs)： 以下の3つの非ガウス誤差分布を用いて、元の個体数スケールに直接適合：
   • 対数正規分布
   • ガンマ分布
   • 逆ガウス分布

すべてのモデルは、温度、平均深度、およびそれらの相互作用を含む線形予測子を利用した。モデル選択は、対数変換されたガウスモデルに対して、元の応答スケールでの妥当な比較を確保するためのヤコビアン補正を適用した上で、一般化赤池情報量基準 (GAIC) を用いて行われた。

最適な支持を得た分布の選択後、著者らは分散を正規化されたビン幅（最大サンプリング深度と最小サンプリング深度の差）の関数として明示的にパラメータ化するようにモデルを拡張した。この「分散モデル」により、サンプリング区間のサイズに基づいて分散を変化させることが可能となった。最終的なモデルの結論の堅牢性を評価するために、感度分析（一回抜きによるケース削除）が行われた。

主な結果

• 分布の適合度： 逆ガウスGLMが、すべての候補モデルの中で最も強い支持を得た（GAIC = -3068.25）。これは、ガンマモデル、対数正規モデル、および対数変換ガウスモデルよりも優れた結果であった。対数変換ガウスモデルと対数正規モデルは同一の係数を示したが、モデルの適合性の観点からは大幅に劣っていた。
• 分散モデリングの影響： ビン幅に依存する分散項を組み込むことで、モデルの適合性がさらに向上した（GAICは -3084.61 に減少）。分析の結果、ビン幅と残差変動の間に正の関係があることが明らかになり、広いサンプリング区間が推定精度の低下を招くことが確認された。
• 推論の変化： 偏差（分散）を考慮することで、共変量の統計的有意性が変化した。具体的には、定数分散の逆ガウスモデルにおいて統計的に有意であった温度と深度の相互作用（$p = 0.0023$）は、精度の不均一性を組み込んだ際に非有意となった（$p = 0.0657$）。一方で、深度と個体数の間の負の関連性は、すべてのモデルにおいて堅牢な予測因子であり続けた。
• 残差診断： 逆ガウスモデル（分散あり）は候補の中で最良の適合を示したが、残差プロットは持続的な下限側の圧縮と可視的な構造を示しており、線形な平均構造が生物学的不均一性（例：種特異的な反応や非線形な深度効果）を完全には捉えきれていないことを示唆している。
• 感度： 深度効果およびビン幅による分散効果に関する主要な結論は、一回抜きによる再適合においても安定していた。しかし、温度と深度の相互作用に関する証拠は、少数の影響力の強い観測値の除去に対して敏感であった。

意義および主張
本論文は、右裾の長い、厳密に正の値をとる動物プランクトンの個体数データに対して、日常的な対数変換ガウスモデルの使用は最適ではない可能性があると主張している。逆ガウス分布がより良い適合を示すことを実証することで、本研究は、元のスケールでの直接的なモデリングが、より正確な係数推定と有意水準をもたらすことを示唆している。

決定的な点として、著者らは、サンプリング設計における不均一性（具体的にはビン幅の変動）を無視することは、過度に自信に満ちた推論を招く可能性があると主張している。ビン幅に依存する分散を明示的にモデル化することで、本研究は、特定の生態学的相互作用（温度と深度の相互作用など）に対する証拠の強さが大きく変わり得ることを示している。著者らは、将来の動物プランクトン個体数の解析において、日常的な対数変換を直接的なスケールでのモデリングに置き換えるべきであり、また、生態学的な結論の堅牢性を高めるために、サンプリング設計に由来する観測レベルの不均一性を明示的に考慮すべきであると結論付けている。