3

在一本旧的统计教科书中，我找到了一个国家人口年龄分布的表格：

        百分比
 年龄人口
------------------
 0-5 8
 5-14 18
14-18 8
18-21 5
21-25 6
25-35 12
35-45 11
45-55 11
55-65 9
65-75 6
75-85 4

我想将此分布绘制为 R 中的直方图，其中年龄范围为中断，人口百分比为密度，但似乎没有一种直接的方法来做到这一点。R 的hist()函数希望您提供单个数据点，而不是像这样的预先计算的分布。

这就是我的做法。

# Copy original textbook table into two data structures
ageRanges <- list(0:5, 5:14, 14:18, 18:21, 21:25, 25:35, 35:45, 45:55, 55:65, 65:75, 75:85)
pcPop <- c(8, 18, 8, 5, 6, 12, 11, 11, 9, 6, 4)
# Make up "fake" age data points from the distribution described by the table
ages <- lapply(1:length(ageRanges), function(i) {
    ageRange <- ageRanges[[i]]
    round(runif(pcPop[i] * 100, min=ageRange[1], max=ageRange[length(ageRange)-1]), 0)
})
ages <- unlist(ages)
# Use the endpoints of the age class intervals as breaks for the histogram
breaks <- append(0, sapply(ageRanges, function(x) x[length(x)]))
hist(ages, breaks=breaks)

似乎必须有一种不那么冗长/hacky的方式来处理它。

编辑：FWIW，这是生成的直方图的样子：

直方图

4

2 回答 2

7

这应该得到你想要的：

test <- read.table(textConnection("age popperc
0-5 8
5-14 18
14-18 8
18-21 5
21-25 6
25-35 12
35-45 11
45-55 11
55-65 9
65-75 6
75-85 4"),header=TRUE,stringsAsFactors=FALSE)

midval <- sapply(strsplit(test$age,"-"),function(x) mean(as.numeric(x)))
breakval <- strsplit(test$age,"-")
breakval <- as.numeric(c(sapply(breakval,head,1),tail(unlist(breakval),1)))
hist(rep(midval,test$popperc),breaks=breakval)

在此处输入图像描述

您还可以定义自己的class直方图对象，然后plot如果您只想绘制频率而不是密度：

# define the histogram object and plot it
histres <- list(
breaks=breakval,
counts=test$popperc,
mids=midval,
xname="ages",
equidist = TRUE
)
class(histres) <- "histogram"
plot(histres)

在此处输入图像描述

于 2013-02-06T04:47:38.253 回答
2

正如评论中所说，使用barplot. 您可以在 a 中指定宽度barplot

barplot(pcPop,  width = seq(0,85,5),space=0)
于 2013-02-06T04:19:55.543 回答