求助:我的winbugs程序有问题,哪位大侠帮我看看-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 会计>>

会计库

>>

求助:我的winbugs程序有问题,哪位大侠帮我看看

求助:我的winbugs程序有问题,哪位大侠帮我看看

发布:dailina1975 | 分类:会计库

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

经管之家新媒体交易平台

提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

提供微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

小女子在作一个研究,具体x~dmnorm(mu[],T[,]),mu是一个6维的列向量omega~逆wish(R[,],7),R是一个6维对角矩阵,并且其对角矩阵上的6个元素都服从dgamma(0.001,0.001),现在我想估计出omega及R,我编的程序如下:m ...
坛友互助群


扫码加入各岗位、行业、专业交流群


小女子在作一个研究,具体x~ dmnorm(mu[], T[,]),mu是一个6维的列向量omega~ 逆wish(R[,], 7),R是一个6维对角矩阵,并且其对角矩阵上的6个元素都服从dgamma(0.001, 0.001),现在我想估计出omega及R,我编的程序如下:
model
{
for(i in 1:N){x[i,1:6] ~ dmnorm(mu[], T[,])}
T[1:6, 1:6]~ dwish(R[,], 7)
R[1,1]~dgamma(0.001, 0.001)
R[2,2]~dgamma(0.001, 0.001)
R[3,3]~dgamma(0.001, 0.001)
R[4,4]~dgamma(0.001, 0.001)
R[5,5]~dgamma(0.001, 0.001)
R[6,6]~dgamma(0.001, 0.001)
omega[1:6,1:6]<-inverse(T[,])
}
list(N=15,mu = c(0, 0, 0, 0, 0, 0), x=structure(
.Data=c(-1.1108, -0.84162, 0.84162, -0.43073, 0.43073, 2.3263,
0.62293,-0.43073, -0.083652, -0.25335, -1.1108, 0.84162,
-0.43073, 0.083652, -0.62293, 0.84162, -1.5011, 0.84162,
0.25335, 2.3263, 3.0902, -0.43073, -0.25335, 0.84162,
0.083652, 0.25335, -0.62293, -0.43073, -0.84162, 1.5011,
-0.84162, -0.083652, 1.1108, -0.43073, 0.62293, 0.84162,
0.43073, -0.62293, 0.083652, 0.43073, -0.43073, 0.84162,
0.84162, 1.5011,0.43073, 1.5011, -0.62293, 0.84162,
-0.083652, 0.43073, -0.43073, 0.083652, 0.083652, 0.84162,
-0.62293, 1.1108, 0.25335, 1.1108, 2.3263, 0.84162,
2.3263, 0.84162, 1.5011, 3.0902, 1.5011, 0.84162,
-1.5011, 0.62293, -0.25335, -0.43073, 0.25335, 0.84162,
1.1108, -1.5011, -0.62293, -0.083652, 1.1108, 0.84162,
1.5011, -0.25335, -0.62293, 0.62293, -0.083652, 1.1108,
-0.25335, -1.1108, 0.62293, 0.25335, 0.84162, 0.84162),
.Dim=c(15,6)),
R(1,2:6)=c(0,0,0,0,0), R(2,1)=0, R(2,3:6)=c(0,0,0,0), R(3,1:2)=c(0,0), R(3,4:6)=c(0,0,0),
R(4,1:3)=c(0,0,0), R(4,5:6)=(0,0), R(5,1:4)=C(0,0,0,0), R(5,6)=0, R(6,1:5)=C(0,0,0,0,0))

list(R[1,1] = 0.01, R[2,2] = 0.01, R[3,3] = 0.01, R[4,4] = 0.01, R[5,5] = 0.01, R[6,6] = 0.01,T = structure(
.Data = c(0.1, 0.3, 0.1, 0.2, 0.1, 0.1,

0.1, 0.1, 0.7, 0.1, 0.1, 0.3,

0.1, 0.3, 0.1, 0.1, 0.2, 0.1,

0.4, 0.1, 0.1, 0.5, 0.1, 0.1,

0.1, 0.6, 0.1, 0.1, 0.1, 0.7,

0.3, 0.1, 0.1, 0.1, 0.2, 0.1),

.Dim=c(6, 6)))

黑色字体没有问题,但在load初始值时,显示:expected an equals sign,我也不知道是哪里错了。请高人帮我指正。
扫码或添加微信号:坛友素质互助


「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-868593-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。