2009/10/26

次回データを用いてMixed Model

前回(2009/10/26)のエビ研は、Mixed Model初級編(2)ということで、初級編(1)の復習と、SPSS(今はPASWですね。)を用いてデータの再構成(SPSS用語か?データのreshaping?)を行いました。

データの再構成自体は、mixed modelと直接の関係はありませんが、3時点以上のデータがある場合にはmixed modelを用いると良いといった話題から、複数の時点で取っているデータをどのように入力したらmixed modelで解析できるのかという現実的な話が出たこともあり、実際にSPSS上で行ってみました。

さらに、土屋政雄さんから、縦断解析を行う際にお薦めの本をご紹介いただきました。Singer & Willett の Applied Longitudinal Data Analysis: Modeling Change and Event Occurrenceです。すばらしい本とのこと。(読んでみたいような、誰かにじっくり読んでもらってそれをあとで教えてもらいたいような。。。)

この日は反復測定など、複数時点でデータ収集を行うような研究について主に考えましたが、複数の地域や施設から横断的にデータを収集するような研究でmixed modelを用いて解析する予定の方もいらっしゃいます。
ということで、次回は縦断、横断両方のデータでmixed modelに触れてみよう!ということになりました。サンプルデータを用いてSPSS(もしくはSASも?)で行う予定です。

Mixed Modelを知りたい人、使用予定の人、どのような方もどうぞいらしてください。
日 時:2009年11月19日(木)11:00~
テーマ:Mixed Model 初級編(3)
場 所:医学部3号館 2階S210(計算機室)

SPSS用のサンプルデータは土屋さんが探してくれるとのことです。ありがとうございます!
ご自分のデータがある方は、それを使っていただいてもかまいません。
また、Mixed Modelを用いた研究で解析方法について丁寧に記述してある論文をお持ちの方は、是非ご紹介ください。

以下は、今回行ったデータの再構成についてです。


1サンプルに関して測定したデータは1行で横長(wide)に入力することが多いと思いますが、SPSSのデータの再構成(Restructure)を行うと、たとえば反復測定を行った場合、同じ人のtime1, time2, time3をそれぞれ1行にして縦長(long)のデータセットにすることができます。その逆もできます。

参加者の多くは(私も含めて)、今までこういったlong形式のデータ入力を行う場合にはエクセル上で加工していたので、SPSSにこんな機能があったなんて!便利!と思った次第です。
SAS等、他のソフトでもできるとのこと。便利ですね。

実際にエビ研で行ったのはwide -> longの例です。
unit1~4に所属する人々(ID1~6)を2つのgroupに割り付け、3時点で何かの値を測定したとします。

wide
↑wide形式のデータ

これを、メニューバーのデータ→再構成から操作していくとlongの形式に変換できます。シンタックスは下記。
VARSTOCASES /MAKE y FROM time1 time2 time3
/INDEX = t(3)
/KEEP = ID unit group
/NULL = KEEP.

long
↑long形式のデータ


UCLAのサイト↓がわかりやすいです!
http://www.ats.ucla.edu/stat/spss/modules/reshapel115.htm

ご関心のある方は次回も是非いらしてください! (文責:宮本有紀)

2009/10/01

次もMixed Model

前回(2009/10/1)のエビ研はMixed Modelを取り上げました。

Mixed Modelとは、そもそも何のMixという意味なのか、いったい何をしているのか、どんなときに使うのか、どんな解析ソフトを用いれば良いのか、といった素朴な疑問から入っていったこの回、参加者全員で確認(納得?)できたことは、
・Mixed ModelはRandom EffectとFixed EffectをMixして解析するもの
・multilevel model, hierarchical linear model, random-effects model等いろいろな呼び方あり
ということです。
まぁ、Mixed Modelを知るのにrandom effect(変量効果)とfixed effect(固定効果)は外せないということですね、きっと。

そしてこのたび、皆でMixed Modelを学ぶにあたり、土屋政雄さんがとっても素敵な参考文献リストをエビ研にご提供くださいました!土屋さん(つっちー)の許可を得て、末尾に掲載させていただきます。

この回では、つっちーのリストにある文献を基に、皆の知識を持ち寄りながら、数学的な話から、MLwiN(multilevel model解析ソフト)の値段といった実用的な話まで、幅広い話をできたように思います。

また、自分のデータのどれをrandom effectと考えたらよいのか?という疑問に対しては、米倉佑貴さんが、Crawleyの書籍STATISTICS: AN INTRODUCTION USING R と、久保拓弥さんの「生態学のデータ解析 - ランダム効果とは?」というページをご紹介くださいました。
こういった知識・情報を共有してくれる皆様がいらっしゃるということは有り難いことですね。

上記のように盛りだくさんだったMixed Model初級編、次回も続けて行うことにしました。
Mixed Modelを知りたい人、使用予定の人などなど、ご関心のある方は、できれば、ご自分の関心領域でMixed Modelを用いている文献を1つご持参いただくと良いかと思います。
日 時:2009年10月26日(月)10:30~
テーマ:Mixed Model 初級編(2)
場 所:医学部3号館 3階S308



