作者clsmbstu ()
看板Statistics
標題[程式] R MCMCglmm prior參數設定
時間Sun Dec 25 23:08:44 2016
[軟體程式類別]:
R
[程式問題]:
Markov chain Monte Carlo (MCMC) multivariate generalized linear mixed model
(MCMCglmm)
[軟體熟悉度]:
熟悉
[問題敘述]:
我有一個data frame如下連結：
https://dl.dropboxusercontent.com/u/63339292/PttHelp.RData
Stimuli：受試者在該trial接受之刺激圖檔名
Dimension.1~6：受試者回答六個問題，每個問題都是0~8九點量表
Subject_Num：受試者編號
Subject_Age：受試者年齡組別
Subject_Sex：受試者生理性別
我希望：
1. 以Dimension.1~6做為六個依變數
2. 以受試者年齡與生理性別為fixed effects並考慮其交互作用
3. 以刺激圖檔與受試者編號為random effects
主要問題：
雖然正在讀原套件作者的論文，
Jarrod D Hadfield (2010). MCMC Methods for Multi-Response Generalized Linear
Mixed Models: The MCMCglmm R Package. Journal of Statistical Software, 33(2),
1-22. URL
http://www.jstatsoft.org/v33/i02/.
但還是完全不知道"prior"這個參數內容應該要放什麼、怎麼寫，
另外"rcov"也不確定寫得對不對。
事實上我的統計知識要理解這篇有點吃不消，
如果板友有提供入門的資料推薦給我也很好。
最後出現的錯誤為：
Error in MCMCglmm(cbind(Dimension.1, Dimension.2, Dimension.3, Dimension.4,
:
Mixed model equations singular: use a (stronger) prior
（也出現過其它幾種不同的錯誤，但看起來都跟G-structure或prior有關。）
[程式範例]:
https://gist.github.com/anonymous/54bf8d9d4153900803c4fd29e30a1be2
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 114.44.100.249
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1482678530.A.63C.html
1F：→ andrew43: 除非你把所有變數的細節交待清楚，否則不能回答。 12/26 00:18
2F：→ andrew43: 主要原因是你不了解G-structure到底要怎麼決定，而這和 12/26 00:19
3F：→ andrew43: 整個模型中參與的變數數量有關。 12/26 00:19
4F：→ andrew43: sorry看到你附資料了 12/26 00:23