作者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/cn.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