SAS第二次作业分析

第二次作业

习题7-3

3、下表数据是退火温度x对黄铜延性效应y的实验结果，y是以延长度计的。

x 300 400 500 600 700 800 y 40 50 55 60 67 70

(1)画出散点图

(2)求y对x的线性回归方程。

解：

（1）

（2）

所以y=24.6+0.06x

（1)作散点图

（2)以模型y=b0+b1x+b2x2+ε,ε~N(0,σ2)拟合数据，其中b0,b1,b2,σ2与x无关，求回归方程y=b0+b1x+b2x2

解：

（1）

所以y=19.03+1.01x-0.02x2

3、钢包容积y和使用次数x的侵蚀数据如表所示：

(1)作散点图

(2)试作变量替换，化非线性回归模型为线性回归模型并讨论回归方程的显著性。解：

（1）

（2）程序

data data123;

input x y @@;

u=log(y);v=1/x;

cards;

2 106.42

3 108.20

4 109.58

5 109.50 7

110.00 8 109.93 10 110.49

11 110.59 14 110.60 15 110.90 16 110.76

18 111.00 19 111.20

;

proc reg;

model u=v;

run;

The REG Procedure

Model: MODEL1

Dependent Variable: u

Analysis of Variance

Sum of Mean

Source DF Squares Square F Value Pr > F

Model 1 0.00174 0.00174 410.17 <.0001 Error 11 0.00004672 0.00000425

Corrected Total 12 0.00179

Root MSE 0.00206 R-Square 0.9739

Dependent Mean 4.69983 Adj R-Sq 0.9715

Coeff Var 0.04385

Parameter Estimates

Parameter Standard

Variable DF Estimate Error t Value Pr > |t|

Intercept 1 4.71408 0.00090629 5201.49 <.0001

v 1 -0.09029 0.00446 -20.25 <.0001

所以lny=4.71+0.09(-1/x)

Pr<0.0001<0.01所以回归方程高度显著。

8.1

例一.有20例肝病患者的四项肝功能指标x1（转氨酶量SGPT)，x2(肝大指数），x3(硫酸锌浊度ZnT)及x4(胎甲球AFP)的观察数据如表所示，试作这四项指标的主成分分析。

3 120 3.0 13 50

4 250 4.

5 18 0

5 120 3.5 9 50

6 10 1.5 12 50

7 40 1.0 19 40

8 270 4.0 13 60

9 280 3.5 11 60

10 170 3.0 9 60

11 180 3.5 14 40

12 130 2.0 30 50

13 220 1.5 17 20

14 160 1.5 35 60

15 220 2.5 14 30

16 140 2.0 20 20

17 220 2.0 14 10

18 40 1.0 10 0

19 20 1.0 12 60

20 120 2.0 20 0 解：

y1=0.99964z1+0.6898z2+0.0879z3+0.1628z4

y2=0.095z1-0.2836z2+0.9042z3+0.305z4

y3=-0.240049z1+0.0585z2-0.2703z3+0.9305z4

y4=-0.6659z1+0.6636z2+0.3189z3-0.1208z4

习题8-2

选拔职员对应聘人员测试6门科目：词汇，阅读，同义词，算数，代数，微积分，分别记为x1,x2,x3,x4,x5,x6,将所有应聘者的考试成绩做计算机处理，得样本相关阵，试对这六科成绩作因子分析，样本相关阵为

1 0.7

2 0.6

3 0.09 0.09 0.00

0.72 1 0.57 0.15 0.16 0.09

0.63 0.57 1 0.14 0.15 0.09

0.09 0.15 0.14 1 0.57 0.72

0.09 0.16 0.15 0.57 1 0.72

0.00 0.09 0.09 0.63 0.72 1

解：程序

data score(type=corr);

input name $ x1 x2 x3 x4 x5 x6;

cards;

x1 1 0.72 0.63 0.09 0.09 0.00

x2 0.72 1 0.57 0.15 0.16 0.09

x3 0.63 0.57 1 0.14 0.15 0.09

x4 0.09 0.15 0.14 1 0.57 0.63

x5 0.09 0.16 0.15 0.57 1 0.72

x6 0.00 0.09 0.09 0.63 0.72 1

;

proc factor data=score;

var x1 x2 x3 x4 x5 x6;

run;

结果：

因子1反映了考生的平均综合

能力，因子2反映了语文能力和

数学能力的差异。第一因子解差

的方差是2.601939，占信息量的

43.37%,是主要因子。→→→→

习题8-3

1、现有6个铅弹头，用“中子活化”方法测得7种微量元素的含量数据如下表所示，试用系统聚类法对6个弹头进行分类。

元素样品号Ag(银)

X1

Al(铝)

X2

Cu(铜）

X3

Ca(钙)

X4

Sb(锑)

X5

Bi(鉍)

X5

Sn(锡)

X7

1 0.05798 5.5150 347.10 21.910 8586 174

2 61.69

2 0.08441 3.9700 347.20 19.710 7947 2000 2440

3 0.07217 1.1530 54.58 3.052 3860 1445 9497

4 0.15010 1.7020 307.50 15.030 12290 1461 6480

5 5.7440 2.8540 229.60 9.657 8099 126

6 12520

6 0.21300 0.7058 240.30 13.910 8980 2820 4135

解：程序

data A;

input no x1-x7;

cards;

1 0.05789 5.5150 347.10 21.910 8586 174

2 61.69

2 0.08441 3.9700 347.20 19.710 7947 2000 2440

3 0.07217 1.1530 54.85 3.052 3860 1445 9497

4 0.15010 1.7020 307.50 15.030 12290 1461 6380

5 5.7440 2.8540 229.60 9.657 8099 126

6 12520

6 0.21300 0.7058 240.30 13.910 8980 2820 4135

;

proc cluster data = A method=sin outtree=out1; var x1-x7;

id no;

proc tree horizontal data=out1;

title'minimum distance';

run; 结果：

分为五类：{1}{2,6}{4}{3}{5}

分成四类：{1,2,6}{4}{3}{5}

分成三类：{1,2,6,4}{3}{5}

分成两类：{1,2,6,4,3}{5}

↓↓↓↓↓↓↓↓↓↓↓↓↓↓