文档建议反馈控制台
首页
学习
活动
专区
圈层
工具
MCP广场
文章/答案/技术大牛
发布
社区首页 >问答首页 >R中的群wise

问R中的群wise
EN

Stack Overflow用户
提问于 2020-05-23 07:25:10
回答 3查看 405关注 0票数 2

我想为不同的车站准备xtabs。但它给我提供了跨车站的整体桌子。我使用了以下代码

代码语言：javascript
复制
library(tidyverse)

df %>% group_by(Station) %>% 
  xtabs( ~ Observed + Forecasted, data = .)

这给了我

代码语言：javascript
复制
#>        Forecasted
#>Observed   1   3   4   5   9
#>       1 132   5  31  31   3
#>       3  16   0   6  13   7
#>       4   6   0  13  23   8
#>       5   4   0  16  33  15
#>       9   0   0   0   2   0

但是我想要像这样的输出

代码语言：javascript
复制
Aizawl
#>        Forecasted
#>Observed   1   3   4   5   9
#>       1 132   5  31  31   3
#>       3  16   0   6  13   7
#>       4   6   0  13  23   8
#>       5   4   0  16  33  15
#>       9   0   0   0   2   0

Serchhip
#>        Forecasted
#>Observed   1   3   4   5   9
#>       1 132   5  31  31   3
#>       3  16   0   6  13   7
#>       4   6   0  13  23   8
#>       5   4   0  16  33  15
#>       9   0   0   0   2   0

然后，我想导出.csv文件中的输出。

数据

代码语言：javascript
复制
df = structure(list(Station = c("Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", "Aizawl", 
"Aizawl", "Aizawl", "Aizawl", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip", 
"Serchhip", "Serchhip", "Serchhip", "Serchhip", "Serchhip"), 
    Observed = c(1, 1, 1, 5, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 
    1, 1, 1, 1, 3, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 3, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 3, 4, 1, 1, 4, 1, 3, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 4, 4, 4, 3, 4, 1, 1, 1, 1, 1, 3, 5, 
    5, 3, 5, 3, 1, 1, 3, 1, 1, 1, 1, 1, 5, 3, 4, 1, 1, 1, 1, 
    1, 3, 1, 4, 1, 1, 1, 1, 1, 4, 4, 5, 1, 5, 4, 5, 5, 5, 5, 
    1, 5, 1, 4, 5, 4, 4, 5, 4, 5, 5, 3, 1, 5, 3, 4, 3, 4, 5, 
    5, 5, 5, 4, 4, 4, 5, 5, 5, 5, 5, 5, 4, 5, 3, 4, 4, 5, 3, 
    5, 4, 4, 5, 5, 5, 5, 4, 5, 5, 5, 5, 5, 5, 5, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 3, 5, 5, 1, 1, 3, 4, 1, 1, 1, 1, 1, 1, 1, 1, 3, 3, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 
    3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 4, 4, 1, 3, 4, 1, 1, 1, 1, 
    1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 4, 3, 9, 5, 5, 4, 1, 5, 1, 1, 1, 1, 4, 5, 5, 5, 
    5, 5, 5, 1, 1, 4, 1, 4, 4, 4, 5, 1, 1, 4, 3, 5, 1, 1, 4, 
    3, 5, 3, 4, 5, 3, 4, 4, 5, 5, 3, 4, 5, 5, 5, 5, 5, 4, 4, 
    4, 4, 5, 1, 9, 5, 5), Forecasted = c(1, 1, 1, 5, 5, 1, 1, 
    1, 5, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 5, 5, 
    5, 9, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 5, 5, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 4, 1, 4, 4, 1, 1, 5, 3, 1, 1, 1, 4, 5, 5, 5, 5, 
    1, 1, 1, 5, 5, 1, 5, 5, 5, 9, 4, 5, 4, 4, 4, 3, 4, 4, 1, 
    1, 5, 5, 4, 4, 4, 1, 1, 1, 4, 4, 4, 4, 4, 4, 1, 1, 5, 4, 
    4, 5, 4, 4, 4, 4, 5, 4, 5, 5, 5, 5, 5, 4, 5, 5, 4, 1, 1, 
    4, 4, 5, 5, 5, 5, 1, 4, 5, 5, 1, 4, 4, 9, 9, 9, 9, 9, 9, 
    9, 9, 9, 9, 9, 9, 9, 9, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 
    5, 5, 5, 9, 1, 1, 1, 5, 4, 1, 1, 1, 5, 4, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 5, 5, 5, 9, 5, 5, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 5, 5, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 5, 5, 4, 1, 1, 1, 1, 1, 4, 1, 1, 1, 
    4, 4, 4, 4, 1, 4, 1, 3, 1, 1, 1, 4, 4, 4, 4, 4, 4, 1, 1, 
    1, 4, 4, 3, 5, 5, 5, 4, 3, 5, 5, 5, 5, 5, 4, 5, 5, 5, 4, 
    5, 4, 4, 5, 5, 4, 4, 5, 4, 1, 4, 4, 5, 5, 4, 5, 4, 5, 4, 
    5, 5, 5, 1, 4, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 
    5, 5, 5, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 9)), row.names = c(NA, 
364L), class = "data.frame")
r
tidyverse
xtabs
EN

回答 3

Stack Overflow用户

回答已采纳

发布于 2020-05-23 08:14:50

