Question

我有一个如下数据框，我想创建一个矩阵，其中平均睡眠持续时间（SLP）根据3个招募网站（网站）和5个招聘年（年）显示。

SLP site year 
8.6  1   2008  
7.2  1   2005  
6.4  2   2006  
9.5  3   2007  
6.1  2   2009  
3.6  1   2005 
8.6  1   2008  
7.2  1   2005  
6.4  2   2006  
9.5  3   2007  
6.1  2   2009  
5.1  3   2008 
2.1  2   2006

我想要的输出是：

       1      2      3 
2005  6.00    -      -
2006   -     4.97    -
2007   -      -     9.5
2008  8.60    -     5.1
2009   -     6.10    -

列名称是站点的变量，行名称是年份的变量，每个单元格中的值是SLP的平均值。我该怎么做？

Answer 1

以下是一些不使用软件包的不同解决方案：

1）tapply 这不使用任何包。它产生一个"matrix"输出，其中包含空单元格的NA值：

tapply(DF$SLP, DF[c("year", "site")], mean)

，并提供：

      site
year     1        2   3
  2005 6.0       NA  NA
  2006  NA 4.966667  NA
  2007  NA       NA 9.5
  2008 8.6       NA 5.1
  2009  NA 6.100000  NA

2）aggregate / xtabs 这会使用aggregate + xtabs。这将创建一个类c("xtabs", "table")的对象，其中空单元格的值为零：

fo <- SLP ~ year + site
xtabs(fo, aggregate(fo, DF, mean))

给予;

      site
year          1        2        3
  2005 6.000000 0.000000 0.000000
  2006 0.000000 4.966667 0.000000
  2007 0.000000 0.000000 9.500000
  2008 8.600000 0.000000 5.100000
  2009 0.000000 6.100000 0.000000

3）聚合/重塑这也使用aggregate但使用reshape而不是xtabs。它为空单元格提供了带有NA的数据帧r。最后一行使列名与先前的解决方案一致，如果这不重要，可以省略。

 ag <- aggregate(SLP ~ site + year, DF, mean)
 r <- reshape(ag, dir = "wide", idvar = "year", timevar = "site")
 names(r) <- sub(".*[.]", "", names(r))

，并提供：

> r
  year   1        2   3
1 2005 6.0       NA  NA
3 2006  NA 4.966667  NA
5 2007  NA       NA 9.5
2 2008 8.6       NA 5.1
4 2009  NA 6.100000  NA

注意：使用的可重现形式的输入DF是：

DF <- structure(list(SLP = c(8.6, 7.2, 6.4, 9.5, 6.1, 3.6, 8.6, 7.2, 
6.4, 9.5, 6.1, 5.1, 2.1), site = c(1L, 1L, 2L, 3L, 2L, 1L, 1L, 
1L, 2L, 3L, 2L, 3L, 2L), year = c(2008L, 2005L, 2006L, 2007L, 
2009L, 2005L, 2008L, 2005L, 2006L, 2007L, 2009L, 2008L, 2006L
)), .Names = c("SLP", "site", "year"), class = "data.frame", row.names = c(NA, 
-13L))

Answer 2

我们可以使用socklen_t client_sin_len; sockaddr_in *client_sin = new sockaddr_in; client_sin_len = sizeof(sockaddr_in ); std::cout << "New User ! " << std::endl; if ((cs = accept(fd, reinterpret_cast<struct sockaddr *>(client_sin), &client_sin_len)) == -1)

acast

或使用library(reshape2) acast(df1, year~site, value.var="SLP", mean)

中的tapply

base R

Answer 3

另一种解决方案

library(tidyr)
library(dplyr)

df%>% 
  group_by(year, site) %>%
    summarise(m=mean(SLP)) %>%
  spread(site, m )%>%
as.matrix()

Answer 4

在@ g-grothendieck使用xtabs的基础上，我们可以将其与table和ifelse结合使用，以返回相同的结果。

# get a count of the number of observations per matrix cell (filling 0s with 1)
tempTab <- ifelse(with(df, table(year, + site)) == 0, 1, with(df, table(year, + site)))

tempTab

year   1 2 3
  2005 3 1 1
  2006 1 3 1
  2007 1 1 2
  2008 2 1 1
  2009 1 2 1

现在使用xtabs，它在单元格中多个观察值时返回值的总和，并除以tempTab得到均值。

xtabs(SLP ~ year + site, df) / tempTab
      site
year          1        2        3
  2005 6.000000 0.000000 0.000000
  2006 0.000000 4.966667 0.000000
  2007 0.000000 0.000000 9.500000
  2008 8.600000 0.000000 5.100000
  2009 0.000000 6.100000 0.000000

如何从数据框创建矩阵？

4 个答案: