我试图获取箱线图的值,但是当我使用波浪号公式时,我无法使其工作。
我觉得这只是一个简单的问题,但是我用错误的方法来找到答案......
以下是我的尝试:
numbers <- c(6, 1, 1, 1, 3, 3, 2, 3, 13, 13, 4, 7, 38, 1, 12, 9, 4, 65, 9, 23)
groups <- c("A", "B", "D", "E", "F", "G", "B", "C", "D", "E", "F", "A", "B", "C", "D", "E", "F", "A", "B", "D")
boxplot(numbers~groups)
summary(numbers~groups)
我只是得到一个输出
Length Class Mode
3 formula call
如何获取boxplot的值(如Min。,1st.Qu.Median等)
输出应该如下:
A B (...)
Min. : ? Min. : ? (...)
1st Qu.: ? 1st Qu.: ? (...)
Median : ? Median : ? (...)
Mean : ? Mean : ? (...)
3rd Qu.: ? 3rd Qu.: ? (...)
Max. : ? Max. : ? (...)
答案 0 :(得分:2)
尝试使用聚合函数
aggregate(numbers,list(groups),summary)
Group.1 x.Min. x.1st Qu. x.Median x.Mean x.3rd Qu. x.Max.
1 A 6.000 6.500 7.000 26.000 36.000 65.000
2 B 1.000 1.750 5.500 12.500 16.250 38.000
3 C 1.000 1.500 2.000 2.000 2.500 3.000
4 D 1.000 9.250 12.500 12.250 15.500 23.000
5 E 1.000 5.000 9.000 7.667 11.000 13.000
6 F 3.000 3.500 4.000 3.667 4.000 4.000
7 G 3.000 3.000 3.000 3.000 3.000 3.000
答案 1 :(得分:2)
数据在boxplot函数中。如果将其保存为对象,则它包含矩阵中的所有统计信息。
bp <- boxplot(numbers~groups)
bp$stats
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 6.0 1.0 1 1.0 1 3.0 3
[2,] 6.5 1.5 1 6.5 5 3.5 3
[3,] 7.0 5.5 2 12.5 9 4.0 3
[4,] 36.0 23.5 3 18.0 11 4.0 3
[5,] 65.0 38.0 3 23.0 13 4.0 3
在boxplot的帮助下,您可以看到它包含一个组件列表,其中一个组件包含您需要的信息。
stats一个矩阵,每列包含较低的极值 晶须,下铰链,中间,上铰链和极端 一组/小区的上晶须。如果所有输入都有 同一个类属性,这个组件也是如此。
相同的统计数据由函数fivenum
计算答案 2 :(得分:1)
我实际上不知道你是否可以使用summary
命令中的波浪号运算符来执行此操作。
但也有其他选择,例如tapply
:
tapply(numbers,groups,summary)
或describeBy
包中的Psych
函数:
# install.packages("psych")
library("psych")
describeBy(numbers,groups,mat=TRUE)
输出:
item group1 vars n mean sd median trimmed mad min max range skew
11 1 A 1 3 26.000000 33.7786915 7.0 26.000000 1.4826 6 65 59 0.38452071
12 2 B 1 4 12.500000 17.3685540 5.5 12.500000 5.9304 1 38 37 0.66132091
13 3 C 1 2 2.000000 1.4142136 2.0 2.000000 1.4826 1 3 2 0.00000000
14 4 D 1 4 12.250000 8.9953692 12.5 12.250000 8.1543 1 23 22 -0.06221018
15 5 E 1 3 7.666667 6.1101009 9.0 7.666667 5.9304 1 13 12 -0.20782656
16 6 F 1 3 3.666667 0.5773503 4.0 3.666667 0.0000 3 4 1 -0.38490018
17 7 G 1 1 3.000000 NA 3.0 3.000000 0.0000 3 3 0 NA
kurtosis se
11 -2.333333 19.5021366
12 -1.756576 8.6842770
13 -2.750000 1.0000000
14 -1.878463 4.4976846
15 -2.333333 3.5276684
16 -2.333333 0.3333333
17 NA NA