R语言：作业五（数据处理均值、方差、标准差等；绘制数据图）
（gdut 3119007134 卢佳康）
3.1 :
首先先在桌面新建一个txt文档，将以下数据复制进去：
74.3 78.8 68.8 78.0 70.4 80.5 80.5 69.7 71.2 73.5
79.5 75.6 75.0 78.8 72.0 72.0 72.0 74.3 71.2 72.0
75.0 73.5 78.8 74.3 75.8 65.0 74.3 71.2 69.7 68.0
73.5 75.0 72.0 64.3 75.8 80.3 69.7 74.3 73.5 73.5
75.8 75.8 68.8 76.5 70.4 71.2 81.2 75.0 70.4 68.0
70.4 72.0 76.5 74.3 76.5 77.6 67.3 72.0 75.0 74.3
73.5 79.5 73.5 74.7 65.0 76.5 81.6 75.4 72.7 72.7
67.2 76.5 72.7 70.4 77.2 68.8 67.3 67.3 67.3 72.7
75.8 73.5 75.0 73.5 73.5 73.5 72.7 81.6 70.3 74.3
73.5 79.5 70.4 76.5 72.7 77.2 84.3 75.0 76.5 70.4
然后在R语言输入：("/Users/fxalll/Desktop/Data.txt"根据自己的电脑更改)
> x <- scan("/Users/fxalll/Desktop/Data.txt");x
Read 100 items
[1] 74.3 78.8 68.8 78.0 70.4 80.5 80.5 69.7 71.2 73.5 79.5 75.6 75.0 78.8 72.0
[16] 72.0 72.0 74.3 71.2 72.0 75.0 73.5 78.8 74.3 75.8 65.0 74.3 71.2 69.7 68.0
[31] 73.5 75.0 72.0 64.3 75.8 80.3 69.7 74.3 73.5 73.5 75.8 75.8 68.8 76.5 70.4
[46] 71.2 81.2 75.0 70.4 68.0 70.4 72.0 76.5 74.3 76.5 77.6 67.3 72.0 75.0 74.3
[61] 73.5 79.5 73.5 74.7 65.0 76.5 81.6 75.4 72.7 72.7 67.2 76.5 72.7 70.4 77.2
[76] 68.8 67.3 67.3 67.3 72.7 75.8 73.5 75.0 73.5 73.5 73.5 72.7 81.6 70.3 74.3
[91] 73.5 79.5 70.4 76.5 72.7 77.2 84.3 75.0 76.5 70.4
- 1、求均值：
> meanX <- mean(x);meanX
[1] 73.668
- 2、求方差：
> varX <- var(x);varX
[1] 15.51513
- 3、求标准差：
> sdX <- sd(x);sdX
[1] 3.938925
- 4、求极差：
方法一：
> diffX <- max(x) - min(x);diffX
[1] 20
方法二：
> diffX <- diff(range(x));diffX
[1] 20
- 5、标准误：
> sqX <- sdX/sqrt(length(x));sqX
[1] 0.3938925
- 6、变异系数：
> changeX <- length(x) * sdX / meanX;changeX
[1] 5.34686
- 7、偏度：
> lenX <- length(x)
> distX <- lenX / ((lenX - 1) * (lenX - 2)) * sum((x - meanX)^3) / sdX^3;distX
[1] 0.0540593
- 8、峰度：
> hillX <- ((lenX * (lenX + 1)) / (lenX - 1) * (lenX - 2) * (lenX - 3) * sum((x - meanX)^4 / sdX^4 - (3 * (lenX - 1)^2) / (lenX - 2) * (lenX - 3)));hillX
[1] -2.822135e+12
- 总结输出：
> data.frame('均值'=meanX,'方差'=varX,'标准差'=sdX,'极差'=diffX,'标准误'=sqX,'变异系数'=changeX,'偏度'=distX,'峰度'=hillX)
3.3 :
- 茎叶图：
> stemX <- stem(x)
The decimal point is at the |
64 | 300
66 | 23333
68 | 00888777
70 | 344444442222
72 | 0000000777777555555555555
74 | 33333333700000004688888
76 | 5555555226
78 | 0888555
80 | 355266
82 |
84 | 3
- 箱线图：
> boxplotX <- boxplot(x)
- 五数总括：
何为五数总括？即为中位数、下四分位数、上四分位数、最小值、最大值。
- 中位数：
> medX <- median(x);medX
[1] 73.5
- 上、下四分位数（即0～100%中处于75%与25%的数值）：
> quantileX <- quantile(x);quantileX
0% 25% 50% 75% 100%
64.3 71.2 73.5 75.8 84.3
上四分位数：75.8
下四分位数：71.2
- 最大值：
> maxX <- max(x);maxX
[1] 84.3
- 最小值：
> minX <- min(x);minX
[1] 64.3
由此五个数据全部出来，但是，其实有一个更简便的命令。
- 五数总括:
> fiveNumX <- fivenum(x);fiveNumX
[1] 64.3 71.2 73.5 75.8 84.3
从左至右分别为：最小值、下四分位数、中位数、上四分位数、最大值。