多元统计分析及R语言建模（王斌会）第二章答案

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

>R=matrix(c(1,0.8,0.26,0.67,0.34,0.8,1,0.33,0.59,0.34,0.26,0.33,1,0.37,0.21,0.67,0.59,0.37,1,0.35,0.34,0.34,0.21,0.35,1),nrow = 5,ncol = 5);

#矩阵求逆运用函数solve()

> solve(R)

#求矩阵的特征值和特征向量运用函数eigen()

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

>R.e $ vectors%*%diag(R.e $ values)%*%t(R.e $ vectors)

在RStudio运行代码如下

图1.用R语言计算函数

1. 某厂对50个计件工人某月份工资进行登记，获得以下原始资料(单位:元):

1465,1405,1355,1225,1000,1760,1755,1710,1605,1535，
1985,1965,1910,1845,1810,2270,2240,2190,2040,2010,
2980,2820,2600,2430,2290,1375,1295,1265,1175,1125,
1735,1645,1625,1595,1575,1940,1880,1865,1835,1815,
2220,2110,2095,2030,2030,2670,2550,2520,2370,2320,
试按组距为300编制频数表,计算频数、频率和累积频率，并绘制直方图。
(1)写出R语言程序;

 

> library(openxlsx)

> d2.2 = read.xlsx('E:/mvexer5.xlsx',sheet ='E2.2')

> head(d2.2) 

> hist(d2.2 $ X,breaks = seq(0,3000,by = 300),col = 1:7,xlab = "工资(元)",ylab = "频数")

  

图2-1.频数直方图

> m <- seq(0,3000,by = 300)

> hist(d2.2 $ X ,m,freq = F,col = 1:7,xlab = "工资(元)",ylab = "频率")

图2-2.频率直方图

> Cumsum <- cumsum(d2.2 $ X)

> M <- seq(0,96000,by = 3000)

> hist(Cumsum,M,freq = F,col = 1:12,las = 3,xlab = "工资(元)",ylab = "累积频率")