MENU

SPSS で線形混合モデルを実行する方法

線形混合モデルをSPSSで実行する方法の解説。

反復測定データで、欠損値があってもできる限りデータを使用できる点で、反復測定分散分析よりも適切な方法だ。

>>もう統計で悩むのは終わりにしませんか? 

↑期間・数量限定で無料プレゼント中!

目次

線形混合モデルをSPSSで実行するために使用するサンプルデータ

EZR公式マニュアル付属のデータ FCZ_CSA.rda を用いる。

EZRでやさしく学ぶ統計学 改訂3版 〜EBMの実践から臨床研究まで〜

3群で3回の反復測定連続データである。

平均値と95%信頼区間のグラフを描くとこんな感じ。

平均値と95%信頼区間の推移のグラフ

推定周辺平均が計算されているので、単純な平均値ではないが、いわゆる平均値の折れ線グラフである。

線形混合モデルの前に反復測定分散分析を行ってみる

SPSSで反復測定分散分析を行うには、「分析」→「一般線型モデル」→「反復測定」を選択。

まず、反復測定の因子の定義を行う必要がある。

反復測定の因子の定義の設定方法は、過去記事参照。

解析結果まで、省略する。

被験者内効果の結果は以下の通りである。

球面性の仮定は置けない結果になっているので、Greenhouse-Geisserの行の結果を見ることになる。

被験者間効果の結果は、以下の通り。

>>もう統計で悩むのは終わりにしませんか? 

↑1万人以上の医療従事者が購読中

SPSS で線形混合モデルを実施してみる

線形混合モデルは、被験者のばらつきを変量効果として考慮したモデルで、少しの欠測値があってもそのまま解析して問題ないところが利点。

サンプルデータは実際ちらほら欠測値がある。

赤矢印部分が欠測値。

まず、データを一人1行ではなく、全目的変数を一列に並べる作業を行う。

反復測定分散分析は、一人1行のデータにしていたが、それを縦に積み重ねる。

これはEZRで実施したほうが簡単なので、SPSSに移植してくる前に縦に積み重ねてから、SPSSに持ってきた。

CDratioという目的変数が縦方向に一列に並んでいる。

それが、被験者と時点が他の変数で区別されるようになっている。

縦持ちに変更したデータ

ここまで準備が終わったら、いよいよ解析する。

「分析」→「混合モデル」→「線型」と選択すると、線形混合モデルが選択できる。

ちなみに、線型と線形は同じものを指しているので、漢字の違いは気にしない。

被験者に id(被験者ID)と FCZ(群別変数)、反復測定に反復測定を区別する変数 Weeks(カテゴリ変数の週数)を投入する。

反復測定共分散は、複合シンメトリ(Compound Symmetry)を選択する。

反復測定データの誤差の分散共分散が、どのような構造になっているかを指定する。

複合シンメトリはそろっていること(等分散・等共分散)を想定している。

あとで、そろっていない「無構造」もやってみる。

続行をクリックした後の窓で、従属変数にCDratio、因子にFCZとWeeksを投入する。

「固定」ボタンをクリックして、モデルを追加する。

左のFCZとWeeksを同時に選択して、追加をクリックすると、以下のようになる。

続行で戻る。

「統計量」ボタンをクリックして、固定効果のパラメータ推定値にチェックを入れる。

続行で戻る。

OKをクリックすると結果が出力される。

線形混合モデルを実施した後の結果の読み方

固定効果のANOVA の結果は以下の通り。

反復測定分散分析とは少し異なるが、傾向は似ている。

固定効果の推定値は、以下の通り。

反復測定分散分析はこういう結果は出せないので、回帰分析グループの線形混合モデルならではの出力である。

欠測値の影響を受けにくく、欠測含め171例中、反復測定分散分析では163例が計算に使われたが、線形混合モデルでは170例が計算に使われていた。

48+87+28=163
被験者数 170

反復測定の時間の項を連続量で投入した場合

反復測定の時間を連続量で投入して、直線的に変化すると仮定するとどうなるか?

反復測定にWeeksNum(数値データ)を投入する。

因子のWeeksを戻し、共変量にWeeksNumを入れる

固定効果は、FCZとWeeksNumの主効果と交互作用をともに設定する。

固定効果の推定の結果は以下の通り。

反復測定共分散を無構造にした場合

無構造というのは、特に想定なしということである。

固定効果の推定結果は、以下の通り。

複合シンメトリと少し異なるがほぼ同じ結果だ。

共分散パラメータの推定値は以下のようになる。

まとめ

線形混合モデルをSPSSで実施する方法を解説した。

欠測値が少々あっても計算に使ってくれるという点でよいモデルだ。

誤差の分散共分散構造は、複合シンメトリにしておけばよいだろう。

何の前提もなしとするなら、無構造がよいかもしれない。

反復測定データは、ぜひ線形混合モデルで解析することをおすすめしたい。

おすすめ書籍

EZRでやさしく学ぶ統計学 改訂3版 〜EBMの実践から臨床研究まで〜
SPSSによる分散分析・混合モデル・多重比較の手順
よかったらシェアしてね!
  • URLをコピーしました!
  • URLをコピーしました!

リサーチクエスチョン探し?データ分析?論文投稿?、、、で、もう悩まない!

第1章臨床研究ではなぜ統計が必要なのか?計画することの重要性
  • 推定ってどんなことをしているの?
  • 臨床研究を計画するってどういうこと?
  • どうにかして標本平均を母平均に近づけられないか?
第2章:研究目的をどれだけ明確にできるのかが重要
  • データさえあれば解析でどうにかなる、という考え方は間違い
  • 何を明らかにしたいのか? という研究目的が重要
  • 研究目的は4種類に分けられる
  • 統計専門家に相談する上でも研究目的とPICOを明確化しておく
第3章:p値で結果が左右される時代は終わりました
  • アメリカ統計協会(ASA)のp値に関する声明で指摘されていること
  • そうは言っても、本当に有意差がなくてもいいの…?
  • なぜ統計専門家はp値を重要視していないのか
  • 有意差がない時に「有意な傾向があった」といってもいい?
  • 統計を放置してしまうと非常にまずい
第4章:多くの人が統計を苦手にする理由
  • 残念ながら、セミナー受講だけで統計は使えません。
  • インプットだけで統計が使えない理由
  • どうやったら統計の判断力が鍛えられるか?
  • 統計は手段なので正解がないため、最適解を判断する力が必要
第5章:統計を使えるようになるために今日から何をすれば良いか?
  • 論文を読んで統計が使えるようになるための5ステップ
第6章:統計を学ぶために重要な環境
  • 統計の3つの力をバランスよく構築する環境

以下のボタンをクリックして、画面に出てくる指示に従って、必要事項を記入してください。

この記事を書いた人

統計 ER ブログ執筆者

元疫学研究者

統計解析が趣味

コメント

コメントする

目次