Agri521'Blog分享 http://blog.sciencenet.cn/u/agri521 研究方向:作物模型、统计计算

博文

SAS与通径分析

已有 12734 次阅读 2010-7-8 14:35 |个人分类:统计计算|系统分类:科研笔记|关键词:学者| SAS, 通径分析

通径分析是数量遗传学家Sewall Wright于1921年提出来,经遗传育种学者不断改进和完善形成的一种多元统计技术。它通过对自变量和因变量之间的相关分解来研究因变量(性状)的相 对重要性,已在众多领域广泛应用。

通径系数是因变量Y关于自变量X1、X2、X3的标准的偏回归系数,即通径系数=自变量的回 归系数×(自变量的标准差/因变量的标准差)。

计算步骤:
1、相关系数计算
2、因变量对自变量的多元回归分析
3、通径系数的计 算
4、间接通径系数的计算

计算涉及的公式:
1、通径系数=自变量的回归系数×(自变量的标准差/因变量的标准差)
2、 通径系数的显著性测验与多元回归分析中对偏回归系数的显著性测验等价
3、回归方程的误差e的通径系数=sqrt(1-R2)
4、x1 通过x2的间接效应(间接通径系数)=X1与的X2相关系数*x2的通径系数
同理x2通过x1的间接效应=X1与的X2相关系数*x1的通径系 数

SAS代码:(本想写成宏,可以直接给出终计算结果,限于能力和时间,本代码只能输出计算通径系数所需的数据 集)
/*系统设置*/
options nodate nonumber nocenter linesize=80 pagesize=60 ;
/*数据读取*/
data xiaomai;
input y x1-x4;
cards;
15.7 10 23 3.6 113
14.5 9 20 3.6 106
17.5 10 22 3.7 111
22.5 13 21 3.7 109
15 5 10 22 3.6 110
16.9 10 23 3.5 103
8.6 8 23 3.3 100
17.0 10 24 3.4 114
13.7 10 20 3.4 104
13.4 10 21 3.4 110
20.3 10 23 3.9 104
10 2 8 21 3.5 109
7.4 6 23 3.2 114
11.6 8 21 3.7 113
12.3 9 22 3 6 105
run;
ods html file=’E:Doctor Degreesasdatamydata.html’ style=science;

ods output SimpleStats=SimStats(keep= variable stddev); /*导出标准差数据集*/
/*相关系数计算*/
proc corr data=xiaomai;
var y x1-x4;
run;
ods trace on ;
ods output FitStatistics=FitStat(keep=label2 cvalue2);
/*导出R2*/
ods output ParameterEstimates=PaterE( keep= variable standardizedest); /*导出标准回归系数*/
/*多元回归分析*/
proc reg data=xiaomai;
model y=x1-x4/stb; /*stb选项:标准回归参数估计*/
run;
ods trace off;

/*合并标准差和标准回归系数数据集*/
data result;
merge PaterE SimStats FitStat ;
label standardizedest=’标准回归系数’ stdDev=’标准差’ variable=’自变量’;
run;
proc print data= result;
run;
ods html close;
run;

参考文献:
1、通径分析在Excel和SPSS中的实现,张琪,农业网络信息,2007/03
2、 偏相关与通径分析的EXCEL VBA程序设计,陈庭木,农业网络信息,2007/03
3、EXCEL在通径分析中的应用,任红松,农业网络信 息,2006/03
4、Excel在通径分析中的应用,何风华,中国卫生统计,2005/05
5、通径分析在EXCEL上的实现,张 天伦,农业网络信息,2004/08



https://m.sciencenet.cn/blog-81938-342036.html

上一篇:R语言AIS软件包详解--look函数
下一篇:使用的是哪个作物生长模型?

1 万润兰

发表评论 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-29 18:31

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部