[論文レビュー] A New Test for One-Way ANOVA with Functional Data and Application to Ischemic Heart Screening

Jin‐Ting Zhang, Ming−Yen Cheng+2 more|arXiv (Cornell University)|Sep 27, 2013

Statistical Methods and Inference参考文献 21被引用数 5

ひとこと要約

本稿では、観測区間全体で点별F統計量の最大値をとる新しいグローバル検定法である$F_{\text{max}}$-testを提案する。ノンパラメトリックブートストラップを用いて帰無分布を近似することで、従来の手法と比較して、特に中程度から高相関の下で良好なサイズ制御と検出力が得られることを示した。心電図（ECG）信号の安静状態での測定に応用することで、ストレス試験を要せず虚血性心筋症を安全かつ低コストでスクリーニング可能であることが示された。

ABSTRACT

We propose and study a new global test, namely the $F_{\max}$-test, for the one-way ANOVA problem in functional data analysis. The test statistic is taken as the maximum value of the usual pointwise $F$-test statistics over the interval the functional responses are observed. A nonparametric bootstrap method is employed to approximate the null distribution of the test statistic and to obtain an estimated critical value for the test. The asymptotic random expression of the test statistic is derived and the asymptotic power is studied. In particular, under mild conditions, the $F_{\max}$-test asymptotically has the correct level and is root-$n$ consistent in detecting local alternatives. Via some simulation studies, it is found that in terms of both level accuracy and power, the $F_{\max}$-test outperforms the Globalized Pointwise F (GPF) test of \cite{Zhang_Liang:2013} when the functional data are highly or moderately correlated, and its performance is comparable with the latter otherwise. An application to an ischemic heart real dataset suggests that, after proper manipulation, resting electrocardiogram (ECG) signals can be used as an effective tool in clinical ischemic heart screening, without the need of further stress tests as in the current standard procedure.

研究の動機と目的

• 一般の自己相関構造下でも正しいサイズと高い検出力を維持できる、関数データに対する頑健なグローバル検定法の開発を目的とする。
• リスクが高く高コストなストレス試験を要しない状態性虚血性心筋症の検出という臨床的課題に応えること。
• 適切に処理された安静状態の心電図信号が、虚血性心疾患の信頼できる非侵襲的スクリーニングツールとして機能しうるかを検討すること。
• $F_{\text{max}}$-検定とグローバライズド・ポイントワイズF（GPF）検定の間で、有意水準の正確性と統計的検出力の面で性能を比較すること。
• 実心電図データを用いた検証を通じて、外来スクリーニングにおける臨床的応用の可能性を示すこと。

提案手法

• $F_{\text{max}}$-検定は、関数的領域$\mathcal{T} = [a,b]$全域における点別F検定統計量の最大値を検定統計量として構築する。
• 帰無分布と臨界値の推定にノンパラメトリックブートストラップ法を用いることで、弱い自己相関仮定下でも正確なサイズ制御が保証される。
• 検定統計量の漸近的分布を導出し、局所代替仮説下での漸近的検出力を検討した。これにより、ルート-nの一貫性が示された。
• 点別F検定を活用することで、一元分散分析を超えた多様な関数データ設定への応用が可能となり、広範な適用可能性が得られる。
• 呼吸アーチファクトの除去と虚血関連信号の強化を目的として、心電図データは適応的リード、パワースペクトル、累積パワースペクトルを用いて事前処理された。
• 実データ応用における有意水準の評価のため、10,000回のブートストラップ再サンプリングを用いてp値を算出している。

実験結果

リサーチクエスチョン

• RQ1 点別F統計量の最大値に基づくグローバル検定は、関数的一元分散分析において、従来手法に比べてより優れたサイズ制御と検出力を達成できるか？
• RQ2 $F_{\text{max}}$-検定は、ややきつい正則性条件の下で、正しいサイズを維持し、局所代替仮説を検出する際にルート-nの一貫性を示すか？
• RQ3 適応的リード信号をパワースペクトルまたは累積パワースペクトルに変換することで、正常群と虚血群の間に検出可能な差が現れるか？
• RQ4 $F_{\text{max}}$-検定は、心電図信号に及ぼされる呼吸の影響といった生理的交絡要因に対して頑健か？
• RQ5 $F_{\text{max}}$-検定は、状態性虚血性心筋症のスクリーニングにおいて、ストレス試験の代替として安全かつ非侵襲的であるか？

主な発見

• 関数的データに中程度から高相関がある場合、$F_{\text{max}}$-検定はGPF検定に比べてサイズ制御と検出力の両面で優れている。
• ややきつい正則性条件の下で、$F_{\text{max}}$-検定は漸近的有意水準の正確性と局所代替仮説の検出におけるルート-nの一貫性を達成している。
• 適応的リード波形に対して、$F_{\text{max}}$-検定はp値0.0006を示し、正常群と虚血群の平均関数の差に強い証拠があることを示している。
• 非正規化された適応的リード信号に適用した場合、$F_{\text{max}}$-検定は有意でない（p = 0.0764）ため、呼吸アーチファクトを除去するための事前処理の重要性が浮き彫りになった。
• 累積パワースペクトルではp値が0.0029となり、低周波数帯域エネルギーの群間シフトが強化されたことにより検出力が向上しており、$F_{\text{max}}$-検定がこれを効果的に捉えていることが確認された。
• 結果から、安静状態の心電図信号のスペクトル特徴に、臨床的に意味のある虚血性心疾患スクリーニング情報が含まれており、ストレス試験の必要性を排除できる可能性があるという仮説が支持された。