リケラボ論文検索は、全国の大学リポジトリにある学位論文・教授論文を一括検索できる論文検索サービスです。

リケラボ 全国の大学リポジトリにある学位論文・教授論文を一括検索するならリケラボ論文検索大学・研究所にある論文を検索できる

リケラボ 全国の大学リポジトリにある学位論文・教授論文を一括検索するならリケラボ論文検索大学・研究所にある論文を検索できる

大学・研究所にある論文を検索できる 「Effect-Size Estimation Using Semiparametric Hierarchical Mixture Models in Disease-Association Studies with Neuroimaging Data」の論文概要。リケラボ論文検索は、全国の大学リポジトリにある学位論文・教授論文を一括検索できる論文検索サービスです。

コピーが完了しました

URLをコピーしました

論文の公開元へ論文の公開元へ
書き出し

Effect-Size Estimation Using Semiparametric Hierarchical Mixture Models in Disease-Association Studies with Neuroimaging Data

江本, 遼 名古屋大学

2021.07.26

概要

【緒言】
認知症などの精神疾患は日本でも五大疾患に指定されており、社会の急速な高齢化に伴い、精神疾患のメカニズム解明や有効な診断法の開発が喫緊の課題となっている。これらの課題の解決に資する研究の一つとして、患者と健常人の脳画像データを比較することで疾患に関連する部位のスクリーニングを行う疾患関連解析がある。この解析における統計学的な課題は、仮説検定による疾患に関連する部位の検出と、各部位の関連の大きさ(効果サイズ)の推定である。

従来の統計的手法の多くは疾患関連部位の検出に重点を置いており、近年では脳機 能の局在性を考慮するために確率場理論の確率モデルを用いるなどの提案がみられる。一方、効果サイズ推定に関しては、検出された部位の単純な推定量が多くの研究で報 告されているが、このような推定量には効果サイズを実際より大きく推定してしまう という推定バイアスの問題が指摘されている。しかし、この問題に対しては、検証用 の集団を別途用意するなどの簡易的な対処法が提案されるに留まっている。そこで、 本研究では、脳画像データの特徴である空間的な従属構造を考慮した、より推定バイ アスの小さな効果サイズの推定法を提案する。

【方法】
脳画像の構成単位である各ボクセルと疾患との関連の大きさ、すなわち、効果サイズとして、ボクセル画素値に関する疾患群と健常群での平均値の差を群内の標準偏差で割った標準化平均差 𝛿 を用いるのが自然である。なお、実際には、𝛿 の推定量として、データから求めた平均値と標準誤差を用いた推定量 𝑌 を用いる。

統計量 𝑌 の統計学的モデルを考える際には、隣接するボクセル間には強い相関が存在するという従属構造、及び、効果サイズ分布の形状に関する情報が欠如していることをモデルに反映させる必要がある。本研究では、ボクセル間の従属構造をモデルに組み込む隠れマルコフ確率場モデルと、効果サイズの分布の柔軟なモデリングを可能にするセミパラメトリック階層混合モデルを統合する。

具体的には、任意のボクセル 𝑠 に対して、これが疾患に関連しているかを表す潜在変数 Θ𝑠( Θ𝑠 = 1 ならばボクセル 𝑠 は疾患と関連あり、Θ𝑠 = 0 ならば関連なし)を導入し、全てのボクセルに関しての Θ𝑠 の同時分布がマフコフ確率場モデルの一つであるイジングモデルに従うと仮定する。サンプルサイズが十分大きい場合、漸近正規性により、疾患と関連しないボクセルの推定量 𝑌𝑠 は平均 δ𝑠 = 0 、分散 𝑐2 の正規分布に従う。すなわち、𝑌𝑠 ∼ 𝑁(0, 𝑐𝑛) である。ここで 𝑐𝑛 はサンプルサイズによって定まる定数である。一方、疾患と関連するボクセルでは、特定の分布形を指定しないノンパラメトリック分布 𝑔 を用いて、𝑌𝑠 の分布は以下のセミパラメトリック階層混合モデルで表されると仮定する。

なお、以上の分布の仮定には正規近似(漸近正規性)を用いているが、サンプルサイズが小さい場合は、正規近似を用いない 𝑡 分布を用いたモデルも考慮する。

