de la trame de données au tableau 4D

Question

Je dois remodeler un cadre de données ( "d" ci-dessous dans le code fictif) à un tableau 4D. Dans ce cas précis, je voudrais que la variable "count" ci-dessous pour être inclus dans un tableau 4D avec des dimensions site par année par date par méthode.

d <- data.frame(site=c(rep("aaa", 4), rep("bbb", 4)), date=c(114:117, 220:223), 
                year=c(2005:2008, 2009:2012), count=rpois(8, 34),
                method=c(rep(1, 2), rep(2, 2), rep(1, 1), rep(2, 3)))

Answer 1

Nous pourrions utiliser xtabs() .

xtabs(count ~ ., d)
# , , year = 2005, method = 1
# 
# date
# site  114 115 116 117 220 221 222 223
# aaa  45   0   0   0   0   0   0   0
# bbb   0   0   0   0   0   0   0   0
# 
# , , year = 2006, method = 1
# 
# date
# site  114 115 116 117 220 221 222 223
# aaa   0  38   0   0   0   0   0   0
# bbb   0   0   0   0   0   0   0   0
# 
# [...]
# , , year = 2005, method = 2
# 
# date
# site  114 115 116 117 220 221 222 223
# aaa   0   0   0   0   0   0   0   0
# bbb   0   0   0   0   0   0   0   0
# 
# , , year = 2006, method = 2
# 
# date
# site  114 115 116 117 220 221 222 223
# aaa   0   0   0   0   0   0   0   0
# bbb   0   0   0   0   0   0   0   0
# 
# [...]

Données

d <- structure(list(site = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 
2L), .Label = c("aaa", "bbb"), class = "factor"), date = c(114L, 
115L, 116L, 117L, 220L, 221L, 222L, 223L), year = 2005:2012, 
    count = c(45L, 38L, 35L, 36L, 30L, 42L, 25L, 31L), method = c(1, 
    1, 2, 2, 1, 2, 2, 2)), class = "data.frame", row.names = c(NA, 
-8L))