多元统计分析及R语言建模答案(王斌会).pdfVIP

《多元统计分析及R语言建模》第2章

王斌会

2020.2.1

rm(list=ls())#清理内存

options(digits=4)#输出结位数

par(mr=c(4,4,2,1))#设置图片输出位置

library(openxlsx)

library(knitr)

2.1

对下面的相关系数矩阵，试用R语言求其逆矩阵、特征根和特征向量。要求

写出R语言计算函数。

1.000.800.260.670.34

0.801.000.330.590.34

=0.260.331.000.370.21

0.670.590.371.000.35

[0.340.340.210.351.00]

R=matrix(c(1.00,0.80,0.26,0.67,0.34,0.80,1.00,0.33,0.59,0.34,0.26,0.33,

1.00,0.37,0.21,0.67,0.59,0.37,1.00,0.35,0.34,0.34,0.21,0.35,1.00),nrow=

5,ncol=5);R#生成矩阵R

[,1][,2][,3][,4][,5]

[1,]1.000.800.260.670.34

[2,]0.801.000.330.590.34

[3,]0.260.331.000.370.21

[4,]0.670.590.371.000.35

[5,]0.340.340.210.351.00

R.=solve(R);R.

[,1][,2][,3][,4][,5]

[1,]3.3881-2.12220.23706-1.0685-0.10623

[2,]-2.12222.9421-0.33593-0.1331-0.16164

[3,]0.2371-0.33591.20699-0.3764-0.08812

[4,]-1.0685-0.1331-0.376372.0091-0.21562

[5,]-0.1062-0.1616-0.08812-0.21561.18505

R.e=eigen(R,symmetric=T);R.e

eigen()decomposition

$vlues

[1]2.79230.82630.77910.42060.1818

$vectors

[,1][,2][,3][,4][,5]

[1,]-0.52550.34022-0.16650.159380.74494

[2,]-0.51870.23435-0.17780.50823-0.62142

[3,]-0.3131-0.90308-0.22870.149430.10844

[4,]-0.49660.03869-0.1186-0.83116-0.21673

[5,]-0.3318-0.110840.93500.056160.01355

2.2

某厂对50个计件工人某月份工资进行登记，获得以下原始资料（单位：元）。

试按组距为300编制频数表，计算频数、频率和累计频率，并绘制直方图。

14651760198522702980137517351940222026701405175519652240

28201295164518802110255013551710191021902600126516251865

20952520122516051845204024301175183520302370100015351810

2010229011251575181520302320

2.2.1

写出R语言程序。

d2.2=read.xlsx(mvexer5.xlsx,sheet=E2.2);tail(d2.2)#读取数并显示

后六行

X