[程式] R GEE的correlation matrix疑問
在做GEE模型評估的時候
比較的模型有3個
1.只放所有解釋變數的主效用，Correlation Structure為independence
2.只放所有解釋變數的主效用，Correlation Structure為exchangeble
3.只放所有解釋變數的主效用，Correlation Structure為AR1
在執行R的結果
得到的output為
[,1] [,2] [,3]
QIC 801.3225 801.3225 823.66768
CIC 41.0889 41.0889 39.22805
在1和2的估計值也都是一樣的
程式的部分是參考下列網址
http://www.unc.edu/courses/2010spring/ecol/562/001/docs/lectures/lecture14.htm
請問是哪裡有問題嗎?
另外還有就是將working correlation matrix 設成 unstructured
得到的估計working correlation matrix裡面的參數有的估計值大於1
這個是因為甚麼原因呢?
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.116.152.115
※ 文章網址: http://www.ptt.cc/bbs/Statistics/M.1418276841.A.F97.html
※ 編輯: tokyo291 (140.116.152.115), 12/11/2014 13:47:41
※ 編輯: tokyo291 (140.116.152.115), 12/11/2014 13:56:25
※ 編輯: tokyo291 (140.116.152.115), 12/11/2014 13:56:30
※ 編輯: tokyo291 (140.116.52.59), 12/11/2014 17:46:56