R2winbugs到出数据时出现错误

[code] we might assume that F{θ; xi (t)} follows a marginalized transition
model (MTM) ， with marginal mean g{μ (t)} = xi (t)β,文献中就是这样提了下，具体程序中就是那样了

文雅13 发表于 2012-12-17 16:11
是在文献中找到的程序，程序就是那样的

那你就向文献的作者要要看
看他是否愿意提供
否则你就要依需求,自行编程

# MTM1Delta is a function written in the WinBUGS development         
# interface for computing the intercept delta[i,j] to satisfy
# marginal mean and serial dependence simultaneously

# Mv is a function written in the WinBUGS development
# interface for computing smooth function estimates
# with low-rank thin-plate bases

epoh 发表于 2012-12-17 16:27
那你就向文献的作者要要看
看他是否愿意提供
否则你就要依需求,自行编程

老师，在检验Winbugs模型时出现loop name already used in outer loop是什么问题？

文雅13 发表于 2013-1-17 09:09
老师，在检验Winbugs模型时出现loop name already used in outer loop是什么问题？

for (i in 1:n) // outer loop
{
    for (i in 1:p) // inner loop
    {
    //...do something
    }
//...do something else
}

loop name already used in outer loop

epoh 发表于 2013-1-17 15:36
for (i in 1:n) // outer loop
{
    for (i in 1:p) // inner loop

老师，我没太看明白这些code是什么意思？

就是for (i in 1:n) 地方有错误

epoh 发表于 2013-1-17 15:36
for (i in 1:n) // outer loop
{
    for (i in 1:p) // inner loop

老师，没明白您的意思，能说的详细点吗？

文雅13 发表于 2013-1-17 19:46
老师，没明白您的意思，能说的详细点吗？

也就是outer loop 用i,inner loop 就要用j,不能用相同的i

epoh 发表于 2013-1-17 20:42
也就是outer loop 用i,inner loop 就要用j,不能用相同的i

用的就是i ，但鼠标还指着这里面的ifor (i in 1:n)，在Winbugs中 gamma as a matrix (two indices)and a
vector (one index)
是指什么？

epoh 发表于 2013-1-17 20:42
也就是outer loop 用i,inner loop 就要用j,不能用相同的i

1.   y[i,j]~dbern(p[i,j])
   
2.                  logit(p[i,j])<-delta[i,j]+gamma[i,j]*y[i,j-1]
   
3.                  etaM[i,j]<-beta0all[i]+beta1all[i]*t[(sumn[i]-nobs[i]+j)]
   
4.            
   
5.                
   
6.                   delta[i,j]<- -log((-b[i,j]+sqrt(b[i,j]*b[i,j]-4*a[i,j]*c[i,j]))/(2*a[i,j]))
   
7.                      
   
8.                   
   
9.                   a[i,j]<-exp(-gamma[i,j])*exp(etaM[i,j])/(1+exp(etaM[i,j]))
   
10. 
    
11.                   b[i,j]<-(exp(etaM[i,j])/(1+exp(etaM[i,j]))+exp(etaM[i,j-1])/(1+exp(etaM[i,j-1]))-1)*exp(-gamma[i,j])+exp(etaM[i,j])/(1+exp(etaM[i,j]))-exp( etaM[i,j-1])/(1+exp(etaM[i,j-1]))
    
12. 
    
13.                   c[i,j]<- exp(etaM[i,j])/(1+exp(etaM[i,j]))-1

复制代码

比如这个code。gamman有矢量也有矩阵吗？不会看，哎