作者warex14 (航)
看板R_Language
標題Package 選擇 (Logit Regression)
時間Mon May 25 16:21:19 2015
[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
請把以下不需要的部份刪除
入門(寫過其他程式，只是對語法不熟悉)
[問題敘述]:
請簡略描述你所要做的事情，或是這個程式的目的
大家好，目前小弟要複製某篇論文的結果，(只有改變樣本期間)
所使用的方法是用Maximum Likelihood Estimates 來估計Logit Model
但是上CRAN 找 Package的時候發現有很多可能的Packages
以下是我找到的：
mlogit: multinomial logit model
mlogitBMA: Bayesian Model Averaging for Multinomial Logit Models
mnlogit: Multinomial Logit Model
pln: Polytomous logit-normit (graded logistic) model estimation
gmnl: Multinomial Logit Models with Random Parameters
logitnorm: Functions for the logitnormal distribution
目前我是覺得mnlogit: Multinomial Logit Model
比較有可能是我需要的工具，以下是他的說明：
(Time and memory efficient estimation of multinomial logit models using
maximum likelihood method. Numerical optimization performed by Newton-Raphson
method using an optimized, parallel C++ library to achieve fast computation
to achieve fast computation of Hessian matrices. Motivated by large scale
multiclass classification problems in econometrics and machine learning.)
不知道有沒有人使用過類似的Package？
他會是我所需要的工具嗎？
(目前打算載來先用用看)
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.108.38
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1432542081.A.439.html
1F：→ andrew43: 內建的glm就可以了不是嗎 05/25 16:28
2F：→ andrew43: 還是你要做multinomial logit model? 05/25 16:30
驚！原來有內建好的？ 我的Logit Model 的dependent variable變數是dummy(0或1)
所以用glm就可以了？
(所以Y變數是分成多類的才是multinomial logit model?)
※ 編輯: warex14 (140.112.108.38), 05/25/2015 17:11:05
3F：→ andrew43: 我不覺得有什麼不能進行的, 貼出你的code才能說是否正確 05/25 17:57
4F：→ andrew43: multinomial logit model 是應變數為三類以上的模型 05/25 17:58
5F：→ andrew43: 另外，如果你只有一個類別的自變數，那做獨立性檢驗就好 05/25 17:59
6F：→ andrew43: 看錯了，上一列應該不是你的情況 05/25 18:02
突然想到我的data好像不能用一般的Logit Model直接下去做
因為他可能有autocorrelation的問題：
年分 公司名稱 Y變數 X1 X2 X3...
1991 AAA 1 12 5.3 789
1991 AAB 0 14 4.8 800
1991 AAC 1 9 2.2 436
...
...
1991 ZZZ 1
1992 AAA 0
1992 AAC 1
...
...
2014 ZZZ 0
不知道這樣還能不能用glm做？
※ 編輯: warex14 (140.112.108.38), 05/26/2015 17:44:32
7F：→ andrew43: glm不適合。存在不獨立樣本。可能可以做mixed model。 05/26 18:57