在同一张图上绘制两个因素

Question

假设我有两个因素，我想将它们绘制在同一个图上，这两个因素具有相同的水平。

s1 <- c(rep("male",20), rep("female", 30))
s2 <- c(rep("male",10), rep("female", 40))
s1 <- factor(s1, levels=c("male", "female"))
s2 <- factor(s2, levels=c("male", "female"))

我本以为使用表函数会产生正确的绘图结果，但它弹出了。

table(s1, s2)
        s2 
s1       male female   
male     10     10   
female    0     30

所以实际上有两个问题，表函数是做什么来获得这个结果的，以及我可以使用什么其他函数来使用具有相同级别的函数创建具有 2 个系列的图形？

另外，如果这是一个因素，我会使用 gplots 包中的 barplot2 来绘制它。

Answer 1

您可以使用lattice包获得稍微更详细的结果：

s1 <- factor(c(rep("male",20), rep("female", 30)))
s2 <- factor(c(rep("male",10), rep("female", 40)))
D <- data.frame(s1, s2)

library(lattice)
histogram(~s1+s2, D, col = c("pink", "lightblue"))

或者如果您想要男性/女性并排以便于比较:

t1 <- table(s1)
t2 <- table(s2)
barchart(cbind(t1, t2), stack = F, horizontal = F)

Answer 2

来自 ?table：

“table”使用交叉分类因素来构建意外事件
每个因子水平组合的计数表。

当您执行 table(s1,s2) 时，函数会将 s1 和 s2 视为配对结果。实际上，它告诉您，如果您采用cbind(s1,s2)，那么将会有 10 行男性-男性，10 行男性-女性，依此类推。

要理解这一点，请考虑一个非常简单的例子：

a <- c("M","M","F","F")
b <- c("F","F","M","M")
table(a,b)

  b
a   F M
  F 0 2
  M 2 0

您应该做的是：

t1 <- table(s1)
t2 <- table(s2)
barplot(cbind(t1,t2), beside=TRUE, col=c("lightblue", "salmon"))

Answer 3

产生稍微不同形式的绘图的两个选项是

plot(s1, s2)

和

plot(table(s1,s2))

前者是 spinplot，马赛克图的一种特殊情况，由 table 的 plot 方法生成（第二个示例）。有关更多详细信息，请参阅 ?spineplot 和 ?mosaicplot，如果您愿意，您可以直接使用这些函数，而不是通用的 plot()。

另请查看 Meyer 等人在 CRAN 上的 vcd 包中的 mosaic() 函数 （CRAN 上的 vcd 链接）

table() 正在生成 这两个因素的列联表。

Answer 4

嗯..我不认为卡梅伦想要创建列联表。如果我理解正确的话，我认为他想创建一个包含两个变量的数据框，其中 s1 和 s2 似乎是相同大小的向量。 （长度（s1）==长度（s2））。

在这种情况下，他只需要使用以下方法创建一个“表”（我认为他指的是 data.frame）：

df = data.frame(s1=s1, s2=s2);

然后在同一个图中绘制 2 个系列。

至于绘制这些东西的第二个问题，我会使用 matplot。例如：

matplot(1:10, data.frame(a=rnorm(10), b=rnorm(10)), type="l", lty=1, lwd=1, col=c("blue","red"))

鉴于他将 2 个向量的数据组织在名为“df”的单个 data.frame 中，他可以执行以下操作：

matplot(df, type="l", lty=1, lwd=1, col=c("blue","red"))

希望这会有所帮助。