狠狠撸

狠狠撸Share a Scribd company logo
構造方程式モデリング基本の「き」 2010 年 6 月 14 日 @ Twittcher 勉強会
Agenda 超基本编 構造方程式モデリングって何ですか? 構造方程式モデリング基本の「き」 論よりrun!~実践編 どうやって分析するの? 分析結果を解釈してみよう 事例绍介
今日の獲得目標 構造方程式モデリングを勉強する下地を作る 一体、何ができるの? 勉強すると何がいいの? どうやって分析?解釈するの? 「ちょっといじってみるかな?」という気分になること 入門書を読む前のウォーミングアップ ※ 簡単のために数式をできるだけ使わず、若干な乱暴な説明をします。  正しい表現は入門書を読みながら身につけましょう!
超基本编
【構造方程式モデリングって何ですか?】 構造 方程式モデリング 構造方程式モデリング(Structural Equation Modeling, SEM) 共分散構造分析(Covariance Structure Analysis)ともいう でも???構造方程式モデリングのほうが良い 共分散だけでなく、平均を構造化する分析もあるため
【構造方程式モデリングって何ですか?】 どういう特徴があるの? ざくっと言えば??? 変数間の関係を統計的に検証できる 潜在変数(構成概念)を導入した分析ができる わかりやすいパス図を使ったビジュアル表現 SEM を勉強すると色々な分析ができる
【構造方程式モデリングって何ですか?】 変数の関係を統計的に検証する 変数の間の関係を記述して、データとのあてはまりを評価する、統計的な分析方法 B A C データに あてはまってる? Data
例えば??? 病気罹患と食品摂取の関係 ???? 総摂取量 乳製品 摂取量 大腸がん 罹患率 誤差 誤差 データへの あてはまり評価 Data 影響力を推定する 変数の関係を描いて 同時に行う! + + -
【構造方程式モデリングって何ですか?】 潜在変数を導入した分析ができる 複数の変数で潜在変数(構成概念)を規定することができる 構成概念って? とりあえずその存在を仮定することによって、複雑に込み入った現象を比較的単純に理解することを目的に構成する概念 「知能」、「抑うつ傾向」、「ブランドエクイティ」など?? 数理力 数学 科学 物理 例えば??? ○○ 傾向 設問 1 設問 3 設問 2
【構造方程式モデリングって何ですか?】 パス図を使ったビジュアル表現 変数間の関係をパス図をつかってビジュアル化 矢印は「もと」から「さき」へ向かっての影響を表現 長方形や楕円は変数の特徴を表現 関係性がわかりやすい 直感的に変数の関係がわかる! e e e e e e ξ
【構造方程式モデリングって何ですか?】 SEM を勉強すると色々な分析ができる 伝統的な統計モデルを下位モデルとして含んでいる 第 2 世代の多変量解析モデル 回帰分析 因子分析 分散分析 共分散分析 時系列分析 パス解析 不完全データ の解析 潜在曲線分析 全て SEM で 表現可能!
【構造方程式モデリングって何ですか?】 よく因子分析と回帰分析の組み合わせとか説明される??? 「因子のパス解析」とか 間違いではないけれど、SEMの一部を表現しているに過ぎない もっといろんなことができるんです! 今日はあんまり触れない???。
【構造方程式モデリング基本の「き」】 SEMを勉強するにあたり 以下の表現を覚えておきましょう。 観測変数と潜在変数 測定方程式と構造方程式 モデルの適合
【構造方程式モデリング基本の「き」】 観測変数と潜在変数 直接に観測される変数を「観測変数」とよぶ 身長、体重、テストの得点など??? 構成概念を「潜在変数」とよぶ 知能、ブランドエクイティ??? ? SEM はこの 2 種類の変数の関係を記述する パス図では??? 観測変数 潜在変数 観測変数は四角形で、潜在変数は円で表現する
【構造方程式モデリング基本の「き」】 構造方程式と測定方程式 2種類の方程式を解くことで、変数間の関係を明らかにする 構造方程式:因果関係を記述する方程式 測定方程式:潜在変数とそれを構成する観測変数を記述する方程式 V1 V2 e V2=α * V1+e 1 α f1 V1 V2 V3 e1 e2 e3 1 1 1 α1 α2 α3 V1=α1 * f1+e1 V2=α2 * f1+e2 V3=α3 * f1+e3 構造方程式 測定方程式 計算は機械がやってくれます???
【構造方程式モデリング基本の「き」】 モデルの適合 解釈を行うためには、構成したモデルがデータに当てはまっている必要がある データとモデルとのあてはまりを 「適合度指標」 を用いて評価する Data model 各種適合度指標であてはまりの良さを評価! (実践しながら説明します!)
【構造方程式モデリング基本の「き」】 分析の流れ 仮説を構築->モデルで表現->母数(パラメータの推定)->適合度評価->モデル修正->???? Plan Do Check Action 仮説を構築する モデル(パス図)を作ってみる 分析用ソフトウェアを使ってモデルを記述 影響力(母数)を推定する 適合度指標をみて、データの当てはまりの良さを確認する 【データ Fit 良い】 結果の解釈を行う 【データ Fit 悪い】 モデル修正の検討-> Plan へ
論より RUN !~実践編
まずは分析環境 いろいろあって一長一短 AMOS がベスト。無料ならば AMOS の student 版、もしくは R で始める AMOS R ( Package sem ) Mplus EQS SAS (CALIS) Mx ???????? グラフィカル グラフィカル コマンド入力 コマンド入力 コマンド入力 ? 分析領域 ほぼカバー 一部未対応 一部未対応 ほぼカバー ほぼカバー ? 価格 有料 (無料版あり) 無料 有料 有料 有料 ? オススメ! オススメ!
AMOSのStudent版を利用 ここからダウンロード可能 http://www.amosdevelopment.com/download/ 8 つの観測変数、 54 の母数に使用が制限されている
まずはパス解析をやってみよう 共分散構造分析[AMOS編]より データ概要:結婚相談所において、女性が男性を自分の結婚相手としてどうかという観点から評価した架空のデータ4変数、210オブザベーション 変数 「学歴」:男性の教育水準 「年収」:男性の年収(百万円) 「職業威信」:職業に対する主観的な地位の格付け 「結婚相手としての評価」:男性の結婚相手としての評価
仮説構築 こんな感じで考えてみる 通常は先行研究やコンセプトワークを行った結果が反映される いきなり仮説作りはじめると???道に迷うことも多々あり 学歴 年収 評価 職業威信 Plan
じゃぁ???AMOSを起動して実行 まずはパス図を描く 誤差はしっかりつける その後???分析設定をする 推定法( estimation )はとりあえず最尤法( maximum likelihood) 1 2 3 標準化解( Standard estimates) 間接効果、総合効果( Indirect ,total effects) 修正指標( Modification Indice) Do
適合度を見てみよう GFI、AGFI、CFI 0.90以上が当てはまりの目安 RMSEA 0.05以下は良適合 0.10未満はグレーゾーン 0.10以上は適合× Check
解釈してみよう~パス係数 非標準化係数:独立変数が1単位変化したときの従属変数の変動量 標準化係数:標準化して相対的な影響力を表現したもの 今回は変数間の相対的な関係に興味があるので、標準化係数を中心に解釈します Action
解釈してみよう~直接効果と間接効果 直接効果:変数->変数の直接的な影響 間接効果:変数を経由した間接的な影響 総合効果:直接効果+間接効果 b a <職業威信から評価への効果> ■ 直接効果 a=0.30 ■ 間接効果 b=0.32*0.49=0.16 ■ 職業威信から評価への 総合効果 a+b=0.30+0.16=0.46 Action
結果のまとめ 結婚相手への評価に最も影響を与えるのは年収 職業威信の直接効果は年収の直接効果より低い しかし、年収を経由した間接効果も合せて考慮するとほぼ年収と同等の影響力がある 単純な回帰分析では導き出せない考察! 学歴->年収->評価の間接効果は0.26。 職業威信の直接効果より低い。
因子を組み合わせた分析 共分散構造分析 [AMOS 編 ] :パソコンセミナーの受講後アンケートデータ( 5 件法) 目的一致:内容は目的と一致していたか 理解度:内容はどの程度理解できたか テキスト:テキストについての評価 ペース:講義内容のペースについての評価 満足度:セミナー全体の満足度 プレゼン:講師のプレゼンについての印象 講師対処:質問の対処についての印象
仮説構築 「テキスト」、「プレゼン」、「ペース」、「講師対処」の背後に「講師の質」が存在すると想定 「満足度」、「理解度」、「目的一致」の背後に「セミナー評価」が存在すると想定 「講師の質」が「セミナー評価」に影響を与えていると想定 講師の質 テキスト プレゼン ペース 満足度 理解度 目的一致 セミナー評価 講師対処 1 2 3 Plan
適合度 GFI、AGFI、CFI、RMSEA全てで適合良し。 Check
結果の解釈 標準化解 講師の質は「プレゼン」の影響が最も高い セミナー評価は「満足度」の影響が最も高い 講師の質->セミナー評価のパス係数も大きい * Action
事例绍介
指導者が評価するサッカータレントの構造 サッカーコーチが評価するユース年代のサッカータレント構造を検証することを目的。 デルファイ法で仮説構築->調査票作成 243名の 有資格サッカーコーチが,プロレベル156名,全国レベル175名,地域レベル173名,計504名のサッカー選手のユース年代におけるサッカータレント を評価 SEMで検証を行う  http://www.taiiku.tsukuba.ac.jp/sc/1_1/09/index.html
まさに、王道! 仮説構築 SEM による検証 シンプレックス構造分析
発信する顧客は優良顧客か? コミュニティユーザーがオンラインコミュニティのどの要因を評価してリピート訪問をしているのかを特定する @ cosme ユーザーに対してアンケートを実施 併設オンラインショップ購買履歴データとのシングルソースデータ分析を行っている http://www.journalarchive.jst.go.jp/english/ jnlabstract_en.php?cdjournal=acs1993&cdvol=11&noissue=1-2&startpage=35
先行研究より、サイト訪問目的の構成概念を図 1 のように整理 その上で、以下の仮説を立てている サイトロイヤルティはサイト訪問目的の構成概念の関数 サイトロイヤルティと優良顧客度との間に有意で正の関数がある
コミュニティ上で書き込むユーザー ?ロイヤルティの源泉は購買前後サーチ、コンサマトリー、 B2C インタラクション(発信行為自体はロイヤルティへの影響は少ない) ?サイトロイヤルティから優良顧客へは正の影響 -0.018 0.040 0.232 0.240 0.604 0.449 0.215 購買データ
ブランドジャパン 日経 BP コンサルティングが主体となり実施 ブランドエクイティの測定 顧客から、一般社会からの評価を測定 国内ブランドを一般消費者、ビジネスパーソンが評価 2001 年から年 1 回のペースで実施し、経年の変化を追う 2003 年までは WEB 調査と郵送調査のハイブリッド 2004 年以降は WEB 調査のみとなり、傾向スコア補正している 構造方程式モデリングにより 2 つの観点からブランドエクティ(構成概念)をスコアリングしている http://www.nikkeibpm.co.jp/bz/chosa/brand_j/index.html
BtoC モデル 1 側面では測定できないブランドエクイティを 4 つの下位因子を構成して規定 下位因子をまとめる g を総合指標としてブランドエクイティを規定 ?????? - ????????? ??? ????????? ???????? g
BtoB モデル 5 つの下位因子と 5 つの観測変数で総合力「 g 」を規定 親和力 活力 先見力 信用力 人材力 g
スコアの推定 出典: http://consult.nikkeibp.co.jp/consult/release/bj_100409.pdf
2010 年のコメント( BtoC 編) ブランド総合力 92.3 ポイントを獲得したユニクロが、今回はじめて首位となった。 過去 5 年間の結果をみると、確実にブランド力を上昇させてきたことが分かる。(中略)しか し最大の理由は、イノベーティブ因子のポイントが驚異的に伸びたことだ。(中略)第 2 位のグーグルとは、 18.2 ポイントという大差がついている。 ブランド間のブランドエクイティ(構成概念スコア)を比較できる 時系列でのエクイティの変遷を追うことができる 下位因子を構成していることで、エクイティの要因分解を行うことができる
今後の勉強のために??? 初級 入門 共分散構造分析の実際(朝野?鈴木?小島) 講談社サイエンティフィック 共分散構造分析 [AMOS 編 ]  (豊田編著) 東京図書 はじめての共分散構造分析― Amos によるパス解析 (小塩) 東京図書 中級 グラフィカル多変量解析― AMOS 、 EQS 、 CALIS による 目で見る共分散構造分析  (狩野?三浦) 現代数学社 共分散構造分析 [ 入門編](豊田) 朝倉書店 共分散構造分析[事例編](豊田編著) 北大路書店
おまけ??? 質問に回答できなかった時系列分析 ラグつき変数の分散共分散行列(トープリッツ行列)を算出 し、そこから分析します。 自己回帰 AR ( 1 )モデルの方程式 V_t=α_t-1 * V_t-1 +e_t AR ( 2 )モデルの方程式 V_t=α_t-1 * V_t-1 + α_t-2 * V_t-2 +e_t 移動平均 MA ( 1 )モデルの方程式 V_t=F_t – β_t-1 * F_t-1 ARMA モデルの方程式 V_t=Σα_t-r * V_t-r  + Σβ_t-r’ * F_t-r’ *t が時系列をあらわす記号 詳しくは???豊田( 2000 )共分散構造分析[応用編]に記載があります

More Related Content

100614 構造方程式モデリング基本の「き」