有阳性对照，正常对照及模型对照的两因素试验如何编写SAS程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我的试验为考察两种（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

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：sas程序 normal 不同水平 组间差异 Norm normal 程序 模型 如何 动物

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;

非常感谢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命令，所以不知它的用法及作用，能否知之一二？万分感谢！！