以上のモデルのもとでの疾患との関連の有無の仮説検定 は local index of significance(LIS)統計量に基づく偽発見率(false discovery rate)の調整により行う。また、ノンパラメトリック分布 𝑔 を離散分布で近似する手法を用いて効果サイズの経験ベイズ推定を考えると、効果サイズの推定量は事後期待値 𝑑𝑠 = 𝐸[𝛿𝑠|𝑦𝑠, Θ𝑠 = 1; 𝝋̂] と事後確率 𝓁𝑠 = Pr( Θ𝑠 = 1|𝒚; 𝝋̂) の積の形で 𝛿̂ = 𝑑𝑠𝓁𝑠 として導かれる。ここで、𝝋 はモデルに含まれるパラメータであり、推定量 𝝋̂ は一般化 EM アルゴリズムを用いて求めることができる。

【結果】
提案する効果サイズ推定量 𝛿̂の真の値 𝛿𝑠からの乖離を、様々な設定のもと数値実験により確認した。サンプルサイズ及びイジングモデルのパラメータの値によらず、提案する推定量の推定バイアスは、従来の推定法よりも小さかった(Figure 1)。また、特にサンプルサイズが小さい場合に 𝑡 分布を用いる提案法の推定量は、正規近似を用いる推定量よりバイアスが小さかった。

実データへの適用として、Open Access Series of Imaging Studies (OASIS)における OASIS-2 という公開データ(Marcus et al., 2010)を用いて、経度認知障害(MCI)からアルツハイマー病(AD)への進行と脳の部位の関連解析を提案手法により行った。画像データの前処理として OASIS-2 で公開されている T1 強調画像をソフトウェア SPM により標準化、正規化し、各部位の容積を表す画像へ変換した後、提案手法を適用した(Figure 2)。Figure 2 に示すように、従来の仮説検定に重点を置いた方法と異なり、関連の有無だけでなく効果サイズの大きさを得ることができる。また、推定された効果サイズの大きい部位の多くは、従来の研究で関連が指摘されている領域であり、従来の研究とも一貫した結果であった(Table 1)。

【考察】
提案法により、多重検定を通じた疾患関連ボクセルの検出だけでなく、個々のボクセルのより正確な効果サイズの推定が可能となる。Figure 2(b)に示すように、個々のボクセルの効果サイズの推定により、検出された脳領域と疾患との関連に関する詳細で偏りのない情報が提供でき、検出された領域の生物学的または臨床的検討に役立つことが期待される。また、検定統計量ではなく効果サイズの統計モデルを用いたことにより、サンプルサイズに依存しない評価が可能となっており、サンプルサイズの異なる研究間での効果サイズの比較に役立つことも期待される。さらに、効果サイズは患者背景などの共変量が存在する場合にも拡張可能であり、提案する推定の枠組みは、線形モデルを用いる fMRI や DTI データなど多くの脳画像を用いた関連解析に適用可能である。

提案した推定法では、正規近似を用いたモデルと 𝑡 分布を用いたモデルそれぞれに対して効果サイズの推定量を得ることができる。数値実験による両者の推定性能の比較では、サンプルサイズが小さい場合、正規近似を用いた推定量は𝑡 分布を用いた推定量に劣るが、サンプルサイズが増えるにつれて同程度になる傾向が認められた。正規近似の利用は計算時間短縮の観点から有用であるため、サンプルサイズが十分に大きい場合(例えば 100 例以上)は正規近似を用いた方法を使用することを推奨する。

最後に、本研究の重要な側面として、脳画像データに含まれる多数のボクセル全体を通しての効果サイズの分布を柔軟にモデリングしていることがあげられる。これにより、ある脳領域もしくは隣接するボクセルといったごく小さな領域に対しても正確な推論を行うことが可能となる。また、提案したモデルの拡張により脳の空間的構造に加えて機能的構造を考慮するといった発展も可能である。

【結語】
本研究は脳画像を用いた疾患関連解析において、効果サイズの過大推定バイアスに対する有効な補正法の開発に取り組んだ最初の試みである。提案手法は、脳画像の従属構造を組み込んだ柔軟な階層セミパラメトリックモデルに基づいて、各ボクセルの効果サイズの正確な推定、および疾患に関連するボクセルの検出を可能にする。

参考文献

[1] R. A. Poldrack, J. A. Mumford, and T. E. Nichols, Handbook of functional MRI data analysis. Cambridge University Press, 2011.

