四分数据的所有排列

Question

我有一个数据集，包含8年的季度数据。如果我从理论上可以从理论上构建“新”年的年份中随机选择每个季度。例如：新年= Q1（2009），Q2（2012），Q3（2010），Q4（2015）。

我遇到的问题是我想构建一个包含所有这些排列的数据集。 8年和4个季度，会给我4^8 = 65536“新”年。这是最好的一个嵌套循环处理，还是那里的功能可以更好地工作？

Answer 1

我们可以使用expand.grid创造一切可能的组合的矩阵：

nrow(do.call('expand.grid', replicate(8, 1:4, simplify=FALSE))) 
[1] 65536 

Answer 2

您可能要稍等一下，看看是否有人给你一少“janky”的答案，但这个例子需要时间系列，每年都会有所有排列并且没有重复的季度，并且将旧年和季度信息作为列返回这些新的年份值。

set.seed(1234) 

# Make some fake data 
q_dat <- data.frame(year = c(rep(2011,4), 
          rep(2012,4), 
          rep(2013,4)), 
        quarters = rep(c("Q1","Q2","Q3","Q4"),3), 
        x = rnorm(12)) 
q_dat 

    year quarters   x 
1 2011  Q1 -1.2070657 
2 2011  Q2 0.2774292 
3 2011  Q3 1.0844412 
4 2011  Q4 -2.3456977 
5 2012  Q1 0.4291247 
6 2012  Q2 0.5060559 
7 2012  Q3 -0.5747400 
8 2012  Q4 -0.5466319 
9 2013  Q1 -0.5644520 
10 2013  Q2 -0.8900378 
11 2013  Q3 -0.4771927 
12 2013  Q4 -0.9983864 

那么要做的就是

1，以时间序列

2的所有可能的组合，删除所有重复所以每个由一年没有同季的它。

# Expand out all possible combinations of our three years 
q_perms <- expand.grid(q1 = 1:nrow(q_dat), q2 = 1:nrow(q_dat) , 
         q3 = 1:nrow(q_dat), q4 = 1:nrow(q_dat)) 

# remove any duplicate combinations 
# EX: So we don't get c(2011Q1,2011Q1,2011Q1,2011Q1) as a year 
q_perms <- q_perms[apply(q_perms,1,function(x) !any(duplicated(x))),] 

# Transpose the grid, remake it as a data frame, and lapply over it 
l_rand_dat <- lapply(data.frame(t(q_perms)),function(x) q_dat[x,]) 

# returns one unique year per list 
l_rand_dat[[30]] 
    year quarters   x 
5 2012  Q1 0.4291247 
6 2012  Q2 0.5060559 
2 2011  Q2 0.2774292 
1 2011  Q1 -1.2070657 


# bind all of those together 
rand_bind <- do.call(rbind,l_rand_dat) 

head(rand_bind) 

     year quarters   x 
X172.4 2011  Q4 -2.3456977 
X172.3 2011  Q3 1.0844412 
X172.2 2011  Q2 0.2774292 
X172.1 2011  Q1 -1.2070657 
X173.5 2012  Q1 0.4291247 
X173.3 2011  Q3 1.0844412 

这是一个相当内存密集的答案。如果有人可以跳过“制定所有可能的组合”的步骤，那么这将是一个重大的改进。

Answer 3

我想你想在4个季度的8年多的组合，这样的组合数为8^4 = 4096：

> x <- years <- 2008:2015 
> length(x) 
[1] 8 
> comb <- expand.grid(x, x, x, x) 
> head(comb) 
    Var1 Var2 Var3 Var4 
1 2008 2008 2008 2008 
2 2009 2008 2008 2008 
3 2010 2008 2008 2008 
4 2011 2008 2008 2008 
5 2012 2008 2008 2008 
6 2013 2008 2008 2008 
> tail(comb) 
    Var1 Var2 Var3 Var4 
4091 2010 2015 2015 2015 
4092 2011 2015 2015 2015 
4093 2012 2015 2015 2015 
4094 2013 2015 2015 2015 
4095 2014 2015 2015 2015 
4096 2015 2015 2015 2015 
> nrow(comb) 
[1] 4096 

每一行都是一年VAR1，VAR2，VAR3，VAR4是4个季度。