如果您想要操作一个xtabs(三维数组)而不是list或data.frame，那么将分层变量(即Station )放在xtabs()中公式的最后位置。

代码语言：javascript
复制
res <- xtabs( ~ Observed + Forecasted + Station, df)
res

# , , Station = Aizawl
# 
#         Forecasted
# Observed  1  3  4  5  9
#        1 56  2 13 13  1
#        3 12  0  4  8  3
#        4  4  0  9 11  4
#        5  3  0  7 23  9
#        9  0  0  0  0  0
# 
# , , Station = Serchhip
# 
#         Forecasted
# Observed  1  3  4  5  9
#        1 76  3 18 18  2
#        3  4  0  2  5  4
#        4  2  0  4 12  4
#        5  1  0  9 10  6
#        9  0  0  0  2  0

class(res)
# [1] "xtabs" "table"

若要将数组打印到两个单独的csv文件，请执行以下操作：

代码语言：javascript
复制
lapply(dimnames(res)$Station, function(x) write.csv(res[,,x], paste0("table_", x, ".csv")))

若要将数组转换为data.frame格式：

代码语言：javascript
复制
library(tidyr)

res %>%
  as.data.frame() %>%
  pivot_wider(names_from = Forecasted, values_from = Freq)

# # A tibble: 10 x 7
#    Observed Station    `1`   `3`   `4`   `5`   `9`
#    <fct>    <fct>    <int> <int> <int> <int> <int>
#  1 1        Aizawl      56     2    13    13     1
#  2 3        Aizawl      12     0     4     8     3
#  3 4        Aizawl       4     0     9    11     4
#  4 5        Aizawl       3     0     7    23     9
#  5 9        Aizawl       0     0     0     0     0
#  6 1        Serchhip    76     3    18    18     2
#  7 3        Serchhip     4     0     2     5     4
#  8 4        Serchhip     2     0     4    12     4
#  9 5        Serchhip     1     0     9    10     6
# 10 9        Serchhip     0     0     0     2     0
票数 3
EN

Stack Overflow用户

发布于 2020-05-23 07:45:51

如果您想要将其写入CSV，最好将数据保持在CSV友好的格式中。这个输出

代码语言：javascript
复制
Aizawl
#>        Forecasted
#>Observed   1   3   4   5   9
#>       1 132   5  31  31   3
#>       3  16   0   6  13   7
#>       4   6   0  13  23   8
#>       5   4   0  16  33  15
#>       9   0   0   0   2   0

Serchhip
#>        Forecasted
#>Observed   1   3   4   5   9
#>       1 132   5  31  31   3
#>       3  16   0   6  13   7
#>       4   6   0  13  23   8
#>       5   4   0  16  33  15
#>       9   0   0   0   2   0

可能会导致问题，而且CSV格式友好的输出很容易用dplyr和tidyr编写.

代码语言：javascript
复制
library(dplyr)
library(tidyr)

df %>% 
  group_by(Station, Observed, Forecasted) %>% 
  tally() %>% 
  pivot_wider(names_from = Forecasted, values_from = n) %>% 
  replace(is.na(.), 0)
代码语言：javascript
复制
# A tibble: 9 x 7
# Groups:   Station, Observed [9]
  Station  Observed   `1`   `3`   `4`   `5`   `9`
  <chr>       <dbl> <int> <int> <int> <int> <int>
1 Aizawl          1    56     2    13    13     1
2 Aizawl          3    12     0     4     8     3
3 Aizawl          4     4     0     9    11     4
4 Aizawl          5     3     0     7    23     9
5 Serchhip        1    76     3    18    18     2
6 Serchhip        3     4     0     2     5     4
7 Serchhip        4     2     0     4    12     4
8 Serchhip        5     1     0     9    10     6
9 Serchhip        9     0     0     0     2     0

你可以很容易地导出到CSV在这一点上或直接链接您的管道。

代码语言：javascript
复制
df %>% 
  group_by(Station, Observed, Forecasted) %>% 
  tally() %>% 
  pivot_wider(names_from = Forecasted, values_from = n) %>% 
  replace(is.na(.), 0) %>%
  readr::write_csv("my_output.csv")
票数 2
EN

Stack Overflow用户

发布于 2020-05-23 07:55:30

使用by。

代码语言：javascript
复制
res <- by(df, df$Station, xtabs, f=~Observed + Forecasted)
res
# df$Station: Aizawl
#          Forecasted
# Observed  1  3  4  5  9
#           1 56  2 13 13  1
#           3 12  0  4  8  3
#           4  4  0  9 11  4
#           5  3  0  7 23  9
# --------------------------------------------------------- 
#   df$Station: Serchhip
#         Forecasted
# Observed  1  3  4  5  9
#           1 76  3 18 18  2
#           3  4  0  2  5  4
#           4  2  0  4 12  4
#           5  1  0  9 10  6
#           9  0  0  0  2  0

## Output .csv
lapply(names(res), function(x) write.csv(res[[x]], paste0("myfile_", x,".csv")))
票数 2
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：

https://stackoverflow.com/questions/61968797

复制
相关文章

相似问题

领券

腾讯云开发者

扫码关注腾讯云开发者

扫码关注腾讯云开发者

领取腾讯云代金券

热门产品

热门推荐

更多推荐

Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有 

深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 粤公网安备44030502008569号

腾讯云计算（北京）有限责任公司 京ICP证150476号 |  京ICP备11018762号

问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档