Question

我正在做一个蒙特卡洛模拟，它会产生一个包含10000个观测值的矩阵，其中包含8个数字变量。我使用dplyr总结了8个变量，如下所示：

# A tibble: 2 x 8
  V1     V2     V3     V4     V5    V6     V7    V8
 <dbl>  <dbl>  <dbl>  <dbl>  <dbl> <dbl>  <dbl> <dbl>
1 29196. 12470. 6821.  5958.  22375. 6512. 10931. 2732.
2  1675.   419.   59.1   15.5  1636.  408.   858.  312.

其中第一行是每个变量的平均值，第二行是每个变量的标准偏差。我将如何使用此摘要统计信息创建一个8个柱形图，其高度为平均值，其误差线为标准偏差？我主要不确定如何填写ggplot的“ aes”部分。

谢谢。

Answer 1

正如评论中的@alistaire所述，使用ggplot2绘制数据时，您的数据实际上并不处于良好状态...所以下面是一个示例，其中我以构造数据的方式设置了一些数据，使用收集以将其拉入列，然后将其重新备份。然后，我使用此处的示例进行绘制：http://www.sthda.com/english/wiki/ggplot2-error-bars-quick-start-guide-r-software-and-data-visualization

我希望这对您有帮助...

library(tidyverse)                                         

df <- data.frame(V1=c(100, 20), V2=c(200,30), V3=c(150,15))
df                                                         
#>    V1  V2  V3
#> 1 100 200 150
#> 2  20  30  15

means <- df[1,]                                            
sds <-   df[2,]                                            

means_long <- gather(means, key='var', value='mean')       
means_long                                                 
#>   var mean
#> 1  V1  100
#> 2  V2  200
#> 3  V3  150

sds_long <- gather(sds, key='var', value='sd')             
sds_long                                                   
#>   var sd
#> 1  V1 20
#> 2  V2 30
#> 3  V3 15

sds_long %>%                                               
inner_join(means_long) ->                                  
combined_long                                              
#> Joining, by = "var"

combined_long                                              
#>   var sd mean
#> 1  V1 20  100
#> 2  V2 30  200
#> 3  V3 15  150

p <- ggplot(combined_long, aes(x=var, y=mean)) +           
geom_bar(stat="identity") +                                
geom_errorbar(aes(ymin=mean-sd, ymax=mean+sd), width=.2)   
p

如何在ggplot2中使用每个变量创建条形图？

1 个答案: