R2winbugs到出数据时出现错误

文雅13 发表于 2012-11-21 19:45
我所给出的是原始数据的纵向格式，您说的是winbugs格式的数据吗？

这个模型比较复杂,
需要提供完整的数据及R SCRIPT比对,
既然是R SCRIPT,当然是纵向的原始数据.

epoh 发表于 2012-11-21 20:53
哈哈!抱歉,zhangtao兄,这个有难度.

epoh老师，您看有什么难度？

zhangtao 发表于 2012-11-21 20:55
epoh老师，您看有什么难度？

哈哈你不是说了,
自体抽样扩大样本数据量,
总该有个依据?

epoh 发表于 2012-11-21 20:53
这个模型比较复杂,
需要提供完整的数据及R SCRIPT比对,
既然是R SCRIPT,当然是纵向的原始数据.

2.txt (1.97 KB)

2012-11-21 21:17:50 上传

这就是原数据了  经R2winbugs到处时出现问题

文雅13 发表于 2012-11-21 21:20
这就是原数据了  经R2winbugs到处时出现问题

初步看了下,你的数据及数据格式很多要修改.
给你重点,你再慢慢修改.
n是样本量为150，n1是14，跟你的2.txt根本不符
y[i,1]~dbern(p[i,1])
y要是一个  n1 x 1 matrix

for (i in (n1+1):n)
{
y[i,1]~dbern(p[i,1])
logit(p[i,1])<-etaM[i,1]
etaM[i,1]<-beta0all+beta1all*t[(sumn-nobs+1)]
  当i=150,nobs早超出范围

BZ[i,1:9]
  BZ要是一个n x 9 matrix
....
....
其余请自行更正

epoh 发表于 2012-11-21 21:57
初步看了下,你的数据及数据格式很多要修改.
给你重点,你再慢慢修改.
n是样本量为150，n1是14，跟你的2. ...

谢谢epoh老师

epoh 发表于 2012-11-21 21:57
初步看了下,你的数据及数据格式很多要修改.
给你重点,你再慢慢修改.
n是样本量为150，n1是14，跟你的 ...

epoh老师，您根据我的winbugs code，能举个例子吗？说明数据格式应该是什么样的格式才能通过R2导出为winbugs所需的数据，谢谢！

文雅13 发表于 2012-11-22 14:00
epoh老师，您根据我的winbugs code，能举个例子吗？说明数据格式应该是什么样的格式才能通过R2导出为winb ...

需要你完整的正确数据,
才能转换处理.

另这个模型我不熟,
需要你自行订出initial values
利用gen inits可能会有问题

epoh 发表于 2012-11-22 14:25
需要你完整的正确数据,
才能转换处理.

epoh老师
1、2.txt就是我数据的一小部分，格式和分析的变量都在里面，您看看那个纵向格式不正确吗？
2、还有我之前传的R2转换数据的程序是不是不能用于处理这样的数据，如果需要改，改哪个地方好呢？

文雅13 发表于 2012-11-22 14:38
epoh老师
1、2.txt就是我数据的一小部分，格式和分析的变量都在里面，您看看那个纵向格式不正确吗？
...

2.txt数据格式有误,无法读入,修改为ecs.csv
   binary_CES.rar (692 Bytes) 本附件包括：

• ces.csv

2012-11-22 19:56:07 上传

######
data=read.csv("ces.csv",header=T)
data
dim(data)  # 18x18
y=matrix(data[,16],18,1)
t=data$t
nobs=data$nobs
cd4=data$SEX
race=data$AGE
BZ=data[,6:14]
变量名称必须与winbugs code一致,
尚缺许多请自行补全
#####
底下的似乎都没定义
knotD[]
OMEGAu[]
del[]
nn1[]
nn2[]
nn3[]
nn4[]