Re: [問題] [R]跑PANEL DATA

看板Statistics作者 (微笑台灣)時間17年前 (2009/03/31 10:10), 編輯推噓0(009)
留言9則, 2人參與, 最新討論串2/2 (看更多)
有人建議我使用 model = "pooling" 但是我用的模型是固定效果模型 plm這個package 只有model = "within" 是固定效果模型 可是我怎麼試都是沒有截距項 想請問 還有別的package 是可以跑panel data的嗎 謝謝 ※ 引述《paiyeu (微笑台灣)》之銘言: : 我現在使用R跑panel data : 使用 plm這個package : 指令為 : fe <- plm(y~x+z, data, model = "within") : 可是我希望有截距項 : 因為跑出來的結果 只有x z這兩個變數的估計直 : 請問要怎麼設才會使式子有截距項呢 : 因為我看這個package的說明文件 : 跑random effect時 有截距項 但是式子中也沒有特別設 : 想請問板上的人 是否有解決之道 : 謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.119.52.177

03/31 10:14, , 1F
啥是 "固定效果模型"?
03/31 10:14, 1F

03/31 10:14, , 2F
自己加一個變數 w, w <- rep(1, nrow(data))
03/31 10:14, 2F

03/31 13:37, , 3F
就是panel-data下的固定效果模型,可是你建議的方法我試了
03/31 13:37, 3F

03/31 13:37, , 4F
還是不行 不過謝謝你
03/31 13:37, 4F

03/31 15:08, , 5F
如果你說的"固定"指得是individual effect. 那你要的截
03/31 15:08, 5F

03/31 15:09, , 6F
距項, 即為individual effect的平均.
03/31 15:09, 6F

03/31 15:17, , 7F
sorry,上面是錯的. 你要得可能是這個.
03/31 15:17, 7F

03/31 15:17, , 8F
z <- fixef(fe); attr(z, "intercept")
03/31 15:17, 8F

03/31 15:19, , 9F
基本上, lme 會比 plm 好用. 參考文件也比較多.
03/31 15:19, 9F
文章代碼(AID): #19qNkaD9 (Statistics)
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 2 篇):
文章代碼(AID): #19qNkaD9 (Statistics)