全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
1798 2
2011-09-28
我的试验为考察两种(A与B)药物(两药物间有互作)治疗动物腹泻时血液指标的影响,每种药物都设三个水平,试验中另设正常组,腹泻组及药物治疗组(对治疗腹泻有效果的药物),现在我想比较A药物不同水平与正常组,腹泻组及药物治疗组间的差异,以及B药物不同水平与正常组,腹泻组及药物治疗组间差异,请教各位高手SAS程序应该怎么编?谢谢!数据如下:

NO

LHC对照

diarrhea对照

normal对照

2.60

6.61

2.64

LHC

3.84

7.66

1.51

LHC

2.71

6.77

1.71

LHC

4.62

6.16

1.13

LHC

3.15

6.09

1.91

B药物

NO

100

300

500

A药物

100

5.89

3.14

7.16

100

3.82

2.14

4.47

100

4.09

2.01

4.94

100

3.52

2.24

3.47

100

4.68

3.35

6.72

300

13.92

9.48

4.09

300

10.80

4.95

3.12

300

10.20

7.04

4.43

300

14.64

9.82

5.66

300

10.49

17.18

5.02

500

3.47

2.73

3.10

500

4.04

3.88

2.15

500

4.43

4.38

4.86

500

4.21

3.71

2.20

500

4.91

4.57

2.63

二维码

扫码加我 拉你入群

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

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

全部回复
2011-9-29 07:37:57
data one;
input group $ a b value;
cards;
normal 100 500 7.16
normal 100 500 4.47
normal 100 500 4.94
normal 100 500 3.47
normal 100 500 6.72
normal 300 500 4.09
normal 300 500 3.12
normal 300 500 4.43
normal 300 500 5.66
normal 300 500 5.02
normal 500 500 3.10
normal 500 500 2.15
normal 500 500 4.86
normal 500 500 2.20
normal 500 500 2.63
normal 0   0   2.64
normal 0   0   1.51
normal 0   0   1.71
normal 0   0   1.13
normal 0   0   1.91
diarrh 100 300 3.14
diarrh 100 300 2.14
diarrh 100 300 2.01
diarrh 100 300 2.24
diarrh 100 300 3.35
diarrh 300 300 9.48
diarrh 300 300 4.95
diarrh 300 300 7.04
diarrh 300 300 9.82
diarrh 300 300 17.18
diarrh 500 300 2.73
diarrh 500 300 3.88
diarrh 500 300 4.38
diarrh 500 300 3.71
diarrh 500 300 4.57
diarrh 0   0   6.61
diarrh 0   0   7.66
diarrh 0   0   6.77
diarrh 0   0   6.16
diarrh 0   0   6.09
lhccon 100 100 5.89
lhccon 100 100 3.82
lhccon 100 100 4.09
lhccon 100 100 3.52
lhccon 100 100 4.68
lhccon 300 100 13.92
lhccon 300 100 10.80
lhccon 300 100 10.20
lhccon 300 100 14.64
lhccon 300 100 10.49
lhccon 500 100 3.47
lhccon 500 100 4.04
lhccon 500 100 4.43
lhccon 500 100 4.21
lhccon 500 100 4.91
lhccon 0   0   2.60
lhccon 0   0   3.84
lhccon 0   0   2.71
lhccon 0   0   4.62
lhccon 0   0   3.15
;
run;

proc glm data=one;
class a b;
model value=a b;
where group="normal";
means a/tukey cldiff;
run;
proc glm data=one;
class a b;
model value=a b;
where group="diarrh";
means a/tukey cldiff;
run;
proc glm data=one;
class a b;
model value=a b;
where group="lhccon";
means a/tukey cldiff;
run;

proc glm data=one;
class a b;
model value=a b;
where group="normal";
means b/tukey cldiff;
run;
proc glm data=one;
class a b;
model value=a b;
where group="diarrh";
means b/tukey cldiff;
run;
proc glm data=one;
class a b;
model value=a b;
where group="lhccon";
means b/tukey cldiff;
run;
二维码

扫码加我 拉你入群

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

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

2011-9-30 11:33:13
非常感谢lchw001,给我很大帮助,但程序中有的地方我还是不太明白,能否再帮我解释下?一是从列出的数据来看正常组(normal)中只有b一个水平(500)与a药物三个水平(100,300,500)的组合,在lhccon和diarrh中也是如此。不知是不是我没把试验设计与过程说清造成的。我试验时是先将所有动物分成两组,一组正常饲喂,另一组通过注射药物造模,然后将造模后的动物分成11组,一组不采取治疗(即diarrh组),一组用已知有治疗效果的药物治疗(Ihccon组),其余9组按A药物与B药物分别为3水平的3*3因子设计进行。最后想了解a药物、b药物总体及不同水平在治疗效果上与正常组,ihccon及diarrh组的差异。二是从程序数据部分来看,好象跟我的试验设计有点不符。三是过程步中where语句有什么作用,我以前一直没看到过glm中可以用where命令,所以不知它的用法及作用,能否知之一二?万分感谢!!
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群