经管之家App
让优质教育人人可得
立即打开
全部版块
我的主页
›
论坛
›
数据科学与人工智能
›
数据分析与数据科学
›
SAS专版
SAS中多水平的方差分析
楼主
慧(会)幸福
5229
4
收藏
2013-09-03
SAS中多水平的方差分析,怎么做,比较三个水平的差异,多谢了
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
全部回复
沙发
妖帝东皇
2013-9-3 19:17:19
不懂,但是帮你顶贴!
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
藤椅
Eternal0601
2013-9-3 19:36:50
/*For reference.*/
proc glm data=sasuser.MGGarlic_Block ;
class Fertilizer Sector;
model BulbWt=Fertilizer Sector;
title "ANOVA for Randomized Block Design";
run;
quit;
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
板凳
cjhong111
2013-9-5 23:08:13
看了你的提问,觉得和我前段时间做了个医学统计的数据分析很像,一个例子,供你参考:
研究不同的氧浓度(A)和不同的药物(B)对人白血病细胞K562抑制效果,因素A分为A1(含氧3%)和A2(含氧20%),因素B分为B1(表阿霉素)、B2(自制中药)、B3(132Ge)、B4(B1+ B3)、B5(B1+ B2)、B6(B2+ B3)、B7(B1+ B2+ B3);进行2*7析因设计并收集到数据如下,分析A、B两因素对K562细胞抑制的效果。
B药物
A氧浓度 重复编号 B1 B2 B3 B4 B5 B6 B7
A1((含氧3%)) 1 0.31 0.46 0.29 0.49 0.72 0.45 0.19
2 0.18 0.39 0.18 0.51 0.49 0.42 0.20
3 0.12 0.40 0.12 0.62 0.55 0.44 0.18
4 0.13 0.34 0.13 0.53 0.37 0.42 0.17
A2(含氧20%) 1 0.29 0.65 0.87 0.74 1.09 1.04 0.81
2 0.27 0.84 0.39 0.78 0.73 0.63 1.01
3 0.29 0.45 0.57 1.45 0.81 1.18 1.18
4 0.28 0.63 0.64 1.41 0.77 1.45 0.94
解题思路:先对数据进行交互作用分析,若无交互作用,那么分析其主效应就可以了,若主效应有多个水平,那么可以采用means a b/snk进行各个主效应的多水平的比较就可以了;
但是如果存在交互效应的话,分析主效应其意义不大,就需要把A或B控制在不同的水平下,对B或A的各个水平之间进行两两比较。
首先对交互效应的分析:
Sas程序:
data h1;
do a=1 to 2;
do r=1 to 4;
do b=1 to 7;
input x@@;
output;
end;
end;
end;
list;
cards;
0.31 0.46 0.29 0.49 0.72 0.45 0.19
0.18 0.39 0.18 0.51 0.49 0.42 0.20
0.12 0.40 0.12 0.62 0.55 0.44 0.18
0.13 0.34 0.13 0.53 0.37 0.42 0.17
0.29 0.65 0.87 0.74 1.09 1.04 0.81
0.27 0.84 0.39 0.78 0.73 0.63 1.01
0.29 0.45 0.57 1.45 0.81 1.18 1.18
0.28 0.63 0.64 1.41 0.77 1.45 0.94
;
proc anova;
class a b;
model x=a b a*b;
means a b a*b;
run;
经过数据分析,发现存在交互效应,那么采用contrast语句分析一个因素的各个水平在另外一个因素的某一水平下的两两比较;
程序如下:
data h1;
do a=1 to 2;
do r=1 to 4;
do b=1 to 7;
input x@@;
output;
end;
end;
end;
list;
cards;
0.31 0.46 0.29 0.49 0.72 0.45 0.19
0.18 0.39 0.18 0.51 0.49 0.42 0.20
0.12 0.40 0.12 0.62 0.55 0.44 0.18
0.13 0.34 0.13 0.53 0.37 0.42 0.17
0.29 0.65 0.87 0.74 1.09 1.04 0.81
0.27 0.84 0.39 0.78 0.73 0.63 1.01
0.29 0.45 0.57 1.45 0.81 1.18 1.18
0.28 0.63 0.64 1.41 0.77 1.45 0.94
;
proc glm;
class a b;
model x=a b a*b;
contrast'(a1- a2)/b1'
a 1 -1 b*a 1 0 0 0 0 0 0 -1 0 0 0 0 0 0;
contrast'(a1 vs a2)/b2'
a 1 -1 b*a 0 1 0 0 0 0 0 0 -1 0 0 0 0 0;
contrast'(a1 vs a2)/b3'
a 1 -1 b*a 0 0 1 0 0 0 0 0 0 -1 0 0 0 0;
contrast'(a1 vs a2)/b4'
a 1 -1 b*a 0 0 0 1 0 0 0 0 0 0 -1 0 0 0;
contrast'(a1 vs a2)/b5'
a 1 -1 b*a 0 0 0 0 1 0 0 0 0 0 0 -1 0 0;
contrast'(a1 vs a2)/b6'
a 1 -1 b*a 0 0 0 0 0 1 0 0 0 0 0 0 -1 0;
contrast'(a1 vs a2)/b7'
a 1 -1 b*a 0 0 0 0 0 0 1 0 0 0 0 0 0 -1;
contrast'(b1-b7)/a1'
b 1 -1 0 0 0 0 0 a*b 1 -1 0 0 0 0 0 0 0 0 0 0 0 0,
b 1 0 -1 0 0 0 0 a*b 1 0 -1 0 0 0 0 0 0 0 0 0 0 0,
b 1 0 0 -1 0 0 0 a*b 1 0 0 -1 0 0 0 0 0 0 0 0 0 0,
b 1 0 0 0 -1 0 0 a*b 1 0 0 0 -1 0 0 0 0 0 0 0 0 0,
b 1 0 0 0 0 -1 0 a*b 1 0 0 0 0 -1 0 0 0 0 0 0 0 0,
b 1 0 0 0 0 0 -1 a*b 1 0 0 0 0 0 -1 0 0 0 0 0 0 0;
contrast'(b1-b7)/a2'
b 1 -1 0 0 0 0 0 a*b 0 0 0 0 0 0 0 1 -1 0 0 0 0 0,
b 1 0 -1 0 0 0 0 a*b 0 0 0 0 0 0 0 1 0 -1 0 0 0 0,
b 1 0 0 -1 0 0 0 a*b 0 0 0 0 0 0 0 1 0 0 -1 0 0 0,
b 1 0 0 0 -1 0 0 a*b 0 0 0 0 0 0 0 1 0 0 0 -1 0 0,
b 1 0 0 0 0 -1 0 a*b 0 0 0 0 0 0 0 1 0 0 0 0 -1 0,
b 1 0 0 0 0 0 -1 a*b 0 0 0 0 0 0 0 1 0 0 0 0 0 -1;
contrast'(b1-b2)/a1'
b 1 -1 0 0 0 0 0 a*b 1 -1 0 0 0 0 0;
contrast'(b1-b3)/a1'
b 1 0 -1 0 0 0 0 a*b 1 0 -1 0 0 0 0;
contrast'(b1-b4)/a1'
b 1 0 0 -1 0 0 0 a*b 1 0 0 -1 0 0 0;
contrast'(b1-b5)/a1'
b 1 0 0 0 -1 0 0 a*b 1 0 0 0 -1 0 0;
contrast'(b1-b6)/a1'
b 1 0 0 0 0 -1 0 a*b 1 0 0 0 0 -1 0;
contrast'(b1-b7)/a1'
b 1 0 0 0 0 0 -1 a*b 1 0 0 0 0 0 -1;
contrast'(b2-b3)/a1'
b 0 1 -1 0 0 0 0 a*b 0 1 -1 0 0 0 0;
contrast'(b2-b4)/a1'
b 0 1 0 -1 0 0 0 a*b 0 1 0 -1 0 0 0;
contrast'(b2-b5)/a1'
b 0 1 0 0 -1 0 0 a*b 0 1 0 0 -1 0 0;
contrast'(b2-b6)/a1'
b 0 1 0 0 0 -1 0 a*b 0 1 0 0 0 -1 0;
contrast'(b2-b7)/a1'
b 0 1 0 0 0 0 -1 a*b 0 1 0 0 0 0 -1;
contrast'(b3-b4)/a1'
b 0 0 1 -1 0 0 0 a*b 0 0 1 -1 0 0 0;
contrast'(b3-b5)/a1'
b 0 0 1 0 -1 0 0 a*b 0 0 1 0 -1 0 0;
contrast'(b3-b6)/a1'
b 0 0 1 0 0 -1 0 a*b 0 0 1 0 0 -1 0;
contrast'(b3-b7)/a1'
b 0 0 1 0 0 0 -1 a*b 0 0 1 0 0 0 -1;
contrast'(b4-b5)/a1'
b 0 0 0 1 -1 0 0 a*b 0 0 0 1 -1 0 0;
contrast'(b4-b6)/a1'
b 0 0 0 1 0 -1 0 a*b 0 0 0 1 0 -1 0;
contrast'(b4-b7)/a1'
b 0 0 0 1 0 0 -1 a*b 0 0 0 1 0 0 -1;
contrast'(b5-b6)/a1'
b 0 0 0 0 1 -1 0 a*b 0 0 0 0 1 -1 0;
contrast'(b5-b7)/a1'
b 0 0 0 0 1 0 -1 a*b 0 0 0 0 1 0 -1;
contrast'(b6-b7)/a1'
b 0 0 0 0 0 1 -1 a*b 0 0 0 0 0 1 -1;
contrast'(b1-b2)/a2'
b 1 -1 0 0 0 0 0 a*b 0 0 0 0 0 0 0 1 -1 0 0 0 0 0;
contrast'(b1-b3)/a2'
b 1 0 -1 0 0 0 0 a*b 0 0 0 0 0 0 0 1 0 -1 0 0 0 0;
contrast'(b1-b4)/a2'
b 1 0 0 -1 0 0 0 a*b 0 0 0 0 0 0 0 1 0 0 -1 0 0 0;
contrast'(b1-b5)/a2'
b 1 0 0 0 -1 0 0 a*b 0 0 0 0 0 0 0 1 0 0 0 -1 0 0;
contrast'(b1-b6)/a2'
b 1 0 0 0 0 -1 0 a*b 0 0 0 0 0 0 0 1 0 0 0 0 -1 0;
contrast'(b1-b7)/a2'
b 1 0 0 0 0 0 -1 a*b 0 0 0 0 0 0 0 1 0 0 0 0 0 -1;
contrast'(b2-b3)/a2'
b 0 1 -1 0 0 0 0 a*b 0 0 0 0 0 0 0 0 1 -1 0 0 0 0;
contrast'(b2-b4)/a2'
b 0 1 0 -1 0 0 0 a*b 0 0 0 0 0 0 0 0 1 0 -1 0 0 0;
contrast'(b2-b5)/a2'
b 0 1 0 0 -1 0 0 a*b 0 0 0 0 0 0 0 0 1 0 0 -1 0 0;
contrast'(b2-b6)/a2'
b 0 1 0 0 0 -1 0 a*b 0 0 0 0 0 0 0 0 1 0 0 0 -1 0;
contrast'(b2-b7)/a2'
b 0 1 0 0 0 0 -1 a*b 0 0 0 0 0 0 0 0 1 0 0 0 0 -1;
contrast'(b3-b4)/a2'
b 0 0 1 -1 0 0 0 a*b 0 0 0 0 0 0 0 0 0 1 -1 0 0 0;
contrast'(b3-b5)/a2'
b 0 0 1 0 -1 0 0 a*b 0 0 0 0 0 0 0 0 0 1 0 -1 0 0;
contrast'(b3-b6)/a2'
b 0 0 1 0 0 -1 0 a*b 0 0 0 0 0 0 0 0 0 1 0 0 -1 0;
contrast'(b3-b7)/a2'
b 0 0 1 0 0 0 -1 a*b 0 0 0 0 0 0 0 0 0 1 0 0 0 -1;
contrast'(b4-b5)/a2'
b 0 0 0 1 -1 0 0 a*b 0 0 0 0 0 0 0 0 0 0 1 -1 0 0;
contrast'(b4-b6)/a2'
b 0 0 0 1 0 -1 0 a*b 0 0 0 0 0 0 0 0 0 0 1 0 -1 0;
contrast'(b4-b7)/a2'
b 0 0 0 1 0 0 -1 a*b 0 0 0 0 0 0 0 0 0 0 1 0 0 -1;
contrast'(b5-b6)/a2'
b 0 0 0 0 1 -1 0 a*b 0 0 0 0 0 0 0 0 0 0 0 1 -1 0;
contrast'(b5-b7)/a2'
b 0 0 0 0 1 0 -1 a*b 0 0 0 0 0 0 0 0 0 0 0 1 0 -1;
contrast'(b6-b7)/a2'
b 0 0 0 0 0 1 -1 a*b 0 0 0 0 0 0 0 0 0 0 0 0 1 -1;
run;
注意:若a和b之间没有交互作用的话,那么可以分别比较a和b因素下的各个水平之间的差异,程序如下:
data h1;
do a=1 to 2;
do r=1 to 4;
do b=1 to 7;
input x@@;
output;
end;
end;
end;
list;
cards;
0.31 0.46 0.29 0.49 0.72 0.45 0.19
0.18 0.39 0.18 0.51 0.49 0.42 0.20
0.12 0.40 0.12 0.62 0.55 0.44 0.18
0.13 0.34 0.13 0.53 0.37 0.42 0.17
0.29 0.65 0.87 0.74 1.09 1.04 0.81
0.27 0.84 0.39 0.78 0.73 0.63 1.01
0.29 0.45 0.57 1.45 0.81 1.18 1.18
0.28 0.63 0.64 1.41 0.77 1.45 0.94
;
proc glm;
class a b;
model x=a b a*b;
lsmeans a b a*b;
contrast'(a1-a2)'
a 1 -1;
contrast'(b1-b2)'
b 1 -1 0 0 0 0 0;
contrast'(b1-b3)'
b 1 0 -1 0 0 0 0;
contrast'(b1-b4)'
b 1 0 0 -1 0 0 0;
contrast'(b1-b5)'
b 1 0 0 0 -1 0 0;
contrast'(b1-b6)'
b 1 0 0 0 0 -1 0;
contrast'(b1-b7)'
b 1 0 0 0 0 0 -1;
contrast'(b2-b3)'
b 0 1 -1 0 0 0 0;
contrast'(b2-b4)'
b 0 1 0 -1 0 0 0;
contrast'(b2-b5)'
b 0 1 0 0 -1 0 0;
contrast'(b2-b6)'
b 0 1 0 0 0 -1 0;
contrast'(b2-b7)'
b 0 1 0 0 0 0 -1;
contrast'(b3-b4)'
b 0 0 1 -1 0 0 0;
contrast'(b3-b5)'
b 0 0 1 0 -1 0 0;
contrast'(b3-b6)'
b 0 0 1 0 0 -1 0;
contrast'(b3-b7)'
b 0 0 1 0 0 0 -1;
contrast'(b4-b5)'
b 0 0 0 1 -1 0 0;
contrast'(b4-b6)'
b 0 0 0 1 0 -1 0;
contrast'(b4-b7)'
b 0 0 0 1 0 0 -1;
contrast'(b5-b6)'
b 0 0 0 0 1 -1 0;
contrast'(b5-b7)'
b 0 0 0 0 1 0 -1;
contrast'(b6-b7)'
b 0 0 0 0 0 1 -1;
run;
但是显然a和b是有交互作用的,这里写的这种情况是没有交互作用的。
程序看起来很吓人,其实也就是个体力活。
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
报纸
cjhong111
2013-9-5 23:11:45
刚才又看了遍楼主的提问,怀疑自己是不是把问题搞复杂了。
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
相关推荐
双因素方差分析,样本量为10可以做吗?怎么做?
求助方差分析
求助方差分析
怎么做直观分析 方差分析
怎么做双因素方差分析
方差分析,单因变量多因素,多水平
如下 题 新手sas做方差分析怎么做,
4乘4的数据 每个中间加入性别因素 怎么做方差分析
二分变量方差分析
怎么做 多组的方差分析
栏目导航
SAS专版
真实世界经济学(含财经时事)
求助成功区
休闲灌水
市场行情分析
管理信息系统
热门文章
中国金融生成式AI多模态内容鉴伪与安全防御 ...
全球能源转型展望2025—全球和区域预测至20 ...
海外资管机构赴上海投资指南(2025版)
understanding climate change perceptions ...
【全美经典】离散数学
全球能源转型展望2025—全球和区域预测至20 ...
世界机器人2025年报告 World Robotics 2025
甲子光年_2025甲子Cool Vendor人形机器人大 ...
AOM:The Boundaries of Trust in a New Era
气象学-山东大学
推荐文章
AI狂潮席卷学术圈,不会编程也能打造专属智 ...
10月重磅来袭|《打造Coze/Dify专属学术智能 ...
最快1年拿证,学费不足5W!热门美国人工智能 ...
关于如何利用文献的若干建议
关于学术研究和论文发表的一些建议
关于科研中如何学习基础知识的一些建议 (一 ...
一个自编的经济学建模小案例 --写给授课本科 ...
AI智能体赋能教学改革: 全国AI教育教学应用 ...
2025中国AIoT产业全景图谱报告-406页
关于文献求助的一些建议
说点什么
分享
微信
QQ空间
QQ
微博
扫码加好友,拉您进群
各岗位、行业、专业交流群