多个样本多个时间点重复观察计量资料的统计

5664

收藏 2012-03-21

下表数据，想计算一下A、B不同组间药后1h、2h的数据分别是否存在差别。

组别	编号	服药前	药后1h	药后2h
A	1	52	49	42
A	2	51	50	46
A	3	50	49	41
A	4	51	49	44
A	5	49	47	40
B	1	51	54	53
B	2	49	47	46
B	3	50	47	44
B	4	49	48	41
B	5	52	50	48

我查了赵耐青的《临床医学研究设计和数据分析》，其中要把数据写为下表格式

数据y	组别g	编号no	t1	t2
52	0	1	0	0
49	0	1	1	0
42	0	1	0	1
……
50	1	10	1	0
48	1	10	0	1

再用xtreg运算。但我的数据比这个要多很多，转换很麻烦，想请问一下，是否有什么方法，可以直接用表1中的数据进行计算，谢谢。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

aolei

2012-3-21 22:37:42

g、t1、t2是什么含义？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jj88625

2012-3-21 23:08:49

g是组别，0代表A组，1代表B组。t1为1，表示药后1h，t2为1，表示药后2h。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

aolei

2012-3-21 23:31:37

clear all

input str10 b n y1 y2 y3
A 1 52 49 42
A 2 51 50 46
A 3 50 49 41
A 4 51 49 44
A 5 49 47 40
B 1 51 54 53
B 2 49 47 46
B 3 50 47 44
B 4 49 48 41
B 5 52 50 48
end

gen t=_n
reshape long y,i(t) j(s)

tostring s,replace
gen t1=strpos(s, "2")!=0
gen t2=strpos(s, "3")!=0
egen g=group(b)
replace g=g-1
keep y g n t1 t2
order y g n t1 t2
l

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jj88625

2012-3-21 23:37:22

谢谢，翻了翻help，好像也只好这么做了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

蓝色

2012-3-22 08:21:49

clear

input ///
str10 group  no y1       y2       y3
A       1       52       49       42
A       2       51       50       46
A       3       50       49       41
A       4       51       49       44
A       5       49       47       40
B       1       51       54       53
B       2       49       47       46
B       3       50       47       44
B       4       49       48       41
B       5       52       50       48
end

reshape long y,i(group no) j(s)

egen g=group(group)
replace g=g-1
gen  t1=(s==2)
gen  t2=(s==3)
drop group s

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

蓝色

2012-3-22 08:22:53

clear

input ///
str10 group  no y1       y2       y3
A       1       52       49       42
A       2       51       50       46
A       3       50       49       41
A       4       51       49       44
A       5       49       47       40
B       1       51       54       53
B       2       49       47       46
B       3       50       47       44
B       4       49       48       41
B       5       52       50       48
end

reshape long y,i(group no) j(s)
egen g=group(group)
replace g=g-1
tab    s,gen(t)

drop group s