|
data a;input x1-xp y1-yp @@;cards;原始数据集数据;
proc standard mean=0 std= 1 out=b;
proc print;
proc pls data= a nfac=2 cv= one cvtest varss censcale details;
model yl-yp=x1-xp/solution;run;
其中:程序standard对数据进行标准化;程序print 要求打印输出标准化数据;程序PLS 进行偏最小二乘回归分析,其语句说明如下: 本例nfac=2表示要求提取成分为2个,可根据实际情况任选成分个数如1,2,…,缺省时提取成分个数与自变量个数相同; cv = one 采用"舍一交叉验证法"估计提取成分数的残差平方和,如果把one 换作split、block或random,表示分别采用分裂样本交叉验证法、分批交叉验证法、随机样本交叉验证法,且采用cv 语句所输出的偏最小二乘回归方程仅包含第一对成分,此时nfac 规定的提取成分数无效; cvtest 则是检验因为提取成分不同导致的残差平方和之间差异的显著性;varss输出提取成分对总的自变量与因变量变差及对各个原始变量变差的解释程度,缺省时只给出所提取成分总的自变量与因变量变差的解释比例; censcale 输出标准化变换时的均数和标准差;details 输出模型效应的载荷量、模型效应权重和因变量权重;solution 输出偏最小二乘回归方程原始变量和标准化变量的参数估计。
|