R 中同一图中 2 个因子的值的绘图频率

Question

我想绘制针对 2 个因子水平编码的变量颜色的频率，例如蓝色条应该是水平 A 的历史记录，绿色条应该是水平 B 的历史记录，两者都在同一个图表中？这可以用 hist 命令实现吗？ hist 的帮助不允许有任何因素。还有其他办法吗？

我设法通过 barplots 手动完成此操作，但我想问是否有更自动的方法

非常感谢 欧共体

PS。我不需要密度图

Answer 1

这是很有可能的。

我没有可以使用的数据，但这里有一个带有不同颜色条的直方图示例。从这里开始，您需要使用我的代码并弄清楚如何使其适用于因子而不是尾部。

基本设置
直方图 <- hist(scale(vector))、breaks= 、plot=FALSE)
绘图（直方图，col=ifelse（abs（直方图$breaks）<#of SD，颜色1，颜色2））

#EXAMPLE
x<-rnorm(1000)
histogram <- hist(scale(x), breaks=20 , plot=FALSE)
plot(histogram, col=ifelse(abs(histogram$breaks) < 2, "red", "green"))

Answer 2

万一其他人没有回答，这是一种令人满意的方式。我最近不得不处理堆积直方图，这就是我所做的：

data_sub <- subset(data, data$V1 == "Yes") #only samples that have V1 as "yes" in my dataset #are added to the subset

hist(data$HL)
hist(data_sub$HL, col="red", add=T)

希望这就是您的意思？

Answer 3

我认为你不能用条形直方图轻松做到这一点，因为你必须“交错”两个因子水平的条形......它需要对现在连续的 x 轴进行某种“离散化”（即它必须分成“类别”，并且在每个类别中，对于每个因子水平，您将有 2 个条形...

但是，如果您擅长绘制密度线函数，那么这很容易并且没有问题：

y <- rnorm(1000, 0, 1)
x <- rnorm(1000, 0.5, 2)
dx <- density(x)
dy <- density(y)
plot(dx, xlim = range(dx$x, dy$x), ylim = range(dx$y, dy$y), 
     type = "l", col = "red")
lines(dy, col = "blue")

Answer 4

我同意其他人的观点，即密度图比合并直方图的彩色条更有用，特别是当组的值混合时。这会非常困难，而且不会告诉你太多信息。您从其他人那里得到了关于密度图的一些很好的建议，这是我有时使用的密度图的 2 美分：

y <- rnorm(1000, 0, 1) 
x <- rnorm(1000, 0.5, 2) 
DF <- data.frame("Group"=c(rep(c("y","x"), each=1000)), "Value"=c(y,x))

library(sm)

with(DF, sm.density.compare(Value, Group, xlab="Grouping"))
title(main="Comparative Density Graph")
legend(-9, .4, levels(DF$Group), fill=c("red", "darkgreen")) 

Answer 5

目前尚不清楚您的数据布局是什么。直方图要求您有一个有序或连续的变量，以便可以创建中断。如果您还有一个单独的分组因素，您可以根据该因素绘制直方图。 lattice 包中的 histogram 函数的帮助页面上的第二个示例提供了这种分组和叠加密度曲线的一个很好的示例。

学习 R 博客是学习lattice和ggplot2绘图相对优点的一个很好的资源。 这是两个绘图系统并排比较的多部分系列的第一篇：

library(lattice)
 library(ggplot2)
 data(Chem97, package = "mlmRev")
#The lattice method:
pl <- histogram(~gcsescore | factor(score), data = Chem97)
 print(pl)

# The ggplot method:
 pg <- ggplot(Chem97, aes(gcsescore)) + geom_histogram(binwidth = 0.5) +
     facet_wrap(~score)
 print(pg)