[2] S. M. Smith and T. E. Nichols, “Threshold-free cluster enhancement: Addressing problems of smoothing, threshold dependence and localisation in cluster inference,” NeuroImage, vol. 44, no. 1, pp. 83 – 98, 2009.

[3] M. Smith and L. Fahrmeir, “Spatial Bayesian variable selection with application to functional magnetic resonance imaging,” Journal of the American Statistical Association, vol. 102, no. 478, pp. 417–431, 2007.

[4] H. Shu, B. Nan, and R. Koeppe, “Multiple testing for neuroimaging via hidden Markov random field,” Biometrics, vol. 71, no. 3, pp. 741–750, 2015.

[5] W. Sun and T. T. Cai, “Large-scale multiple testing under dependence,” Journal of the Royal Statistical Society: Series B (Statistical Methodology), vol. 71, no. 2, pp. 393–424, 2009.

[6] D. A. Brown, N. A. Lazar, G. S. Datta, W. Jang, and J. E. McDowell, “Incorporating spatial dependence into Bayesian multiple testing of statistical parametric maps in functional neuroimaging,” NeuroImage, vol. 84, pp. 97–112, 2014.

[7] M. C. Reddan, M. A. Lindquist, and T. D. Wager, “Effect size estimation in neuroimaging,” JAMA Psychiatry, vol. 74, no. 3, pp. 207–208, 2017.

[8] M. A. Lindquist and A. Mejia, “Zen and the art of multiple comparisons,” Psychosomatic Medicine, vol. 77, no. 2, pp. 114–125, 2015.

[9] B. Efron, “Empirical Bayes estimates for large-scale prediction problems,” Journal of the American Statistical Association, vol. 104, no. 487, pp. 1015–1028, 2009.

[10] S. Matsui and H. Noma, “Estimating effect sizes of differentially expressed genes for power and sample-size assessments in microarray experiments,” Biometrics, vol. 67, no. 4, pp. 1225–1235, 2011.

[11] B. Efron, “Microarrays, empirical bayes and the two-groups model,” Statistical Science, vol. 23, no. 1, pp. 1–22, 2008.

[12] L. D. Brown, “Admissible estimators, recurrent diffusions, and insoluble boundary value problems,” The Annals of Mathematical Statistics, vol. 42, no. 3, pp. 855–903, 1971.

[13] C. M. Stein, “Estimation of the mean of a multivariate normal distribution,” The annals of Statistics, vol. 9, no. 6, pp. 1135–1151, 1981.

[14] W. Shen and T. A. Louis, “Empirical Bayes estimation via the smoothing by roughening approach,” Journal of Computational and Graphical Statistics, vol. 8, no. 4, pp. 800–823, 1999.

[15] S. Matsui and H. Noma, “Estimation and selection in high-dimensional genomic studies for developing molecular diagnostics,” Biostatistics, vol. 12, no. 2, pp. 223– 233, 2011.

[16] P. Mu ̈ller and R. Mitra, “Bayesian nonparametric inference–why and how,” Bayesian analysis, vol. 8, no. 2, 2013.

[17] B. Efron, Large-Scale Inference: Empirical Bayes Methods for Estimation, Testing, and Prediction. Institute of Mathematical Statistics Monographs, Leiden: Cambridge University Press, 2010.

[18] N. Tzourio-Mazoyer, B. Landeau, D. Papathanassiou, F. Crivello, O. Etard, N. Delcroix, B. Mazoyer, and M. Joliot, “Automated Anatomical Labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain,” NeuroImage, vol. 15, no. 1, pp. 273 – 289, 2002.

[19] Z. Wei, W. Sun, K. Wang, and H. Hakonarson, “Multiple testing in genome-wide association studies via hidden Markov models,” Bioinformatics, vol. 25, no. 21, pp. 2802–2808, 2009.

[20] D. S. Marcus, A. F. Fotenos, J. G. Csernansky, J. C. Morris, and R. L. Buckner, “Open Access Series of Imaging Studies: Longitudinal MRI data in nondemented and demented older adults,” Journal of Cognitive Neuroscience, vol. 22, no. 12, pp. 2677–2684, 2010.

[21] J. Ashburner, “A fast diffeomorphic image registration algorithm,” NeuroImage, vol. 38, no. 1, pp. 95 – 113, 2007.

