[問題] Package 選擇 (Logit Regression)

看板R_Language作者 (航)時間9年前 (2015/05/25 16:21), 9年前編輯推噓0(007)
留言7則, 1人參與, 最新討論串1/1
[問題類型]: 程式諮詢(我想用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

05/25 16:28, , 1F
內建的glm就可以了不是嗎
05/25 16:28, 1F

05/25 16:30, , 2F
還是你要做multinomial logit model?
05/25 16:30, 2F
驚!原來有內建好的? 我的Logit Model 的dependent variable變數是dummy(0或1) 所以用glm就可以了? (所以Y變數是分成多類的才是multinomial logit model?) ※ 編輯: warex14 (140.112.108.38), 05/25/2015 17:11:05

05/25 17:57, , 3F
我不覺得有什麼不能進行的, 貼出你的code才能說是否正確
05/25 17:57, 3F

05/25 17:58, , 4F
multinomial logit model 是應變數為三類以上的模型
05/25 17:58, 4F

05/25 17:59, , 5F
另外,如果你只有一個類別的自變數,那做獨立性檢驗就好
05/25 17:59, 5F

05/25 18:02, , 6F
看錯了,上一列應該不是你的情況
05/25 18:02, 6F
突然想到我的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

05/26 18:57, , 7F
glm不適合。存在不獨立樣本。可能可以做mixed model。
05/26 18:57, 7F
文章代碼(AID): #1LOjk1Gv (R_Language)
文章代碼(AID): #1LOjk1Gv (R_Language)