2016-05-03 26 views
1

訓練生にとって、私はメンバーがイニシアチブを近隣の人たちに宣伝して参加を促す環境取り組みを評価しています。マルチレベルロジスティック回帰R

Iは、予測子の集合から(バイナリ)のグループメンバーシップを予測する:

  1. 人口統計(4カテゴリ、2つの連続変数)
  2. 心理的変数(9つの連続変数)近隣の
  3. インプレッション(8つの連続変数)イニシアチブの
  4. インプレッション(2カテゴリ、2連続)

の30の異なる地域のデータがで収集されたため、私の上司が分析で近所の管理を提案しました。私は今、何をしたいか

がある:

モデルの適合が上がるかどうかを確認するために、別の後に予測因子の1の4つのセットを追加します - もしそうなら、そして、原因その中の変数へ近傍についてはを制御している間にset_となる。

このようなモデルを構築するためのドキュメントがあるRパッケージはありますか?この種のモデリング、前提条件、結果の解釈などに関する文献のヒントも歓迎すべきものです。

以前は非常に基本的なマルチレベルモデルでしか作業していなかったので、この状況ではかなり迷っています。事前に多くの感謝!

+1

ロジスティック回帰を使用するのは正しい考えがあるようです。私はGoogleのロジスティック回帰をRにするだけで、ビデオ、記事、Stack Overflowなどのサンプルをたくさん見つけることができると確信しています。 –

答えて

3

質問にはコードがありませんので、Cross Validatedに移動することを検討してください。 Rでモデリングする方法を学習すると、Googleはmany suggestionsを与えるでしょう。しかし、あなたが、私はD1 + ...を入れているあなたの応答は、あなたの予測因子が行く「グループ」は、基本的な出発点

glm_out <- glm(group ~ D1 + ... + Neighborhood, family = binomial(), data = D) 

を与え、あなたにもモデルに近所を含めます。あなたが記述したものはネストされていませんので、これはNeighborhoodエフェクトを制御するために機能するはずです。次に、MASSパッケージのstepAICを使用してモデルの選択を実行できます。

glm_aic <- MASS::stepAIC(glm_out) 
summary(glm_aic) 

他にも複雑な手順があります。

+1

お返事ありがとうございます! イニシアチブ参加者と非参加者の割合がかなり異なるため、Neighborhoodのダミーを含めるだけでは十分ではない可能性があります。 lme4パッケージを使用して、各近隣に対してランダムな代行受信を伴うロジスティック回帰モデルを作成しましたが、それは同様に実行可能な方法でしょうか? さらに、単一のプレディクタを順に処理するのではなく、段階的にプレディクタのセット全体を含める方法がありますか? –

+0

あなたは正しい軌道に乗っているように聞こえます。もしあなたが私の基本的な例を過ぎてしまっているのであれば、正しい軌道にはいります。データの例や特定の質問がなければ、より多くのヘルプを提供することは難しいです。私は、あなたが「予測器のセット全体を段階的に含める」ということをどういう意味か分かりません。あなたはモデル式を意味しますか?もしそうなら、あなたはそれを行う必要があります...これは分析の人間の部分です。 – JohnSG