[问答] R2winbugs到出数据时出现错误 [推广有奖]

101楼

文雅13 发表于 2013-2-6 20:55:26

epoh 发表于 2013-2-6 15:53
请在这个帖子回答问题,
这样才有连贯性.
请先确认底下数据:

现在加上y=matrix(data = y, nrow = 317, ncol = 4, byrow = FALSE)，在R中出现
错误于FUN(X[[3L]], ...) : 不能改变因子的存储方式，更没法进入Winbugs中。

102楼

epoh 发表于 2013-2-6 21:16:05

文雅13 发表于 2013-2-6 20:55
现在加上y=matrix(data = y, nrow = 317, ncol = 4, byrow = FALSE)，在R中出现
错误于FUN(X[[3L]], ... ...

那是因为你的错误一箩筐,
你只改那行是没用的,
我已经跑到gen.inits(),才出现大问题.
需要你提供正确的数据(n1,sumn),再想办法

display(log)
check(c:/Bugs/ces_d/ces_d.txt)
model is syntactically correct
data(c:/Bugs/ces_d/data.txt)
data loaded
compile(1)
model compiled
gen.inits()

103楼

文雅13 发表于 2013-2-6 21:41:04

epoh 发表于 2013-2-6 21:16
那是因为你的错误一箩筐,
你只改那行是没用的,
我已经跑到gen.inits(),才出现大问题.

1111111111

1.rar
下载链接: https://bbs.pinggu.org/a-1268106.html

758 Bytes

本附件包括：

1.sav

104楼

epoh 发表于 2013-2-6 22:02:35

文雅13 发表于 2013-2-6 21:41
1111111111

这个数据还是不对
还是会产生负值
i=1,sumn[1]-nobs[1]+1= -2,

105楼

epoh 发表于 2013-2-7 12:14:05

文雅13 发表于 2013-2-6 20:55
现在加上y=matrix(data = y, nrow = 317, ncol = 4, byrow = FALSE)，在R中出现
错误于FUN(X[[3L]], ... ...

由程序看出数据尚存矛盾之处:

1.sample size n=317,但t的长度是1268
显然不对,需要重上传
2.only one observation,n1=2,

for (i in 1:n1)
{
y[i,1]~dbern(p[i,1])
...
etaM[i,1]<-beta0all[i]+beta1all[i]*t[(sumn[i]-nobs[i]+1)]
}

复制代码

当i=1,2时
sumn=[ 1 2 5 7...1202]
nobs=[ 4 4 2 4.............4 4]

这个相对应的数据对吗?

106楼

epoh 发表于 2013-2-7 18:14:04

文雅13 发表于 2013-2-6 21:41
1111111111

Data are sorted by the race and baseline CD4 groups,
and consist of the binary outcome variable (y[,])

107楼

文雅13 发表于 2013-2-9 13:30:32

epoh 发表于 2013-2-6 21:16
那是因为你的错误一箩筐,
你只改那行是没用的,
我已经跑到gen.inits(),才出现大问题.

方便的话发送消息看看您的数据和程序？

108楼

epoh 发表于 2013-2-9 22:45:30

文雅13 发表于 2013-2-9 13:30
方便的话发送消息看看您的数据和程序？

数据格式变更的只有y,       #2013.02.06,我已回帖告诉你
数据要更正的仅是n1=2,    #n1是由你短信息提供的
其余的通通没变.

y<-response$y
t<-response$t
sumn=response$sumnobs

sqbh<-cov$sq    #1 x 317
sex<-cov$sex    #1 x 317
nobs<-cov$nobs    #1 x 317
n<-length(nobs)
n1=2

knotD
OMEGA_all
...
nn3
nn4

但是t,sumn,nobs,数据有误,2013.02.07,我已告诉你要检验
这会造成gen.inits()后,Trap:undefined real result
只是你一直没回复,
所以程序是为了错误数据能运行所做出的非常措施,
严格说你的程序我并没改变.
但底下数据你终究要面对,检验,解释,尤其是sort,才不会原地踏步

1.sample size n=317,但t的长度是1268, 显然不对,

2.only one observation,n1=2,
  for (i in 1:n1)

   {
   y[i,1]~dbern(p[i,1])
   ...
   etaM[i,1]<-beta0all+beta1all*t[(sumn-nobs+1)]
   }
当i=1,2时
sumn=[ 1 2 5 7...1202]
nobs=[ 4 4 2 4.............4 4]