[程式] R MCMCglmm prior參數設定

[clsmbstu ()]

[軟體程式類別]:
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://www.ptt.cc/bbs/Statistics/M.1482678530.A.63C.html
→ andrew43: 除非你把所有變數的細節交待清楚，否則不能回答。
→ andrew43: 主要原因是你不了解G-structure到底要怎麼決定，而這和
→ andrew43: 整個模型中參與的變數數量有關。
→ andrew43: sorry看到你附資料了