Question

任何人都可以帮我计算一个新变量，它会从某些因素中得到一个独特的组合吗？

假设主题因素中有4个（A，B，C，D），对于10个主题中的任何一个，每个组合重复8次，这就是如何我的数据看起来似乎代表了它的实际结构：

library(AlgDesign) #for generating a factorial design)
df <-gen.factorial(c(2,2,2,2,8,10), factors = "all", 
                      varNames = c("A", "B", "C", "D", "replication", "Subject"))
> head(df)
  A B C D replication Subject
1 1 1 1 1           1       1
2 2 1 1 1           1       1
3 1 2 1 1           1       1
4 2 2 1 1           1       1
5 1 1 2 1           1       1
6 2 1 2 1           1       1
> tail(df)
     A B C D replication Subject
1275 1 2 1 2           8      10
1276 2 2 1 2           8      10
1277 1 1 2 2           8      10
1278 2 1 2 2           8      10
1279 1 2 2 2           8      10
1280 2 2 2 2           8      10

在这个例子中，replication被简单地生成以强制执行8个代表，但它并没有＃34;代码＆＃34;结合本身。

我的原始数据只包含变量A，B，C，D和Subject，我想计算 replication以一种具有不同价值的方式但，适用于A，B，C，D

的每个组合

Answer 1

library(AlgDesign) 
library(dplyr)

df <-gen.factorial(c(2,2,2,2,8,10), factors = "all", 
                   varNames = c("A", "B", "C", "D", "replication", "Subject"))

df %>%
  rowwise() %>%                                             # for each row
  mutate(factors = paste0(c(A,B,C,D), collapse = "_")) %>%  # create a combination of your factors
  ungroup() %>%                                             # forget the row grouping
  mutate(replication_upd = as.numeric(factor(factors)))     # create a number based on the combination you have

# # A tibble: 1,280 x 8
#   A     B     C     D     replication Subject factors replication_upd
#   <fct> <fct> <fct> <fct> <fct>       <fct>   <chr>             <dbl>
# 1 1     1     1     1     1           1       1_1_1_1               1
# 2 2     1     1     1     1           1       2_1_1_1               9
# 3 1     2     1     1     1           1       1_2_1_1               5
# 4 2     2     1     1     1           1       2_2_1_1              13
# 5 1     1     2     1     1           1       1_1_2_1               3
# 6 2     1     2     1     1           1       2_1_2_1              11
# 7 1     2     2     1     1           1       1_2_2_1               7
# 8 2     2     2     1     1           1       2_2_2_1              15
# 9 1     1     1     2     1           1       1_1_1_2               2
#10 2     1     1     2     1           1       2_1_1_2              10
# # ... with 1,270 more rows

您可以删除任何不必要的变量。我把它们留在那里，这样你就可以看到这个过程是如何运作的。

另一种选择是

# create a look up table based on unique combinations and assign them a number
df %>% distinct(A,B,C,D) %>% mutate(replication_upd = row_number()) -> look_up

# join back to original dataset
df %>% inner_join(look_up, by=c("A","B","C","D")) %>% tbl_df()

# # A tibble: 1,280 x 7
#   A     B     C     D     replication Subject replication_upd
#   <fct> <fct> <fct> <fct> <fct>       <fct>             <int>
# 1 1     1     1     1     1           1                     1
# 2 2     1     1     1     1           1                     2
# 3 1     2     1     1     1           1                     3
# 4 2     2     1     1     1           1                     4
# 5 1     1     2     1     1           1                     5
# 6 2     1     2     1     1           1                     6
# 7 1     2     2     1     1           1                     7
# 8 2     2     2     1     1           1                     8
# 9 1     1     1     2     1           1                     9
# 10 2     1     1     2     1           1                    10
# # ... with 1,270 more rows

请注意，第一种方法是根据我们创建的新变量（即订单A，B，C，D）选择数字，第二种方法使用数据集的初始顺序来选择每个唯一组合的数字。

根据R中的一组因子创建组合ID号

1 个答案: