Package 選擇 (Logit Regression)

[warex14 (航)]

[問題類型]:

程式諮詢(我想用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://www.ptt.cc/bbs/R_Language/M.1432542081.A.439.html
→ andrew43: 內建的glm就可以了不是嗎
→ andrew43: 還是你要做multinomial logit model?
驚！原來有內建好的？ 我的Logit Model 的dependent variable變數是dummy(0或1)
所以用glm就可以了？
(所以Y變數是分成多類的才是multinomial logit model?)
※ 編輯: warex14 (140.112.108.38), 05/25/2015 17:11:05
→ andrew43: 我不覺得有什麼不能進行的, 貼出你的code才能說是否正確
→ andrew43: multinomial logit model 是應變數為三類以上的模型
→ andrew43: 另外，如果你只有一個類別的自變數，那做獨立性檢驗就好
→ andrew43: 看錯了，上一列應該不是你的情況
突然想到我的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
→ andrew43: glm不適合。存在不獨立樣本。可能可以做mixed model。