一般化線形混合モデルのソースを表示

{{otheruses|一般化線形混合モデル|一般線形混合モデル|混合モデル}}

'''一般化線形混合モデル'''（いっぱんかせんけいこんごうモデル、{{lang-en-short|Generalized linear mixed model, GLMM}}）とは、統計学において[[一般化線形モデル]]を拡張した統計解析モデルである。さらにこの一般化線形混合モデルを拡張し、事前分布に含まれる母数の事前分布を導入する場合には、[[階層ベイズモデル]]とみなされる。

* 固定効果に加えて変量効果を考慮している。変量効果は通常、[[正規分布]]を仮定される。変量効果の積分を含めて最尤法に基づき解を求めるが、一般に解析形式で式を表現することができない。そのため、[[数値積分]]や[[マルコフ連鎖モンテカルロ法]]等を用いて計算機により[[数値解析]]的に解を求める。コンピュータの高速化により実用化されてきている。

* [[シグモイド関数|逆ロジット関数]]のようなリンク関数のほかにも、多くのリンク関数を用いて、[[ベルヌーイ分布|ベルヌイ分布]]以外の確率分布に回帰式をリンクできることが利点である。

==モデル==
一般化線形混合モデルは、一般に、ランダム効果 <math>u</math> を条件とする従属変数 <math>y</math> が[[指数型分布族]]に従って分布し、その期待値がリンク関数 <math>g</math> を介して線形予測子 <math>X \beta + Zu</math> に関連する、として定義される。

:<math>g\left(E\left[y \mid u\right]\right) = X \beta + Zu</math>

ここで、<math>X</math> は固定効果 <math>\beta</math> のデザイン行列、<math>Z</math> はランダム効果 <math>u</math> のデザイン行列である。
この簡潔な定義を理解するためには、[[一般化線形モデル]]と[[混合モデル]]についての説明を参照のこと。
一般化線形混合モデルは、[[階層型一般化線形モデル]]のうち、ランダム効果の正規分布を仮定したものと考えることもできる。

完全尤度 <math>\ln{p(y)}</math> は次式で表され、多くの場合、膨大な計算量が必要になる<ref>{{cite book |last1=Pawitan |first1=Yudi |title=In All Likelihood: Statistical Modelling and Inference Using Likelihood |publisher=OUP Oxford |isbn=978-0199671229 |pages=459 |edition=Paperbackition}}</ref>。

:<math>\ln{p(y)} = \ln{\int p(y \mid u)\,p(u)\,du}</math>

Gauss-Hermite 求積法などによる数値積分のほか、ラプラス近似を用いた方法も提案されている
<ref>{{cite journal |last1=Breslow |first1=N. E. |last2=Clayton |first2=D. G. |title=Approximate Inference in Generalized Linear Mixed Models |journal=Journal of the American Statistical Association |date=20 December 2012 |volume=88 |issue=421 |pages=9–25 |doi=10.1080/01621459.1993.10594284}}</ref>。
罰則つき疑似尤度法（penalized quasi-likelihood method）は、重み付き正規混合モデルを作業変数で繰り返しフィッティングするものであり、多くの商用ないしオープンソースの統計プログラムで実装されている<ref>{{cite journal |last1=Wolfinger |first1=Russ |last2=O'connell |first2=Michael |title=Generalized linear mixed models a pseudo-likelihood approach |journal=Journal of Statistical Computation and Simulation |date=December 1993 |volume=48 |issue=3-4 |pages=233–243 |doi=10.1080/00949659308811554}}</ref>。

==モデルのフィッティング==
最尤法による一般化線形混合モデルのフィッティングにはランダム効果の積分が必要だが、一般的に、解析的な形では表現できない。様々な近似法が開発されたが、どんなモデルとデータセットにも上手く適用できるような方法はない。このため、計算能力の向上と手法の進歩により、[[数値積分]]や[[マルコフ連鎖モンテカルロ法]]が広く用いられるようになった。

モデル選択の基準としては、[[赤池情報量規準]]（AIC）が一般的である。最近、ある[[指数型分布族]]の分布に基づく、一般化線形混合モデルの AIC の推定値が得られた<ref name="Saefken">{{Citation | last1= Saefken | first1= B. | last2= Kneib | first2= T. | last3= van Waveren | first3= C.-S. | last4= Greven | first4= S. | year= 2014 | title= A unifying approach to the estimation of the conditional Akaike information in generalized linear mixed models | journal= [[Electronic Journal of Statistics]] | volume= 8 | pages= 201–225 | doi= 10.1214/14-EJS881| url= http://goedoc.uni-goettingen.de/goescholar/bitstream/handle/1/10640/euclid.ejs.1393510264.pdf?sequence=1 | doi-access= free }}</ref>

==ソフトウェア==
* [[R言語]]<ref>{{Citation | first1= J. C. | last1= Pinheiro| first2= D. M. | last2= Bates| title= Mixed-effects models in S and S-PLUS | year= 2000 | publisher= [[Springer Publishing|Springer, New York]]}}</ref><ref>{{Citation | first1= D. M. | last1= Berridge | first2= R. | last2= Crouchley | title= Multivariate Generalized Linear Mixed Models Using R | year= 2011 | publisher= [[CRC Press]]}}</ref>
* [[SAS]] および [[SPSS]]<ref>{{cite web|title=IBM Knowledge Center|url=https://www.ibm.com/support/knowledgecenter/en/SSLVMB_22.0.0/com.ibm.spss.statistics.help/spss/advanced/idh_idd_genlin_typeofmodel.htm|website=www.ibm.com|accessdate=6 December 2017}}</ref>
* [[MATLAB]] – 関数 <code>fitglme</code> で一般化線形混合モデルをフィッティング
* [[Python]] – Statsmodels パッケージ<ref>{{cite web|title=Statsmodels Documentation|url=https://www.statsmodels.org/stable/mixed_glm.html|website=www.statsmodels.org|accessdate=17 March 2021}}</ref>

==関連項目==
* [[一般化線形モデル]]
* [[混合モデル]]
* [[階層ベイズモデル]]

== 脚注 ==
{{脚注ヘルプ}}
=== 出典 ===
{{Reflist}}

==参考文献==
* {{cite journal| last=Breslow|first=N.E. |last2=Clayton |first2=D.G. |authorlink2=David Clayton |title=[[Approximate Inference]] in Generalized Linear Mixed Models|journal=[[Journal of the American Statistical Association]]| volume=88| issue=421| pages=9–25| year=1993| doi=10.2307/2290687| jstor=2290687 }}

* {{cite book|title=Applied longitudinal analysis|first= Garrett M. | last=Fitzmaurice |coauthors= Laird, Nan M.; Ware, James H. |publisher= Wiley-Interscience|location= Hoboken, NJ|year= 2004| isbn= 0-471-21487-6}}

{{統計学}}
{{Statistics-stub}}

{{デフォルトソート:いつはんかせんけいこんこうもてる}}
[[Category:統計モデル]]
[[Category:分散分析]]
[[Category:数学に関する記事]]