Mixed model おすすめ文献集 (by 土屋政雄さん)

<マルチレベル一般>
・Diez Roux AV. A glossary for multilevel analysis. J Epidemiol Community Health. 2002 Aug;56(8):588-94.
・土居淳子. 社会科学データ分析のための線形混合モデル入門. 京都光華女子大学研究紀要. 45:211-232, 2007
・マルチレベルモデリング. 科学技術振興機構. 尾崎幸謙 http://www010.upp.so-net.ne.jp/koken/multilevel.ppt

[JECHのシリーズ] 
Merlo J, Chaix B, Yang M, Lynch J, Råstam L. A brief conceptual tutorial of multilevel analysis in social epidemiology: linking the statistical concept of clustering to the idea of contextual phenomenon. J Epidemiol Community Health.2005 Jun;59(6):443-9. Review.
Merlo J, Yang M, Chaix B, Lynch J, Råstam L. A brief conceptual tutorial on multilevel analysis in social epidemiology: investigating contextual phenomena indifferent groups of people. J Epidemiol Community Health. 2005 Sep;59(9):729-36.
Merlo J, Chaix B, Yang M, Lynch J, Råstam L. A brief conceptual tutorial onmultilevel analysis in social epidemiology: interpreting neighbourhood differences and the effect of neighbourhood characteristics on individual health.J Epidemiol Community Health. 2005 Dec;59(12):1022-8. Review.
Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, Råstam L, Larsen K. A brief conceptual tutorial of multilevel analysis in social epidemiology:using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health. 2006 Apr;60(4):290-7.


<介入研究用>
Gueorguieva R, Krystal JH. Move over ANOVA: progress in analyzing repeated-measures data and its reflection in papers published in the Archives of General Psychiatry. Arch Gen Psychiatry. 2004 Mar;61(3):310-7.
Hamer RM, Simpson PM. Last observation carried forward versus mixed models in the analysis of psychiatric clinical trials. Am J Psychiatry. 2009Jun;166(6):639-41.
Vickers AJ. Analysis of variance is easily misapplied in the analysis of randomized trials: a critique and discussion of alternative statistical approaches. Psychosom Med. 2005 Jul-Aug;67(4):652-5.
※これの補足として
Vickers AJ, Altman DG. Statistics notes: Analysing controlled trials with baseline and follow up measurements. BMJ. 2001 Nov 10;323(7321):1123-4.


<看護系雑誌での初心者向け特集>
Wu YW, Wooldridge PJ. The impact of centering first-level predictors on individual and contextual effects in multilevel data analysis. Nurs Res. 2005 May-Jun;54(3):212-6.
Adewale AJ, Hayduk L, Estabrooks CA, Cummings GG, Midodzi WK, Derksen L. Understanding hierarchical linear models: applications in nursing research. Nurs Res. 2007 Jul-Aug;56(4 Suppl):S40-6.
Estabrooks CA, Midodzi WK, Cummings GG, Wallin L. Predicting research use in nursing organizations: a multilevel analysis. Nurs Res. 2007 Jul-Aug;56(4Suppl):S7-23.
Shin JH. Application of repeated-measures analysis of variance and hierarchical linear model in nursing research. Nurs Res. 2009 May-Jun;58(3):211-7.


<統計ソフト>
[SPSS]
・Peugh JL. Using the SPSS Mixed Procedure to Fit Cross-Sectional and Longitudinal Multilevel Models. Educational and Psychological Measurement. 65:717-741 (2005)
・Hayes AF. A Primer on Multilevel Modeling. Human Communication Research. 32:385-410 (2006).
・Painter J. Designing Multilevel Models Using SPSS. http://www.unc.edu/~painter/SPSSMixed/SPSSMixedModels.PDF


[SAS]
・Singer JD. Using SAS PROC MIXED to Fit Multilevel Models, Hierarchical Models, and Individual Growth Models. Journal of Educational and Behavioral Statistics. 23:323-355 (1998)

<個人ページ>
John B. Willett
http://gseacademic.harvard.edu/~willetjo/
Judith D. Singer
http://gseweb.harvard.edu/~faculty/singer/
Tom Snijders
http://stat.gamma.rug.nl/multilevel.htm
Don Hedeker
http://tigger.uic.edu/~hedeker/ml.html
兼清さんの「まぜなれ」とか
http://bm.hus.osaka-u.ac.jp/~kanekiyo/study/index.html


[機関のページ]
The Centre for Multilevel Modelling (CMM)
http://www.cmm.bristol.ac.uk/index.shtml

UCLA Stat Computing Portal Multilevel Modeling Portal
http://statcomp.ats.ucla.edu/mlm/default.htm


つっちーリスト以上

ご関心のある方は次回も是非いらしてください!(文責:宮本)