r中自己的设计矩阵答案

【问题标题】：own design matrix in rr中自己的设计矩阵
【发布时间】：2012-07-08 06:08:34
【问题描述】：

这是我的数据：

sub <- paste ("s", 1:6, sep = "")
mark1a <- c("A", "A", "B", "d1", "A", 2)
mark1b <- c("A", "B", "d1", 2, "d1", "A")
myd <- data.frame (sub, mark1a, mark1b)
myd 
     sub mark1a mark1b
1  s1      A      A
2  s2      A      B
3  s3      B     d1
4  s4     d1      2
5  s5      A     d1
6  s6      2      A

我想创建一个变量对（列）的设计矩阵 - mark1a 和 mark1b。设计矩阵将由每个唯一 (c(mark1a, mark1b)) 的长度 (unique (c(mark1a, mark1b))) 组成。然后是 1 或 2，具体取决于列中是否存在特定数字一次或两次，否则为 0 . 以下是预期输出（不是数字）：

我可以理解如何做到这一点：

【问题讨论】：

对了，你可能对paste0函数感兴趣。

标签： r loops matrix

【解决方案1】：

你可以试试这样的：

cbind(myd, t(apply(myd, 1, function(x) sapply(unique(unlist(myd[, 2:3])), function(y) sum(x==y)))))
1  s1      A      A 2 0  0 0
2  s2      A      B 1 1  0 0
3  s3      B     d1 0 1  1 0
4  s4     d1      2 0 0  1 1
5  s5      A     d1 1 0  1 0
6  s6      2      A 1 0  0 1

【讨论】：

我如何知道级别的顺序（是否按字母顺序显示，因为这里没有显示

【解决方案2】：

我想说@jmsigner 的解决方案是单线的解决方案，但我通常会被那些嵌套的apply（及其亲属）解决方案弄糊涂。

这是一个类似的解决方案：

# Identify all the levels in `mark1a` and `mark1b`
mydLevels = unique(c(levels(myd$mark1a), levels(myd$mark1b)))
# Use these levels and an anonymous function with `lapply`
temp = data.frame(lapply(mydLevels, 
                         function(x) rowSums(myd[-1] == x)+0))
colnames(temp) = mydLevels
# This gives you the correct output, but not in the order
# that you have in your original question.
cbind(myd, temp)
#   sub mark1a mark1b 2 A B d1
# 1  s1      A      A 0 2 0  0
# 2  s2      A      B 0 1 1  0
# 3  s3      B     d1 0 0 1  1
# 4  s4     d1      2 1 0 0  1
# 5  s5      A     d1 0 1 0  1
# 6  s6      2      A 1 1 0  0

【讨论】：

【解决方案3】：

首先，确保mark1a 和mark1b 列共享相同的级别：

all.levels <- levels(myd["mark1a", "mark1b"])
levels(myd$mark1a) <- all.levels
levels(myd$mark1b) <- all.levels

然后你可以计算两个频率表的总和并将其绑定到myd：

library(plyr)
cbind(myd, ddply(myd, "sub", function(x)table(x$mark1a) + table(x$mark1b))[,-1])
#   sub mark1a mark1b 2 A B d1
# 1  s1      A      A 0 2 0  0
# 2  s2      A      B 0 1 1  0
# 3  s3      B     d1 0 0 1  1
# 4  s4     d1      2 1 0 0  1
# 5  s5      A     d1 0 1 0  1
# 6  s6      2      A 1 1 0  0

【讨论】：

矩阵加法与“表”结果。我喜欢它。