我试图在直方图上获得最大频率。我有一个值列表。然后,我执行以下操作:
hist(list, breaks=length(list), freq=TRUE)
它会自动生成 x 和 y 轴的范围。y 轴是频率,x 轴是列表中的值。
那么,我怎样才能找到将显示在此图表上的最大频率?
我正在尝试在图表的右上角制作图例,因此我需要获得最大频率值。或者有没有办法告诉 R 在图表的右上角放置一个图例框?
直方图的值可以存储为 R 中的数据框。以 OP 的示例数据框“列表”为例,您可以:
list_histo <- hist(list, breaks=length(list), freq=TRUE)
只需键入
list_histo
回到 R 将显示新的“元”数据框,其中包含有关直方图的信息(此处显示的数据是任意的,仅用于说明目的):
$breaks
[1] 0.40 0.42 0.44 0.46 0.48 0.50 0.52 0.54 0.56 0.58 0.60 0.62 0.64 0.66 0.68
[16] 0.70 0.72 0.74 0.76
$counts
[1] 1 15 112 878 4734 17995 51094 110146 178855 216454
[11] 194536 130591 64218 23017 6117 1070 144 23
$intensities
[1] 0.00005 0.00075 0.00560 0.04390 0.23670 0.89975 2.55470 5.50730
[9] 8.94275 10.82270 9.72680 6.52955 3.21090 1.15085 0.30585 0.05350
[17] 0.00720 0.00115
$density
[1] 0.00005 0.00075 0.00560 0.04390 0.23670 0.89975 2.55470 5.50730
[9] 8.94275 10.82270 9.72680 6.52955 3.21090 1.15085 0.30585 0.05350
[17] 0.00720 0.00115
$mids
[1] 0.41 0.43 0.45 0.47 0.49 0.51 0.53 0.55 0.57 0.59 0.61 0.63 0.65 0.67 0.69
[16] 0.71 0.73 0.75
$xname
[1] "list_histo"
$equidist
[1] TRUE
attr(,"class")
[1] "histogram"
调用最大值现在很简单——只需使用
max(list_histo$counts)
将返回最大值。
set.seed(100)
x = rnorm(100, mean = 5, sd = 2)
res = hist(x)
res$mids[which.max(res$counts)]
[1] 4.5
根据休息时间,条的宽度会发生变化,但中点将为您提供您正在寻找的条的中点。这将找到具有最高计数(最大频率)的条形的中点
而不是 x, y 参数到图例,您可以使用legend('topright',...)
你也可以使用table(list)
它将返回值列表和它们重复的次数:
> list<-c(0.2, 0.6, 0.4, 0.5, 0.1, 0.5, 0.6, 0.6, 0.6, 0.1, 0.1, 0.6, 0.6, 0.6, 0.6)
> table(list)
list
0.1 0.2 0.4 0.5 0.6
3 1 1 2 8
> max(table(list))
[1] 8