- 对下面相关系数矩阵,试用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语言计算函数
- 某厂对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