更新

Question

这是我的问题，我无法解决所有问题。

假设我们有如下代码：

## A data frame named a    
a <- data.frame(A = c(0,0,1,1,1), B = c(1,0,1,0,0), C = c(0,0,1,1,0), D = c(0,0,1,1,0), E = c(0,1,1,0,1))
## 1st function calculates all the combinaisons of colnames of a and the output is a character vector named item2
items2 <- c()
countI <- 1 
while(countI <= ncol(a)){
        for(i in countI){
                countJ <- countI + 1
                while(countJ <= ncol(a)){
                        for(j in countJ){
                                items2 <- c(items2, paste(colnames(a[i]), colnames(a[j]), collapse = '', sep = ""))
                        }
                        countJ <- countJ + 1
                }
                countI <- countI + 1
        }
}

这是我尝试解决的代码（输出是一个名为 count_1 的数字向量）：

## 2nd function
colnames(a) <- NULL ## just for facilitating the calculation
count_1 <- numeric(ncol(a)*2)
countI <- 1
while(countI <= ncol(a)){
        for(i in countI){
                countJ <- countI + 1
                while(countJ <= ncol(a)){
                        for(j in countJ){
                                s <- a[, i]
                                p <- a[, j]
                                count_1[i*2] <- as.integer(s[i] == p[j] & s[i] == 1)
                        }
                        countJ <- countJ + 1
                }
                countI <- countI + 1
        }
}

但是当我在RStudio Console中执行此代码时，会返回非期望结果！：

 count_1
 [1] 0 0 0 0 0 1 0 1 0 0

但是，我期待以下结果：

count_1
[1] 1 2 2 2 1 1 1 1 2 1

您可以看到访问以下网址，您可以在Dropbox上找到图片以获取详细说明。 https://www.dropbox.com/s/5ylt8h8wx3zrvy7/IMAG1074.jpg?dl=0

我会尝试解释一下，我发布了第一个功能（代码）只是为了向您展示我正在寻找的东西，这是一个例子。我试图从第二个函数（代码）中得到的是计算每行中出现的数字 1 （首先我们放counter = 0）的次数（每行两列（例如AB）必须等于两列中的一列，以表示counter = counter + 1）我们通过将所有其他列（AC，AD，AE，BC，BD，BE，CD， CE，然后是DE），组合是n!/2!(n-2)!，这意味着例如我是否有以下数据框：

a =

A B C D E
0 1 0 0 0

0 0 0 0 1

1 1 1 1 1

1 0 0 1 0

1 0 1 0 1

然后，通过组合两个第一列，每行的数字 1 的出现次数如下：（请注意，我将colnames(a) <- NULL仅用于方便工作并且更加清晰）

0 1 0 0 0

0 0 0 0 1

1 1 1 1 1

1 0 0 1 0

1 0 1 0 1

### Example 1: #####################################################

所以我从这里开始（对于A和B列（AB））

s <- a[, i]
## s is equal to
## [1] 0 0 1 1 1
p <- a[, j]
## p is equal to
## [1] 1 0 1 0 0

然后我会在两个向量中查找数字 1 的出现，条件必须相同，即a[, i] == 1 && a[, j] == 1 && a[, i] == a[, j]，对于此示例，数字向量将是[1] 1

### Example 2: #####################################################

从这里开始（对于A和D（AD）列）

s <- a[, i]
## s is equal to
## [1] 0 0 1 1 1
p <- a[, j]
## p is equal to
## [1] 0 0 1 1 0

然后我会在两个向量中查找数字 1 的出现，条件是它必须相同，即a[, i] == 1 && a[, j] == 1 && a[, i] == a[, j]，对于此示例，数字向量将是[1] 2

等等，我的名为count_1的数字向量等于：

[1] 1 2 2 2 1 1 1 1 2 1

虽然count_1的每个索引都是其他列的组合（没有数据框的名称）

AB AC AD AE BC BD BE CD CE DE

1 2 2 2 1 1 1 1 2 1

Answer 1

根本不清楚你完成了什么。

对于第一个代码块，这是一些丑陋的R编码，涉及一大堆不必要的while / for循环。

您可以在一行中获得相同的结果items2。

items2 <- sort(toupper(unlist(sapply(1:4, function(i)
    sapply(5:(i+1), function(j)
        paste(letters[i], letters[j], sep = ""))))));
items2;
# [1] "AB" "AC" "AD" "AE" "BC" "BD" "BE" "CD" "CE" "DE"

关于第二个代码块，请解释您尝试计算的内容。这些while / for循环可能与第一种情况一样不必要。

更新

请注意，这是基于帖子开头定义的a。您的预期输出基于不同的a，您在帖子中进一步更改。

不需要for / while循环，这两种＆＃34;功能＆＃34;可以用两个单行写成。

# Your sample dataframe a
a <- data.frame(A = c(0,0,1,1,1), B = c(1,0,1,0,0), C = c(0,0,1,1,0), D = c(0,0,1,1,0), E = c(0,1,1,0,1))

# Function 1
items2 <- toupper(unlist(sapply(1:(ncol(a) - 1), function(i) sapply(ncol(a):(i+1), function(j)
        paste(letters[i], letters[j], sep = "")))));
# Function 2
count_1 <- unlist(sapply(1:(ncol(a) - 1), function(i) sapply(ncol(a):(i+1), function(j)
        sum(a[, i] + a[, j] == 2))));

# Add names and sort
names(count_1) <- items2;
count_1 <- count_1[order(names(count_1))];
# Output
count_1;
#AB AC AD AE BC BD BE CD CE DE
# 1  2  2  2  1  1  1  2  1  1

如何计算R

1 个答案:

更新