R语言 同一图上的两个箱线图

llew8vvj  于 2023-06-19  发布在  其他
关注(0)|答案(6)|浏览(203)

我有两个不同的数据集,它们有不同的观察次数。我想在同一个图表上绘制两个箱线图,这样比较起来就更容易了。我可以绘制一个箱线图,但如果没有它们并排,就很难发现任何差异。
我有一些假数据。

Group A
V1    V2   V3    V4     V5
6.5   2    11    0.5    6
7     1    8     0.34   8
5.4   4    7.8   0.45   5
3.4   6    9.1   0.72   5

Group B
V1    V2   V3    V4     V5
5.0   5    9     0.4    7
2     7    5.2   0.69   5
3.2   2    2.9   0.79   2
6.8   9    6.5   0.43   6
4.7   3    3.8   0.49   4
5.5   4    7.4   0.94   3

我不知道怎么画这个,所以我没有一个例子。我将尽力描述情节。我想将A组和B组的变量1绘制在同一张图上。所以在一个图中,我会有A组的箱形图和B组的另一个箱形图填充V1的数据。所以这两个箱线图是并排的。如果有5个变量,我将有5个图表,每个图表都有2个并排的箱线图。如果我说得不够清楚请告诉我。谢谢你。

kuhbmx9i

kuhbmx9i1#

ggplot最适合于“长格式”数据(例如,值、变量和组中的每一个都有一列)。您可以按如下方式重新排列数据:

A <- read.table(text='V1    V2   V3    V4     V5
6.5   2    11    0.5    6
7     1    8     0.34   8
5.4   4    7.8   0.45   5
3.4   6    9.1   0.72   5', header=TRUE)

B <- read.table(text='V1    V2   V3    V4     V5
5.0   5    9     0.4    7
2     7    5.2   0.69   5
3.2   2    2.9   0.79   2
6.8   9    6.5   0.43   6
4.7   3    3.8   0.49   4
5.5   4    7.4   0.94   3', header=TRUE)

d <- rbind(cbind(stack(A), group='A'), cbind(stack(B), group='B'))

前几行看起来像这样:

head(d)

##   values ind group
## 1    6.5  V1     A
## 2    7.0  V1     A
## 3    5.4  V1     A
## 4    3.4  V1     A
## 5    2.0  V2     A
## 6    1.0  V2     A

现在我们可以这样做:

library(ggplot2)
ggplot(d, aes(group, values)) + 
  geom_boxplot() +
  facet_wrap(~ind, scales='free_y')

f5emj3cl

f5emj3cl2#

我想出的解决方案是将两个data.frame和一个变量组合起来,指示观察值属于哪个组。然后,您可以使用reshape2中的melt函数将数据转换为data.frame,以便打印。可以使用facet_gridfacet_wrap为不同变量创建单独的图。这是一种方法:

library(ggplot2)
library(reshape2)

# Combine two data.frame
df <- rbind(GroupA, GroupB)

# Create variable Group
df$Group <- rep(c("A", "B"), c(dim(GroupA)[1], dim(GroupB)[1]))

# Transform to long format
df <- melt(df, "Group")

ggplot(df, aes(x=Group, y=value)) + geom_boxplot() + facet_grid(~ variable)

ymdaylpp

ymdaylpp3#

假设数据集的名称是grpa(Group A)和grpb(Group B)。首先,给它们每个添加一个变量Group
grpa$Group <-"A"
grpb$Group <-"B"
然后将它们组合成一个数据框架
combined <- rbind(grpa,grpb)
然后使用ggplot绘制,如下所示:
ggplot(combined,aes(x= factor(Group), y=V1))+geom_boxplot()

根据需要贴标签。

zhte4eai

zhte4eai4#

# Adding a variable to the dataframes Group_A & Group_B as done from pervious users
Group_A$fac <- "A"
Group_B$fac <- "B"
Group_c <- rbind(Group_A,Group_B)
df <- melt(Group_c)

#You can plot the same in bwplot from library(lattice) 

bwplot(value~fac|variable,data=df,scales=list(relation="free"),as.table=T)

50pmv0ei

50pmv0ei5#

group_a<-data.frame(
  V1 = c(6.5, 7, 5.4, 3.4),
  V2 = c(2, 1, 4, 6),
  V3 = c(11, 8, 7.8, 9.1),
  V4 = c(0.5, 0.34, 0.45, 0.72),
  V5 = c(6, 8, 5, 5)
  )
group_a
group_b<-data.frame(
  V1 =c(5.0,2,3.2,6.8,4.7,5.5),
  V2 =c(5,7,2,9,3,4),
  V3=c(9,5.2,2.9,6.5,3.8,7.4),
  V4=c(0.4,0.69,0.79,0.43,0.49,0.94),
  V5=c(7,5,2,6,4,3)
  )
group_b
# Create a new empty data frame with header
empty_df <- data.frame(Value = numeric(), var = character(), group = 
character(), stringsAsFactors = FALSE)
# Loop through each column in group_a
for (col in names(group_a)) {
  # Extracting column values
  column_values <- group_a[[col]]
  # Create a temporary data frame to hold the extracted column
  temp_df <- data.frame(Value = column_values, var = col, group = "A")
  # Append the temporary data frame to empty_df
  empty_df <- rbind(empty_df, temp_df)
  }

for (colb in names(group_b)) {
 column_values<-group_b[,colb]
 tempp_df<- data.frame(Value=column_values, var = colb, group = "B")
 empty_df<-rbind(empty_df,tempp_df)
}
ggplot(empty_df, aes(x = group, y = Value)) +
geom_boxplot() +
facet_wrap(~ var, scales = "free_y")

solution

bd1hkmkf

bd1hkmkf6#

par(mfrow=c(1,2))
   summary(A)
   summary(B)
   boxplot(A,ylim=summary(A)[[1]][1]) ##not sure about this just find where y is min
   boxplot(B,ylim=summary(B)[[1]][1]) ## still not sure
    ## adjusts the ylims in a way so that they are easy to compare you can also use boxplot(A,B) but that would make the graph look weird

相关问题