2016-05-03 70 views
1

對於實習生來說,我正在評估一個支持環保的倡議,其中成員可以向鄰居宣傳倡議以說服他們加入。R中的多級Logistic迴歸

我想從一組預測的預測組成員(二進制):

  1. 人口統計學(4分類,2個連續變量)
  2. 心理變量(9個連續變量)鄰里的
  3. 曝光(8個連續的變量)的主動
  4. 曝光(2分類,2連續)

數據收集在30個不同社區所以我的主管建議我在我的分析中控制鄰里。

我想現在要做的是:

添加四套預測此起彼伏,看看模型擬合度上升 - 而且,如果是這樣,由於在該變量設置 - ,而控制鄰里。

是否有任何R軟件包的文檔可以建立這樣的模型?任何有關這種建模,假設,結果解釋等的文獻提示都非常值得歡迎!

我以前只使用過非常基本的多級模型,所以在這種情況下我迷失了方向。提前謝謝了!

+1

這聽起來像你有正確的想法使用邏輯迴歸。我只是谷歌邏輯迴歸在R',我敢肯定你會發現大量的視頻,文章,堆棧溢出的例子等。 –

回答

3

您的問題中沒有任何代碼,因此可以考慮將其移至交叉驗證。至於學習如何在R模型,谷歌將給many suggestions。但是,爲了給你一個基本的起點

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

其中「組」是你的迴應,你的預測去哪裏我已經把D1 + ...和你在模型中包括鄰居爲好。您所描述的任何內容都是嵌套的,所以這應該可以控制鄰居效應。然後,您可以使用MASS包中的stepAIC執行模型選擇。

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

當然還有其他更復雜的程序可以使用。

+1

嘿,謝謝你的回答! 我認爲只包括鄰居的假人可能不夠,因爲主動中的分享者和非分享者的比例相差很大。 我現在使用lme4軟件包爲每個鄰域做一個隨機截取的邏輯迴歸模型,你會說這是一個同樣可行的方法嗎? 此外,是否有辦法逐步包含整套預測變量,而不僅僅是單個預測變量? –

+0

這聽起來像你在正確的軌道上 - 如果你使用lme,你已經超過了我的基本範例,並且肯定是在正確的軌道上。沒有示例數據和特定問題很難提供更多幫助。我不確定你的意思是「逐步包括整套預測變量」。你的意思是在模型公式中?如果是這樣,你只需要做到這一點......這是分析的人的一部分。 – JohnSG