[22] R. T. Shinohara, E. M. Sweeney, J. Goldsmith, N. Shiee, F. J. Mateen, P. A. Calabresi, S. Jarso, D. L. Pham, D. S. Reich, C. M. Crainiceanu, et al., “Statistical normalization techniques for magnetic resonance imaging,” NeuroImage: Clinical, vol. 6, pp. 9–19, 2014.

[23] Y. Deng, K. Liu, L. Shi, Y. Lei, P. Liang, K. Li, W. C. W. Chu, and D. Wang for the Alzheimer’s Disease Neuroimaging Initiative, “Identifying the alteration patterns of brain functional connectivity in progressive mild cognitive impairment patients: A longitudinal whole-brain voxel-wise degree analysis,” Frontiers in Aging Neuroscience, vol. 8, p. 195, 2016.

[24] P. J. Modrego, N. Fayed, and M. Sarasa, “Magnetic resonance spectroscopy in the prediction of early conversion from amnestic mild cognitive impairment to dementia: a prospective cohort study,” BMJ open, vol. 1, no. 1, 2011.

[25] E. Pravatà, J. Tavernier, R. Parker, H. Vavro, J. E. Mintzer, and M. V. Spampinato, “The neural correlates of anomia in the conversion from mild cognitive impairment to Alzheimer’s disease,” Neuroradiology, vol. 58, no. 1, pp. 59–67, 2016.

[26] S. F. Eskildsen, P. Coupé, V. S. Fonov, J. C. Pruessner, and D. L. Collins, “Structural imaging biomarkers of Alzheimer’s disease: predicting disease progression,” Neurobiology of Aging, vol. 36, pp. S23–S31, 2015.

[27] Y. Liu, T. Paajanen, Y. Zhang, E. Westman, L.-O. Wahlund, A. Simmons, C. Tunnard, T. Sobow, P. Mecocci, M. Tsolaki, B. Vellas, S. Muehlboeck, A. Evans, C. Spenger, S. Lovestone, and H. Soininen, “Analysis of regional MRI volumes and thicknesses as predictors of conversion from mild cognitive impairment to Alzheimer’s disease,” Neurobiology of Aging, vol. 31, no. 8, pp. 1375–1385, 2010.

[28] X. Tang, D. Holland, A. M. Dale, L. Younes, M. I. Miller, and A. D. N. Initiative, “Shape abnormalities of subcortical and ventricular structures in mild cognitive impairment and Alzheimer’s disease: detecting, quantifying, and predicting,” Human Brain Mapping, vol. 35, no. 8, pp. 3701–3725, 2014.

[29] H.-A. Yi, C. Möller, N. Dieleman, F. H. Bouwman, F. Barkhof, P. Scheltens, W. M. van der Flier, and H. Vrenken, “Relation between subcortical grey matter atrophy and conversion from mild cognitive impairment to Alzheimer’s disease,” Journal of Neurology, Neurosurgery & Psychiatry, vol. 87, no. 4, pp. 425–432, 2016.

[30] G. Karas, J. Sluimer, R. Goekoop, W. van der Flier, S. Rombouts, H. Vrenken, P. Scheltens, N. Fox, and F. Barkhof, “Amnestic mild cognitive impairment: Structural MR imaging findings predictive of conversion to Alzheimer disease,” American Journal of Neuroradiology, vol. 29, no. 5, pp. 944–949, 2008.

[31] T. Kato, Y. Inui, A. Nakamura, and K. Ito, “Brain fluorodeoxyglucose (FDG) PET in dementia,” Ageing research reviews, vol. 30, pp. 73–84, 2016.

[32] C. Salvatore, A. Cerasa, and I. Castiglioni, “MRI characterizes the progressive course of AD and predicts conversion to Alzheimer’s dementia 24 months before probable diagnosis,” Frontiers in Aging Neuroscience, vol. 10, p. 135, 2018.

[33] J. E. Desmond and G. H. Glover, “Estimating sample size in functional MRI (fMRI) neuroimaging studies: statistical power analyses,” Journal of Neuroscience Methods, vol. 118, no. 2, pp. 115 – 128, 2002.

[34] B. Efron, “Large-scale simultaneous hypothesis testing: the choice of a null hypoth- esis,” Journal of the American Statistical Association, vol. 99, no. 465, pp. 96–104, 2004.

参考文献をもっと見る

全国の大学の
卒論・修論・学位論文

一発検索!

この論文の関連論文を見る