无法将数据帧转换为r中的设计矩阵
我几天前开始使用R,并且确实可以使用一些帮助:D 我目前有一个包含200个观察值和12个变量的数据框(它们表示点击,所以他们被命名为c1-c12)。无法将数据帧转换为r中的设计矩阵
背景: 用户在字段上单击(最多12次)。这些字段下令6行和4列和所述用户的点击次数被记录为电话号码。:
[,1] [,2] [,3] [,4]
[1,] "11" "21" "31" "41"
[2,] "12" "22" "32" "42"
[3,] "13" "23" "33" "43"
[4,] "14" "24" "34" "44"
[5,] "15" "25" "35" "45"
[6,] "16" "26" "36" "46"
每个值的第一位数字表示行,第二个列。这就是为什么他们被命名为11-16,21-26 ...... 41-46。 值范围从c(11:16,21:26,31:36,41:46),因此有24个可能的值。
我能够将数据读入R和第10组的意见是这样的:
我的数据:
c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 c12
1 33 43 63 23 34 32 31 41 61 21 NA NA
2 24 23 22 21 31 61 41 NA NA NA NA NA
3 61 62 63 64 31 32 33 34 41 42 43 44
4 31 32 33 34 21 22 23 24 41 NA NA NA
5 11 12 13 14 22 23 32 33 62 63 42 52
6 51 52 53 54 61 62 63 64 31 32 33 34
7 31 21 61 62 63 64 33 23 NA NA NA NA
8 41 42 43 44 32 33 62 63 52 53 61 64
9 61 62 63 64 21 22 24 23 34 31 41 44
10 51 52 53 54 24 34 21 31 33 23 61 63
我想我的数据转换成某种设计矩阵,这是基于在后台给出的位置表,但: 作为所需的输出我需要一个矩阵每个观察,再次包含6行4列和1时,他们点击该领域(即位置的值是ob的一部分为所有其他职位0)。 先观察会再看看这样的:
0 1 1 1 0 1
0 0 1 0 0 0
0 1 1 1 0 1
0 0 1 0 0 0
能给我哪些软件包等等。我应该看看,也许有技巧作为一个准则或类似这样的东西的提示?
我的想法是为位置创建一个矩阵,然后将它们与观察值相乘,但我目前真的很挣扎,不知道从哪里开始。
R中我的位置矩阵现在看起来是这样的:
#construct design matrix
alt1 <- c(paste0(1, 1:6))
alt2 <- c(paste0(2, 1:6))
alt3 <- c(paste0(3, 1:6))
alt4 <- c(paste0(4, 1:6))
positions <- matrix(c(alt1, alt2, alt3, alt4), nrow = 6)
(背景提供)
**有可能是一种方式更容易和更聪明的解决方案,因为我的价值观包含ROWNUMBER(第一位数)和columnnr(第二位),但我不知道如何通过位置并将其写入到矩阵**
dput(head(data))
structure(list(c1 = c("33", "24", "61", "31", "11", "51"), c2 = c("43",
"23", "62", "32", "12", "52"), c3 = c("63", "22", "63", "33",
"13", "53"), c4 = c("23", "21", "64", "34", "14", "54"), c5 = c("34",
"31", "31", "21", "22", "61"), c6 = c("32", "61", "32", "22",
"23", "62"), c7 = c("31", "41", "33", "23", "32", "63"), c8 = c("41",
"", "34", "24", "33", "64"), c9 = c("61", NA, "41", "41", "62",
"31"), c10 = c("21", NA, "42", "", "63", "32"), c11 = c("", NA,
"43", NA, "42", "33"), c12 = c(NA, NA, "44", NA, "52", "34")), .Names = c("c1",
"c2", "c3", "c4", "c5", "c6", "c7", "c8", "c9", "c10", "c11",
"c12"), row.names = c(NA, 6L), class = "data.frame")
欢呼和感谢, Sidebob
您可以通过一个简单的循环获取各条线并将它们转换为矩阵 。
for(row in 1:nrow(data)) {
x <- as.numeric(data[row,])
i <- x %% 10
j <- x %/% 10
mat <- matrix(0, 6, 6)
mat[cbind(i,j)] <- 1
print(mat)
}
这给
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0 1 1 1 0 1
[2,] 0 0 1 0 0 0
[3,] 0 1 1 1 0 1
[4,] 0 0 1 0 0 0
[5,] 0 0 0 0 0 0
[6,] 0 0 0 0 0 0
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0 1 1 1 0 1
[2,] 0 1 0 0 0 0
[3,] 0 1 0 0 0 0
[4,] 0 1 0 0 0 0
[5,] 0 0 0 0 0 0
[6,] 0 0 0 0 0 0
...
我不知道你期望哪种输出。您可以打印 (并输出到文件的sink()
),您可以将write()
更改为 文件,也可以将write.table()
更改为类似csv的格式。
顺便说一句,你可能会遇到声称循环很慢,应该 避免R,在某种程度上,它是真实的,但在这种情况下,它工作正常,并 简单易懂。
您能更具体地了解您的投入和预期产出? –
当然:用户在字段上点击(最多12次)。这些字段分为6行和4列。这就是为什么他们被命名为11-16,21-26 ...... 41-46。 我现在输入的是200个观察值(用户)及其12个点击,记录为数字。 作为输出,我需要每个用户一个矩阵,包含6行和4列,当他们点击字段时为1,当他们没有时为0。 – sidebob
请提供易于复制和粘贴表格的数据。 –