楼主: 一诺9257
1985 4

[问答] WinBUGS纵向数据怎么输入? [推广有奖]

  • 10关注
  • 7粉丝

教授

66%

还不是VIP/贵宾

-

威望
0
论坛币
23 个
通用积分
3.8021
学术水平
28 点
热心指数
34 点
信用等级
17 点
经验
301 点
帖子
949
精华
0
在线时间
1321 小时
注册时间
2010-3-21
最后登录
2020-2-22

楼主
一诺9257 发表于 2013-6-20 10:09:28 |AI写论文
5论坛币
请问WinBUGS怎么做随机截距和斜率、纵向数据?下面数据怎么输入到输入到WinGUGS中?每一个个体观察八次)谢谢                                  week
Subject Group 1 2 3 4 5 6 7 8
1 A 45 45 45 45 80 80 80 90
2 A 20 25 25 25 30 35 30 50
3 A 50 50 55 70 70 75 90 90
4 A 25 25 35 40 60 60 70 80
5 A 100 100 100 100 100 100 100 100
6 A 20 20 30 50 50 60 85 95
7 A 30 35 35 40 50 60 75 85
8 A 30 35 45 50 55 65 65 70
9 B 40 55 60 70 80 85 90 90
10 B 65 65 70 70 80 80 80 80
11 B 30 30 40 45 65 85 85 85
12 B 25 35 35 35 40 45 45 45
13 B 45 45 80 80 80 80 80 80
14 B 15 15 10 10 10 20 20 20
15 B 35 35 35 45 45 45 50 50
16 B 40 40 40 55 55 55 60 65
17 C 20 20 30 30 30 30 30 30
18 C 35 35 35 40 40 40 40 40
19 C 35 35 35 40 40 40 45 45
20 C 45 65 65 65 80 85 95 100
21 C 45 65 70 90 90 95 95 100
22 C 25 30 30 35 40 40 40 40
23 C 25 25 30 30 30 30 35 40
24 C 15 35 35 35 40 50 65 65


关键词:winbugs WINBUG 纵向数据 BUGS bug

沙发
一诺9257 发表于 2013-6-20 17:17:05
model
{
# likelihood
for(subject in 1:N){ # loop in subject
intercept[subject]~dnorm(0,tau.subject);
for(time in 1:T){ # loop in time
mu[subject,time]<-intercept[subject] +
alpha.c[group[subject]] + (beta[group[subject]]*time);
ability[subject,time]~dnorm(mu[subject,time],tau.resid);
} # end of time loop
} # end of subject loop
# priors
var.subject~dunif(0,1.0E4);
var.resid~dunif(0,1.0E4);
beta[1]~dnorm(0,1.0E-4); # Linear effect of time (group=A)
beta[2]~dnorm(0,1.0E-4); # Linear effect of time (group=B)
beta[3]~dnorm(0,1.0E-4); # Linear effect of time (group=C)
alpha.c[1]~dnorm(0,1.0E-4); # Centered intercept (group=A)
alpha.c[2]~dnorm(0,1.0E-4); # Centered intercept (group=B)
alpha.c[3]~dnorm(0,1.0E-4); # Centered intercept (group=C)
# scalars
tau.subject<-1/var.subject;
tau.resid<-1/var.resid;
rho<-var.subject/(var.resid+var.subject); # Within-sub corr.
b.diff[1]<-beta[2]-beta[1];
b.diff[2]<-beta[3]-beta[1];
alpha[1]<-alpha.c[1]+50; # re-scale intercepts
alpha[2]<-alpha.c[2]+50;
alpha[3]<-alpha.c[3]+50;
a.diff[1]<-alpha[2]-alpha[1];
a.diff[2]<-alpha[3]-alpha[1];
}
这是model代码!不知道怎么载入数据

藤椅
liuqi99 发表于 2013-8-31 19:36:08
你把数据换成TXT,用WINBUGS打开

板凳
liuqi99 发表于 2013-8-31 20:00:30
你的数据是哪里来的?你还需要初值数据。

报纸
一诺9257 发表于 2013-9-2 07:16:13
是的!谢谢

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-27 00:42